[BioC] Target file for Limma
Caimiao Wei
caimiaow at u.washington.edu
Tue Nov 16 23:12:59 CET 2004
I used the following code lines to read in Agilent data within limma:
targets=readTargets("C:/caimiao/two_channel_one_channel/target1.txt")
> targets
SlideNumber FileName Cy3 Cy5
5ugDay1s1 1 5ugDay1s1.txt Sample1 Uref
5ugDay1s2 2 5ugDay1s2.txt Sample1 Uref
5ugDay1s3 3 5ugDay1s3.txt Sample1 Uref
5ugDay1s4 4 5ugDay1s4.txt Uref Sample1
5ugDay1s5 5 5ugDay1s5.txt Uref Sample1
5ugDay1s6 6 5ugDay1s6.txt Uref Sample1
RG=read.maimages(files=targets$FileName, source="agilent", path="C:/caimiao/day1",
names=NULL, columns=list(gene="ProbeName", col="Col", row="Row",
Rf="rMeanSignal", Rb="gBGMeanSignal",Gf="gMeanSignal", Gb="gBGMeanSignal"),
other.columns=list( col="Col", row= "Row", controltype="ControlType"),
annotation=list (gene="ProbeName"))
> RG$targets
FileName
5ugDay1s1 5ugDay1s1.txt
5ugDay1s2 5ugDay1s2.txt
5ugDay1s3 5ugDay1s3.txt
5ugDay1s4 5ugDay1s4.txt
5ugDay1s5 5ugDay1s5.txt
5ugDay1s6 5ugDay1s6.txt
What I don't understand is why the output for the command "RG$targets" does
not have the Cy3 and Cy5 information?
Thanks,
Caimiao
More information about the Bioconductor
mailing list