使用R包RandomFields
时,您可以通过添加不同的组件来指定描述数据空间分布的模型。
例如:您的数据中较近的点彼此比较熟悉,因此它们遵循指数变异函数。但是仍然存在相当多的噪音,导致你的变异函数具有金块效应。你会写:
RMexp(var = partial sill, scale = range) + RMnugget(var = nugget)
单独添加金块。
想象一下,您的数据值另外增加为x坐标的线性函数。我很确定您仍然会在上面的模型中添加+ RMtrend()
,指定趋势
现在这就是我陷入困境的地方。我真的没有得到help(RMtrend)
- 页面 - 特别是这些例子让我陷入困境(这个R.p
- 函数做了什么?它的帮助页面对我来说完全是胡言乱语,不知道任何C
- 语言)。
无论如何,RMtrend
只接受一个参数:mean
,它不是公式,只是数字向量或其他RMmodel
。在我看来似乎没有RMmodel
来定义线性或多项式趋势,因此我可以将其交给RMtrend
。我还尝试通过
RMtrend
给它一个公式
RMtrend(RFformula(response Variable ~ x))
稍后会出错,说RFformula
被误用了。并且似乎不可能给RMtrend()
一个矢量具有模拟网格的网格节点数量的长度,指定"均值" (=模拟趋势)每个人
例如。我的网格应该有62250点。如果我在没有趋势的情况下进行模拟,我确实得到了一个包含249行和250列的矩阵,但如果我将带有62250个数字条目的向量移交给RMtrend()
RFsimulate()
,仍然会给出错误(或结果)在R的致命错误中:
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)
那么如何赋予RMtrend
线性或多项式趋势呢?