Influence analysis for Poisson inverse Gaussian regression models based on the EM algorithm

For Poisson inverse Gaussian regression models, it is very complicated to obtain the influence measures based on the traditional method, because the associated likelihood function involves intractable expressions, such as the modified Bessel function. In this paper, the EM algorithm is employed as a basis to derive diagnostic measures for the models by treating them as a mixed Poisson regression with the weights from the inverse Gaussian distributions. Several diagnostic measures are obtained in both case-deletion model and local influence analysis, based on the conditional expectation of the complete-data log-likelihood function in the EM algorithm. Two numerical examples are given to illustrate the results.     

GENERALIZED AUTOREGRESSIVE SCORE MODELS WITH APPLICATIONS

We propose a class of observation‐driven time series models referred to as generalized autoregressive score (GAS) models. The mechanism to update the parameters over time is the scaled score of the likelihood function. This new approach provides a unified and consistent framework for introducing time‐varying parameters in a wide class of nonlinear models. The GAS model encompasses other well‐known models such as the generalized autoregressive conditional heteroskedasticity, autoregressive conditional duration, autoregressive conditional intensity, and Poisson count models with time‐varying mean. In addition, our approach can lead to new formulations of observation‐driven models. We illustrate our framework by introducing new model specifications for time‐varying copula functions and for multivariate point processes with time‐varying parameters. We study the models in detail and provide simulation and empirical evidence. Copyright © 2012 John Wiley & Sons, Ltd.     

Comovements and asymmetric tail dependence in state housing prices in the USA: A nonparametric approach

We reexamine the methods used in estimating comovements among US regional home prices and find that there are insufficient moments to ensure a normal limit necessary for employing the quasi‐maximum likelihood estimator. Hence we propose applying the self‐weighted quasi‐maximum exponential likelihood estimator and a bootstrap method to test and account for the asymmetry of comovements as well as different magnitudes across state pairs. Our results reveal interstate asymmetric tail dependence based on observed house price indices rather than residuals from fitting autoregressive–generalized autoregressive conditional heteroskedasticity (AR‐GARCH) models.     

Improved likelihood ratio tests for cointegration rank in the VAR model

We suggest improved tests for cointegration rank in the vector autoregressive (VAR) model and develop asymptotic distribution theory and local power results. The tests are (quasi-)likelihood ratio tests based on a Gaussian likelihood, but as usual the asymptotic results do not require normally distributed innovations. Our tests differ from existing tests in two respects. First, instead of basing our tests on the conditional (with respect to the initial observations) likelihood, we follow the recent unit root literature and base our tests on the full likelihood as in, e.g., Elliott et al. (1996). Second, our tests incorporate a “sign” restriction which generalizes the one-sided unit root test. We show that the asymptotic local power of the proposed tests dominates that of existing cointegration rank tests.     

On Maximum Likelihood Estimation of Dynamic Panel Data Models

We analyse the finite sample properties of maximum likelihood estimators for dynamic panel data models. In particular, we consider transformed maximum likelihood (TML) and random effects maximum likelihood (RML) estimation. We show that TML and RML estimators are solutions to a cubic first‐order condition in the autoregressive parameter. Furthermore, in finite samples both likelihood estimators might lead to a negative estimate of the variance of the individual‐specific effects. We consider different approaches taking into account the non‐negativity restriction for the variance. We show that these approaches may lead to a solution different from the unique global unconstrained maximum. In an extensive Monte Carlo study we find that this issue is non‐negligible for small values of T and that different approaches might lead to different finite sample properties. Furthermore, we find that the Likelihood Ratio statistic provides size control in small samples, albeit with low power due to the flatness of the log‐likelihood function. We illustrate these issues modelling US state level unemployment dynamics.     

Heteroscedasticity diagnostics for <Emphasis Type="Italic">t</Emphasis> linear regression models

