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.

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

Generating contingency tables with fixed marginal probabilities and dependence structures described by loglinear models

IPF for contingency tables with missing data via the ECM algorithm

Limiting behavior of recursive M-estimators in multivariate linear regression models and their asymptotic efficiencies

Fast maximum likelihood estimation for general hierarchical models

Mixed‐Effect Hybrid Models for Longitudinal Data with Nonignorable Dropout

Decorrelated Empirical Likelihood for Generalized Linear Models with High-Dimensional Longitudinal Data

Recursive Least-Squares Method with Membership Functions

The Best Fit Bayesian Hierarchical Generalized Linear Model Selection Using Information Complexity Criteria in the MCMC Approach

A High-dimensional M-estimator Framework for Bi-level Variable Selection

On the correspondence of deviances and maximum likelihood and interval estimates from log-linear to logistic regression modelling

A Constructive Procedure for Modeling Categorical Variables: Log-Linear and Logit Models

On the maximum likelihood estimation in general log-linear models

Supervised Dimensionality Reduction Methods via Recursive Regression

LOCALLY D-OPTIMAL DESIGNS FOR HIERARCHICAL RESPONSE EXPERIMENTS

A Recursive Partitioning Approach for Dynamic Discrete Choice Modeling in High Dimensional Settings

Fitting Fully Observed Recursive Mixed-process Models with cmp

Residual Normality Assumption and the Estimation of Multiple Membership Random Effects Models

Learning Structure in Nested Logit Models

Recursive estimation in large panel data models: Theory and practice