我需要对多元数据集进行移动窗口线性回归。我在包“robfilter”中使用过 lqd.filter() 和类似的东西,这是我所追求的,但我需要包含 x 变量,因为它不是常量。我是否错过了另一个允许您执行此操作的包?
示例数据框,其中 x=Year 和 y=Value。如果它可以包含错误(对于 y),那将是理想的,但它并不重要。
structure(list(Year = c(1950L, 1923L, 1920L, 1907L, 1896L, 1865L,
1845L, 1822L, 1796L, 1754L, 1732L, 1708L, 1654L, 1613L, 1548L
), Value = c(-0.8, -0.5, -0.3, 0.4, 0.3, -0.2, -0.6, -0.7, -0.9,
-1, -0.7, -0.8, -1.2, -1.4, -1.6), Error = c(0.2, 0.3, 0.2, 0.5,
0.2, 0.2, 0.4, 0.2, 0.1, 0.2, 0.3, 0.2, 0.4, 0.2, 0.2)), .Names = c("Year",
"Value", "Error"), class = "data.frame", row.names = c(NA, -15L))
我的目标是像一个 50 年的移动窗口;以及每个移动窗口的输出回归值,然后我可以将其绘制为曲线(针对其他类似数据集)。
提前谢谢了