HAREstimate | R Documentation |
HAR estimation
HAREstimate(RM, BPV = NULL, RQ = NULL, periods = c(1,5,22), periodsJ = NULL, periodsRQ = NULL, type = "HAR", insanityFilter = TRUE, h = 1)
RM |
A |
BPV |
A |
RQ |
A |
periods |
A |
periodsJ |
A |
periodsRQ |
A |
type |
A |
insanityFilter |
A |
h |
A |
The estimates for the HARQ and HARQ-J models differ slightly from the results of BPQ (2016). This is due to a small difference in the demeaning approach for the realized quarticity. Here, the demeaning is done with mean(RQ) over all periods.
A HARModel
object
Emil Sjoerup
Corsi, F. 2009, A Simple Approximate Long-Memory Model
of Realized Volatility, Journal of Financial Econometrics, 174–196.
Bollerslev, T., Patton, A., Quaedvlieg, R. 2016, Exploiting the errors: A simple approach for improved volatility forecasting, Journal of Econometrics , vol.192, issue 1, 1-18.
#Vanilla HAR from Corsi(2009) #load data data("SP500RM") SP500rv = SP500RM$RV #Estimate the HAR model: FitHAR = HAREstimate(RM = SP500rv, periods = c(1,5,22)) #extract the estimated coefficients: coef(FitHAR) #plot the fitted values plot(FitHAR) #calculate the Q-like loss-function: mean(qlike(FitHAR)) #HAR-J: #load data data("SP500RM") SP500rv = SP500RM$RV SP500bpv = SP500RM$BPV #Estimate the HAR-J model: FitHARJ = HAREstimate(RM = SP500rv, BPV = SP500bpv, periods = c(1,5,22), periodsJ = c(1,5,22), type = "HARJ" ) #Calculate the Q-like loss-function: mean(qlike(FitHARJ)) #HAR-Q of BPQ(2016) with weekly aggregation #load data data("SP500RM") SP500rv = SP500RM$RV SP500rq = SP500RM$RQ #Estimate the HAR-Q model: FitHARQ = HAREstimate(RM = SP500rv, RQ = SP500rq, periods = c(1,5,22), periodsRQ = c(1,5,22), type = "HARQ", h = 5) #Show the model: show(FitHARQ) #Extract the coefficients: HARQcoef = coef(FitHARQ) #HARQ-J of BPQ(2016) with monthly aggregation #load data data("SP500RM") SP500rv = SP500RM$RV SP500rq = SP500RM$RQ SP500bpv = SP500RM$BPV #Estimate the HARQ-J model: FitHARQJ = HAREstimate(RM = SP500rv, BPV = SP500bpv, RQ = SP500rq, periods = c(1,5,22), periodsJ = c(1), periodsRQ = c(1), type = "HARQ-J", h = 22) #show the model: show(FitHARQJ)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.