The t regression models provide a useful extension of the normal regression models for datasets involving errors with longer-than-normal tails. Homogeneity of variances (if they exist) is a standard assumption in t regression models. However, this assumption is not necessarily appropriate. This paper is devoted to tests for heteroscedasticity in general t linear regression models. The asymptotic properties, including asymptotic Chi-square and approximate powers under local alternatives of the score tests, are studied. Based on the modified profile likelihood (Cox and Reid in J R Stat Soc Ser B 49(1):1–39, 1987), an adjusted score test for heteroscedasticity is developed. The properties of the score test and its adjustment are investigated through Monte Carlo simulations. The test methods are illustrated with land rent data (Weisberg in Applied linear regression. Wiley, New York, 1985). The project supported by NSFC 10671032, China, and a grant (HKBU2030/07P) from the Grant Council of Hong Kong, Hong Kong, China.     

Two are better than one: Volatility forecasting using multiplicative component GARCH-MIDAS models

We examine the properties and forecast performance of multiplicative volatility specifications that belong to the class of generalized autoregressive conditional heteroskedasticity–mixed-data sampling (GARCH-MIDAS) models suggested in Engle, Ghysels, and Sohn (Review of Economics and Statistics, 2013, 95, 776–797). In those models volatility is decomposed into a short-term GARCH component and a long-term component that is driven by an explanatory variable. We derive the kurtosis of returns, the autocorrelation function of squared returns, and the R² of a Mincer–Zarnowitz regression and evaluate the QMLE and forecast performance of these models in a Monte Carlo simulation. For S&P 500 data, we compare the forecast performance of GARCH-MIDAS models with a wide range of competitor models such as HAR (heterogeneous autoregression), realized GARCH, HEAVY (high-frequency-based volatility) and Markov-switching GARCH. Our results show that the GARCH-MIDAS based on housing starts as an explanatory variable significantly outperforms all competitor models at forecast horizons of 2 and 3 months ahead.     

Stochastic volatility models with ARMA innovations: An application to G7 inflation forecasts

We introduce a new class of stochastic volatility models with autoregressive moving average (ARMA) innovations. The conditional mean process has a flexible form that can accommodate both a state space representation and a conventional dynamic regression. The ARMA component introduces serial dependence, which results in standard Kalman filter techniques not being directly applicable. To overcome this hurdle, we develop an efficient posterior simulator that builds on recently developed precision-based algorithms. We assess the usefulness of these new models in an inflation forecasting exercise across all G7 economies. We find that the new models generally provide competitive point and density forecasts compared to standard benchmarks, and are especially useful for Canada, France, Italy, and the U.S.     

Modeling and forecasting of stock index volatility with APARCH models under ordered restriction

This article examines volatility models for modeling and forecasting the Standard & Poor 500 (S&P 500) daily stock index returns, including the autoregressive moving average, the Taylor and Schwert generalized autoregressive conditional heteroscedasticity (GARCH), the Glosten, Jagannathan and Runkle GARCH and asymmetric power ARCH (APARCH) with the following conditional distributions: normal, Student's t and skewed Student's t‐distributions. In addition, we undertake unit root (augmented Dickey–Fuller and Phillip–Perron) tests, co‐integration test and error correction model. We study the stationary APARCH (p) model with parameters, and the uniform convergence, strong consistency and asymptotic normality are prove under simple ordered restriction. In fitting these models to S&P 500 daily stock index return data over the period 1 January 2002 to 31 December 2012, we found that the APARCH model using a skewed Student's t‐distribution is the most effective and successful for modeling and forecasting the daily stock index returns series. The results of this study would be of great value to policy makers and investors in managing risk in stock markets trading.     

Two-stage non Gaussian QML estimation of GARCH models and testing the efficiency of the Gaussian QMLE

