AgiMicroRna package error message
1
0
Entering edit mode
@wang-yonghong-nihnci-1287
Last seen 8.3 years ago
Dear Bioconductor users: I am trying to use AgiMicroRna to perform the microRNA array analysis (Agilent platform) these days, always have some problems. The main problem is the error when I used the function "filterMicroRna()". Here is the error message: FILTERING BY ControlType Error in data.frame(as.character(PROBE_ID), as.character(GENE_ID), as.character(probe.chr), : arguments imply differing number of rows: 795, 0 When I loaded in the data, I used the "read.maimages()" so I can specify which columns to import. I have searched the mailing list and Pedro has suggested to use "readMicroRnaAFE()" to load the data. The first time I tried, for the "genes" information, it include "sequence", "ProbeUID", "ControlType" rather than the required "ControlType, ProbeName, and GeneName". It seems to me the package picked up the column information based on the position rather than column header. After I remove the "sequence, ProbeUID" columns from the raw data files, checked the loaded data, the "genes" information is correct now, but still when I ran the "filterMicroRna() function, same error appear. I have run the sample data set come from the package, everything works fine. I have the following several questions: First: how can I solve above problem? Second: I doubt my file format is not exact same to the sample one (probably with more other column information will interfere the correct loading data), if this is the case, who can show me where to download the sample FE file used in the AgiMicroRna package (one is enough, just for the purpose of comparison). third: there are two functions, one is the tgsNormalization() with the parameter "quantile", another is rmaMicroRna(). I checked the output, they are different. Since RMA is also quantile normalization, I think I miss something here, if anyone can tell me the difference of these two functions, that would be very helpful. Fourth: This is about the agilent FE file. How the "gProcessedSignal" file being generated, I mean how this is different to the other intensity information? Fifth: Since this is my first time to work on the microRNA dataset, any suggests about the data analysis would be very helpful Thanks in advance Best YH
Normalization microRNA AgiMicroRna Normalization microRNA AgiMicroRna • 588 views
ADD COMMENT
0
Entering edit mode
@wang-yonghong-nihnci-1287
Last seen 8.3 years ago
Thanks for Richard's previous email to Neel. The problem I have to run "filterMicroRna()" is that Agilent changed the header for "chr_coord" to "probe_mapping" in its new version. After I changed the name back. everything works fine so far. Best YH
ADD COMMENT

Login before adding your answer.

Traffic: 215 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6