Question: Getting information about samples(.CEL files) from GEO using R
gravatar for osieman52
13 months ago by
osieman5220 wrote:

Hello all,

Using the GEOquery library in R I have managed to download 850 samples I am interested in. I used the following code for this:

x <- as.list(scan("fileWithGEOsampleAccesionCodes.txt", what ="" , sep = "" ))

data= lapply(x,getGEOSuppFiles)

I have read the CEL files into an affyBatch using - > r <- ReadAffy(celfile.path = pathToFiles,compress = TRUE)

I want to look which micro-array platform is being for each sample. With micro-array platform I mean --> AffyA , AffyPlus2 , also know as hgu133a ,hgu133plus2 .

How can I do this using the data I already have or using the accessioncodes? Basically I want to know the annotation for each sample (.CEL file) .

Thank you in advance!

micro-array R • 598 views
ADD COMMENTlink modified 13 months ago by Kevin Blighe39k • written 13 months ago by osieman5220
gravatar for Kevin Blighe
13 months ago by
Kevin Blighe39k
Republic of Ireland
Kevin Blighe39k wrote:

If you just have a list of the CEL files, then you should be able to get this via the CDF entry in the CEL header with:

 [1] "SampleFiles//AC31_(BovGene-1_0-st).CEL"          
 [2] "SampleFiles//AC32+AC33_(BovGene-1_0-st).CEL"     
 [3] "SampleFiles//AC33+AC32_(BovGene-1_0-st).CEL"     
 [4] "SampleFiles//AC33+OA23_(BovGene-1_0-st).CEL"     
 [5] "SampleFiles//AC34_(BovGene-1_0-st).CEL"          
 [6] "SampleFiles//NL13_(BovGene-1_0-st).CEL"          
 [7] "SampleFiles//NL14_(BovGene-1_0-st).CEL"          
 [8] "SampleFiles//NL16_(BovGene-1_0-st).CEL"          
 [9] "SampleFiles//NL17_(BovGene-1_0-st).CEL"          
[10] "SampleFiles//NM1_(BovGene-1_0-st).CEL"           
[11] "SampleFiles//NM2_(BovGene-1_0-st).CEL"           
[12] "SampleFiles//NM3_(BovGene-1_0-st).CEL"           
[13] "SampleFiles//NM41_(BovGene-1_0-st).CEL"          
[14] "SampleFiles//NM42_(BovGene-1_0-st).CEL"          
[15] "SampleFiles//NM43_(BovGene-1_0-st).CEL"          
[16] "SampleFiles//NM44_(BovGene-1_0-st).CEL"          
[17] "SampleFiles//NM4_(BovGene-1_0-st).CEL"           
[18] "SampleFiles//OA21_(BovGene-1_0-st).CEL"          
[19] "SampleFiles//OA22+OA23+OA24_(BovGene-1_0-st).CEL"
[20] "SampleFiles//OA22+OA24_(BovGene-1_0-st-v1).CEL"  
[21] "SampleFiles//OA23+AC32_(BovGene-1_0-st-v1).CEL"  
[22] "SampleFiles//OA24+OA22_(BovGene-1_0-st-v1).CEL"  

data.frame(CELFiles, unlist(lapply(CELFiles, function(x) affyio::read.celfile.header(x)[["cdfName"]])))
                                           CELFiles  Type
1            SampleFiles//AC31_(BovGene-1_0-st).CEL  BovGene-1_0-st
2       SampleFiles//AC32+AC33_(BovGene-1_0-st).CEL  BovGene-1_0-st
3       SampleFiles//AC33+AC32_(BovGene-1_0-st).CEL  BovGene-1_0-st
4       SampleFiles//AC33+OA23_(BovGene-1_0-st).CEL  BovGene-1_0-st
5            SampleFiles//AC34_(BovGene-1_0-st).CEL  BovGene-1_0-st
6            SampleFiles//NL13_(BovGene-1_0-st).CEL  BovGene-1_0-st
7            SampleFiles//NL14_(BovGene-1_0-st).CEL  BovGene-1_0-st
8            SampleFiles//NL16_(BovGene-1_0-st).CEL  BovGene-1_0-st
9            SampleFiles//NL17_(BovGene-1_0-st).CEL  BovGene-1_0-st
10            SampleFiles//NM1_(BovGene-1_0-st).CEL  BovGene-1_0-st
11            SampleFiles//NM2_(BovGene-1_0-st).CEL  BovGene-1_0-st
12            SampleFiles//NM3_(BovGene-1_0-st).CEL  BovGene-1_0-st
13           SampleFiles//NM41_(BovGene-1_0-st).CEL  BovGene-1_0-st
14           SampleFiles//NM42_(BovGene-1_0-st).CEL  BovGene-1_0-st
15           SampleFiles//NM43_(BovGene-1_0-st).CEL  BovGene-1_0-st
16           SampleFiles//NM44_(BovGene-1_0-st).CEL  BovGene-1_0-st
17            SampleFiles//NM4_(BovGene-1_0-st).CEL  BovGene-1_0-st
18           SampleFiles//OA21_(BovGene-1_0-st).CEL  BovGene-1_0-st
19 SampleFiles//OA22+OA23+OA24_(BovGene-1_0-st).CEL  BovGene-1_0-st
20   SampleFiles//OA22+OA24_(BovGene-1_0-st-v1).CEL  BovGene-1_0-st
21   SampleFiles//OA23+AC32_(BovGene-1_0-st-v1).CEL  BovGene-1_0-st
22   SampleFiles//OA24+OA22_(BovGene-1_0-st-v1).CEL  BovGene-1_0-st
ADD COMMENTlink written 13 months ago by Kevin Blighe39k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2233 users visited in the last hour