0

我目前的代码是这样的:

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回归方程,但我不断收到错误。任何帮助将非常感激。

4

1 回答 1

3

目前尚不清楚您要做什么,但您可以像这样分解公式:

form <- y ~ x1 + x2
FUN <- function(x) fitted(lm(form, data = as.data.frame(x)))
rollapplyr(z, 131, FUN, by.column = FALSE)
于 2013-03-26T00:17:53.060 回答