我想用克里格污染样本进行插值,但是当我拟合我的变异函数时,范围是负的,它必须是正的。我运行了这个:
head(muestras)
Index SAMPLE ZONE EASTING NORTHING Latitude Longitude Mg Bal Al
1 1 Al01 19K 535234 7528392 -22.35000 -68.65781 20155.623 523711.3 14897.99
2 2 Al02 19K 538091 7533652 -22.30242 -68.63019 18240.623 527488.0 30534.27
3 3 Al03 19K 540850 7549211 -22.16180 -68.60380 12797.157 533882.2 45389.78
4 4 Al04 19K 543208 7561915 -22.04697 -68.58127 7748.070 539969.0 60457.23
5 5 Al05 19K 552494 7579475 -21.88807 -68.49184 6646.537 542206.5 61427.67
6 6 Al06 19K 553310 7571159 -21.96318 -68.48367 6764.483 545600.1 62231.17
coordinates(muestras) <- ~ Latitude + Longitude
#Model
mod <- lm(Ca ~ Latitude + Longitude, as.data.frame(muestras))
#Variogram
variog_mod <- variogram(Ca ~ 1, data = muestras)
plot(variog_mod)
fit.vgm <- fit.variogram(variog_mod, vgm("Sph", nugget = 1060365658,
psill = 151716381,
range = -1.776595))
尝试运行后fit.variogram()
,R 抛出以下错误:
Error in vgm("Sph", nugget = 1060365658, psill = 151716381, range = -1.776595) :
range should be positive
你能帮我设定正的范围吗?我会很感激。
问候。