[R-sig-Geo] Masking interpolations

2009-02-15 Thread Wesley Roberts
Dear R-sig-geo'ers I am currently running some interpolations using automap written by Paul Hiemstra. So far my interpolations have been producing suitable results except for one problem. From the code you will see that the boundaries of the spatial grid are determined using the range of the X

Re: [R-sig-Geo] Point pattern analysis

2009-02-15 Thread Virgilio Gomez Rubio
Dear Michel, > I'm new to Spatial Data Analysis and have just begun working through > "Applied Spatial Data Analysis wit R" by Bivand et al. For my research I > would like to use SDA to be able to tell more about my restaurant data set > than just pinpointing them on a google map. So far, from rea

Re: [R-sig-Geo] Vector to raster?

2009-02-15 Thread Jonathan Greenberg
Roger: Thanks, I'll try to incorporate some of your suggestions and post the code + some sample imagery next time I ask a question. I'm trying to build up some R-based code that replicates some of what my "starspan" program does (starspan.casil.ucdavis.edu) -- essentially taking vector da

Re: [R-sig-Geo] SpatialFiltering in spdep: How to create an 'nb' object from an imported spatial weights matrix

2009-02-15 Thread Roger Bivand
On Sun, 15 Feb 2009, Timo Mitze wrote: Dear Roger and List, thanks for your advice. I tried it with the code as (as well as some related with different style parameter etc.): x <- mat2listw(mat) lw <- nb2listw(x$neighbours, style="B", zero.policy=TRUE) migfilt

Re: [R-sig-Geo] SpatialFiltering in spdep: How to create an 'nb' object from an imported spatial weights matrix

2009-02-15 Thread Timo Mitze
Dear Roger and List, thanks for your advice. I tried it with the code as (as well as some related with different style parameter etc.): x <- mat2listw(mat) lw <- nb2listw(x$neighbours, style="B", zero.policy=TRUE) migfilt <- SpatialFiltering(lnmr_i ~ 1, data=migr

Re: [R-sig-Geo] localmoran p-values with/without permutation

2009-02-15 Thread Roger Bivand
On Sun, 15 Feb 2009, Valerio Bartolino wrote: Dear list, I've the objective to identify hotspot areas from a model prediction over a high resolution grid. After calculating a spatial weight object I easily applied the localmoran function from the spdep library. It's not really clear to me the me

Re: [R-sig-Geo] localmoran p-values with/without permutation

2009-02-15 Thread Danlin Yu
Valerio: For a local Moran's I, a randomization approximation is probably not a good idea since for each permutation, the mean and variance (which are used to calculated the p-values) of a specific location will change, which causes the randomized p-values not poolable (in global Moran's case

[R-sig-Geo] Point pattern analysis

2009-02-15 Thread Michel Barbosa
Dear all, I'm new to Spatial Data Analysis and have just begun working through "Applied Spatial Data Analysis wit R" by Bivand et al. For my research I would like to use SDA to be able to tell more about my restaurant data set than just pinpointing them on a google map. So far, from reading the li

[R-sig-Geo] Point pattern analysis

2009-02-15 Thread Michel Barbosa
Dear all, I'm new to Spatial Data Analysis and have just begun working through "Applied Spatial Data Analysis wit R" by Bivand et al. For my research I would like to use SDA to be able to tell more about my restaurant data set than just pinpointing them on a google map. So far, from reading the li

Re: [R-sig-Geo] Vector to raster?

2009-02-15 Thread Roger Bivand
On Sat, 14 Feb 2009, Jonathan Greenberg wrote: Thanks -- another few questions along these lines: Perhaps you could simplify by providing code examples with a built-in data set and so on? Verbatim copies do let others see what is going on. So I'll start by trying to reconstruct what you seem

Re: [R-sig-Geo] SpatialFiltering in spdep: How to create an 'nb' object from an imported spatial weights matrix

2009-02-15 Thread Roger Bivand
On Sat, 14 Feb 2009, Timo Mitze wrote: Dear List, I have a complete (distance/border based binary) spatial weights matrix from Stata (e.g. as a txt-file). Try using read.dta() in foreign of the exported Stata matrix. This reads a data.frame, which you coerce to an R matrix with as.matrix().

[R-sig-Geo] localmoran p-values with/without permutation

2009-02-15 Thread Valerio Bartolino
Dear list, I've the objective to identify hotspot areas from a model prediction over a high resolution grid. After calculating a spatial weight object I easily applied the localmoran function from the spdep library. It's not really clear to me the meaning of the p-values associated to the localmora