2

Eu tenho o seguinte processo:yt=2+0.7yt-1+Et. Quero fazer pseudo previsões fora da amostra, através de rolling regressions e com janela móvel contendo 50 observações. A previsão deve ser realizada utilizando-se os verdadeiros coeficientes:y^t+1=2+0.7yt e comparar com a previsão para passeio aleatório, ou seja,y^t+1=yt. Se fosse um AR (1) qualquer, comparado com um random walk, eu consigo fazer da seguinte forma:

set.seed(1) 
T=100             
e=rnorm(T)
rt(100,5)                      #Para gerar distribuição Student com 10 graus de liberdade                             
mu=2
phi=0.7 
y=matrix(0,T,1)                   
y[1]=mu+e[1]                   
year_start=1
year_end=100
first_date=c(year_start)
final_date=c(year_end)
x=ts(e,start=first_date,frequency=1)  
T=length(x)
R=round(T/2)    

f=ts(0*1:(T-R-1),start=c(year_start+R+1),frequency=1)                  
frw=ts(0*1:(T-R-1),start=c(year_start+R+1),end=final_date,frequency=1)

Rolling Forecast:

ar=1
for (i in 1:length(f)) {
    xt=window(x,start=c(year_start+i),end=c(year_start,R+i))
        xt1=window(x,start=c(year_start+i-1),end=c(year_start,R+i-1))
    for(p in 1:ar){xt1=cbind(xt1,lag(xt,-p))}
    y=na.omit(xt1)
    reg=lm(y[,1]~y[,2])
    cond.x=rev(tail(xt,ar))
    far1[i]=reg$coefficients%*%c(1,cond.x)                  frw[i]=tail(xt,1)               
}
far1    #valores previstos pelo AR(1)   
frw     #valores previstos pelo RW

Mas não estou conseguindo realizar esse procedimento incorporando os coeficientes verdadeiros (mu=2 e phi=0.7), utilizando por exemplo for (i in 2:R){y[i+1]=mu+phi*y[i]}.

2

0

Você deve fazer log-in para responder a esta pergunta.

Pesquise outras perguntas com a tag .