Question: problems in reading a file in AgiMicroRna package
0
18 months ago by
National Institute of Technology
neelakandanbt0 wrote:

## the

>dd.micro1=readMicroRnaAFE(targets1.micro,verbose=TRUE)

reading file  2  -  GSM1511527_US10353831_252929710061_S01_miRNA_1010_Sep10_1_2.txt 
reading file  3  -  GSM1511528_US10353831_252929710061_S01_miRNA_1010_Sep10_1_3.txt
reading file  6  -  GSM1511531_US10353831_252929710061_S01_miRNA_1010_Sep10_2_3.txt 
reading file  7  -  GSM1511532_US10353831_252929710061_S01_miRNA_1010_Sep10_2_4.txt
Error in Newagi[[a]][, n] <- obj[, columns[[a]]] :
number of items to replace is not a multiple of replacement length 2  -  

## the command and the output is given. i dont know what went wrong kindly help me in this.

uRNAList:

dd$TGS: 'gTotalGeneSignal' dd$TPS: 'gTotalProbeSignal'

dd$meanS: 'gMeanSignal' dd$procS: 'gProcessedSignal'

Warning message:

In scan(file = file, what = what, sep = sep, quote = quote, dec = dec,  :

embedded nul(s) found in input

modified 18 months ago by Axel Klenk920 • written 18 months ago by neelakandanbt0

Hi Neela,

I don't know the answer to your problem I'm afraid. The only thing I can think of is when I ran these arrays a few years back there was an issue initially with the feature extraction on the scanner. It essentially defaulted to saving a minimal reduced output rather than a full output. That might be where your problems are? Have a dig around the software and make sure you're saving the full results.

Good luck

0
18 months ago by
Axel Klenk920
Switzerland
Axel Klenk920 wrote:

Dear Neela,

Agilent's reduced output is lacking the 'gMeanSignal' column and this will cause an issue. I'm not sure it is the source of your problem, though, since I'd expect it to affect *all* files rather than just the 13th.

You can inspect the source code of AgiMicroRna::readMicroRnaAFE by typing it like that, i.e., without parentheses,

and you will see that it essentially passes your target data.frame on to a function read.agiMicroRna but is also hardcoding the column names.

To find out if the missing 'gMeanSignal' column is causing your problem, try:

dd <- AgiMicroRna:::read.agiMicroRna(targets1.micro, columns = list(TGS = "gTotalGeneSignal",
TPS = "gTotalProbeSignal", meanS = "gMedianSignal", procS = "gProcessedSignal"),
other.columns = list(IsGeneDetected = "gIsGeneDetected",
IsSaturated = "gIsSaturated", IsFeatNonUnifOF = "gIsFeatNonUnifOL",
IsFeatPopnOL = "gIsFeatPopnOL", BGKmd = "gBGMedianSignal",
BGKus = "gBGUsed"), annotation = c("ControlType",
"ProbeName", "GeneName"), verbose = TRUE)

Note the three colons and the change from 'gMeanSignal' to 'gMedianSignal'. Also note that it's a hack. No warranty. :-)

Good luck,

- axel