Displaying 1 result from an estimated 1 matches for "gsm286039".
Did you mean:
gsm286031
2010 Mar 29
1
stuck with affy / limma
...outcomes.
I did this based on the affy and limma manuals like:
targets file:
Name FileName Target
0h1 GSM286031.CEL control
0h2 GSM286032.CEL control
0h3 GSM286033.CEL control
3h1 GSM286034.CEL three
3h2 GSM286035.CEL three
3h3 GSM286036.CEL three
6h1 GSM286037.CEL six
6h2 GSM286038.CEL six
6h3 GSM286039.CEL six
library(affy)
library(limma)
library(vsn)
pd <- read.AnnotatedDataFrame("er_for_affy.txt", header = TRUE, row.names =
2)
pData(pd)
#### load
a <- ReadAffy(filenames = rownames(pData(pd)), phenoData = pd, verbose =
TRUE)
#### normalize
x <- expresso(a, bg.correct = FALS...