[R-sig-Geo] How to use R scripts to estimate two-weight matrix SAR models

2015-10-17 Thread Chen Lee
Seek expert help: How to use R scripts to estimate two-weight matrix SAR models? The two-weight matrix SAR model is specified as: Y=ρ1W1Y+ρ2W2Y+Xβ+ε. [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@r-project.org

Re: [R-sig-Geo] Convert rasters to data frame with time stamp

2015-10-17 Thread Thiago V. dos Santos
Thank you Dominik, Vijay and Robers for all the valuable inputs. Robert, the function you devised is really convenient. However, it fails when the observed and modelled raster objects have different number of layers. Please notice that the obs and model objects on the example from the qmap

Re: [R-sig-Geo] Convert rasters to data frame with time stamp

2015-10-17 Thread Robert J. Hijmans
Thiago, > Therefore, I need to load data as rasters and iterate through all individual > gridcells to create a data frame containing: > date1, cell1, cell2, cell3, ..., cell3600 > date2, cell1, cell2, cell3, ..., cell3600 > date3, cell1, cell2, cell3, ..., cell3600... > date408, cell1, cell2,

[R-sig-Geo] Predicted kriging variances don't match errors in cross-validation

2015-10-17 Thread Roelof Coster
Dear list, I'm working on a local (nmax=100) space-time kriging model. I did a cross-validation in which I made my model predict the values at 2000 randomly selected data points, based on the rest of the observations. The results are quite good (the average error is very small, the errors are