[BioC] Problems with affycomp - quantile.default

Markus Schmidberger schmidb at ibe.med.uni-muenchen.de
Mon Aug 13 10:26:38 CEST 2007


Hi Jim,

thanks a lot, now it works.
But what is wrong by reading the phenodata-File

 > pd95<-read.AnnotatedDataFrame(filename="spikein_pdata.txt", 
header=TRUE, row.names="filename", sep=" ")

 > colnames(pData(pd95))
 [1] "X37777_at" "X684_at"   "X1597_at"  "X38734_at" "X39058_at" "X36311_at"
 [7] "X36889_at" "X1024_at"  "X36202_at" "X36085_at" "X40322_at" "X407_at" 
[13] "X1091_at"  "X1708_at"  "X33818_at" "X546_at" 

Attached the pdata-File from the affycompII website 
(http://affycomp.biostat.jhsph.edu/)
(set="\t" also does not work)

Thanks
Markus


James W. MacDonald schrieb:
> Are you sure you have the correct phenoData? What do you get from
>
> colnames(pData(eset))?
>
> You should get this:
>
> > colnames(pData(eset))
>  [1] "37777_at" "684_at"   "1597_at"  "38734_at" "39058_at" "36311_at"
>  [7] "36889_at" "1024_at"  "36202_at" "36085_at" "40322_at" "407_at"
> [13] "1091_at"  "1708_at"  "33818_at" "546_at"
>
> If you don't, you have the incorrect phenoData, and you should use
>
> pd <- as(data(spikein.phenodata), "AnnotatedDataFrame")
>
> Best,
>
> Jim
>
>
>
> Markus Schmidberger wrote:
>> Hello,
>>
>> Phenodata do not help. I changed my code to this:
>>
>> library(affy)
>> library(affycomp)
>> pathRawData <- "~/Microarray/hgu95a-spikein/rawdata"
>> pathPhenoData <- "~/Microarray/hgu95a-spikein/spikein_pdata.txt"
>> celFile <- list.celfiles(path=pathRawData,full.names=TRUE);
>> pd<-read.AnnotatedDataFrame(filename=pathPhenoData, header=TRUE, 
>> row.names="filename", sep=" ")
>> affyBatch <- ReadAffy(filenames=celFile, phenoData=pd);
>> eset <- rma(affyBatch)
>> assessFC(eset,method.name="RMA")
>>
>> and I still have the same error:
>> Fehler in quantile.default(m[-spikeindex], prob = probs) :
>>         missing values and NaN's not allowed if 'na.rm' is FALSE
>>
>> (same error with the hgu133 dataset)
>>
>> What I have to change?
>>
>> Thanks,
>> Markus
>>
>> James MacDonald schrieb:
>>> Hi Markus,
>>>
>>> Markus Schmidberger wrote:
>>>> Hello,
>>>>
>>>> I want to do some comparison with affycomp and use the affycompII 
>>>> data. This is my code
>>>>
>>>> library(affy)
>>>> library("affycomp")
>>>> path <- "/Microarray/hgu95a-spikein/rawdata"
>>>> celFile <- list.celfiles(path=path,full.names=TRUE);
>>>> affyBatch <- ReadAffy(filenames=celFile);
>>>> eset1 <- rma(affyBatch)
>>>> assessFC(eset1,method.name="RMA.2")
>>>>
>>>> and I get the ERROR:
>>>> Fehler in quantile.default(m[-spikeindex], prob = probs) :
>>>>         missing values and NaN's not allowed if 'na.rm' is FALSE
>>> If you look at the code for assessFC(), you will see that the 
>>> spikeindex comes from the phenoData slot of the ExpressionSet object 
>>> (in particular the genenames column, which I assume lists the 
>>> probesets that were spiked in). If you don't have this information 
>>> in your phenoData slot, then assessFC() won't work, and I assume 
>>> this to be true of the other assessXX() functions as well.
>>>
>>> The phenoData can be accessed using data(spikein.phenodata). 
>>> Unfortunately this is an old-style phenoData object instead of an 
>>> AnnotatedDataFrame, so you have to update first.
>>>
>>> pd <- data(spikein.phenodata)
>>> pd <- as(pd, "AnnotatedDataFrame")
>>> ## you will see warnings here that I think you can ignore
>>> affyBatch <- ReadAffy(filenames=celFile, phenoData=pd)
>>>
>>> should get you going.
>>>
>>> Best,
>>>
>>> Jim
>>>
>>>
>>>> How to fix this?
>>>>
>>>> R version 2.5.0 (2007-04-23)
>>>> i686-pc-linux-gnu
>>>> locale:
>>>> LC_CTYPE=de_DE.UTF-8;LC_NUMERIC=C;LC_TIME=de_DE.UTF-8;LC_COLLATE=de_DE.UTF-8;LC_MONETARY=de_DE.UTF-8;LC_MESSAGES=de_DE.UTF-8;LC_PAPER=de_DE.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=de_DE.UTF-8;LC_IDENTIFICATION=C 
>>>>
>>>> attached base packages:
>>>> [1] "tools"     "stats"     "graphics"  "grDevices" "utils"     
>>>> "datasets"
>>>> [7] "methods"   "base"    other attached packages:
>>>> hgu95acdf  affycomp      affy    affyio   Biobase
>>>>  "1.16.0"  "1.12.0"  "1.14.1"   "1.4.0"  "1.14.0"
>>>>
>>>> Thanks,
>>>> Markus
>>>>
>>
>>
>


-- 
Dipl.-Tech. Math. Markus Schmidberger

Ludwig-Maximilians-Universität München
IBE - Institut für medizinische Informationsverarbeitung,
Biometrie und Epidemiologie
Marchioninistr. 15, D-81377 Muenchen
URL: http://ibe.web.med.uni-muenchen.de 
Mail: Markus.Schmidberger [at] ibe.med.uni-muenchen.de
Tel: +49 (089) 7095 - 4599

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: spikein_pdata.txt
Url: https://stat.ethz.ch/pipermail/bioconductor/attachments/20070813/b92a8a61/attachment.txt 


More information about the Bioconductor mailing list