1

使用 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线性或多项式趋势呢?

4

0 回答 0