我想使用可分离模型来拟合我的变差函数。 Graler 说我应该使用“具有单独块效应和(联合)基台值为 1 的标准化空间和时间变异函数”
我的密码是:
spatialVgm = vgm(psill=7,sp, range=575, nugget=0)
temporalVgm = vgm(psill=7,t, range=700, nugget=0)
separable <- vgmST("separable", space = spatialVgm, time = temporalVgm,
sill = 14)
fit.ssum.os = fit.StVariogram(mod1,separable,fit.method = 8)
其中 mod1 是经验变异函数。它是否正确?或者我是否必须对变异函数做一些事情来标准化它们并使联合门槛 = 1?
我没有收到任何错误 - 我只是想知道我是否做对了?
我的变异函数代码如下
#dataframe includes the observed data for O18.res and dat is a data.frame with all my other data.
COOR.formatted=SpatialPoints(COOR,proj4string=CRS("+proj=longlat +datum=WGS84"))
DATE.formatted <- unique(as.Date(dat$date))
datST = STSDF(COOR.formatted, DATE.formatted, dataframe,
index)
mod1=variogramST(O18.res~ 1,
datST,tlags=0:23,
na.omit=TRUE)