拟合变异函数R.

Question

我创建了一个像这样的数据框

 head(df)
DISTBin   MSDBin
1       0  4.301515
2      20  5.599879
3      40  5.883236
4      60  6.430584
5      80  7.396374
6     100  7.561652

其中第一列是距离矢量（分箱），第二列是给定数量的均方差矢量。

我希望使用高斯模型拟合df $ MSDBin~df $ DISTBin。

我不能使用gstat或类似物，因为我没有变异图对象（我必须自己计算变异函数以丢弃来自同一ID的对）。

任何帮助？

Answer 1

我自己找到了答案。

Gau <- function(h,tau,sigma.sq,range){
    return(tau + sigma.sq * (1-exp(-(h*h)/(range*range))))
}

Gau_Fit <- nls(MSDBin ~ Gau(DISTBin,tau,sigma.sq,range),
                 data = df,
                 start = list(tau = 0, sigma.sq = 10, range = 1000))