尝试从1979-2015年的{https://www.northwestknowledge.net/metdata/data/monthly/pr_gridMET.nc)中删除美国每月降水量栅格堆栈中的线性趋势(趋势)。这些数据足够大,以这里的数据为例,在这里有些不守规矩,因此出于效率考虑,我将使用栅格数据包中的数据。我目前的工作模型是在线性模型上使用`raster“ :: calc``并提取残差。我的理解是那些残差是去趋势序列,但是我不是100%确信那是正确的。我正在使用如下:
super.x["foo"]
我觉得我在这里想念什么。谁能确认我在这里的路线正确?如果我对如何适当降低趋势数据没有任何建议,将很有帮助!谢谢!
答案 0 :(得分:1)
残差去除斜率和截距,您将得到异常。也许您只想去除斜坡?在这种情况下,您可以将截距添加到get_residuals
q <- residuals(m) + coefficients(m)[1]
或更佳:
q <- residuals(m) + predict(m)[1]
因此您以1年(而不是0年)为基准,并且如果时间为2000:2004
您还可以将去年,年中或平均值作为基准。