2013-03-26 83 views
0

我当前的代码是这样的:rollapply嵌入式参考

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

+0

您好!请通过看看[**如何使一个伟大的可重复的例子**],使您的文章重现性(http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-例如)为我们提供帮助。谢谢。 – Arun 2013-03-26 00:06:07

+0

另外,请编辑您遇到的错误。 – Arun 2013-03-26 00:06:57

回答

3

目前尚不清楚,你想做什么做什么,但你可以分解出这样的公式:

form <- y ~ x1 + x2 
FUN <- function(x) fitted(lm(form, data = as.data.frame(x))) 
rollapplyr(z, 131, FUN, by.column = FALSE) 
+0

谢谢大家!我会看看这个并相应地编辑/改进我的帖子... – user2209633 2013-03-26 01:41:56