On Fri, 15 Jan 2010, rick reeves wrote:
Hello List:
I have been reviewing techniques within R for performing the GIS-like
operation 'polygon algebra' (computing the union and intersection of
polygons within two SpatialPolygonDataFrame objects). The goal is to
combine two polygon data frames,
Hello List:
I have been reviewing techniques within R for performing the GIS-like
operation
'polygon algebra' (computing the union and intersection of polygons within
two SpatialPolygonDataFrame objects). The goal is to combine two polygon
data frames, to produce a new polygon data frames that
Hello,
If I generate a variogram with an outside program on the x-y plane, would this
variogram be valid for ordinary kriging in gstat? I am giving gstat the x,y,z
coordinates and I am trying to predict at the elevations of a digital elevation
model.
Thank you in advance for taking the time to
mykrige.cv = gstat:::krige.cv.spatial
mykrige.cv = edit(mykrige.cv)
then, you will find out that you will have to do the same with gstat.cv,
and that gstat.cv is rather cryptic as it deals with cokriging of which
ordinary kriging is a special case.
Maybe it wasn't such a good idea after all i
Hello Prof. Pebesma,
I would like to use leave-one-out cross validation to have a fair comparison
with my other methods.
Could you explain further how to copy krige.cv to my own function and modify?
Can I open up the krige.cv code somewhere? I am still new to R.
Thank you in advance,
Cara
__
Dear Cara,
one way would be to copy krige.cv to your own function and modify this;
I would recommend this if you want leave-one-out cross validation.
If n-fold cross validation is good enough for a few folds, e.g. because
you have plenty of data, then I would do the whole procedure "manually"