現在、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% 信頼帯を追加したいと思います。これを行う機能はありますか？精度を示す別の方法はありますか?