Question on histograms from Raster and RasterVIS packages

Nuno,

Does rasterVIS use base or grid graphics?

The place I usually see the error that plot.new has not been called yet is
when trying to use a base graphics function to augment a plot created using
grid graphics.  I suspect that the plot created used grid graphics and that
is why abline is not working for you.

There is a 'grid.abline' function in the grid package, though it does not
look like it would work easily for a vertical line, so you may need to dive
into the details of using grid in more detail, probably using the
grid.lines function.  There is also the gridBase package which helps grid
and base graphics play nicely together, though it could be a bit of
overkill for just adding a vertical line.

On Tue, Nov 18, 2014 at 11:21 AM, Nuno S? <nunocesardesa at gmail.com> wrote:

Hey Thomas!

Yes thanks, I know it works for the normal hist function which uses 3% of
all data (mind I am talking about rasters) which I am not sure that is a
good enough sample, I leave that question in the air

The problem is it does not work in rasterVIS (which does use all the
raster). The error it gives me is the following:

abline(v=m.rst,col="blue",lwd=2)

Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) :
  plot.new has not been called yet

Regards!!
Thank you once again!
Nuno

On 18 November 2014 18:04, Thomas Adams <tea3rd at gmail.com> wrote:

Nuno,

This works for me:

require(stats)
set.seed(14)
x <- rchisq(100, df = 4)


hist(x)
meanx<-mean(x)
meanx

[1] 3.734162

abline(v=meanx, col="red")

Cheers!
Tom

On Tue, Nov 18, 2014 at 10: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]]

--

Nuno C?sar de S?
+351 91 961 90 37

        [[alternative HTML version deleted]]

--
Gregory (Greg) L. Snow Ph.D.
538280 at gmail.com