So, it seems you have one raster for your dependent variable: species
richness and one for each explanatory variable.
I have never accomplished this but I think you should make a raster stack of
all these variables using raster package: the dependent one + the
explanatory ones then fit a model,
Solved: it was caused by having upgraded to raster 1.8.40 but using
raster objects created under 1.8.38
Agus
2011/7/6 Agustin Lobo alobolis...@gmail.com:
A function of mine outputs several raster objects in a list l:
class(l$lisaobsv)
[1] RasterLayer
attr(,package)
[1] raster
But by some
Dear Robert,
Thank you for sorting this out. And thanks to all the others for their
valuable comments.
Christian
--
View this message in context:
http://r-sig-geo.2731867.n2.nabble.com/Subsetting-RasterBrick-very-slow-tp6537890p6554009.html
Sent from the R-sig-geo mailing list archive at
I summarize here the different solutions I got to my question regarding
how to plot several raster layers on the same page and with the same
color scale.
Sorry it took a long time, I was interrupted by something else for weeks.
Hope this is useful.
Many thanks to all.
require(raster)
Hi,
Just to add one more solution. It uses the rasterVis package (at
R-Forge now and should be at CRAN today).
library(raster)
library(rasterVis)
levelplot(s, par.settings=RdBuTheme)
More examples here: http://rastervis.r-forge.r-project.org/
Best,
Oscar.
---
Oscar
Horacio
See this thread from December:
http://www.mail-archive.com/r-sig-geo@r-project.org/msg00154.html
I haven't made much progress or further tested this approach since then. It
seems to work for me and my needs.
The example with the North Carolina map relies on pretty() to determine where
I'm using sp::overlay between a SpPolDF and a SpPointsDF:
delme - overlay(rndp,q1km, fn = mean)
class(q1km)
[1] SpatialPolygonsDataFrame
attr(,package)
[1] sp
class(rndp)
[1] SpatialPointsDataFrame
attr(,package)
[1] sp
class(delme)
[1] data.frame
but I'm confused by the fact that the
I'm not too familiar with GIS terminology, but would you agree that
sp::over does (or allows for) a spatial table join?
On 07/06/2011 06:10 PM, Agustin Lobo wrote:
Yes!
I'm doing:
delme2 - over(q1km, rndp,fn = mean)
And now not only the row.names are correct and the order of the 2
