Reparameterized inverse gamma regression models with varying precision

In this article, we propose a mean linear regression model where the response variable is inverse gamma distributed using a new parameterization of this distribution that is indexed by mean and precision parameters. The main advantage of our new parametrization is the straightforward interpretation of the regression coefficients in terms of the expectation of the positive response variable, as usual in the context of generalized linear models. The variance function of the proposed model has a quadratic form. The inverse gamma distribution is a member of the exponential family of distributions and has some distributions commonly used for parametric models in survival analysis as special cases. We compare the proposed model to several alternatives and illustrate its advantages and usefulness. With a generalized linear model approach that takes advantage of exponential family properties, we discuss model estimation (by maximum likelihood), black further inferential quantities and diagnostic tools. A Monte Carlo experiment is conducted to evaluate the performances of these estimators in finite samples with a discussion of the obtained results. A real application using minerals data set collected by Department of Mines of the University of Atacama, Chile, is considered to demonstrate the practical potential of the proposed model.     

On solving bias‐corrected non‐linear estimation equations with an application to the dynamic linear model

In a seminal paper, Mak, Journal of the Royal Statistical Society B, 55, 1993, 945, derived an efficient algorithm for solving non‐linear unbiased estimation equations. In this paper, we show that when Mak's algorithm is applied to biased estimation equations, it results in the estimates that would come from solving a bias‐corrected estimation equation, making it a consistent estimator if regularity conditions hold. In addition, the properties that Mak established for his algorithm also apply in the case of biased estimation equations but for estimates from the bias‐corrected equations. The marginal likelihood estimator is obtained when the approach is applied to both maximum likelihood and least squares estimation of the covariance matrix parameters in the general linear regression model. The new approach results in two new estimators when applied to the profile and marginal likelihood functions for estimating the lagged dependent variable coefficient in the dynamic linear regression model. Monte Carlo simulation results show the new approach leads to a better estimator when applied to the standard profile likelihood. It is therefore recommended for situations in which standard estimators are known to be biased.     

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.     

A Simple Improvement of the IV‐estimator for the Classical Errors‐in‐Variables Problem

Two measures of an error‐ridden variable make it possible to solve the classical errors‐in‐Variable problem by using one measure as an instrument for the other. It is well known that a second IV‐estimate can be obtained by reversing the roles of the two measures. We explore the optimal linear combination of these two estimates. In a Monte Carlo study, we show that the gain in precision is significant. The proposed estimator also compares well with full information maximum likelihood under normality. We illustrate the method by estimating the capital elasticity in the Norwegian ICT‐industry.     

Estimation and testing for a partially linear single-index spatial regression model

ABSTRACT

Observations recorded on ‘locations’ usually exhibit spatial dependence. In an effort to take into account both the spatial dependence and the possible underlying non-linear relationship, a partially linear single-index spatial regression model is proposed. This paper establishes the estimators of the unknowns. Moreover, it builds a generalized F-test to determine whether or not the data provide evidence on using linear settings in empirical studies. Their asymptotic properties are derived. Monte Carlo simulations indicate that the estimators and test statistic perform well. The analysis of Chinese house price data shows the existence of both spatial dependence and a non-linear relationship.     

Sequential Monte Carlo Methods for Estimating Dynamic Microeconomic Models

This paper develops estimators for dynamic microeconomic models with serially correlated unobserved state variables using sequential Monte Carlo methods to estimate the parameters and the distribution of the unobservables. If persistent unobservables are ignored, the estimates can be subject to a dynamic form of sample selection bias. We focus on single‐agent dynamic discrete‐choice models and dynamic games of incomplete information. We propose a full‐solution maximum likelihood procedure and a two‐step method and use them to estimate an extended version of the capital replacement model of Rust with the original data and in a Monte Carlo study. Copyright © 2015 John Wiley & Sons, Ltd.     

MCMC maximum likelihood for latent state models

This paper develops a pure simulation-based approach for computing maximum likelihood estimates in latent state variable models using Markov Chain Monte Carlo methods (MCMC). Our MCMC algorithm simultaneously evaluates and optimizes the likelihood function without resorting to gradient methods. The approach relies on data augmentation, with insights similar to simulated annealing and evolutionary Monte Carlo algorithms. We prove a limit theorem in the degree of data augmentation and use this to provide standard errors and convergence diagnostics. The resulting estimator inherits the sampling asymptotic properties of maximum likelihood. We demonstrate the approach on two latent state models central to financial econometrics: a stochastic volatility and a multivariate jump-diffusion models. We find that convergence to the MLE is fast, requiring only a small degree of augmentation.     

Simultaneous dependence between firm-level stock returns

