OCK in gstat and co-located [SEC=UNCLASSIFIED]

Hi Jin, you're not the first to run into this one; read the fit.lmc 
documentation.

fit.lmc fits partial sill matrices that may have perfect correlation 
among variables. To avoid so, you could instead of

...

mud.ock.fit=fit.lmc(x, mud.ock)

try:
mud.ock.fit = fit.lmc(x, mud.ock, correct.diagonal = 1.01)

or a larger value. This gets the correlation down to say .99 instead of 
1.0; 1.0 is allowed (I believe) for a LMC, but then you run into trouble 
when the secondary variable is collocated with the first.

There's also ways of avoiding predict.gstat to do this check, but that 
might lead to trouble further down the pipeline.
--
Edzer

mud.pred<-predict(mud.ock.fit, newdata=data.file.pred)

#plot(x, model=mud.ock.fit)

mud.pred<- as.data.frame(mud.pred)

mud.pred$bt.pred<-(mud.pred$md.pred)^2

mud.pred$bt.pred[mud.pred$bt.pred>=100]=100

 

Regards,

 

Jin

 

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.springer.com/978-0-387-78170-9 e.pebesma at wwu.de