我需要对多变量数据集进行移动窗口线性回归。我使用了lqd.filter()和类似的“robfilter”包,这是我追求的东西,但我需要包含x变量,因为它不是常量。我错过了另一个允许你这样做的套餐吗?
示例数据框,其中x =年,y =值。如果它可以包含那些理想的错误(对于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年移动的窗户;以及每个移动窗口的输出回归值,然后我可以将其绘制为曲线(与其他类似数据集对比)。
非常感谢提前