如何在R中计算Variogram时更改滞后距离

Question

我正在尝试在不同的滞后距离下计算实验-半方差值，所以我正在使用variogramm命令。

variog1 <- variogram((Copper)~1,ds)

但是我不知道如何指定所需的滞后距离。 例如，我想在h=(15, 30, 45, 60)处得到一个值。

- Ahmed Emam

---

欢迎来到Stackoverflow！请提供更多信息 - 您是否在使用这个铜数据集？ - B--rian

Answer 1

您可以在参数"width"中修复延迟距离，或者您可以使用参数"boundary"来准确固定距离。请参见下一个示例：

library(sp)
library(gstat)
cord <- data.frame(x=rnorm(50,-20,6),y=rnorm(50,40,10))
predictors <- SpatialPointsDataFrame(coords = cord, data = data.frame(val=cord$x**2+cord$y+rnorm(50,sd=0.4)),
                                     proj4string = CRS("+proj=longlat +datum=WGS84"))
vario <- variogram(object = val ~ 1,data =  predictors,width=50)
plot(vario)

vario <- variogram(object = val ~ 1,data =  predictors,width=350)
plot(vario)
vario <- variogram(object = val ~ 1,data =  predictors,boundaries=c(50,100,500,700,800)) 
plot(vario)

- Santiago I. Hurtado