Add Loess Fit to Hexplot
Fit a loess line using the hexagon centers of mass as the x and y coordinates and the cell counts as weights.
hexMA.loess(pMA, span = 0.4, col = "red", n = 200, ...) hexVP.loess(hbin, hvp = NULL, span = 0.4, col = "red", n = 200, ...)
hbin |
an object of class |
hvp |
A |
pMA |
the list returned by |
span |
the parameter alpha which controls the degree of smoothing. |
col |
line color for the loess fit. |
n |
number of points at which the fit should be evaluated. |
... |
Additional graphical parameter settings for the |
Returns invisibly the object associated with the loess fit.
Nicholas Lewin-Koh
if(require(marray)){
data(swirl)
hb <- plotMAhex(swirl[,1], main = "M vs A plot with hexagons", legend=0)
hexVP.abline(hb$plot, h=0, col= gray(.6))
hexMA.loess(hb)
}
dat <- data.frame(x=rnorm(1000), y=rnorm(1000))
bin <- hexbin(dat$x, dat$y)
hb <- plot(bin)
hexVP.loess(bin, hvp = hb$plot.vp, span = 0.4, n = 200, col = "blue", lwd = 3, lty = "dashed")Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.