Maximum Likelihood Estimation of a Labour Demand System. An Application of a Model of Seemingly Unrelated Regression Equations with the Regression Errors Composed of Two Components

Abstract  In this paper maximum likelihood estimation is applied to a system of labour demand equations. In view of the data used for estimation, the model of seemingly unrelated regression equations containing errors composed of two components which has been treated by A very [1] is applied. A very's results are presented and slightly generalized.     

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.     

Autocorrelation and dynamic methodology with an application to wage determination models

A demonstration is provided of rigorous, statistical methodology whereby both the type and order of an error process can be identified in dynamic, single equation econometric models. The paper relies heavily upon maximum likelihood estimation, nested likelihood ratio tests and the overfitting or exponentially weighted procedure for model selection. An application of the methodology to a class of quarterly wage determination models is included.     

Missing observations in ARIMA models: Skipping approach versus additive outlier approach

Optimal estimation of missing values in ARMA models is typically performed by using the Kalman filter for likelihood evaluation, ‘skipping’ in the computations the missing observations, obtaining the maximum likelihood (ML) estimators of the model parameters, and using some smoothing algorithm. The same type of procedure has been extended to nonstationary ARIMA models in Gómez and Maravall (1994). An alternative procedure suggests filling in the holes in the series with arbitrary values and then performing ML estimation of the ARIMA model with additive outliers (AO). When the model parameters are not known the two methods differ, since the AO likelihood is affected by the arbitrary values. We develop the proper likelihood for the AO approach in the general non-stationary case and show the equivalence of this and the skipping method. Finally, the two methods are compared through simulation, and their relative advantages assessed; the comparison also includes the AO method with the uncorrected likelihood.     

The multi-state latent factor intensity model for credit rating transitions

A new empirical reduced-form model for credit rating transitions is introduced. It is a parametric intensity-based duration model with multiple states and driven by exogenous covariates and latent dynamic factors. The model has a generalized semi-Markov structure designed to accommodate many of the stylized facts of credit rating migrations. Parameter estimation is based on Monte Carlo maximum likelihood methods for which the details are discussed in this paper. A simulation experiment is carried out to show the effectiveness of the estimation procedure. An empirical application is presented for transitions in a 7 grade rating system. The model includes a common dynamic component that can be interpreted as the credit cycle. Asymmetric effects of this cycle across rating grades and additional semi-Markov dynamics are found to be statistically significant. Finally, we investigate whether the common factor model suffices to capture systematic risk in rating transition data by introducing multiple factors in the model.     

On the estimation of the parameters for the Littlewood model in software reliability

A very well-known model in software reliability theory is that of Littlewood (1980). The (three) parameters in this model are usually estimated by means of the maximum likelihood method. The system of likelihood equations can have more than one solution. Only one of them will be consistent, however. In this paper we present a different, more analytical approach, exploiting the mathematical properties of the log-likelihood function itself. Our belief is that the ideas and methods developed in this paper could also be of interest for statisticians working on the estimation of the parameters of the generalised Pareto distribution. For those more generally interested in maximum likelihood the paper provides a 'practical case', indicating how complex matters may become when only three parameters are involved. Moreover, readers not familiar with counting process theory and software reliability are given a first introduction.     

Spatial clustering of time series via mixture of autoregressions models and Markov random fields

Time series data arise in many medical and biological imaging scenarios. In such images, a time series is obtained at each of a large number of spatially dependent data units. It is interesting to organize these data into model‐based clusters. A two‐stage procedure is proposed. In stage 1, a mixture of autoregressions (MoAR) model is used to marginally cluster the data. The MoAR model is fitted using maximum marginal likelihood (MMaL) estimation via a minorization–maximization (MM) algorithm. In stage 2, a Markov random field (MRF) model induces a spatial structure onto the stage 1 clustering. The MRF model is fitted using maximum pseudolikelihood (MPL) estimation via an MM algorithm. Both the MMaL and MPL estimators are proved to be consistent. Numerical properties are established for both MM algorithms. A simulation study demonstrates the performance of the two‐stage procedure. An application to the segmentation of a zebrafish brain calcium image is presented.     

Analysis of a mixed model for ordinal data by iterative re-weighted REML

