maxdist for kriging with an external drift
Dear Edzer, Would it be correct to say then that if a neighbourhood is specified in the krige command the result would be that of Kriging with an External Drift (KED), otherwise it would be that of Simple Kriging with varying local means (SKlm)? Apologies for posting on this thread but I was about to post a similiar question. Thanks, Cutberto. 2010/1/26 Edzer Pebesma <edzer.pebesma at uni-muenster.de>:
Yes, that is right. Els Verfaillie wrote:
Dear list, I want to use Kriging with an external drift for a sedimentological dataset of grain-size that has a linear relation with the depth. Am I correct that when I set a 'maxdist' using the krige command, that a trend for the primary variable (grain-size) is calculated as a local linear function of the secondary variable (depth)? Is this function thus different for each interpolation window? d50.ked.dir50 <- krige(D50F~depth, locations=ds50, newdata=Depth, model=d50.fit.var.50, nmin=2, nmax=16, maxdist=9000) Thank you for your help. Best regards, Els Verfaillie
______________________________________________ Dr. Els Verfaillie Carto-GIS cluster Ghent University (UGent) - Department of Geography ______________________________________________ ? ? ? ?[[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
-- Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster Weseler Stra?e 253, 48151 M?nster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 ?http://ifgi.uni-muenster.de http://www.52north.org/geostatistics ?e.pebesma at wwu.de
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo