Factorial kriging in GStat - tips for implementation?

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 at 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,

Tim

----------------------
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 :
http://www.ie.unimelb.edu.au/people/staff.php?person_ID=141135
Research Gate : https://www.researchgate.net/profile/Tim_Peterson7
Google Scholar:
http://scholar.google.com.au/citations?user=kkYJLF4AAAAJ&hl=en&oi=ao