Question: affyPLM boxplots for subset of samples
0
gravatar for Juliet Hannah
3.1 years ago by
Juliet Hannah360
United States
Juliet Hannah360 wrote:

I have hundreds of arrays for which I would like to obtain NUSE and RLE plots. I would like to make one call to affyPLM, but then plot a subset (boxplot), perhaps a few random samples that cover the entire set. Otherwise, there will be too many samples to show effectively in one boxplot. What is a way to go about this?  Example below taken from affyPLM vignette.

require(affydata)
library("affyPLM")
data(Dilution)
Dilution = updateObject(Dilution)
Pset <- fitPLM(Dilution)
boxplot(Pset)

 

 

affyplm • 432 views
ADD COMMENTlink modified 3.1 years ago by James W. MacDonald50k • written 3.1 years ago by Juliet Hannah360
Answer: affyPLM boxplots for subset of samples
1
gravatar for James W. MacDonald
3.1 years ago by
United States
James W. MacDonald50k wrote:

Unfortunately the PLMset class never had a [ function, so it's not as simple as

boxplot(Pset[,1:3])

In addition, if you want NUSE or RLE values, you should use either NUSE or RLE rather than boxplot, because you can return the values.

Let's say you have 200 samples, and you want to do a NUSE plot on a random subset of 25:

nuse.dat <- NUSE(Pset, type = "values")

boxplot(nuse.dat[,sample(1:200, 25)])

You can do the same for RLE.

ADD COMMENTlink written 3.1 years ago by James W. MacDonald50k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 121 users visited in the last hour