Entering edit mode
Alyaa Mahmoud
▴
440
@alyaa-mahmoud-4670
Last seen 4.9 years ago
Good Morning
I am trying to analyze a single channel agilent dataset. But I cannot
seem
to find how to retrieve probe info from pData object, i.e. the matrix
I end
up with contains no probe IDs. Please advise as I am really stuck at
this
point.
gset <- getGEO("GSE48113", GSEMatrix =TRUE)
if (length(gset) > 1) idx <- grep("GPL15331", attr(gset, "names"))
else idx
<- 1
gset <- gset[[idx]]
ex <- exprs(gset)
> head (ex) [1:2,1:5]
GSM1168586 GSM1168587 GSM1168588 GSM1168589 GSM1168590
1 NA 17.60373 18.32132 18.22676 17.5361
2 NA NA NA NA NA
Any help is much appreciated
Thanks a lot
Alyaa
--
Alyaa Mahmoud
"Love all, trust a few, do wrong to none"- Shakespeare
[[alternative HTML version deleted]]