2013-06-01 11 views
5

La prima volta che faccio una domanda qui, farò del mio meglio per essere esplicito - ma fammi sapere se dovrei fornire maggiori informazioni! Secondo, questa è una lunga domanda ... speriamo sia semplice da risolvere per qualcuno;)! Quindi, usando "R", sto modellando modelli GARCH multivariati basati su alcuni documenti (Manera et al., 2012).R - Modellazione multivariato GARCH (rugarch e ccgarch)

I modelli di correlazione condizionale costante (CCC) e di correlazione dinamica condizionata (DCC) con regressori esterni nelle equazioni medie; utilizzando la versione 3.0.1 "R" con il pacchetto "rugarch" versione 1.2-2 per GARCH univariate con regressori esterni e il pacchetto "ccgarch" (versione 0.2.0-2) per i modelli CCC/DCC. (Attualmente sto esaminando il pacchetto "rmgarch", ma sembra essere solo per il DCC e anch'io ho bisogno del modello CCC.)

Ho problemi nelle equazioni medie dei miei modelli. Nel documento che ho menzionato sopra, le stime dei parametri dell'equazione media tra i modelli CCC e DCC cambiano! E non so come lo farei in R ... (attualmente, guardando su Google e nel libro di Tsay "analisi delle serie temporali finanziarie" e nel libro di Engle "Anticipazione delle correlazioni" per trovare il mio errore)

Quello che intendo per "le mie equazioni medie non cambiano tra i modelli CCC e DCC", è il seguente: Specifico il GARCH univariato per le mie serie temporali n = 5 con il pacchetto rugarch. Quindi, utilizzo i parametri delle stime del GARCH (termini ARCH + GARCH) e li utilizzo sia per le funzioni CCC e DCC "eccc.sim()" e "dcc.sim()". Quindi, dalle funzioni eccc.estimation() e dcc.estimation(), posso recuperare le stime per le equazioni di varianza e le matrici di correlazione. Ma non per l'equazione media.

Inserisco il codice R (riproducibile e il mio originale) solo per i modelli univariati e solo per il modello CCC. Grazie già per aver letto il mio post !!!!!

Nota: nel seguente codice, "data.repl" è un oggetto "zoo" di dim 843x22 (9 serie di restituzioni merceologiche giornaliere e serie di variabili esplicative). Il GARCH multivariato è solo per 5 serie.

codice riproducibile:

# libraries: 
library(rugarch) 
library(ccgarch) 
library(quantmod) 
# Creating fake data: 
dataRegr <- matrix(rep(rnorm(3149, 11, 1),1), ncol=1, nrow=3149) 
dataFuelsLag1 <- matrix(rep(rnorm(3149, 24, 8),2), ncol=2, nrow=3149) 
#S&P 500 via quantmod and Yahoo Finance 
T0 <- "2000-06-23" 
T1 <- "2012-12-31" 
getSymbols("^GSPC", src="yahoo", from=T0, to=T1) 
sp500.close <- GSPC[,"GSPC.Close"], 
getSymbols("UBS", src="yahoo", from=T0, to=T1) 
ubs.close <- UBS[,"UBS.Close"] 
dataReplic <- merge(sp500.close, ubs.close, all=TRUE) 
dataReplic[which(is.na(dataReplic[,2])),2] <- 0 #replace NA 

### (G)ARCH modelling ### 
######################### 
# External regressors: macrovariables and all fuels+biofuel Working's T index 
ext.regr.ext <- dataRegr 
regre.fuels <- cbind(dataFuelsLag1, dataRegr) 
### spec of GARCH(1,1) spec with AR(1) ### 
garch11.fuels <- as.list(1:2) 
for(i in 1:2){ 
    garch11.fuels[[i]] <- ugarchspec(mean.model = list(armaOrder=c(1,0), 
                external.regressors = as.matrix(regre.fuels[,-i]))) 
} 

### fit of GARCH(1,1) AR(1) ### 
garch11.fuels.fit <- as.list(1:2) 
for(i in 1:2){ 
    garch11.fuels.fit[[i]] <- ugarchfit(garch11.fuels[[i]], dataReplic[,i]) 
} 
################################################################## 
#### CCC fuels: with external regression in the mean eqaution #### 
################################################################## 
nObs <- length(data.repl[-1,1]) 
coef.unlist <- sapply(garch11.fuels.fit, coef) 
cccFuels.a <- rep(0.1, 2) 
cccFuels.A <- diag(coef.unlist[6,]) 
cccFuels.B <- diag(coef.unlist[7, ]) 
cccFuels.R <- corr.test(data.repl[,fuels.ind], data.repl[,fuels.ind])$r 

# model=extended (Jeantheau (1998)) 
ccc.fuels.sim <- eccc.sim(nobs = nObs, a=cccFuels.a, A=cccFuels.A, 
          B=cccFuels.B, R=cccFuels.R, model="extended") 
ccc.fuels.eps <- ccc.fuels.sim$eps 
ccc.fuels.est <- eccc.estimation(a=cccFuels.a, A=cccFuels.A, 
           B=cccFuels.B, R=cccFuels.R, 
           dvar=ccc.fuels.eps, model="extended") 
ccc.fuels.condCorr <- round(corr.test(ccc.fuels.est$std.resid, 
             ccc.fuels.est$std.resid)$r,digits=3) 

mio codice originale:

### (G)ARCH modelling ### 
######################### 
# External regressors: macrovariables and all fuels+biofuel Working's T index 
ext.regr.ext <- as.matrix(data.repl[-1,c(10:13, 16, 19:22)]) 
regre.fuels <- cbind(fuel.lag1, ext.regr.ext) #fuel.lag1 is the pre-lagged series 
### spec of GARCH(1,1) spec with AR(1) ### 
garch11.fuels <- as.list(1:5) 
for(i in 1:5){ 
    garch11.fuels[[i]] <- ugarchspec(mean.model = list(armaOrder=c(1,0), 
            external.regressors = as.matrix(regre.fuels[,-i]))) 
}# regre.fuels[,-i] => "-i" because I model an AR(1) for each mean equation 

### fit of GARCH(1,1) AR(1) ### 
garch11.fuels.fit <- as.list(1:5) 
for(i in 1:5){ 
    j <- i 
    if(j==5){j <- 7} #because 5th "fuels" is actually column #7 in data.repl 
    garch11.fuels.fit[[i]] <- ugarchfit(garch11.fuels[[i]], as.matrix(data.repl[-1,j]))) 
} 

