AgiMicroRna problem
0
0
Entering edit mode
Paulo Nuin ▴ 200
@paulo-nuin-3012
Last seen 7.7 years ago
Hi We're trying to analyse some Agilent microRNA chips with AgiMicroRna without much success. Our chips are V1 and I saw in the documentation that the package is intended to V2. Anyway, we were able to work around the columns issue, but know we are stuck on another problem. Our dataset is composed of 30 chips, each one from a different sample. Eleven of those 30 are one treatment, while the remainder 19 are in another treatment. We created a targets file as specified by the manual with 4 columns: Filename Treatment GErep Subject fileA A 1 1 fileB A 1 2 fileC B 2 3 fileD B 2 4 .... ... ... 30 We were able to read the images into a variable using dd=read.maimages(files=targets$FileName,source="agilent", columns=list(Rf="gTotalGeneSignal", Gf="gTotalProbeSignal", Rb="gTotalGeneSignal", Gb="gProcessedSignal"), other.columns=list(IsGeneDetected="gIsGeneDetected", IsSaturated="gIsSaturated", IsFeatNonUnifOF="gIsFeatNonUnifOL", IsFeatPopnOL="gIsFeatPopnOL", probe_mappings="probe_mappings", BGKmd="gBGMedianSignal", BGKus="gBGUsed"), annotation = c( "ControlType","ProbeName","GeneName"), verbose=TRUE,sep="\t",quote="") names(dd$others) will output names(dd\$other) [1] "gIsGeneDetected" "gIsSaturated" "gIsFeatNonUnifOL" "gIsFeatPopnOL" "gBGMedianSignal" But when we run the rmaMicroRna with ddTGS.rma = rmaMicroRna(dd, normalize = TRUE, background = FALSE) we get the following error: Error in split.default(0:(length(pNList) - 1), pNList) : Group length is 0 but data length > 0 Any help is very appreciated. Thanks in advance Paulo Nuin
Annotation microRNA Annotation microRNA • 672 views