如何在 R 中计算变差函数时更改滞后距离
how to change the lag distance while calculating a Variogram in R
我正在尝试计算不同滞后距离的实验变差函数值,因此我正在使用 variogramm 命令
variog1 <- variogram((Copper)~1,ds)
但我不知道如何指定所需的滞后距离。
例如我想在 h=(15, 30, 45, 60)
获得一个值
您可以在参数 "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)
我正在尝试计算不同滞后距离的实验变差函数值,因此我正在使用 variogramm 命令
variog1 <- variogram((Copper)~1,ds)
但我不知道如何指定所需的滞后距离。
例如我想在 h=(15, 30, 45, 60)
您可以在参数 "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)