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: samanmonfar...@gmail.com Tel: +98 917 5305167 _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo