[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