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

[Edzer Pebesma]

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

Jin.Li at ga.gov.au wrote:
...
>
> 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