[R-sig-Geo] GWR Analysis

2010-06-11 Thread Pinar Aslantas Bostan
Dear all, I want to make GWR analysis to predict precipitation distribution measured from 225 meteorological stations. I have three independent variables (Z, V1, V2). I want to obtain predictions on the grid dataset which has 31203 number of pixels. I gave some details about my datasets

Re: [R-sig-Geo] How to extract coordinates values from a shapefile?

2010-06-11 Thread Paul Hiemstra
On 06/09/2010 11:01 PM, Dylan Beaudette wrote: On Wednesday 09 June 2010, Hadley Wickham wrote: On Wed, Jun 9, 2010 at 3:24 PM, Edzer Pebesma edzer.pebe...@uni-muenster.de wrote: The example provided by Matt assumes that each polygon consists of a single ring, and doesn't have

Re: [R-sig-Geo] GWR Analysis

2010-06-11 Thread Roger Bivand
On Fri, 11 Jun 2010, Pinar Aslantas Bostan wrote: Dear all, I want to make GWR analysis to predict precipitation distribution measured from 225 meteorological stations. I have three independent variables (Z, V1, V2). I want to obtain predictions on the grid dataset which has 31203 number

Re: [R-sig-Geo] How to locate local minima points on a Krige-fitted surface?

2010-06-11 Thread Paul Hiemstra
Dear Jeff, I think you cannot escape discretizing your kriging surface as there is no mathematical expression that captures the kriging surface. To use optim you could make an objective function that was two inputs, x and y location. The functions calls the fields package to estimate the

Re: [R-sig-Geo] test for CSR

2010-06-11 Thread Paul Hiemstra
On 06/08/2010 12:37 PM, adrian.badde...@csiro.au wrote: Hamidhamid200...@yahoo.com writes: I use the following function to simulate CSR point pattern nsim times. This is a question about the package 'spatstat'. Is there a way to reduce running time (maybe by avoiding the

Re: [R-sig-Geo] GWR Analysis

2010-06-11 Thread Pinar Aslantas Bostan
Dear Roger, I tried with data frames instead of SPDF. I gave some details of names and class types of datasets. Grid data has 31204 rows (first row is column names) and 5 columns (Z,X,Y,V1,V2). Station data has 226 rows (first row is column names) and 6 columns (PREC,Z,X,Y,V1,V2). I gave

[R-sig-Geo] variogram with nominal data

2010-06-11 Thread Julien Beguin
Dear list, I would like to evaluate spatial dependence for about 1 spatial data points. However, the variable associated with those points is not continuous as in the data meuse but is a land cover type (e.g. WATER). What should I consider as variable Z in the variogram() function with

Re: [R-sig-Geo] variogram with nominal data

2010-06-11 Thread ashton
Hi Julian, and list, the variogram calculates the average squared difference between z-values separated by particular lags. Calculating differences with nominal data, and certainly IDs, is wrong. You will get results, but garbage in, garbage out applies here. That said, you can calculate

Re: [R-sig-Geo] GWR Analysis

2010-06-11 Thread Roger Bivand
On Fri, 11 Jun 2010, Pinar Aslantas Bostan wrote: Dear Roger, I tried with data frames instead of SPDF. I gave some details of names and class types of datasets. Grid data has 31204 rows (first row is column names) and 5 columns (Z,X,Y,V1,V2). Station data has 226 rows (first row is