An estimation procedure will be presented for a class of threshold models for ordinal data. These models may include both fixed and random effects with associated components of variance on an underlying scale. The residual error distribution on the underlying scale may be rendered greater flexibility by introducing additional shape parameters, e.g. a kurtosis parameter or parameters to model heterogeneous residual variances as a function of factors and covariates. The estimation procedure is an extension of an iterative re-weighted restricted maximum likelihood procedure, originally developed for generalized linear mixed models. This procedure will be illustrated with a practical problem involving damage to potato tubers and with data from animal breeding and medical research from the literature.     

Large sample estimation and testing procedures for dynamic equation systems

Joint two-step estimation procedures which have the same asymptotic properties as the maximum likelihood (ML) estimator are developed for the final equation, transfer function and structural form of a multivariate dynamic model with normally distributed vector-moving average errors. The ML estimator under fixed and known initial values is obtained by iterating the procedure until convergence. The asymptotic distribution of the two-step estimators is used to construct large sample testing procedures for the different forms of the model.     

Estimating a generalized correlation coefficient for a generalized bivariate probit model

In this paper we consider semiparametric estimation of a generalized correlation coefficient in a generalized bivariate probit model. The generalized correlation coefficient provides a simple summary statistic measuring the relationship between the two binary decision processes in a general framework. Our semiparametric estimation procedure consists of two steps, combining semiparametric estimators for univariate binary choice models with the method of maximum likelihood for the bivariate probit model with nonparametrically generated regressors. The estimator is shown to be consistent and asymptotically normal. The estimator performs well in our simulation study.     

An overview of variance component estimation

Variance components estimation originated with estimating error variance in analysis of variance by equating error mean square to its expected value. This equating procedure was then extended to random effects models, first for balanced data (for which minimum variance properties were subsequently established) and later for unbalanced data. Unfortunately, this ANOVA methodology yields no optimum properties (other than unbiasedness) for estimation from unbalanced data. Today it is being replaced by maximum likelihood (ML) and restricted maximum likelihood (REML) based on normality assumptions and involving nonlinear equations that have to be solved numerically. There is also minimum norm quadratic unbiased estimation (MINQUE) which is closely related to REML but with fewer advantages.An invited paper for the ProbaStat '94 conference, Smolenice, Slovakia, May 30–June 3, 1994 Paper number BU-677 in the Biometrics Unit. Cornell University Ithaca NY     

The stochastic frontier production function and average efficiency: An empirical analysis

This paper provides an empirical estimation of a stochastic frontier Cobb-Douglas production function using micro data from a cross-section of Brazilian manufacturing firms. Following a procedure developed by Aigner, Lovell and Schmidt incorporating both stochastic and efficiency disturbance terms in the estimating model, maximum likelihood techniques are used for the estimation of the stochastic frontier. A measure of mean technical efficiency is also developed and employed with the Brazilian data. Unlike the previous empirical exercises carried out with aggregated data, the efficiency disturbance with the Brazilian micro data estimates is not swamped by the stochastic disturbance.     

Markov-switching models with endogenous explanatory variables II: A two-step MLE procedure

This paper proposes a two-step maximum likelihood estimation (MLE) procedure to deal with the problem of endogeneity in Markov-switching regression models. A joint estimation procedure provides us with an asymptotically most efficient estimator, but it is not always feasible, due to the ‘curse of dimensionality’ in the matrix of transition probabilities. A two-step estimation procedure, which ignores potential correlation between the latent state variables, suffers less from the ‘curse of dimensionality’, and it provides a reasonable alternative to the joint estimation procedure. In addition, our Monte Carlo experiments show that the two-step estimation procedure can be more efficient than the joint estimation procedure in finite samples, when there is zero or low correlation between the latent state variables.     

Bayes estimates of muIti-criteria decision alternatives using Monte Carlo integration

A Bayesian procedure is proposed for the estimation of the weights of the alternatives in a multi-criteria decision model with data that stem from pair-wise comparison of alternatives. The prior information restricts the weights to the unit simplex. The posterior results are computed by Monte Carlo integration procedures based on importance sampling. The Bayesian procedure is applied to a case study concerning the choice of a professor of Operations Research (OR). Results are: (1) according to the Bayesian procedure a different candidate would be chosen as professor of OR than according to the maximum likelihood procedure; (2) given the prior and data information, there exists a substantial probability of taking the wrong decision; (3) there exists a ranking of the candidates with a posterior probability greater than one half.     

