R:需要帮助理解RMtrend(RandomFields)

时间:2016-08-14 13:16:31

标签: r geospatial sp geostatistics

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

0 个答案:

没有答案