[BioC] oligo and pd.aragene.1.1.st annotation

Stephen Rolfe s.rolfe at sheffield.ac.uk
Tue Aug 7 09:27:09 CEST 2012


  Thanks for the advice. I think my confusion arose probeset vs core in 
rma (and an interaction between oligo and affy). The affy files do 
contain AGI indices, gene symbols and descriptions and I've not yet got 
those into the analysis, but am trying the alternatives below.
Cheers
Steve

On 07/08/2012 01:43, Benilton Carvalho wrote:
> and getNetAffx returns whatever Affymetrix provides with their
> annotation files... These are quite new arrays, so it's likely that on
> future releases they'll add more information. b
>
> On 2 August 2012 18:33, Marc Carlson<mcarlson at fhcrc.org>  wrote:
>> Hi Stephen,
>>
>> It sounds like you might need to make an annotation package for your new
>> platform.  You can see a vignette that shows how to do that here:
>>
>> http://www.bioconductor.org/packages/2.10/bioc/vignettes/AnnotationDbi/inst/doc/SQLForge.pdf
>>
>> Please let me know if there are any other issues that crop up.
>>
>>
>>    Marc
>>
>>
>>
>>
>> On 08/02/2012 09:20 AM, Stephen Rolfe wrote:
>>> I have a set of files from the new Aragene-1_1-st microarrays. I am using
>>> oligo to analyse these but have hit a problem with the annotation.
>>>
>>> I have the most up-to-date packages and don't hit any error messages.
>>> However, using getNetAffx doesn't seem to add annotation data to my rma
>>> normalised data set.
>>>
>>> The feature labels are all there and the probsetids have been set. All of
>>> the other data are NA.
>>>
>>>
>>>
>>> Any insight would be much appreciated. Eventually I want to analyse the
>>> data
>>> in limma.
>>>
>>> (I tried building my own pd files using pdinfoBuilder as an alternative,
>>> but
>>> hit a bunch of error messages that I couldn't overcome.)
>>>
>>> Thanks for your help
>>>
>>> Steve Rolfe
>>>
>>>
>>>
>>>
>>>
>>> R version 2.15.1 (2012-06-22) -- "Roasted Marshmallows"
>>>
>>>
>>>
>>>> library("pd.aragene.1.1.st")
>>>> library(oligo)
>>>> library(Biobase)
>>>> celfiles<-list.celfiles()
>>>> affyRaw<-read.celfiles(files=celfiles)
>>> Platform design info loaded.
>>>
>>> Reading in : WT16HypCon1_01_A05.CEL
>>>
>>> Reading in : WT16HypCon2_02_A07.CEL
>>>
>>> Reading in : WT16HypCon3_03_B05.CEL
>>>
>>> Reading in : WT16HypInf1_04_B07.CEL
>>>
>>> Reading in : WT16HypInf2_05_C05.CEL
>>>
>>> Reading in : WT16HypInf3_06_C07.CEL
>>>
>>> Reading in : WT16RootCon1_07_D05.CEL
>>>
>>> Reading in : WT16RootCon2_08_D07.CEL
>>>
>>> Reading in : WT16RootCon3_09_E05.CEL
>>>
>>> Reading in : WT16RootInf1_10_E07.CEL
>>>
>>> Reading in : WT16RootInf2_11_F05.CEL
>>>
>>> Reading in : WT16RootInf3_12_F07.CEL
>>>
>>> Reading in : WT26HypCon1_13_H07.CEL
>>>
>>> Reading in : WT26HypCon2_14_G07.CEL
>>>
>>> Reading in : WT26HypCon3_15_H05.CEL
>>>
>>> Reading in : WT26HypInf1_16_G05.CEL
>>>
>>> Reading in : WT26HypInf2_17_A05.CEL
>>>
>>> Reading in : WT26HypInf3_18_A07.CEL
>>>
>>> Reading in : WT26RootCon1_19_B05.CEL
>>>
>>> Reading in : WT26RootCon2_20_B07.CEL
>>>
>>> Reading in : WT26RootCon3_21_C05.CEL
>>>
>>> Reading in : WT26RootInf1_22_C07.CEL
>>>
>>> Reading in : WT26RootInf2_23_D05.CEL
>>>
>>> Reading in : WT26RootInf3_24_D07.CEL
>>>
>>>> eset<-rma(affyRaw,target="core")
>>> Background correcting
>>>
>>> Normalizing
>>>
>>> Calculating Expression
>>>
>>>> featureData(eset)<-getNetAffx(eset,"transcript")
>>>> eset
>>> ExpressionSet (storageMode: lockedEnvironment)
>>>
>>> assayData: 38408 features, 24 samples
>>>
>>>     element names: exprs
>>>
>>> protocolData
>>>
>>>     rowNames: WT16HypCon1_01_A05.CEL WT16HypCon2_02_A07.CEL ...
>>>
>>>       WT26RootInf3_24_D07.CEL (24 total)
>>>
>>>     varLabels: exprs dates
>>>
>>>     varMetadata: labelDescription channel
>>>
>>> phenoData
>>>
>>>     rowNames: WT16HypCon1_01_A05.CEL WT16HypCon2_02_A07.CEL ...
>>>
>>>       WT26RootInf3_24_D07.CEL (24 total)
>>>
>>>     varLabels: index
>>>
>>>     varMetadata: labelDescription channel
>>>
>>> featureData
>>>
>>>     featureNames: 13320001 13320003 ... 13545870 (38408 total)
>>>
>>>     fvarLabels: transcriptclusterid probesetid ... category (18 total)
>>>
>>>     fvarMetadata: labelDescription
>>>
>>> experimentData: use 'experimentData(object)'
>>>
>>> Annotation: pd.aragene.1.1.st
>>>
>>>> featureData(eset)
>>> An object of class "AnnotatedDataFrame"
>>>
>>>     rowNames: 13320001 13320003 ... 13545870 (38408 total)
>>>
>>>     varLabels: transcriptclusterid probesetid ... category (18 total)
>>>
>>>     varMetadata: labelDescription
>>>
>>>> varLabels(featureData(eset))
>>> [1] "transcriptclusterid" "probesetid"          "seqname"
>>>
>>>    [4] "strand"              "start"               "stop"
>>>
>>>    [7] "totalprobes"         "geneassignment"      "mrnaassignment"
>>>
>>> [10] "swissprot"           "unigene"             "gobiologicalprocess"
>>>
>>> [13] "gocellularcomponent" "gomolecularfunction" "pathway"
>>>
>>> [16] "proteindomains"      "crosshybtype"        "category"
>>>
>>>> pData(featureData(eset))[1000:1002,]
>>>            transcriptclusterid probesetid seqname strand start stop
>>> totalprobes
>>>
>>> 13321999            13321999   13321999<NA>     <NA>      NA   NA
>>> NA
>>>
>>> 13322001            13322001   13322001<NA>     <NA>      NA   NA
>>> NA
>>>
>>> 13322003            13322003   13322003<NA>     <NA>      NA   NA
>>> NA
>>>
>>>                                             geneassignment
>>>
>>> 13321999 --- // --- // intron-exon control // --- // ---
>>>
>>> 13322001 --- // --- // intron-exon control // --- // ---
>>>
>>> 13322003 --- // --- // intron-exon control // --- // ---
>>>
>>>
>>> mrnaassignment
>>>
>>> 13321999 --- // --- // intron-exon control // --- // --- // --- // --- //
>>> --- // --- //
>>>
>>> 13322001 --- // --- // intron-exon control // --- // --- // --- // --- //
>>> --- // --- //
>>>
>>> 13322003 --- // --- // intron-exon control // --- // --- // --- // --- //
>>> --- // --- //
>>>
>>>            swissprot unigene gobiologicalprocess gocellularcomponent
>>>
>>> 13321999        NA      NA                  NA                  NA
>>>
>>> 13322001        NA      NA                  NA                  NA
>>>
>>> 13322003        NA      NA                  NA                  NA
>>>
>>>            gomolecularfunction pathway proteindomains crosshybtype
>>>
>>> 13321999                  NA      NA             NA           NA
>>>
>>> 13322001                  NA      NA             NA           NA
>>>
>>> 13322003                  NA      NA             NA           NA
>>>
>>>                    category
>>>
>>> 13321999 normgene->intron
>>>
>>> 13322001 normgene->intron
>>>
>>> 13322003 normgene->intron
>>>
>>>
>>>          [[alternative HTML version deleted]]
>>>
>>> _______________________________________________
>>> 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
>>
>> _______________________________________________
>> 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
> _______________________________________________
> 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



More information about the Bioconductor mailing list