Question on histograms from Raster and RasterVIS packages
Here is an example, with the mean in dotted red line and the median in
dotted blue line:
library(rasterVis)
f <- system.file("external/test.grd", package="raster")
r <- raster(f)
mn <- cellStats(r, mean)
md <- cellStats(r, median)
histogram(r) + layer(panel.abline(v = mn, col = "red", lty = 2)) +
layer(panel.abline(v = md, col = "blue", lty = 2))
On Wed, Nov 19, 2014 at 2:56 AM, Nuno S? <nunocesardesa at gmail.com> wrote:
Hello!
My aim is to add a "Median" or a "Mean" line to an histogram plot in R.
The problem is the following:
The "hist" function from the raster package uses a maximum of 100 000
values for generating the histogram but allows me to edit the plot, so I
can easily add a line within the plot using abline.
The "histogram" function in rasterVIS uses all my dataset but does not
allow editing of the plot area, so I cannot use the abline function to add
this line (or I do not know how to do it).
If you have an alternative or a solution to work around this, I am all hears
Thank you for any help in advance!
Ciao!
--
Nuno C?sar de S?
+351 91 961 90 37
[[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Pascal Oettli Project Scientist JAMSTEC Yokohama, Japan