📘 MORISTAT: Advanced Econometrics

Chapter 4: Advanced Topics in Econometrics

Version 1.0  ·  July 2026  ·  naglal@linux-mauritius.com
4.1 Limited Dependent Variable Models

What Are Limited Dependent Variables?

Limited dependent variables are outcomes that are not continuous and unbounded. Common types include:

TypeDescriptionExamples
BinaryTwo outcomesBuy/Don't buy, Employed/Unemployed
MultinomialMultiple unordered outcomesMode of transport (Car/Bus/Train)
OrderedMultiple ordered outcomesSatisfaction (Low/Medium/High)
CensoredPartially observedHours worked (≥0), Expenditure
TruncatedSelectively observedIncome (only above poverty line)
CountNon-negative integersNumber of patents, Doctor visits

Why Can't We Use OLS?

ProblemIssueConsequence
Non-normal errorsBinary outcome → heteroskedastic errorsInefficient estimates
Non-linear relationshipEffect of X on P(Y=1) is non-linearOLS mis-specifies the relationship
Predicted probabilities outside [0,1]OLS can predict negative probabilitiesNonsensical predictions
HeteroskedasticityVar(Y|X) = P(1-P) depends on XInvalid inference
4.2 Binary Choice Models: Logit and Probit

The Binary Choice Framework

We want to model:

