[BioC] boxplot of MM data and normalization
James W. MacDonald
jmacdon at med.umich.edu
Tue Oct 3 21:35:21 CEST 2006
Shi, Tao wrote:
> Hi list,
>
> My first question is a follow-up of an old post (
> https://stat.ethz.ch/pipermail/bioconductor/2003-October/002660.html
> ). I don't think boxplot(myDat, with="MM") works. Please see below.
> I'm using R 2.3.0, WinXP, affy 1.10.0.
You are correct; with = "MM" doesn't work because R is case sensitive,
and the correct argument is "mm".
>
> My second question is: do differenet normalization options in
> 'expresso' do anything to the MM probes?
Well, a quick look at ?normalize.loess gives this:
Details:
The type arguement should be one of
'"separate","pmonly","mmonly","together"' which indicates whether
to normalize only one probe type (PM,MM) or both together or
separately.
Looking at the help pages for the other methods may be just as helpful ;-P
Best,
Jim
>
> thanks,
>
> ...Tao
>
>
>
>> data(affybatch.example) par(mfrow=c(1,2))
>> boxplot(affybatch.example, with="PM") boxplot(affybatch.example,
>> with="MM")
>>
>> boxplot(affybatch.example, with="PM")$stat
>
> [,1] [,2] [,3] [1,] 6.175924 5.672425 5.942515 [2,]
> 6.794416 6.375039 6.574404 [3,] 7.179909 6.870365 6.954196 [4,]
> 8.321921 7.927185 8.053383 [5,] 14.176664 13.515946 14.188434
>
>> boxplot(affybatch.example, with="MM")$stat
>
> [,1] [,2] [,3] [1,] 6.175924 5.672425 5.942515 [2,]
> 6.794416 6.375039 6.574404 [3,] 7.179909 6.870365 6.954196 [4,]
> 8.321921 7.927185 8.053383 [5,] 14.176664 13.515946 14.188434
>
>>
>> boxplot(data.frame(log2(pm(affybatch.example))))$stat
>
> [,1] [,2] [,3] [1,] 6.175924 5.700440 5.942515 [2,] 6.768184
> 6.345183 6.549205 [3,] 7.108524 6.791814 6.882643 [4,] 8.012624
> 7.660353 7.749199 [5,] 9.878051 9.612868 9.541097
>
>> boxplot(data.frame(log2(mm(affybatch.example))))$stat
>
> [,1] [,2] [,3] [1,] 6.205549 5.672425 5.995485 [2,]
> 6.820179 6.397461 6.589464 [3,] 7.285402 6.965784 7.056908 [4,]
> 8.657851 8.224002 8.359309 [5,] 11.398904 10.959278 10.990104
>
>
> _______________________________________________ Bioconductor mailing
> list Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor Search the
> archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor
--
James W. MacDonald, M.S.
Biostatistician
Affymetrix and cDNA Microarray Core
University of Michigan Cancer Center
1500 E. Medical Center Drive
7410 CCGC
Ann Arbor MI 48109
734-647-5623
**********************************************************
Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues.
More information about the Bioconductor
mailing list