我目前正在使用 gstat 包来创建西弗吉尼亚州各县的石油和天然气生产的变异函数。
这是我的 SpatialPointsDataFrame 数据集的链接，生成变异函数的代码如下
library(gstat)
library(rgdal)
prod_wv = readOGR(dsn = "~/variogram_wv", layer = "variogram_wv")
wv.vgm = variogram(log(BOE + 0.0001) ~ 1, data = prod_wv,
prod_wv2, cutoff = cf, width = cf/width_num)
wv.fit = fit.variogram(wv.vgm, vgm("Sph"))
plot(wv.vgm, wv.fit)
我想为这个变异函数添加一个 95% 的置信带。有这样做的功能吗？有没有另一种方法来表示精度？