[BioC] superimpose lowess on smoothScatter
Anderson, S. Keith
Anderson.S at mayo.edu
Mon Aug 30 17:13:42 CEST 2010
Hi Andrew,
Another option is to use the package "hexbin"
----------
library(hexbin)
## a vector of your M-values
m.ii
## a vector of your A-values
a.ii
temp.plot <- list(M=m.ii,A=a.ii)
hb <- plotMAhex(temp.plot,newpage=TRUE,main="Main Title Goes
Here",legend=0)
hexVP.abline(hb$plot.vp, h=0, col="black")
hexMA.loess(hb,span=.75)
# if a Median Line is desired
hexVP.abline(hb$plot.vp,v=median(means.log2), col="blue")
---------
Have a great day,
S. Keith Anderson, M.S.
Statistician/Assistant Supervisor
Cancer Center Statistics
Division of Biomedical Statistics and Informatics
Phone: 507-266-9059
Fax: 507-266-2477
E-mail: anderson.s at mayo.edu
_______________________________
Mayo Clinic
200 First Street S.W.
Rochester, MN 55905
www.mayoclinic.org
-----Original Message-----
From: bioconductor-bounces at stat.math.ethz.ch
[mailto:bioconductor-bounces at stat.math.ethz.ch] On Behalf Of Andrew
Einhorn
Sent: Monday, August 30, 2010 10:03 AM
To: bioconductor at stat.math.ethz.ch
Subject: [BioC] superimpose lowess on smoothScatter
I am trying to create an MvA plot using smoothScatter that has a lowess
curve fit superimposed. I can't seem to find documentation on
smoothScatter
that will allow this. Currently I have:
smoothScatter(x,y, main = "My MvA Plot")
Can a lowess curve fit be added as an argument to smoothScatter? If not,
is
there a way to add it externally to the smoothScatter function?
[[alternative HTML version deleted]]
_______________________________________________
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
More information about the Bioconductor
mailing list