[BioC] Loading data into AgiMicroRna

Steve Taylor stephen.taylor at imm.ox.ac.uk
Fri Oct 15 14:56:48 CEST 2010


Hi again,

>
> I am copying the mail from Pedro to the mailing list. I am thinking if
> you can read the files using read.maimages function it should work.
>

I have loaded my data in AgiMicroRna using:

dd=read.maimages(files=targets.micro$FileName,source="agilent",
     columns=list(Rf="gTotalGeneSignal",
     Gf="gTotalProbeSignal",
     Rb="gMedianSignal",
     Gb="gProcessedSignal"),
     other.columns=list(IsGeneDetected="gIsGeneDetected",
     IsSaturated="gIsSaturated",
     IsFeatNonUnifOF="gIsFeatNonUnifOL",
     IsFeatPopnOL="gIsFeatPopnOL",
     BGKmd="gBGMedianSignal"),
     annotation = c( "ControlType", "ProbeName"),
     verbose=TRUE,sep="\t",quote="")


There was no Rb="gMeanSignal", but I assume this should be ok.

Everything looked ok until I tried to normalize the data.

> ddTGS = tgsMicroRna(dd.micro, half = TRUE, makePLOT = FALSE,verbose = FALSE)
> ddNORM = tgsNormalization(ddTGS, "quantile", makePLOTpre = FALSE,makePLOTpost = FALSE, targets.micro, verbose = TRUE)
Error in xy.coords(x, y) : 'x' and 'y' lengths differ

I see this problem has come up before on the list:

http://article.gmane.org/gmane.science.biology.informatics.conductor/31182/match=xy+coords+x+y+lengths+differ+agimicrorna

But I was wondering if there is a solution for it?  Apologies if I have missed something.

Thanks,

Steve



More information about the Bioconductor mailing list