In generalized autoregressive conditional heteroskedastic (GARCH) models, the standard identifiability assumption that the variance of the iid process is equal to 1 can be replaced by an alternative moment assumption. We show that, for estimating the original specification based on the standard identifiability assumption, efficiency gains can be expected from using a quasi-maximum likelihood (QML) estimator based on a non Gaussian density and a reparameterization based on an alternative identifiability assumption. A test allowing to determine whether a reparameterization is needed, that is, whether the more efficient QMLE is obtained with a non Gaussian density, is proposed.     

Temporal aggregation of volatility models

In this paper, we consider temporal aggregation of volatility models. We introduce semiparametric volatility models, termed square-root stochastic autoregressive volatility (SR-SARV), which are characterized by autoregressive dynamics of the stochastic variance. Our class encompasses the usual GARCH models and various asymmetric GARCH models. Moreover, our stochastic volatility models are characterized by multiperiod conditional moment restrictions in terms of observables. The SR-SARV class is a natural extension of the class of weak GARCH models. This extension has four advantages: (i) we do not assume that fourth moments are finite; (ii) we allow for asymmetries (skewness, leverage effect) that are excluded from weak GARCH models; (iii) we derive conditional moment restrictions and (iv) our framework allows us to study temporal aggregation of IGARCH models.     

Copula Gaussian graphical models with penalized ascent Monte Carlo EM algorithm

Typical data that arise from surveys, experiments, and observational studies include continuous and discrete variables. In this article, we study the interdependence among a mixed (continuous, count, ordered categorical, and binary) set of variables via graphical models. We propose an ?₁‐penalized extended rank likelihood with an ascent Monte Carlo expectation maximization approach for the copula Gaussian graphical models and establish near conditional independence relations and zero elements of a precision matrix. In particular, we focus on high‐dimensional inference where the number of observations are in the same order or less than the number of variables under consideration. To illustrate how to infer networks for mixed variables through conditional independence, we consider two datasets: one in the area of sports and the other concerning breast cancer.     

Power of Tests for Unit Roots in the Presence of a Linear Trend*

Dickey and Fuller [Econometrica (1981) Vol. 49, pp. 1057–1072] suggested unit‐root tests for an autoregressive model with a linear trend conditional on an initial observation. TPower of tests for unit roots in the presence of a linear trendightly different model with a random initial value in which nuisance parameters can easily be eliminated by an invariant reduction of the model. We show that invariance arguments can also be used when comparing power within a conditional model. In the context of the conditional model, the Dickey–Fuller test is shown to be more stringent than a number of unit‐root tests motivated by models with random initial value. The power of the Dickey–Fuller test can be improved by making assumptions to the initial value. The practitioner therefore has to trade‐off robustness and power, as assumptions about initial values are hard to test, but can give more power.     

Pseudo maximum likelihood estimation of spatial autoregressive models with increasing dimension

Pseudo maximum likelihood estimates are developed for higher-order spatial autoregressive models with increasingly many parameters, including models with spatial lags in the dependent variables both with and without a linear or nonlinear regression component, and regression models with spatial autoregressive disturbances. Consistency and asymptotic normality of the estimates are established. Monte Carlo experiments examine finite-sample behaviour.     

Testing for Linearity

The problem of testing for linearity and the number of regimes in the context of self‐exciting threshold autoregressive (SETAR) models is reviewed. We describe least‐squares methods of estimation and inference. The primary complication is that the testing problem is non‐standard, due to the presence of parameters which are only defined under the alternative, so the asymptotic distribution of the test statistics is non‐standard. Simulation methods to calculate asymptotic and bootstrap distributions are presented. As the sampling distributions are quite sensitive to conditional heteroskedasticity in the error, careful modeling of the conditional variance is necessary for accurate inference on the conditional mean. We illustrate these methods with two applications — annual sunspot means and monthly U.S. industrial production. We find that annual sunspots and monthly industrial production are SETAR(2) processes.     

Forecasting risk measures using intraday data in a generalized autoregressive score framework