Bayesian estimation and model selection for spatial Durbin error model with finite distributed lags

In this paper we investigate a spatial Durbin error model with finite distributed lags and consider the Bayesian MCMC estimation of the model with a smoothness prior. We study also the corresponding Bayesian model selection procedure for the spatial Durbin error model, the spatial autoregressive model and the matrix exponential spatial specification model. We derive expressions of the marginal likelihood of the three models, which greatly simplify the model selection procedure. Simulation results suggest that the Bayesian estimates of high order spatial distributed lag coefficients are more precise than the maximum likelihood estimates. When the data is generated with a general declining pattern or a unimodal pattern for lag coefficients, the spatial Durbin error model can better capture the pattern than the SAR and the MESS models in most cases. We apply the procedure to study the effect of right to work (RTW) laws on manufacturing employment.     

Modeling microstructure price dynamics with symmetric Hawkes and diffusion model using ultra-high-frequency stock data

This study examines the theoretical and empirical perspectives of the symmetric Hawkes model of the price tick structure. Combined with the maximum likelihood estimation, the model provides a proper method of volatility estimation specialized in ultra-high-frequency analysis. Empirical studies based on the model using the ultra-high-frequency data of stocks in the S&P 500 are performed. The performance of the volatility measure, intraday estimation, and the dynamics of the parameters are discussed. A new approach of diffusion analogy to the symmetric Hawkes model is proposed with the distributional properties very close to the Hawkes model. As a diffusion process, the model provides more analytical simplicity when computing the variance formula, incorporating skewness and examining the probabilistic property. An estimation of the diffusion model is performed using the simulated maximum likelihood method and shows similar patterns to the Hawkes model.     

Semiparametric three‐step estimation methods for simultaneous equation systems

This paper proposes a new method for estimating a structural model of labour supply in which hours of work depend on (log) wages and the wage rate is considered endogenous. The main innovation with respect to other related estimation procedures is that a nonparametric additive structure in the hours of work equation is permitted. Though the focus of the paper is on this particular application, a three‐step methodology for estimating models in the presence of the above econometric problems is described. In the first step the reduced form parameters of the participation equation are estimated by a maximum likelihood procedure adapted for estimation of an additive nonparametric function. In the second step the structural parameters of the wage equation are estimated after obtaining the selection‐corrected conditional mean function. Finally, in the third step the structural parameters of the labour supply equation are estimated using local maximum likelihood estimation techniques. The paper concludes with an application to illustrate the feasibility, performance and possible gain of using this method. Copyright © 2005 John Wiley & Sons, Ltd.     

A stationary point for the stochastic frontier likelihood

The likelihood function for the stochastic frontier model is shown to possess an unusual stationary point which may or may not be a maximum. A condition is given to determine if the point is a maximum, and the result is interpreted in the context of specification and estimation.     

Estimating the natural rate of unemployment and testing the natural rate hypothesis

How should one measure the natural rate of unemployment? This paper proposes a systems procedure as an alternative to NAIRU. The natural rate is treated as an unobserved state variable in a system that includes measurement equations for the unemployment rate, the rate of wage growth and the rate of inflation. The model is derived from a version of the wage bargaining model of Blanchard and embodies a version of the natural rate hypothesis. The model is estimated by embedding the Kalman filter within the full-information maximum likelihood procedure. For US data, the estimated model implies substantial post-war variation in the natural rate and a negative, but weak, effect of inflation surprises on unemployment. Copyright © 1999 John Wiley & Sons, Ltd.     

Serial correlation in latent discrete variable models

We consider the problems of estimation and testing in models with serially correlated discrete latent variables. A particular case of this is the time series regression model in which a discrete explanatory variable is measured with error. Test statistics are derived for detecting serial correlation in such a model. We then show that the likelihood function can be evaluated by a recurrence relation, and thus maximum likelihood estimation is computationally feasible. An illustrative example of these methods is given, followed by a brief discussion of their applicability to a Markov model of switching regressions.