使用 R 包RandomFields
时,您可以通过添加不同的组件来指定描述数据空间分布的模型。
例如:数据中较近的点比较远的点彼此更熟悉,因此它们遵循指数变异函数。但是仍然存在相当多的噪音,导致您的变异函数具有块金效应。你会写:
RMexp(var = partial sill, scale = range) + RMnugget(var = nugget)
分别添加金块。
想象一下,您的数据值作为 x 坐标的线性函数额外增加。我很确定您仍然会添加+ RMtrend()
到上述模型中,指定趋势。
这就是我现在卡住的地方。我真的没有得到help(RMtrend)
-page - 尤其是示例让我陷入困境(这个R.p
-function 有什么作用?它的帮助页面对我来说完全是胡言乱语,不知道任何C
-language)。
无论如何,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
线性或多项式趋势呢?