从栅格堆栈R删除线性趋势

时间:2018-10-29 21:46:35

标签: r raster

尝试从1979-2015年的{https://www.northwestknowledge.net/metdata/data/monthly/pr_gridMET.nc)中删除美国每月降水量栅格堆栈中的线性趋势(趋势)。这些数据足够大,以这里的数据为例,在这里有些不守规矩,因此出于效率考虑,我将使用栅格数据包中的数据。我目前的工作模型是在线性模型上使用`raster“ :: calc``并提取残差。我的理解是那些残差是去趋势序列,但是我不是100%确信那是正确的。我正在使用如下:

super.x["foo"]

我觉得我在这里想念什么。谁能确认我在这里的路线正确?如果我对如何适当降低趋势数据没有任何建议,将很有帮助!谢谢!

1 个答案:

答案 0 :(得分:1)

残差去除斜率截距,您将得到异常。也许您只想去除斜坡?在这种情况下,您可以将截距添加到get_residuals

的残差中
q <- residuals(m) + coefficients(m)[1]

或更佳:

q <- residuals(m) + predict(m)[1]

因此您以1年(而不是0年)为基准,并且如果时间为2000:2004

您还可以将去年,年中或平均值作为基准。