R:理解RMtrend所需的帮助(随机字段)
使用R软件包R:理解RMtrend所需的帮助(随机字段),r,geospatial,sp,geostatistics,R,Geospatial,Sp,Geostatistics,使用R软件包随机字段时,可以通过添加不同的组件来指定描述数据空间分布的模型。 例如:数据中较近的点比较远的点彼此更熟悉,因此它们遵循指数变异函数。但是仍然有相当多的噪音,导致你的变异函数有一个金块效应。你会写: RMexp(var = partial sill, scale = range) + RMnugget(var = nugget) 单独添加熔核。 假设数据值作为x坐标的线性函数额外增加。我很确定您仍然会在上述模型中添加+RMtrend(),指定趋势。 这就是我被卡住的地方。我真的没有
随机字段时,可以通过添加不同的组件来指定描述数据空间分布的模型。
例如:数据中较近的点比较远的点彼此更熟悉,因此它们遵循指数变异函数。但是仍然有相当多的噪音,导致你的变异函数有一个金块效应。你会写:
RMexp(var = partial sill, scale = range) + RMnugget(var = nugget)
单独添加熔核。
假设数据值作为x坐标的线性函数额外增加。我很确定您仍然会在上述模型中添加+RMtrend()
,指定趋势。
这就是我被卡住的地方。我真的没有得到帮助(RMtrend)
-页面——尤其是示例让我束手无策(这个R.p
-函数做什么?它的帮助页面对我来说完全是胡言乱语,不知道任何C
-语言)
无论如何,RMtrend
只接受一个参数:mean
,它不能是公式,只能是数值向量或另一个RMmodel
。在我看来,似乎没有RMmodel
来定义线性或多项式趋势,这样我就可以把它交给RMtrend
。我还试图欺骗RMtrend
,通过
RMtrend(RFformula(response Variable ~ x))
后来出现了一个错误,说RFformula
被误用了。而且似乎不可能给RMtrend()
我的网格应该有62250个点。如果我在没有趋势的情况下进行模拟,我确实会得到一个有249行250列的矩阵,但是如果我将一个有62250个数字项的向量交给RMtrend()
RFsimulate()
仍会给我一个错误(或导致致命错误R):
那么,如何才能给RMtrend
一个线性或多项式趋势?我想要一个标签“地质统计学”我想要一个标签“地质统计学”
Error in rfInit(model = list("Simulate", setseed = eval(parse(text = "quote(set.seed(seed=seed))")), :
'RFsimulate' : Number of data (62250) not a multiple of the number of locations (62567 x 0)