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]]
_______________________________________________ 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]]