[BioC] oligo package attaching phenodata for HuGene2.0 ST array

James W. MacDonald jmacdon at uw.edu
Wed Jul 24 23:07:47 CEST 2013


You have to create the AnnotatedDataFrame first.

pd<-read.AnnotatedDataFrame("covdesc.txt",header=TRUE,row.name="Filename" ,sep="\t")
x <- varMetadata(pd)

Best,

Jim


On 7/24/2013 5:03 PM, Jia, Li (NIH/NCI) [C] wrote:
> Hi Jim,
>
> Thanks for the prompt reply. When I tested on the first command you sent
> to us, it generated an error.
>
>> x<- varMetadata(pd)
> Error in (function (classes, fdef, mtable)  :
>    unable to find an inherited method for function ‘varMetadata’ for
> signature ‘"data.frame"’
>
>
>
> Please advise again. Thanks a lot.
>
> Li
>
> On 7/24/13 4:58 PM, "James W. MacDonald"<jmacdon at uw.edu>  wrote:
>
>>
>> On 7/24/2013 4:40 PM, Guest [guest] wrote:
>>> Hi,
>>>
>>> We are working with affymetrix HuGene2.0 ST arrays and are trying to
>>> read in our CEL files with the oligo package. However, we are having
>>> trouble integrating our phenodata to the expression data. Please see our
>>> code and error message below. The covdesc.txt is the file with our
>>> phenodata.
>>>
>>> code:
>>>
>>>
>>> pd<-read.AnnotatedDataFrame("covdesc.txt",header=TRUE,row.name="Filename"
>>> ,sep="\t")
>>>
>>> OligoRaw<-read.celfiles(filenames=list.celfiles(), phenoData=pd,
>>> verbose = TRUE)
>>>
>>>> pData(pd)
>>>           Filename Treatment Transfection
>>> 1      ACsiNS.CEL   control     ns_siRNA
>>> 2     ACSsi03.CEL   control      siRNA03
>>> 3    ANCSsi03.CEL   Treated      siRNA03
>>> 4  ANCSsiNS_2.CEL   Treated     ns_siRNA
>>> 5      DCsi03.CEL   control     ns_siRNA
>>> 6      DCsiNS.CEL   control      siRNA03
>>> 7    DNCSsi03.CEL   Treated      siRNA03
>>> 8    DNCSsiNS.CEL   Treated     ns_siRNA
>>> 9      ECsi03.CEL   control     ns_siRNA
>>> 10     ECsiNS.CEL   control      siRNA03
>>> 11   ENCSsi03.CEL   Treated      siRNA03
>>> 12   ENCSsiNS.CEL   Treated     ns_siRNA
>>>
>>> Error:
>>>> OligoRaw<-read.celfiles(filenames=list.celfiles(), phenoData=pd1,
>>>> verbose = TRUE)
>>> Platform design info loaded.
>>> Reading in : ACsiNS.CEL
>>> Reading in : ACSsi03.CEL
>>> Reading in : ANCSsi03.CEL
>>> Reading in : ANCSsiNS_2.CEL
>>> Reading in : DCsi03.CEL
>>> Reading in : DCsiNS.CEL
>>> Reading in : DNCSsi03.CEL
>>> Reading in : DNCSsiNS.CEL
>>> Reading in : ECsi03.CEL
>>> Reading in : ECsiNS.CEL
>>> Reading in : ENCSsi03.CEL
>>> Reading in : ENCSsiNS.CEL
>>> Error in validObject(out) : invalid class “GeneFeatureSet” object:
>>>     'NChannelSet' varMetadata must have a 'channel' column
>>>
>>> Any suggestions will be appreciated!
>> There is a varMetadata slot to the AnnotatedDataFrame, and oligo expects
>> that there will be a column called 'channel'. A really cheap way to fix
>> this would be to do something like
>>
>> x<- varMetadata(pd)
>> x<- data.frame(x, channel = "_ALL_")
>> varMetadata(pd)<- x
>>
>> and then proceed with read.celfiles().
>>
>> Best,
>>
>> Jim
>>
>>
>>> Thanks!
>>>
>>>
>>>    -- output of sessionInfo():
>>>
>>>> sessionInfo()
>>> R version 3.0.0 (2013-04-03)
>>> Platform: x86_64-apple-darwin10.8.0 (64-bit)
>>>
>>> locale:
>>> [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
>>>
>>> attached base packages:
>>> [1] parallel  stats     graphics  grDevices utils     datasets  methods
>>>   base
>>>
>>> other attached packages:
>>>    [1] pd.hugene.2.0.st_3.8.0     RSQLite_0.11.4             DBI_0.2-7
>>>    [4] oligo_1.24.1               oligoClasses_1.22.0
>>> affyPLM_1.36.0
>>>    [7] preprocessCore_1.22.0      arrayQualityMetrics_3.16.0
>>> affyQCReport_1.38.0
>>> [10] lattice_0.20-15            simpleaffy_2.36.1          gcrma_2.32.0
>>> [13] genefilter_1.42.0          affy_1.38.1
>>> Biobase_2.20.1
>>> [16] BiocGenerics_0.6.0         BiocInstaller_1.10.2
>>>
>>> loaded via a namespace (and not attached):
>>>    [1] affxparser_1.32.3    affyio_1.28.0        annotate_1.38.0
>>> AnnotationDbi_1.22.6
>>>    [5] beadarray_2.10.0     BeadDataPackR_1.12.0 Biostrings_2.28.0
>>> bit_1.1-10
>>>    [9] Cairo_1.5-2          cluster_1.14.4       codetools_0.2-8
>>> colorspace_1.2-2
>>> [13] ff_2.2-11            foreach_1.4.1        GenomicRanges_1.12.4
>>> grid_3.0.0
>>> [17] Hmisc_3.12-2         hwriter_1.3          IRanges_1.18.2
>>> iterators_1.0.6
>>> [21] KernSmooth_2.23-10   latticeExtra_0.6-24  limma_3.16.6
>>> plyr_1.8
>>> [25] RColorBrewer_1.0-5   reshape2_1.2.2       rpart_4.1-1
>>> setRNG_2011.11-2
>>> [29] splines_3.0.0        stats4_3.0.0         stringr_0.6.2
>>> survival_2.37-4
>>> [33] SVGAnnotation_0.93-1 tools_3.0.0          vsn_3.28.0
>>> XML_3.95-0.2
>>> [37] xtable_1.7-1         zlibbioc_1.6.0
>>>
>>> --
>>> Sent via the guest posting facility at bioconductor.org.
>>>
>>> _______________________________________________
>>> 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
>>

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



More information about the Bioconductor mailing list