We show that use of ordinary least-squares to explore relationships involving firm-level stock returns as the dependent variable in the face of structured dependence between individual firms leads to an endogeneity problem. This in turn leads to biased and inconsistent least-squares estimates. A maximum likelihood estimation procedure that will produce consistent estimates in these situations is illustrated. This is done using methods that have been developed to deal with spatial dependence between regional data observations, which can be applied to situations involving firm-level observations that exhibit a structure of dependence. In addition, we show how to correctly interpret maximum likelihood parameter estimates from these models in the context of firm-level dependence, and provide a Monte Carlo as well as applied illustration of the magnitude of bias that can arise.     

The Spatial Durbin Model and the Common Factor Tests

Abstract

The spatial Durbin model occupies an interesting position in the field of spatial econometrics. It is the reduced form of a model with cross-sectional dependence in the errors and it may be used as the nesting equation in a more general approach of model selection. Specifically, in this equation we obtain the common factor tests (of which the likelihood ratio is the best known) whose objective is to discriminate between substantive and residual dependence in an apparently misspecified equation. Our paper tries to delve deeper into the role of the spatial Durbin model in the problem of specifying a spatial econometric model. We include a Monte Carlo study related to the performance of the common factor tests presented in the paper in small sample sizes.     

Dynamic spatial autoregressive models with autoregressive and heteroskedastic disturbances

We propose a new class of models specifically tailored for spatiotemporal data analysis. To this end, we generalize the spatial autoregressive model with autoregressive and heteroskedastic disturbances, that is, SARAR(1, 1), by exploiting the recent advancements in score‐driven (SD) models typically used in time series econometrics. In particular, we allow for time‐varying spatial autoregressive coefficients as well as time‐varying regressor coefficients and cross‐sectional standard deviations. We report an extensive Monte Carlo simulation study in order to investigate the finite‐sample properties of the maximum likelihood estimator for the new class of models as well as its flexibility in explaining a misspecified dynamic spatial dependence process. The new proposed class of models is found to be economically preferred by rational investors through an application to portfolio optimization.     

Estimating dynamic equilibrium economies: linear versus nonlinear likelihood

This paper compares two methods for undertaking likelihood‐based inference in dynamic equilibrium economies: a sequential Monte Carlo filter and the Kalman filter. The sequential Monte Carlo filter exploits the nonlinear structure of the economy and evaluates the likelihood function of the model by simulation methods. The Kalman filter estimates a linearization of the economy around the steady state. We report two main results. First, both for simulated and for real data, the sequential Monte Carlo filter delivers a substantially better fit of the model to the data as measured by the marginal likelihood. This is true even for a nearly linear case. Second, the differences in terms of point estimates, although relatively small in absolute values, have important effects on the moments of the model. We conclude that the nonlinear filter is a superior procedure for taking models to the data. Copyright © 2005 John Wiley & Sons, Ltd.     

Quasi-likelihood inference in a generalized linear mixed model for balanced data

For modelling the effect of crossed, fixed factors on the response variable in balanced designs with nested stratifications, a generalized linear mixed model is proposed. This model is based on a set of quasi-likelihood assumptions which imply quadratic variance functions. From these variance functions, deviances are obtained to quantify the variation per stratification. The effects of the fixed factors will be tested, an dispersion components will be estimated. The practical use of the model is illustrated by reanalysing a soldering failures problem.     

Heteroskedasticity-consistent covariance matrix estimators for spatial autoregressive models

In the presence of heteroskedasticity, conventional test statistics based on the ordinary least squares (OLS) estimator lead to incorrect inference results for the linear regression model. Given that heteroskedasticity is common in cross-sectional data, the test statistics based on various forms of heteroskedasticity-consistent covariance matrices (HCCMs) have been developed in the literature. In contrast to the standard linear regression model, heteroskedasticity is a more serious problem for spatial econometric models, generally causing inconsistent extremum estimators of model coefficients. This paper investigates the finite sample properties of the heteroskedasticity-robust generalized method of moments estimator (RGMME) for a spatial econometric model with an unknown form of heteroskedasticity. In particular, it develops various HCCM-type corrections to improve the finite sample properties of the RGMME and the conventional Wald test. The Monte Carlo results indicate that the HCCM-type corrections can produce more accurate results for inference on model parameters and the impact effects estimates in small samples.     

DO PEERS AFFECT STUDENT ACHIEVEMENT? EVIDENCE FROM CANADA USING GROUP SIZE VARIATION

We provide the first empirical application of a new approach proposed by Lee (Journal of Econometrics 2007; 140 (2), 333–374) to estimate peer effects in a linear‐in‐means model when individuals interact in groups. Assumingsufficient group size variation, this approach allows to control for correlated effects at the group level and to solve the simultaneity (reflection) problem. We clarify the intuition behind identification of peer effects in the model. We investigate peer effects in student achievement in French, Science, Mathematics and History in secondary schools in the Province of Québec (Canada). We estimate the model using conditional maximum likelihood and instrumental variables methods. We find some evidence of peer effects. The endogenous peer effect is large and significant in Mathematics but imprecisely estimated in the other subjects. Some contextual peer effects are also significant. In particular, for most subjects, the average age of peers has a negative effect on own test score. Using calibrated Monte Carlo simulations, we find that high dispersion in group sizes helps with potential issues of weak identification. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimating and Forecasting with a Dynamic Spatial Panel Data Model*

