Re: [R-sig-Geo] map2SpatialLines and map('world2')

2011-07-03 Thread Roger Bivand
On Sat, 2 Jul 2011, Rob Hodges wrote: Greetings R-sig-geo'ders; I want to use the Pacific-centric world map 'world2' from the maps package for use with hurricane observations. I am also making use of the sp spatial data classes. The following (Atlantic-centric) example works: x = map('world',

Re: [R-sig-Geo] testing linear combinations of coefficients after spatial regression

2011-07-03 Thread Roger Bivand
On Sat, 2 Jul 2011, Jason Sorens wrote: I'm interested in testing linear combinations of coefficients after models like spautolm, errorsarlm, and lagsarlm in the spdep package. I'm used to using "anova" for this, but the anova.sarlm implementation doesn't seem to permit it straightforwardly. For

Re: [R-sig-Geo] testing linear combinations of coefficients after spatial regression

2011-07-03 Thread Jason Sorens
Thanks, Roger. And thanks for alerting me to the possible issues with the coefficient covariance matrix. For my immediate purposes, this yields plausible results. On Sun, Jul 3, 2011 at 10:54 AM, Roger Bivand wrote: > On Sat, 2 Jul 2011, Jason Sorens wrote: > > I'm interested in testing linear

[R-sig-Geo] Assessing and ranking the relationships and contribution of environmental correlates to species richness

2011-07-03 Thread ah3881
I am trying to run some analysis to determine the percentage contribution of various factors (minimum temp in lgm, npp, etc) on determining species richness of 171 species throughout Southeast Asia (on a km by km basis-so over 4 million rows of cells, and about 14 columns). I have read about variou

Re: [R-sig-Geo] Assessing and ranking the relationships and contribution of environmental correlates to species richness

2011-07-03 Thread Mathieu Rajerison
Hi, How is represented your data: a regular grid: raster, irregular samples: points, or slots/parcels: areas? If you want to determine species richness, you can use simpson and shannon indexes with diversity package. On a point pattern object, use spatstat::marktable function as mentionned in h