I am currently using the gstat package to create variograms of oil and gas production in counties in West Virginia.
Here is the link to my dataset of SpatialPointsDataFrame and the code to produce a variogram is below
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)
I would like to add a 95% confidence band to this variogram. Is there a function to do this? Is there another way to give an indication of precision?