[R-sig-Geo] Question on histograms from Raster and RasterVIS packages

[Pascal Oettli]

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