使用 fit.variogram 估计的系数与 model 拟合的 plot 不匹配（gstat R 包）

Question

我正在使用 package gstat 将指数 model 拟合到经验半变异函数。 我的数据集（数据）具有十进制度的纬度和经度值。 使用fit.variogram()估计的范围与我通过查看 model 的 plot 所期望的范围不匹配。 请参阅下面的代码和数字。

这是来自fit.variogram()的 output ：

   model    psill    range
1   Nug 26.66719  0.00000
2   Exp 26.26643 25.76976

从这里开始，我预计与 ~53 (26.66719+26.26643) 的窗台相对应的范围为 25.76976。 然而，当查看 plot 时，对应于 ~53 的门槛的范围大于 60 公里。

library(gstat)  
library(sp)

coordinates(data)= ~ longitude+latitude
proj4string(data)=CRS("+longlat")
TheVariogram=variogram(Svmean~1, data=data, cutoff=90)
plot(TheVariogram) #to get seed parameters
lzn.fit <- fit.variogram(TheVariogram, model = vgm(psill=60, model="Exp", nugget=20, range=70),fit.method=6)
plot(TheVariogram, pl=F, model=lzn.fit, xlab="Distance (km)", ylab=expression(Semivariance (dB^2)))

Answer 1

报告的范围是 model exp(-h/a) 中的范围参数 a，而不是指数 model 达到其（部分）门槛的 95% 时的值。 要获得后一个值，请乘以（大约）3。