Ich habe eine bivariate Zeitreihe, z_tin der z_1tdie Veränderung der monatlichen US-Schatzwechsel (Laufzeit 3 Monate) und z_2tdie prozentuale Inflationsrate des monatlichen US-Verbraucherpreisindex (VPI) angegeben sind. Der verwendete VPI ist der Verbraucherpreisindex für alle städtischen Verbraucher: alle Artikel (CPIAUCSL). Die Originaldaten werden von der Federal Reserve Bank of St. Louis heruntergeladen. Die CPI-Rate beträgt das 100-fache der ersten Differenz des Protokoll-CPI-Index. Ich möchte das angegebene VAR-Modell anpassen und die Anpassung durch einen Befehl R ( refVaraus Paket MTSoder restrictaus Paket vars) mit dem Schwellenwert 1,65 vereinfachen .
Ich fand diese Übung ( pdf ) auf der Website von R. Tsay an der Universität von Chicago. Die Daten sind hier .
Was ich bis jetzt gemacht habe, ist folgendes:
y <- diff(zt[,3])
lot(y, type="l", ylab="tb3m")
# difference
x <- diff(log(zt[,4]))
plot(x, type="l", ylab="CPI rate")
new <- cbind(x, y)
# order selection gives VAR(6)
VARselect(new, lag.max=9, type="const")
data1 <- data[,c("tb3m","cpiaucsl")]
fit <- VAR(data1,p=6)
fit
restrict(fit, method="ser", thresh=1.65, resmat=T)
restrictund VARgib mir nicht die richtigen Ergebnisse oder die gleichen Koeffizienten des Var-Modells in den Antworten im PDF.