Abstract:At the heart of many econometric models are a linear function and a normal error. Examples include the classical small-sample linear regression model and the probit, ordered probit, multinomial probit, tobit, interval regression, and truncated-distribution regression models. Because the normal distribution has a natural multidimensional generalization, such models can be combined into mul-tiequation systems in which the errors share a multivariate normal distribution. The literature has historically focused on multistage procedures for fitting mixed models, which are more efficient computationally, if less so statistically, than maximum likelihood. Direct maximum likelihood estimation has been made more practical by faster computers and simulated likelihood methods for estimating higher-dimensional cumulative normal distributions. Such simulated likelihood methods include the Geweke–Hajivassiliou-Keane algorithm (Geweke, 1989, Econometrica 57: 1317–1339; Hajivassiliou and McFadden, 1998, Econometrica 66: 863–896; Keane, 1994, Econometrica 62: 95–116). Maximum likelihood also facilitates a generalization to switching, selection, and other models in which the number and types of equations vary by observation. The Stata command cmp fits seemingly unrelated regressions models of this broad family. Its estimator is also consistent for recursive systems in which all endogenous variables appear on the right-hand sides as observed. If all the equations are structural, then estimation is full-information maximum likelihood. If only the final stage or stages are structural, then estimation is limited-information maximum likelihood. cmp can mimic a score of built-in and user-written Stata commands. It is also appropriate for a panoply of models that previously were hard to estimate. Heteroskedasticity, however, can render cmp inconsistent. This article explains the theory and implementation of cmp and of a related Mata function, ghk2(), that implements the Geweke–Hajivassiliou–Keane algorithm.

IPF for contingency tables with missing data via the ECM algorithm

Inverse Probability-Weighted Estimation for Dynamic Structural Equation Model with Missing Data

Maximum-Likelihood-Estimation Via the Ecm Algorithm - A General Framework

RECURSIVE PROCEDURES FOR HIERARCHICAL LOGLINEAR MODELS ON HIGH-DIMENSIONAL CONTINGENCY TABLES

Parameter Expansion to Accelerate Em: the Px-Em Algorithm

The ECME Algorithm: A Simple Extension of EM and ECM with Faster Monotone Convergence

A Robust Functional EM Algorithm for Incomplete Panel Count Data

Semiparametric Inference for Estimating Equations with Nonignorably Missing Covariates

An Imputation-Consistency Algorithm for High-Dimensional Missing Data Problems and Beyond

Accelerating L1-penalized expectation maximization algorithm for latent variable selection in multidimensional two-parameter logistic models

A Computational Note on the Application of the Supplemented EM Algorithm to Item Response Models

An MCMC method for the evaluation of the Fisher information matrix for non-linear mixed effect models

Empirical Likelihood for Estimating Equations with Missing Values

Semiparametric Double Robust and Efficient Estimation for Mean Functionals with Response Missing at Random

Evaluating tree-based imputation methods as an alternative to MICE PMM for drawing inference in empirical studies

Fitting Fully Observed Recursive Mixed-process Models with cmp

A Robust and Flexible EM Algorithm for Mixtures of Elliptical Distributions with Missing Data

Sparse high-dimensional linear mixed modeling with a partitioned empirical Bayes ECM algorithm

The Improved EMS Algorithm for Latent Variable Selection in M3PL Model

An empirical comparison of some missing data treatments in PLS-SEM