[R-sig-Geo] Factorial kriging in GStat - tips for implementation?

Tim Peterson timjp at unimelb.edu.au
Thu Aug 27 02:48:51 CEST 2015

Hi all,

I'm looking to undertake factorial colocated cokriging with GStat but 
factorial kriging does not appear to be an available option in GStat (eg 
help.search("factorial kriging") ). There was a post in 2007 requesting 
this feature (http://marc.info/?l=r-sig-geo&m=119810046000747) and a 
proposed approach in 2000 
(https://www.mail-archive.com/gstat-info@geog.uu.nl/msg00045.html) but 
no updates since. As far as I can tell, the required changes would be to 
'gls.c' lines 412-423, specifically the call to GCV0 for calculating C0 
(right-hand side of the kriging matrix) should be edited so that only 
the covariance for the required factor is calculated.

Anyway, before I consider adding this feature to GStat, I'd be grateful 
if anyone could share tips or thoughts. Lastly, any changes I make will 
(of course) be freely shared.

Many thanks,


Dr. Tim Peterson

The Department of Infrastructure Engineering
The University of Melbourne, 3010 Australia
T: +61 3 8344 9950, M: +61 0438 385 937

Dept. profile : 
Research Gate : https://www.researchgate.net/profile/Tim_Peterson7
Google Scholar: 

More information about the R-sig-Geo mailing list