This study focuses on the estimation and predictive performance of several estimators for the dynamic and autoregressive spatial lag panel data model with spatially correlated disturbances. In the spirit of Arellano and Bond (1991) and Mutl (2006) , a dynamic spatial generalized method of moments (GMM) estimator is proposed based on Kapoor, Kelejian and Prucha (2007) for the spatial autoregressive (SAR) error model. The main idea is to mix non‐spatial and spatial instruments to obtain consistent estimates of the parameters. Then, a linear predictor of this spatial dynamic model is derived. Using Monte Carlo simulations, we compare the performance of the GMM spatial estimator to that of spatial and non‐spatial estimators and illustrate our approach with an application to new economic geography.     

Maximum likelihood estimation and Lagrange multiplier tests for panel seemingly unrelated regressions with spatial lag and spatial errors: An application to hedonic housing prices in Paris

This paper proposes maximum likelihood estimators for panel seemingly unrelated regressions with both spatial lag and spatial error components. We study the general case where spatial effects are incorporated via spatial errors terms and via a spatial lag dependent variable and where the heterogeneity in the panel is incorporated via an error component specification. We generalize the approach of Wang and Kockelman (2007) and propose joint and conditional Lagrange multiplier tests for spatial autocorrelation and random effects for this spatial SUR panel model. The small sample performance of the proposed estimators and tests are examined using Monte Carlo experiments. An empirical application to hedonic housing prices in Paris illustrate these methods. The proposed specification uses a system of three SUR equations corresponding to three types of flats within 80 districts of Paris over the period 1990-2003. We test for spatial effects and heterogeneity and find reasonable estimates of the shadow prices for housing characteristics.     

NON‐LINEAR DSGE MODELS AND THE CENTRAL DIFFERENCE KALMAN FILTER

This paper introduces a quasi maximum likelihood approach based on the central difference Kalman filter to estimate non‐linear dynamic stochastic general equilibrium (DSGE) models with potentially non‐Gaussian shocks. We argue that this estimator can be expected to be consistent and asymptotically normal for DSGE models solved up to third order. These properties are verified in a Monte Carlo study for a DSGE model solved to second and third order with structural shocks that are Gaussian, Laplace distributed, or display stochastic volatility. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimation of Markov regime-switching regression models with endogenous switching

Following Hamilton [1989. A new approach to the economic analysis of nonstationary time series and the business cycle. Econometrica 57, 357–384], estimation of Markov regime-switching regressions typically relies on the assumption that the latent state variable controlling regime change is exogenous. We relax this assumption and develop a parsimonious model of endogenous Markov regime-switching. Inference via maximum likelihood estimation is possible with relatively minor modifications to existing recursive filters. The model nests the exogenous switching model, yielding straightforward tests for endogeneity. In Monte Carlo experiments, maximum likelihood estimates of the endogenous switching model parameters were quite accurate, even in the presence of certain model misspecifications. As an application, we extend the volatility feedback model of equity returns given in Turner et al. [1989. A Markov model of heteroskedasticity, risk, and learning in the stock market. Journal of Financial Economics 25, 3–22] to allow for endogenous switching.     

Convenient links for the estimation of hedonic price indexes: the case of unique,infrequently traded assets

Hedonic methods are a prominent approach in the construction of quality‐adjusted price indexes. This paper shows that the process of computing such indexes is substantially simplified if arithmetic (geometric) price indexes are computed based on exponential (log‐linear) hedonic functions estimated by the Poisson pseudo‐maximum likelihood (ordinary least squares) method. A Monte Carlo simulation study based on housing data illustrates the convenience of the links identified and the very attractive properties of the Poisson estimator in the hedonic framework.     

Computing marginal likelihoods from a single MCMC output

In this article, we propose new Monte Carlo methods for computing a single marginal likelihood or several marginal likelihoods for the purpose of Bayesian model comparisons. The methods are motivated by Bayesian variable selection, in which the marginal likelihoods for all subset variable models are required to compute. The proposed estimates use only a single Markov chain Monte Carlo (MCMC) output from the joint posterior distribution and it does not require the specific structure or the form of the MCMC sampling algorithm that is used to generate the MCMC sample to be known. The theoretical properties of the proposed method are examined in detail. The applicability and usefulness of the proposed method are demonstrated via ordinal data probit regression models. A real dataset involving ordinal outcomes is used to further illustrate the proposed methodology.