[BioC] simpleaffy and yeast2cdf problems
Jenny Drnevich
drnevich at uiuc.edu
Mon Jan 8 16:42:07 CET 2007
Hi Dan,
See the very recent post by Greg Alford:
https://stat.ethz.ch/pipermail/bioconductor/2007-January/015514.html
They have modified the code of qc() so that it will work for unknown chip
types. Although he is calling his code soybean.mod.qc.affy(), there is
nothing specific for the soybean array so it can be used for any chip type.
However, it doesn't give all the information that qc() does; the spike-in
probes, other qc probes and bioB hyb controls are specific for each chip
and had to be left out. If you know these for your chip type, you might be
able to modify the original code (commented out in soybean.mod.qc.affy) to
include them.
Cheers,
Jenny
At 05:49 AM 1/8/2007, Dan Swan wrote:
>Hi Bioconductorites,
>
>I've seen references to this problem in 2005
>(https://stat.ethz.ch/pipermail/bioconductor/2005-November/010865.html)
>and 2006 (https://stat.ethz.ch/pipermail/bioconductor/2006-March/012480.html),
>the former suggesting it should be fixed, but I'm still having the
>same issue.
>
>I just want some quick and dirty QC measures on some chips I am pretty
>sure are duff.
>
>following along with the outline in the "QC and Affymetrix data" pdf:
>
> > ampli.data<-read.affy("phenodata")
> > ampli.eset<-call.exprs(ampli.data,"mas5")
>Background correcting
>Retrieving data from AffyBatch...done.
>Computing expression calls...
>..............................done.
>scaling to a TGT of 100 ...
>Scale factor for: Yeast0wt1.CEL 0.1877573661054
>Scale factor for: Yeast0wt2.CEL 0.169107349891647
>Scale factor for: Yeast0wt3.CEL 0.176890842235261
>Scale factor for: Yeast1wt1.CEL 0.114515875199183
>Scale factor for: Yeast1wt2.CEL 0.132867931840117
>Scale factor for: Yeast1wt3.CEL 0.201211860669323
>Scale factor for: Yeast2wt1.CEL 0.120553407334518
>Scale factor for: Yeast2wt2.CEL 0.110981821370179
>Scale factor for: Yeast2wt3.CEL 0.12635084302873
>Scale factor for: Yeast3wt1.CEL 0.133156738139789
>Scale factor for: Yeast3wt2.CEL 0.119446187371786
>Scale factor for: Yeast3wt3.CEL 0.136591273338234
>Scale factor for: Yeast4wt1.CEL 0.143174910484902
>Scale factor for: Yeast4wt2.CEL 0.127534874666620
>Scale factor for: Yeast4wt3.CEL 0.177048335355844
>Scale factor for: Yeast0mut1.CEL 0.126568563373901
>Scale factor for: Yeast0mut2.CEL 0.177828594244113
>Scale factor for: Yeast0mut3.CEL 0.113037148372224
>Scale factor for: Yeast1mut1.CEL 0.125912932405511
>Scale factor for: Yeast1mut2.CEL 0.167309068269027
>Scale factor for: Yeast1mut3.CEL 0.134897445983855
>Scale factor for: Yeast2mut1.CEL 0.136284006808245
>Scale factor for: Yeast2mut2.CEL 0.197169068283275
>Scale factor for: Yeast2mut3.CEL 0.149844757694802
>Scale factor for: Yeast3mut1.CEL 0.119095781239726
>Scale factor for: Yeast3mut2.CEL 0.155941398564195
>Scale factor for: Yeast3mut3.CEL 0.144702217239922
>Scale factor for: Yeast4mut1.CEL 0.183634382620250
>Scale factor for: Yeast4mut2.CEL 0.180182357563251
>Scale factor for: Yeast4mut3.CEL 0.152303319246271
> > qc.data<-qc(ampli.data, ampli.eset)
>Error in qc.affy(unnormalised, ...) : I'm sorry, I do not know about
>chip type: yeast2cdf
>
>which is the same error that has been reported before.
>
>sessionInfo() follows:
>
> > sessionInfo()
>R version 2.4.1 (2006-12-18)
>i386-pc-mingw32
>
>locale:
>LC_COLLATE=English_United Kingdom.1252;LC_CTYPE=English_United
>Kingdom.1252;LC_MONETARY=English_United
>Kingdom.1252;LC_NUMERIC=C;LC_TIME=English_United Kingdom.1252
>
>attached base packages:
>[1] "splines" "tools" "stats" "graphics" "grDevices"
>"utils" "datasets"
>[8] "methods" "base"
>
>other attached packages:
> yeast2cdf simpleaffy genefilter survival affy affyio Biobase
> "1.14.0" "2.8.0" "1.12.0" "2.30" "1.12.2" "1.2.0" "1.12.2"
>
>any help would be appreciated.
>
>regards,
>
>Dan
>
>--
>Research Support Officer, Bioinformatics Support Unit,
>Institute for Cell and Molecular Biosciences,
>Faculty of Medical Sciences, Framlington Place,
>University of Newcastle upon Tyne, Newcastle, NE2 4HH
>Tel: +44 (0)191 222 7253 (Leech: Rooms M.2046/M.2046A - Mon/Wed)
>Tel: +44 (0)191 246 4833 (Devonshire: Rooms G.25/G.26 - Thu/Fri)
>Website: http://bsu.ncl.ac.uk/
>
>_______________________________________________
>Bioconductor mailing list
>Bioconductor at stat.math.ethz.ch
>https://stat.ethz.ch/mailman/listinfo/bioconductor
>Search the archives:
>http://news.gmane.org/gmane.science.biology.informatics.conductor
Jenny Drnevich, Ph.D.
Functional Genomics Bioinformatics Specialist
W.M. Keck Center for Comparative and Functional Genomics
Roy J. Carver Biotechnology Center
University of Illinois, Urbana-Champaign
330 ERML
1201 W. Gregory Dr.
Urbana, IL 61801
USA
ph: 217-244-7355
fax: 217-265-5066
e-mail: drnevich at uiuc.edu
More information about the Bioconductor
mailing list