[BioC] qqplot

Hrishikesh Deshmukh d_hrishikesh at yahoo.com
Sat Mar 19 17:02:17 CET 2005


Hi,

I want to make a qqplot for each CEL file in an
affybatch.
library(affy)
library(hgu95av2cdf)
library(hgu95av2probe)
library(matchprobes)
data(hgu95av2probe)
summary(hgu95av2probe)

getwd()
setwd("C:/Program Files/R/rw2001/Normal")
Data <- ReadAffy()
length(Data)
file.names<-c("01.CEL",  "02.CEL",  "03.CEL", 
"04.CEL",  "05.CEL",
"06.CEL","07.CEL",  "08.CEL",  "09.CEL",  "10.CEL", 
"11.CEL", "12.CEL", "13.CEL",
"14.CEL","15.CEL","16.CEL","17.CEL")
M<-read.affybatch(filenames=file.names,
                 description=NULL,notes="",compress=F,
                
rm.mask=F,rm.outliers=F,rm.extra=F,verbose=T)

But i get this error:
> qqplot(M[,1])
Error in sort(x) : `x' must be atomic

Thanks,
Hrishi



More information about the Bioconductor mailing list