A new framework for the joint estimation and forecasting of dynamic value at risk (VaR) and expected shortfall (ES) is proposed by our incorporating intraday information into a generalized autoregressive score (GAS) model introduced by Patton et al., 2019 to estimate risk measures in a quantile regression set-up. We consider four intraday measures: the realized volatility at 5-min and 10-min sampling frequencies, and the overnight return incorporated into these two realized volatilities. In a forecasting study, the set of newly proposed semiparametric models are applied to four international stock market indices (S&P 500, Dow Jones Industrial Average, Nikkei 225 and FTSE 100) and are compared with a range of parametric, nonparametric and semiparametric models, including historical simulations, generalized autoregressive conditional heteroscedasticity (GARCH) models and the original GAS models. VaR and ES forecasts are backtested individually, and the joint loss function is used for comparisons. Our results show that GAS models, enhanced with the realized volatility measures, outperform the benchmark models consistently across all indices and various probability levels.     

REALIZED BETA GARCH: A MULTIVARIATE GARCH MODEL WITH REALIZED MEASURES OF VOLATILITY

We introduce a multivariate generalized autoregressive conditional heteroskedasticity (GARCH) model that incorporates realized measures of variances and covariances. Realized measures extract information about the current levels of volatilities and correlations from high‐frequency data, which is particularly useful for modeling financial returns during periods of rapid changes in the underlying covariance structure. When applied to market returns in conjunction with returns on an individual asset, the model yields a dynamic model specification of the conditional regression coefficient that is known as the beta. We apply the model to a large set of assets and find the conditional betas to be far more variable than usually found with rolling‐window regressions based exclusively on daily returns. In the empirical part of the paper, we examine the cross‐sectional as well as the time variation of the conditional beta series during the financial crises. Copyright © 2014 John Wiley & Sons, Ltd.     

A Note on the Power of Money-Output Causality Tests

This study suggests that some empirical findings against money-output causality can be the consequence of ignoring autoregressive conditional heteroskedastic (ARCH) errors. Monte Carlo results confirm that ARCH effects drastically reduce the power of the standard causality test. The maximum likelihood approach allowing for ARCH effects, on the other hand, provides a good power performance. Using different specifications and sample period, Friedman and Kuttner (1993) and Thomas (1994) report limited evidence of money causing output. We detect significant ARCH effects in the models considered by these studies. Once ARCH effects are explicitly accounted for, we find that the monetary effect is significant though its magnitude is quite small.     

A new marked point process model for the federal funds rate target: Methodology and forecast evaluation

Forecasts of key interest rates set by central banks are of paramount concern for investors and policy makers. Recently it has been shown that forecasts of the federal funds rate target, the most anticipated indicator of the Federal Reserve Bank's monetary policy stance, can be improved considerably when its evolution is modeled as a marked point process (MPP). This is due to the fact that target changes occur in discrete time with discrete increments, have an autoregressive nature and are usually in the same direction. We propose a model which is able to account for these dynamic features of the data. In particular, we combine Hamilton and Jordà's [2002. A model for the federal funds rate target. Journal of Political Economy 110(5), 1135–1167] autoregressive conditional hazard (ACH) and Russell and Engle's [2005. A discrete-state continuous-time model of financial transactions prices and times: the autoregressive conditional multinomial-autoregressive conditional duration model. Journal of Business and Economic Statistics 23(2), 166 – 180] autoregressive conditional multinomial (ACM) model. The paper also puts forth a methodology to evaluate probability function forecasts of MPP models. By improving goodness of fit and point forecasts of the target, the ACH–ACM qualifies as a sensible modeling framework. Furthermore, our results show that MPP models deliver useful probability function forecasts at short and medium term horizons.     

Geometrical Interpretation of the Mean and the Constant in a Box-Jenkins Time Series Model

Certain manuals and computer programs mistakenly identify the mean with the constant in Box-Jenkins time series models. In this paper, it will be shown that (a) the mean and the constant have different values in autoregressive models, and (b) they have an algebraic and graphical relationship.