[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

regression kriging and x-validation




Greetings mailing list members,

I have two queries about  gstat.

1. I want to perform regression kriging (i.e. krig the residuals of a regression model between two variables thereafter adding the kriged residuals back to the predicted map). Can this be done in 'universal kriging' mode in gsat, using the predicted surface from the regression model as the drift? I would assume if this was the case, the mask would be the predicted variable (Y) at all locations and X=4 in the command line, for example would, be a column containing the predicted (Y) variable at the data locations.

2. Is there any way I can look into the details of the cross-validation statistics? i.e. generate a bivariate plot of the errors, compute other evaluation statistics, etc.?

My very best wishes
Jonathan

---------------------------------------------
Jonathan Seaquist
Dept. of Physical Geography
Remote Sensing Lab.
Lund University
Sölvegatan 13
S--223 62 Lund
Sweden

Phone: +46 +46 222 42 62
Fax:   : +46 +46 222 40 11
email  : jonathan.seaquist@natgeo.lu.se