[R-sig-Geo] Cross variogram fitting

Saman Monfared samanmonfared1 at gmail.com
Wed Dec 28 21:38:11 CET 2011


Dear all
I want fitt a cross variogram with two covariables.

I have two problems.

1)-How can I find the best model for cross variogram??

2)-When I use fit.lmc() with fit.ranges = F or  fit.ranges = T it is
non possituve
 definite.

B is main variable and E,T are covariables.

coordinates(dd) <-~x+y
g = gstat(NULL, "T", T~1, dd)
g = gstat(g, "E", E~1, dd)
g = gstat(g, "B", B~1,dd)
v = variogram(g)
g = gstat(g, model = vgm(1814.2042, "Lin",39061.93, 199.2033), fill.all= T)
g.fit = fit.lmc(v,g,fit.ranges = F)
plot(v, g.fit,type="b",col="black")
but when I use
predict.gstat(g.fit, dd)
I faced with error
""non-positive definite coefficient matrix in structure 1non-positive
definite coefficient
 matrix in structure 2Warning: No Intrinsic Correlation or Linear
Model of Coregionalization found""

-- 

Saman Monfared
Msc Student,
Department of Statistics, Shiraz University,
Shiraz 71454, Iran
Email: Samanmonfared1 at gmail.com

Tel: +98 917 5305167



More information about the R-sig-Geo mailing list