我目前的代码是这样的:
Model <- head(rollapplyr(z, width = 131, function(x)
fitted(lm(y ~x1+ x2, data = as.data.frame(x))), by.column = FALSE)[,131],3429)
我想这样做:
fit1 <- lm(y ~x1+ x2)
Model <- head(rollapplyr(z, width = 131, function(x) fitted(fit1),
data = as.data.frame(x))), by.column = FALSE)[,131],3429)
基本上,我只是想加入fit1
回归方程,但我不断收到错误。任何帮助将非常感激。