[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