question regarding qc function in simpleaffy package
Entering edit mode
Last seen 9.6 years ago
Hello James, I am surprised to see that - apparently - simpleaffy detects that another chip is provided than returned via the 'annotation(' function. Such a message would worry me indeed. Having said this: why would you want to perform QC utilizing the MBNI Custom CDF? The Custom CDF won't 'magically' improve bad quality arrays! In other words: performing QC utilizing the (default) Affymetrix CDF is just fine to assess thequality of your arrays. Regards, Dr. Philip de Groot Ph.D. Bioinformatics Researcher Wageningen University / TIFN Nutrigenomics Consortium Nutrition, Metabolism & Genomics Group Division of Human Nutrition PO Box 8129, 6700 EV Wageningen Visiting Address: Erfelijkheidsleer: De Valk, Building 304 Dreijenweg 2, 6703 HA Wageningen Room: 0052a T: +31-317-485786 F: +31-317-483342 E-mail: Philip.deGroot at <mailto:philip.degroot at=""""> Internet: <http:""/> <http:""/> <https:""/> ________________________________ From: James Anderson [] Sent: Wed 31-3-2010 15:31 To: bioconductor at Subject: [BioC] question regarding qc function in simpleaffy package Hi, I am trying to use the qc function in the simpleaffy package to generate those qc metrics for a bunch of cel files in U133A, but I need to use a different cdf file, not the Affy CDF. Below is the code library(simpleaffy) cel.path <- "C:\Cel" fn <- list.celfiles(cel.path, = T) <- ReadAffy(filenames = fn,cdfname = "hgu133ahsentrezg") x.mas5 <- call.exprs(,"mas5",sc = 100) qcmas5 <- qc(,x.mas5) Everything runs smoothly except the last line, it has the following warning message: Warning message: In qc.affy(unnormalised, ...) : CDF Environment name ' hgu133acdf ' does not match cdfname ' hgu133ahsentrezgcdf ' Is this problem negligible? or it is serious? I've looked into the qc function and it does not seem to have an argument that allows you to specify cdfname. Thanks, -James
cdf affy simpleaffy cdf affy simpleaffy • 935 views

Login before adding your answer.

Traffic: 372 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6