[R] mva.pairs

From: Vijay A Raghavan <rva_ec_at_yahoo.com>
Date: Tue 21 Feb 2006 - 04:38:06 EST


Hello,  

 I am using the following code to plot an MVA plot.  

 library(affy)
 library(Biobase)
 library(limma)
 library(gcrma)
 
     pd<-read.phenoData("Clk.targets.2.txt",header=TRUE,
     row.names=1,as.is=TRUE,sep="\t")
     Data <- ReadAffy(filenames=pData(pd)$FileName,phenoData=pd) 
     Print(Data)
     
     eset <- gcrma(Data) 
     write.exprs(eset, file="clk.6-23-05.txt")
       
         
     bitmap("clk-1.mva.jpg",width=15, height=12)
     mva.pairs(exprs(eset)[,c(19,20,21,22,23,24)],log.it=FALSE)
     dev.off()
 

 but I am getting an subscript out of bound error.  

 Any Idea ?  

 Thanks,  

 Vijay  

        [[alternative HTML version deleted]]



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Tue Feb 21 04:43:27 2006

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:42:38 EST