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

[Greg Snow]

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]]
> >>
> >> _______________________________________________
> >> R-sig-Geo mailing list
> >> R-sig-Geo at r-project.org
> >> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >>
> >
> >
> >
> >
> >
>
>
> --
>
> 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
>



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

	[[alternative HTML version deleted]]