[BioC] boxplot of subset of PLMset

Rebecca tagett tagett at ipsogen.com
Fri Aug 22 17:00:28 CEST 2008


Hello!
I am trying to do a boxplot using a subset of a PLMset (just certain
boxes). Without trying to subset (3 lines of code below), it works just
fine. (It is the same thing as NUSE)

	# THIS IS OK
	fplm <- fitPLM(myAffyBatch)
	ss <- boxplot(fplm)  
	bxp(ss)

How can I get a subset of the PLMset ? Here are 2 things that I tried,
using the subset argument in boxplot, which did not work (where
cels4Plot is a vector containing the CEL file names):

ss <- boxplot(fplm, 
subset = sampleNames(fplm at phenoData)== cels4Plot)

ss <- boxplot(fplm,
subset = which(sampleNames(fplm at phenoData) %in% cels4Plot))

Thanks for your help!



More information about the Bioconductor mailing list