[BioC] edgeR query

[James W. MacDonald]

Hi Shravanthi,

 > library(GEOquery)
 > dat <- getGEO("GSE9936")[[1]]
 > dat
ExpressionSet (storageMode: lockedEnvironment)
assayData: 22283 features, 105 samples
   element names: exprs
protocolData: none
phenoData
   sampleNames: GSM251354 GSM251355 ... GSM251458 (105 total)
   varLabels: title geo_accession ... data_row_count (37 total)
   varMetadata: labelDescription
featureData
   featureNames: 1007_s_at 1053_at ... AFFX-TrpnX-M_at (22283 total)
   fvarLabels: ID GB_ACC ... Gene Ontology Molecular Function (16 total)
   fvarMetadata: Column Description labelDescription
experimentData: use 'experimentData(object)'
Annotation: GPL96

 > table(pData(phenoData(dat))[,11])

                    treated with 300nM genistein
                                              17
treated with 300nM genistein and 3uM ICI182,780
                                              13
                      treated with 300nM S-Equol
                                              13
  treated with 300nM S-Equol and 3uM ICI 182,780
                                               9
               treated with 6nM 17beta-estradiol
                                              14
                      treated with 6nM genistein
                                              18
                    treated with vehicle control
                                              21

Note that the limma package is quite happy to use an ExpressionSet as 
input, and the 11th column of the phenoData slot contains the treatment. 
That should get you well on your way. You should at the very least read 
the majority of the limma User's Guide.

Best,

Jim



On 9/5/2012 8:16 AM, Shravanthi P wrote:
> Thanks a lot for correcting my error. I have downloaded the Limma
> package(it seems to be majorly for spotted microarray as opposed to
> oligonucleotides microarry). I am having trouble even reading my data
> into R. I have the tab delimited format file for a series (eg GSE9936)
> .How do I read and  initialize it to an object in R? And how to I
> define groups within this file for a comparison.Limma user manual
> directed me to use the affy package for normalizing data.Can I not use
> the already available normalized data ?
> I have a time constraint and was hoping to avoid this step.
>
> My aim is to find the up regulated and downregulated genes in breast
> cancer cells when treated with Genistein on a time progression basis
> (with a constant concentration)
>
> On 9/5/12, Steve Lianoglou<mailinglist.honeypot at gmail.com>  wrote:
>> Hi Shravanthi,
>>
>> On Tue, Sep 4, 2012 at 6:33 AM, Shravanthi P<shravanthipsg at gmail.com>
>> wrote:
>>> I am new to using R and my project is on differential gene expression .I
>>> was told to use Edger but I am having trouble reading a tab delimited
>>> format data for affy oligonucleotide microarray data .
>>> How can I read the data into R ? and how can I identify the upregulated
>>> and
>>> downregulated genes  ?
>> The bad news is that you've got a long way to go before you can get
>> from where you are to where you want to be. The good news is that
>> there are a lot of resources available for your self study that can
>> help you get to where you want to be.
>>
>> It would be in your best interest to learn more about R first (start
>> with an intro to R) before you try to do anything more serious.
>>
>> After you do that, and you then read the (very helpful) edgeR user
>> manual (as Wenhuo already suggested) -- you'll probably come to the
>> realization that edgeR is the wrong tool for this job. You have
>> microarray data and edgeR is used for sequencing data.
>>
>> In this case, you'll find that the limma package is probably what
>> you're after. Fortunately for you, there is a very thorough limma user
>> guide that you can read through that gives you some idea of how to
>> process different types of microarray data. The goods at the link
>> below (and the associated packages) are worth some study as well:
>>
>> http://bioconductor.org/help/workflows/oligo-arrays/
>>
>> HTH,
>> -steve
>>
>> --
>> Steve Lianoglou
>> Graduate Student: Computational Systems Biology
>>   | Memorial Sloan-Kettering Cancer Center
>>   | Weill Medical College of Cornell University
>> Contact Info: http://cbio.mskcc.org/~lianos/contact
>>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor

-- 
James W. MacDonald, M.S.
Biostatistician
University of Washington
Environmental and Occupational Health Sciences
4225 Roosevelt Way NE, # 100
Seattle WA 98105-6099