[BioC] simpleaffy 2.14.02, setQCEnvironment

Crispin Miller CMiller at picr.man.ac.uk
Thu Oct 11 13:10:31 CEST 2007


Hi Stefan,
I've have just fixed this - the changes will appear in version 2.14.04...
Cheers,

Crispin


On 11/10/07 10:33, "Stefan Thomsen" <stt26 at cam.ac.uk> wrote:

> Dear all,
> 
> I am trying to perform some simpleaffy analyses on Drosophila Genome 2.0
> genechips. Following the instructions of the simpleaffy 2.14.02 vignette I
> set up the following .qcdef file:
> 
> array drosophila2cdf
> alpha1 0.05
> alpha2 0.065
> spk bioB AFFX-r2-Ec-bioB-3_at
> spk bioC AFFX-r2-Ec-bioC-3_at
> spk bioD AFFX-r2-Ec-bioD-3_at
> spk creX AFFX-r2-P1-cre-3_at
> ratio actin3/actin5 AFFX-Dros-ACTIN_3_at AFFX-Dros-ACTIN_5_at
> ratio actin3/actinM AFFX-Dros-ACTIN_3_at AFFX-Dros-ACTIN_M_f_at
> ratio gapdh3/gapdh5 AFFX-Dros-GAPDH_3_at AFFX-Dros-GAPDH_5_at
> ratio gapdh3/gapdhM AFFX-Dros-GAPDH_3_at AFFX-Dros-GAPDH_M_at
> 
> I saved it as drosophila2cdf.qcdef.txt and removed the file extension .txt
> afterwards.
> 
> When I try to run qc() on my affybatch and affybatch.mas5 (created by
> calling call.exprs(affybatch, "mas5") objects I get the following error
> message:
> 
>> library(simpleaffy)
> Loading required package: affy
> Loading required package: Biobase
> Loading required package: tools
> Loading required package: affyio
> Loading required package: preprocessCore
> Loading required package: genefilter
> Loading required package: survival
> Loading required package: splines
> Loading required package: gcrma
> Loading required package: matchprobes
> 
>> setQCEnvironment("drosophila2cdf", "/Users/Stefan")
>> qc.data=qc(affybatch.raw, affybatch.mas5)
> Error in setQCEnvironment(cdfn) :
>   Could not find array definition file ' drosophila2cdf.qcdef '. Simpleaffy
> does not know the QC parameters for this array type. See the package
> vignette for details about how to specify QC parameters manually.
> 
> 
> I confess, my programming skills are rather limited. However, I tried to
> look at setQCEnvironment() and its subfunction qc.read.file(). If I
> directly assign the .qcdef file including its path to the 'fn' argument,
> qc.read.file seems to be able to open, read and print the lines of my file.
> 
> Thus, why do keep on running into this error?
> 
> Any help on this issue would be highly appreciated.
> 
> 
> My sessionInfo is as follows:
> 
>> sessionInfo()
> R version 2.6.0 (2007-10-03)
> powerpc-apple-darwin8.10.1
> 
> locale:
> en_GB.UTF-8/en_GB.UTF-8/en_GB.UTF-8/C/en_GB.UTF-8/en_GB.UTF-8
> 
> attached base packages: [1] splines tools stats graphics grDevices utils
> datasets methods base
> 
> other attached packages: [1] simpleaffy_2.14.02 gcrma_2.10.0
> matchprobes_1.10.0 genefilter_1.16.0 survival_2.32 affy_1.16.0 [7]
> preprocessCore_1.0.0 affyio_1.6.0 Biobase_1.16.0
> 
> loaded via a namespace (and not attached): [1] AnnotationDbi_1.0.4
> DBI_0.2-3 RSQLite_0.6-3 annotate_1.16.0
>>
 
--------------------------------------------------------


This email is confidential and intended solely for the u...{{dropped:13}}



More information about the Bioconductor mailing list