fit_Exp <- fit.variogram(vgm, Exp_model, fit.method = 2),如何查看这个模型的块金效应等参数
时间: 2024-03-27 09:35:41 浏览: 118
在使用fit.variogram()函数拟合半变异函数模型时,可以设置参数fit.method = 2来查看拟合模型的块金效应等参数。具体步骤如下:
1. 在fit.variogram()函数中设置fit.method = 2。例如:
```
library(gstat)
data(meuse)
coordinates(meuse) <- ~x+y
v <- variogram(log(zinc)~1, meuse)
Exp_model <- vgm("Exp", range = 100, nugget = 0.1, sill = 1)
fit_Exp <- fit.variogram(v, Exp_model, fit.method = 2)
```
2. 运行fit.variogram()函数进行半变异函数模型拟合,程序会输出一些关于拟合结果的信息,例如:
```
[using ordinary least squares fit]
[variogram fit OK]
[model fit: range = 129.6759, sill = 0.875, nugget = 0.125]
```
其中,[using ordinary least squares fit]表示正在使用最小二乘法进行拟合;[variogram fit OK]表示半变异函数模型拟合成功;[model fit: range = 129.6759, sill = 0.875, nugget = 0.125]表示拟合的半变异函数模型的参数,包括范围参数range、平台值sill和截距nugget。
通过查看这些信息,可以了解拟合的半变异函数模型的性质,例如它是否具有块金效应、是否存在趋势等。如果模型具有块金效应,可以进一步进行克里金插值分析。
阅读全文