我目前正在使用gstat包来创建西弗吉尼亚州的石油和天然气生产的变异图。
以下是我的SpatialPointsDataFrame数据集的link，生成变异函数的代码如下
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％置信区间。有这个功能吗？还有另一种方法来表示精确度吗？