2013-03-26 6 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

+0

Также, пожалуйста, отредактируйте с ошибкой, которую вы получите. – Arun

ответ

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

Смежные вопросы