P(Y = 1 | X) = F(X'β)

Where:

The Latent Variable Interpretation

Assume there is an unobserved (latent) variable Y*:

Y* = X'β + ε

We observe:

Y = 1 if Y* > 0 Y = 0 if Y* ≤ 0

Then:

P(Y = 1 | X) = P(Y* > 0 | X) = P(ε > -X'β | X) = 1 - F(-X'β) = F(X'β)

The Two Main Models

ModelDistribution of εLink FunctionF(z)
LogitLogisticlog[P/(1-P)] = X'βΛ(z) = e^z/(1+e^z)
ProbitNormalΦ⁻¹(P) = X'βΦ(z)
Logit Model

Link Function:

log[P/(1-P)] = X'β

Interpretation: The log-odds of Y=1 is a linear function of X.

Odds Ratio Interpretation:

Odds = P/(1-P) OR = exp(βⱼ)

A 1-unit increase in Xⱼ multiplies the odds by exp(βⱼ).

Probit Model

Link Function:

Φ⁻¹(P) = X'β

Interpretation: The inverse normal CDF of P is a linear function of X.

Marginal Effects:

∂P/∂Xⱼ = βⱼ × φ(X'β)

Where φ is the standard normal PDF.

Estimation: Maximum Likelihood

The likelihood for a binary choice model is:

L(β) = ∏ᵢ [F(Xᵢ'β)]^(Yᵢ) [1 - F(Xᵢ'β)]^(1-Yᵢ)

The log-likelihood:

ln L(β) = Σ [Yᵢ ln F(Xᵢ'β) + (1-Yᵢ) ln(1 - F(Xᵢ'β))]

Marginal Effects

For Logit:

∂P/∂Xⱼ = βⱼ × Λ(X'β) × (1 - Λ(X'β))

For Probit:

∂P/∂Xⱼ = βⱼ × φ(X'β)

Interpretation Guidelines

MethodInterpretation
Logit coefficientsChange in log-odds for a 1-unit increase in X
Odds ratiosexp(βⱼ) — multiplicative effect on odds
Marginal effectsChange in probability for a 1-unit increase in X
Probit coefficientsChange in the latent variable index

Comparing Logit and Probit

AspectLogitProbit
DistributionLogisticNormal
TailsHeavierLighter
Coefficient Ratioβ_Logit ≈ 1.6 × β_Probitβ_Probit ≈ 0.625 × β_Logit
InterpretationOdds ratiosLatent index
ComputationEasierSlightly harder
When to UseDefault choiceWhen normality is preferred
4.3 Multinomial and Ordered Models

Multinomial Logit (MNL)

When the outcome has more than two unordered categories.

Model:

P(Y = j | X) = exp(X'βⱼ) / Σₖ exp(X'βₖ)

Normalization: β_J = 0 (base category)

Interpretation: Relative risk ratios.

Conditional Logit

When the choices have attributes that vary by alternative.

Model:

P(Y = j | Xⱼ) = exp(Xⱼ'β) / Σₖ exp(Xₖ'β)

Ordered Logit/Probit

When the outcome has ordered categories.

Model:

P(Y ≤ j | X) = F(αⱼ - X'β)

Where:

Interpretation: Coefficients indicate the direction of the effect on the underlying latent variable.

4.4 Tobit Model (Censored Data)

When to Use Tobit

The Tobit model is used when the dependent variable is censored — we observe the value only within a certain range.

Examples:

The Tobit Model

Latent variable:

Y* = X'β + ε

Observed variable:

Y = Y* if Y* > 0 Y = 0 if Y* ≤ 0

The Likelihood Function

For observations with Y > 0:

L(β, σ) = (1/σ) φ((Y - X'β)/σ)

For observations with Y = 0:

L(β, σ) = Φ(-X'β/σ)

Interpretation of Tobit

Two components:

  1. Probability of being uncensored: P(Y > 0 | X) = Φ(X'β/σ)
  2. Expected value conditional on being uncensored: E(Y | Y > 0, X) = X'β + σ × λ(X'β/σ)

Where λ(z) = φ(z)/Φ(z) is the inverse Mills ratio.

Marginal Effects in Tobit

EffectFormula
On all observations∂E(Y|X)/∂Xⱼ = βⱼ × Φ(X'β/σ)
On uncensored observations∂E(Y|Y>0,X)/∂Xⱼ = βⱼ × [1 - λ(z)(z + λ(z))]
On probability of censoring∂P(Y>0|X)/∂Xⱼ = (βⱼ/σ) × φ(X'β/σ)
4.5 Count Data Models

What Is Count Data?

Count data are non-negative integers measuring the number of times an event occurs.

Examples:

The Poisson Model

Assumption: Y | X ~ Poisson(λ), where λ = exp(X'β)

Density:

P(Y = y | X) = exp(-λ) × λ^y / y!

Properties:

Interpretation:

∂E(Y|X)/∂Xⱼ = βⱼ × exp(X'β)

A 1-unit increase in Xⱼ multiplies the expected count by exp(βⱼ).

The Negative Binomial Model

Problem with Poisson: The assumption that Var(Y) = E(Y) is often violated (overdispersion).

Solution: Add a random term to the mean:

λ = exp(X'β + ν)

Where ν is a Gamma-distributed random effect.

Properties:

Zero-Inflated Models

Problem: Many count data have "too many" zeros.

Zero-Inflated Poisson (ZIP):

P(Y = 0 | X) = p + (1-p) × exp(-λ) P(Y = y | X) = (1-p) × exp(-λ) × λ^y / y!

Zero-Inflated Negative Binomial (ZINB): Similar but with Negative Binomial instead of Poisson.

4.6 Duration (Survival) Models

What Is Duration Data?

Duration data measure the time until an event occurs.

Examples:

Survival Functions

Survival function: S(t) = P(T > t) (probability of surviving past time t)

Hazard function: h(t) = lim_{Δt→0} P(t ≤ T < t+Δt | T ≥ t) / Δt

Relationship: h(t) = -d ln S(t) / dt

The Cox Proportional Hazards Model

Model:

h(t | X) = h₀(t) × exp(X'β)

Where:

Interpretation: exp(βⱼ) is the hazard ratio for a 1-unit increase in Xⱼ.

Parametric Survival Models

ModelHazard FunctionDistribution
Exponentialh(t) = λConstant hazard
Weibullh(t) = λp t^{p-1}Monotonic hazard
Log-normalh(t) = λp t^{p-1}Non-monotonic hazard
Log-logistich(t) = λp t^{p-1}/(1+λt^p)Non-monotonic hazard

Kaplan-Meier Estimator

Non-parametric estimate of the survival function:

Ŝ(t) = ∏_{j: tⱼ ≤ t} (1 - dⱼ / nⱼ)

Where:

4.7 Machine Learning in Econometrics

The ML-Econometrics Interface

Machine LearningEconometricsPurpose
PredictionForecastingOut-of-sample prediction
RegularizationVariable selectionShrinkage (Ridge, Lasso)
Cross-validationModel selectionChoosing model complexity
Ensemble methodsCombining modelsImproving predictions
Feature engineeringSpecification searchFinding functional forms
Ridge Regression

Problem: OLS overfits when k is large relative to n.

Solution: Add a penalty term:

β̂_Ridge = argmin_β {Σ(yᵢ - xᵢ'β)² + λΣβⱼ²}

Properties:

Lasso Regression

Penalty:

β̂_Lasso = argmin_β {Σ(yᵢ - xᵢ'β)² + λΣ|βⱼ|}

Properties:

Elastic Net

Combines Ridge and Lasso:

β̂_EN = argmin_β {Σ(yᵢ - xᵢ'β)² + λ[(1-α)Σβⱼ² + αΣ|βⱼ|]}

Post-Lasso / Double ML

  1. Use Lasso for variable selection
  2. Estimate the selected model with OLS (post-Lasso)
  3. Or use Double ML for causal inference

Machine Learning Diagnostics

ConceptDescription
OverfittingModel fits training data too well, performs poorly on test data
Cross-validationSplitting data to evaluate model performance
Bias-Variance TradeoffSimpler models have higher bias, complex models have higher variance
HyperparametersParameters that control model complexity (λ, α)
4.8 Spatial Econometrics

What Is Spatial Econometrics?

Spatial econometrics accounts for spatial dependence and spatial heterogeneity in data.

Spatial Dependence:
Spatial Heterogeneity:

Spatial Weight Matrix

Definition: W is an N×N matrix where w_ij captures the spatial relationship between observation i and j.

Common Specifications:

Spatial Autoregressive (SAR) Model

Model:

Y = ρWY + Xβ + ε

Where:

Interpretation:

Spatial Error Model (SEM)

Model:

Y = Xβ + u u = λWu + ε

Where:

Spatial Durbin Model (SDM)

Y = ρWY + Xβ + WXθ + ε

Where WXθ captures the effect of neighboring X variables.

Moran's I Test

Testing for spatial autocorrelation:

I = (N/ΣᵢΣⱼ w_ij) × [ΣᵢΣⱼ w_ij (Yᵢ - Ȳ)(Yⱼ - Ȳ)] / [Σᵢ(Yᵢ - Ȳ)²]

Interpretation:

4.9 Mixture Models

What Are Mixture Models?

Mixture models assume the data come from multiple subpopulations (mixture components).

Examples:

Finite Mixture Model

Density:

f(y | X) = Σ_{k=1}^{K} πₖ fₖ(y | X, θₖ)

Where:

Latent Class Analysis

P(Y = y | X) = Σ_{k=1}^{K} πₖ P(Y = y | X, Class = k)

Interpretation: Each class has a different relationship between X and Y.

Mixture of Regressions

Y = Σ_{k=1}^{K} πₖ (X'βₖ + εₖ)

Interpretation: Different regression coefficients for each component.

Model Selection for Mixture Models

CriterionDescription
AICAkaike Information Criterion
BICBayesian Information Criterion
ICLIntegrated Completed Likelihood
EntropyMeasures classification quality
4.10 Implementation in Practice

Checklist for Advanced Models

For Binary Choice (Logit/Probit):

For Multinomial Models:

For Count Data:

For Tobit:

Common Mistakes to Avoid

MistakeWhy It's WrongWhat To Do Instead
Interpreting coefficients as marginal effectsCoefficients are not probabilitiesCompute marginal effects
Ignoring overdispersion in count dataInefficient estimatesUse Negative Binomial
Using Tobit without censoringInconsistent estimatesUse OLS if no censoring
Ignoring spatial dependenceInefficient inferenceUse spatial models
Overfitting with complex modelsPoor out-of-sample performanceUse cross-validation
4.11 Quick Reference Cards

Binary Choice Summary

ModelLink FunctionCDFInterpretation
Logitlog[P/(1-P)] = X'βΛ(z) = e^z/(1+e^z)Log-odds, odds ratios
ProbitΦ⁻¹(P) = X'βΦ(z)Latent index

Count Data Summary

ModelMeanVarianceWhen to Use
PoissonλλEquidispersion
Negative Binomialλλ + αλ²Overdispersion
ZIP(1-p)λ(1-p)λ(1+pλ)Excess zeros

Tobit Summary

Y* = X'β + ε Y = max(0, Y*) ∂E(Y|X)/∂X = β × Φ(X'β/σ)

Spatial Models Summary

ModelEquationInterpretation
SARY = ρWY + Xβ + εSpatial spillovers in Y
SEMY = Xβ + u, u = λWu + εSpatial error correlation
SDMY = ρWY + Xβ + WXθ + εSpillovers in X and Y
4.12 Chapter Summary

Key Takeaways

  1. Binary choice models (Logit/Probit) are used for binary outcomes — interpret via marginal effects
  2. Multinomial/Ordered models extend binary choice to multiple categories
  3. Tobit models handle censored data — two-part interpretation
  4. Count data models (Poisson/NB) handle non-negative integers — check for overdispersion
  5. Duration models analyze time-to-event data — use survival functions and hazard ratios
  6. Machine learning methods (Ridge, Lasso) are useful for prediction and variable selection
  7. Spatial models account for spatial dependence — use Moran's I to test for spatial correlation
  8. Mixture models capture heterogeneity across subpopulations

Before Proceeding

4.13 What's Next?
📘 Part 2: Software Implementation

In Part 2, we will implement these models in MORISTAT:

  1. Logit and Probit — Binary choice models
  2. Multinomial Logit — Multiple categories
  3. Tobit Model — Censored outcomes
  4. Poisson and Negative Binomial — Count data
  5. Cox Proportional Hazards — Duration analysis
  6. Ridge and Lasso — Regularized regression
  7. Spatial Models — SAR, SEM, SDM
  8. Mixture Models — Latent class analysis