[BioC] Questions about MAplot in the oligo package
Benilton Carvalho
beniltoncarvalho at gmail.com
Tue Mar 5 02:26:34 CET 2013
Hi Yiwen,
1. The IQR values refer to the log-ratio (M)
2. In general, there are no requirements for the loess line to be as
close to zero as possible. Note that your call to MAplot() produces an
MA plot comparing each of your samples to a (virtual) average array.
Therefore, if all your samples belong to the same group, then you can
expect the loess line to be very close to zero. For samples that
differ "significantly" from the average array, the line will not be
very close to zero.
3. Even after RMA, I wouldn't expect median+IQR to be zero. How many
samples are you running? Just to confirm (this may sound silly),
you're sure that all of your CEL files correspond to different
hybridizations, right?
b
2013/3/4 He, Yiwen (NIH/CIT) [C] <heyiwen at mail.nih.gov>:
> Hi,
>
> I am using the oligo package for ST gene arrays, and I have some questions about the MAplot() function in the oligo package.
> This is part of my sessionInfo():
>
>> sessionInfo()
> R version 2.15.2 (2012-10-26)
> Platform: i386-w64-mingw32/i386 (32-bit) other attached packages:
> [1] pd.hugene.1.0.st.v1_3.8.0 RSQLite_0.11.2
> [3] DBI_0.2-5 oligo_1.22.0
> [5] Biobase_2.18.0 oligoClasses_1.20.0
> [7] BiocGenerics_0.4.0
>
> And this is the MAplots that I tried to generate:
>
>> myData <- read.celfiles(celFiles)
>> MAplot(myData, ylim=c(-2, 2))
>
>> eset <- rma(myData, target="core")
>> MAplot(myData, ylim=c(-2, 2))
>
> So I have two sets for MA plots one for the raw CEL file values and one for the RMA normalized values. My questions are:
> 1. About the median and IQR values reported on each plots, are they for the M values?
> 2. In the MA plots for the RMA values, should the loess line be very close to the horizontal 0 line? What does it mean if it's not?
> 3. My MAplots for the RMA values all have median and IQR values of 0. I understand that the median should be 0, but how about IQR? Should that all be 0 as well? If one array is not, is that an outlier?
>
> Thank you very much for your help!
>
> Yiwen
> CIT/NIH/HHS
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list