R Прогноз с запаздывающей зависимой переменной

Question

Использовать функцию lm для установки (Pt = a Pt-1 + b Xt + Dummy variable для каждой четверти), чтобы соответствовать данным образца. Как создать прогноз n.ahead = 12? Не удалось выяснить, как настроить манекен и Pt-1 перед итерацией. Любая помощь приветствуется!

Answer 1

Может быть, это может помочь

#store your model 
model<-your_model 

#get the last pt observation 
last<-dato[nrows(dato$pt), c('pt', 'age')] 

years<-12/4 

#create dummy 
t1<-rep(c(1,0,0,0) , years) 
t2<-rep(c(0,1,0,0) , years) 
t3<-rep(c(0,0,1,0) , years) 
t4<-rep(c(0,0,0,1) , years) 

#create pt observation 
pt<-c(last$pt, rep(NA, length(t1)-1)) 

df<-data.frame(t1=t1,t2=t2,t3=t3,t4=t4,lag_pt=pt, age=last$age) 

df$predict<-NA 

for (i in 1:nrow(df)) 
{ 
df$predict[i]<-predict(model, data=df[i,]) 

if (i!=nrow(df)) 
{df$lag_pt[i+1]<-df$predict[i] } 

}