#fuelsLag1.names <- paste(cmdty.names[fuels.ind], "(-1)") 
fuelsLag1.names <- cmdty.names[fuels.ind] 
rowNames.ext <- c("Constant", fuelsLag1.names, "Working's T Gasoline", "Working's T Heating Oil", 
       "Working's T Natural Gas", "Working's T Crude Oil", 
       "Working's T Soybean Oil", "Junk Bond", "T-bill", 
       "SP500", "Exch.Rate") 
ic.n <- c("Akaike", "Bayes") 
garch11.ext.univSpec <- univ.spec(garch11.fuels.fit, ols.fit.ext, rowNames.ext, 
            rowNum=c(1:15), colNames=cmdty.names[fuels.ind], 
            ccc=TRUE) 
################################################################## 
#### CCC fuels: with external regression in the mean eqaution #### 
################################################################## 
# From my GARCH(1,1)-AR(1) model, I extract ARCH and GARCH 
# in order to model a CCC GARCH model: 
nObs <- length(data.repl[-1,1]) 
coef.unlist <- sapply(garch11.fuels.fit, coef) 

cccFuels.a <- rep(0.1, length(fuels.ind)) 
cccFuels.A <- diag(coef.unlist[17,]) 
cccFuels.B <- diag(coef.unlist[18, ]) 
#based on Engle(2009) book, page 31: 
cccFuels.R <- corr.test(data.repl[,fuels.ind], data.repl[,fuels.ind])$r 

# model=extended (Jeantheau (1998)) 
# "allow the squared errors and variances of the series to affect 
# the dynamics of the individual conditional variances 
ccc.fuels.sim <- eccc.sim(nobs = nObs, a=cccFuels.a, A=cccFuels.A, 
            B=cccFuels.B, R=cccFuels.R, model="extended") 
ccc.fuels.eps <- ccc.fuels.sim$eps 
ccc.fuels.est <- eccc.estimation(a=cccFuels.a, A=cccFuels.A, 
              B=cccFuels.B, R=cccFuels.R, 
              dvar=ccc.fuels.eps, model="extended") 
ccc.fuels.condCorr <- round(corr.test(ccc.fuels.est$std.resid, 
             ccc.fuels.est$std.resid)$r,digits=3) 
colnames(ccc.fuels.condCorr) <- cmdty.names[fuels.ind] 
rownames(ccc.fuels.condCorr) <- cmdty.names[fuels.ind] 
lowerTri(ccc.fuels.condCorr, rep=NA) 

risposta

3

Sei consapevole che c'è un intero pacchetto rmgarch per i modelli GARCH multivariati?

Per la sua descrizione, si copre

fattibili modelli GARCH multivariati tra cui DCC, GO-GARCH e Copula-GARCH.

+0

Grazie per la risposta! Ne sono davvero consapevole. Ma sfortunatamente avrei bisogno del modello CCC che non è incluso (almeno, da quello che ho visto nel file di aiuto, vignetta del pacchetto rmgarch). – cmembrez

+0

Userò semplicemente il pacchetto ccgarch per il modello CCC e il rmgarch per il modello DCC. Grazie per il tuo contributo – cmembrez