Financial Analytics & Risk Modeling

Chapter 3
Statistical and Econometric Methods for Finance

Graduate Program in Business Analytics
Financial Analytics & Risk Modeling   •   Graduate

Statistical and Econometric Methods for Finance

Introduction

Financial analytics converts uncertain economic and market information into estimates, forecasts, decisions, and controls. Statistical methods describe uncertainty and extract patterns from data; econometrics adds economic structure, causal reasoning, and explicit treatment of time dependence. Together, these methods support portfolio construction, credit scoring, valuation, stress testing, market-risk measurement, and policy analysis.

Financial data are unusually challenging. Returns can be heavy-tailed, volatility changes over time, observations are often dependent, firms differ systematically, and economic relationships may change after crises or regulatory reforms. Consequently, a statistically significant relationship is not automatically economically meaningful or causally valid. A sound analysis combines theory, data diagnostics, estimation, validation, and governance.

Learning Objectives

By the end of this chapter, students should be able to:

  1. Select probability distributions appropriate for financial variables.
  2. Estimate parameters and conduct hypothesis tests with correct standard errors.
  3. Construct and interpret cross-sectional, time-series, and panel regressions.
  4. Explain maximum likelihood estimation and its use in volatility and default models.
  5. Diagnose heteroskedasticity, autocorrelation, multicollinearity, endogeneity, nonstationarity, and structural breaks.
  6. Use simulation to quantify risk and communicate model limitations.

1. Probability Models for Financial Data

A random variable X represents an uncertain quantity. Its probability density function, f(x), describes the likelihood of values in a continuous model, while its cumulative distribution function is F(x)=P(X≤x). The expected value and variance are

E[X] = ∫ x f(x) dx, Var(X) = E[(X − μ)²] = σ²

where μ=E[X], σ² is variance, and integration is over the support of X. For a discrete variable, the integral is replaced by a sum. Variance has squared units; standard deviation σ=√(σ²) has the same units as the variable.

A simple return is Rt=(Pt−Pt−1)/Pt−1, where Pt is price. A continuously compounded, or log, return is rt=ln(Pt/Pt−1). Log returns are additive over time, whereas simple returns compound multiplicatively. For small changes, they are approximately equal, but they should not be treated as identical for large movements.

Common distributions

DistributionTypical financial useImportant properties
BernoulliDefault/no defaultBinary outcome with probability p
BinomialNumber of defaultsFixed number of independent binary trials
NormalApproximate measurement errorsSymmetric, thin tails; mean μ, variance σ²
Student’s tHeavy-tailed returns or residualsDegrees of freedom control tail thickness
LognormalPositive prices, losses, or claim sizesLogarithm is normally distributed
PoissonEvent counts in a fixed intervalMean equals variance under the basic model
ExponentialWaiting timesMemoryless duration model

If Z∼N(0,1), standardization is Z=(X−μ)/σ. The normal model is convenient but often inadequate for daily returns because empirical returns may exhibit skewness, excess kurtosis, volatility clustering, and occasional jumps. A Student’s t model can represent thicker tails, but distributional choice should be tested rather than assumed.

For two assets, covariance is

Cov(X,Y) = E[(X−μX)(Y−μY)]

and correlation is ρXY=Cov(X,Y)/(σXσY), provided both standard deviations are positive. Correlation is unit-free and lies between −1 and 1. Portfolio return for weights wi is Rp=∑i=1n wiRi. Its variance is

Var(Rp) = w′Σw

where w is the weight vector and Σ is the covariance matrix. The formula assumes finite second moments; it does not require normally distributed returns. Normality becomes relevant when translating variance into probability statements or parametric value-at-risk.

2. Sampling, Estimation, and Inference

A sample is used to learn about a population. The sample mean is \bar{x}=n-1∑i=1n xi, and the unbiased sample variance is

s² = [1/(n−1)] ∑i=1n (xi−x̄)²

The denominator n−1 corrects, under random sampling, for estimation of the population mean. The standard error of the sample mean is SE(x̄)=s/√n when observations are independent; dependence requires an appropriate long-run variance estimate.

An estimator is unbiased if its expected value equals the parameter, consistent if it converges to the parameter as sample size increases, and efficient if it has comparatively small variance within a specified class. In finance, a biased but lower-variance estimator can sometimes have lower mean squared error. Model selection should therefore consider out-of-sample performance, not only unbiasedness.

A confidence interval for a mean, under independent observations and approximate normality of the standardized mean, is

x̄ ± t1−α/2,n−1 [s/√n]

where α is the significance level and t1−α/2,n−1 is a Student’s t critical value. The interval describes a procedure with long-run coverage; it does not mean that the fixed parameter has a stated probability of lying in one realized interval.

Hypothesis testing

A null hypothesis H0 is tested against an alternative H1. For a scalar estimate \hat{θ}, the Wald statistic is commonly

t = (θ̂−θ₀)/SE(θ̂)

where θ0 is the null value. A small p-value indicates that the observed result would be relatively unusual under H0, conditional on the model and assumptions. It does not measure the probability that the null hypothesis is true, nor does statistical significance establish economic importance.

Multiple testing, data snooping, and large samples can produce apparently significant relationships by chance. Pre-specified hypotheses, economic effect sizes, confidence intervals, false-discovery controls, and out-of-sample validation improve credibility.

3. Regression and Econometric Interpretation

The classical linear regression model is

yi = β0 + β1x1i + ⋯ + βkxki + εi

where yi is the dependent variable, xji are explanatory variables, βj are parameters, and εi is an unobserved disturbance. Ordinary least squares (OLS) chooses coefficients minimizing

∑i=1n (yi−β0−∑j=1k βjxji)²

In matrix form, y=Xβ+ε, and, when X′X is invertible,

β̂OLS = (X′X)−1X′y

The estimator is linear and unbiased under exogeneity, E(ε|X)=0, and has minimum variance among linear unbiased estimators under homoskedastic, uncorrelated disturbances. Normality is not required for unbiasedness, but helps produce exact small-sample t and F distributions.

The coefficient βj is the expected change in y associated with a one-unit increase in xj, holding other included variables constant. This conditional interpretation is not automatically causal. Causal interpretation requires credible identification, such as randomized variation, a valid instrument, a natural experiment, or strong assumptions justified by institutional knowledge.

For asset pricing, a single-factor regression may be written

Ri,t−Rf,t = αi + βi(Rm,t−Rf,t) + εi,t

where Ri is asset return, Rf is the risk-free return, and Rm is market return. βi measures systematic market exposure; αi is abnormal average performance conditional on the model. A nonzero estimated alpha may reflect omitted risk, costs, measurement error, or luck rather than managerial skill.

Maximum likelihood estimation

Maximum likelihood chooses the parameter value making the observed sample most plausible. For observations with density or probability f(yi|θ), the likelihood is

L(θ)=∏i=1n f(yi|θ), ℓ(θ)=ln L(θ)=∑i=1n ln f(yi|θ)

The maximum likelihood estimator is \hat{θ}=argmax_θ ℓ(θ). Taking logarithms improves numerical stability without changing the maximizer. Standard large-sample inference generally requires correct specification or robust variance methods, identification, and regularity conditions.

For a binary default model, logistic regression specifies

P(Yi=1|xi)=pi=1/[1+exp(−xi′β)]

where Yi=1 denotes default. The log-likelihood is

ℓ(β)=∑i=1n [Yi ln(pi)+(1−Yi)ln(1−pi)]

A coefficient is not a direct percentage-point effect. The marginal effect of variable xj is ∂pi/∂xji=pi(1−pi)βj, evaluated at a specified observation or averaged across observations.

4. Panel Data and Dependence

Panel data observe units such as firms or borrowers repeatedly over time. A fixed-effects model is

yit=αi+λt+xit′β+εit

where αi captures time-invariant unit characteristics and λt captures common time shocks. The within transformation subtracts each unit’s time average, eliminating αi:

yit−ȳi=(xit−x̄i)′β+(εit−ε̄i)

Fixed effects are useful when unobserved, time-invariant factors correlate with regressors. Time-invariant regressors cannot be separately estimated in a standard fixed-effects model. Random effects require a stronger condition: unit effects must be uncorrelated with all regressors. Standard errors are commonly clustered by unit, and sometimes by both unit and time, to address within-group dependence.

Panel estimators do not by themselves solve reverse causality, time-varying omitted variables, or measurement error. Difference-in-differences, instrumental variables, event studies, and dynamic panel methods may be appropriate when their identifying assumptions are defensible.

5. Diagnosing Econometric Problems

Heteroskedasticity

Heteroskedasticity means Var(εi|Xi) varies across observations. It does not necessarily bias OLS coefficients under exogeneity, but conventional standard errors become unreliable and OLS loses efficiency. Breusch–Pagan and White tests can detect forms of heteroskedasticity, although tests are sample-sensitive. Use heteroskedasticity-robust standard errors, transform variables, model the variance, or use weighted least squares when weights are justified.

Autocorrelation

Autocorrelation occurs when disturbances are correlated across time, for example Cov(εt,εt−1)≠0. It can make standard errors too small and can signal omitted dynamics. Durbin–Watson is limited in dynamic models; Breusch–Godfrey and residual correlograms are more flexible. Remedies include lagged variables, ARIMA-type errors, Newey–West standard errors, or explicit volatility models. Newey–West inference requires a suitable bandwidth and is not a substitute for correcting a badly specified mean model.

Multicollinearity

Multicollinearity is strong linear association among regressors. It inflates coefficient variances without necessarily biasing estimates. The variance inflation factor is

VIFj=1/(1−Rj²)

where Rj² comes from regressing regressor xj on the other regressors. There is no universal cutoff. Remedies include theory-based variable selection, combining related variables, increasing informative variation, regularization, or focusing on prediction rather than individual coefficient interpretation.

Endogeneity

Endogeneity arises when E(ε|X)≠0, commonly through omitted variables, simultaneity, or measurement error. OLS then generally has biased and inconsistent estimates. Instrumental variables use Z satisfying relevance, Cov(Z,X)≠0, and exclusion, meaning Z affects the outcome only through X after controls. The exclusion restriction is substantive and cannot be verified solely by a statistical test. Weak instruments produce unstable estimates and misleading inference.

Nonstationarity and spurious regression

A weakly stationary process has stable moments and dependence structure over time. A unit-root process, such as yt=yt−1+ut, is nonstationary. Regressing unrelated nonstationary series can yield high R² and significant t-statistics despite no meaningful relationship. Examine plots, autocorrelation, and unit-root tests such as augmented Dickey–Fuller, while recognizing their low power and dependence on specification. Differencing may produce stationarity but can remove long-run information. Cointegration methods model stationary linear combinations of nonstationary variables and should be grounded in economic theory.

Structural breaks

A structural break is a change in parameters or distribution, such as a regulatory reform or crisis. Chow tests apply to a known break date under restrictive conditions; unknown-break procedures such as Bai–Perron search over possible dates. Rolling estimates, recursive residuals, and stability tests are useful diagnostics. Regime-switching models, break indicators, time-varying parameters, and separate pre- and post-break validation are possible responses. A model trained across incompatible regimes may be statistically precise but operationally dangerous.

6. Simulation and Risk Quantification

Simulation propagates uncertainty through a model. In Monte Carlo analysis, one specifies a data-generating process, draws random inputs, computes an outcome, and repeats the procedure. If simulated losses are L(1),…,L(N), an empirical quantile estimates value-at-risk:

VaR_α = inf{l:P(L≤l)≥α}

where α is the confidence level and the loss convention must be stated clearly. Expected shortfall is the conditional average loss beyond the VaR threshold, subject to a precise treatment of quantile ties. Simulation results depend on distributions, dependence assumptions, parameter uncertainty, and the number of replications. Scenario analysis and historical stress tests complement, rather than replace, Monte Carlo simulation.

Variance reduction, antithetic draws, common random numbers, and quasi-random sequences can improve computational efficiency. Confidence intervals for simulation outputs should be reported. Reproducible random seeds, version-controlled code, and independent validation are essential in regulated environments.

Applied Case: Credit Portfolio Stress Testing

Background and problem

An illustrative bank has 10,000 small-business loans. Management wants to estimate one-year losses under ordinary conditions and under a recession. Each loan has exposure at default (EAD), probability of default (PD), and loss given default (LGD). The simple expected-loss relationship is

EL=∑i=1n EADi×PDi×LGDi

This is an average-loss measure, not a complete loss-distribution model.

Analysis and approach

The analytics team first estimates PD using logistic regression with leverage, liquidity, interest coverage, industry, and lagged delinquency. Because firms are observed over several years, firm fixed effects and year effects are considered in exploratory models; however, a production PD model must handle new borrowers and time-varying macroeconomic conditions. Standard errors are clustered by firm in estimation samples.

Residual and calibration diagnostics reveal heteroskedasticity and a break during a pandemic period. The team therefore uses robust inference, includes macroeconomic interactions, performs rolling-origin validation, and estimates separate recession sensitivities. Correlated defaults are simulated through a common macroeconomic factor rather than treating loans as independent.

For each simulation draw, the model generates a macroeconomic state, conditional default probabilities, default indicators, recoveries, and portfolio loss. Results are reported as mean loss, selected quantiles, and expected shortfall. Sensitivity analysis varies recovery rates, unemployment, house prices, and concentration by industry.

Outcome and governance

The bank finds that expected loss is moderate, but tail loss is concentrated in construction and hospitality. A model using independent defaults materially understates tail risk. Management increases sector limits, raises capital buffers, and introduces early-warning triggers. Model risk governance requires documented assumptions, challenger models, back-testing, monitoring for drift, approval of material changes, and human review of unusual outputs.

Learning points

The case illustrates that estimation, diagnostics, dependence modeling, and scenario design are inseparable. A sophisticated algorithm cannot compensate for biased data, unstable regimes, weak validation, or unclear loss definitions. Ethical practice also requires checking whether proxy variables produce discriminatory outcomes and ensuring that automated credit decisions remain explainable, reviewable, and compliant with applicable law.

Limitations, Technology, and Best Practice

Statistical models simplify reality. Parameter uncertainty, model uncertainty, regime change, data revisions, survivorship bias, look-ahead bias, and transaction costs can invalidate impressive back-tests. Machine-learning methods may improve prediction but can reduce interpretability and do not eliminate endogeneity or distribution shift. Explainable AI tools are aids to diagnosis, not proof of causal validity.

Best practice includes:

  • begin with an economic mechanism and a clearly defined target;
  • inspect data lineage, missingness, outliers, revisions, and measurement units;
  • separate training, validation, and genuinely out-of-sample test periods;
  • use economically justified transformations and avoid leakage;
  • report effect sizes, uncertainty intervals, diagnostics, and sensitivity analyses;
  • select robust or clustered inference appropriate to dependence;
  • test calibration and stability across time, sectors, and demographic groups;
  • document assumptions, code, versions, and decisions;
  • monitor live performance and establish escalation thresholds;
  • maintain human accountability, privacy protection, fairness review, and auditability.

Chapter Summary

Financial econometrics combines probability, estimation, inference, regression, likelihood methods, panel techniques, and simulation. Distributional assumptions influence tail-risk estimates; OLS provides a useful baseline but depends on exogeneity and appropriate inference. Maximum likelihood supports binary-response and volatility models. Fixed effects control for time-invariant unit heterogeneity, while causal claims require stronger identification. Heteroskedasticity, autocorrelation, multicollinearity, endogeneity, nonstationarity, and structural breaks require distinct diagnostics and remedies. Simulation extends models into portfolio-loss distributions, but its conclusions remain conditional on assumptions. Reliable financial analytics therefore requires theory, careful data management, robust validation, governance, and transparent communication.

Key Terms

Probability distribution: A mathematical description of possible values and their probabilities.

Expected value: The probability-weighted average of a random variable.

Variance: The expected squared deviation from the mean.

Covariance: A measure of joint variation between two variables.

Correlation: Standardized covariance ranging from −1 to 1.

Simple return: Percentage change in an asset price.

Log return: Natural logarithm of the price ratio.

Estimator: A rule for using sample data to estimate a population parameter.

Standard error: Estimated sampling variability of an estimator.

Confidence interval: An interval generated by a procedure with specified long-run coverage.

p-value: Probability, under the null model, of a result at least as extreme as observed.

Exogeneity: A condition such as E(ε|X)=0 linking regressors and disturbances.

Maximum likelihood: Parameter estimation by maximizing the likelihood of observed data.

Fixed effects: Unit-specific effects that may correlate with regressors and remain constant over time.

Heteroskedasticity: Nonconstant conditional error variance.

Autocorrelation: Dependence among errors across observations or time.

Endogeneity: Correlation between an explanatory variable and the disturbance.

Nonstationarity: Time variation in a process’s probabilistic properties.

Structural break: A change in a model’s parameters or data-generating process.

Value-at-risk: A specified quantile of a loss distribution.

Review Questions

Short Answer Questions

  1. Distinguish simple returns from log returns.
  2. What does the condition E(ε|X)=0 imply for OLS interpretation?
  3. Why can heteroskedasticity make conventional standard errors unreliable?
  4. What is removed by the within transformation in a fixed-effects model?
  5. Explain why a high R² can occur in a spurious regression.

Long Answer Questions

  1. Explain the assumptions, interpretation, and limitations of the classical linear regression model in an asset-pricing application.
  2. Compare fixed-effects and random-effects panel models, including their identifying assumptions.
  3. Describe maximum likelihood estimation and derive the log-likelihood for logistic regression.
  4. Discuss how nonstationarity and structural breaks affect financial forecasting and inference.
  5. Explain how a bank should combine statistical estimation, Monte Carlo simulation, stress testing, and model governance.

Analytical / Application Questions

  1. A portfolio contains two assets with weights 0.60 and 0.40, standard deviations 10% and 20%, and correlation 0.25. Derive and calculate the portfolio standard deviation.
  2. A default model has estimated probability 0.08 and coefficient 1.5 for a continuous risk variable. Calculate the marginal effect at that observation and explain why it is not 1.5 percentage points.
  3. A regression of returns on market returns shows significant autocorrelation and a break after a regulatory reform. Design a revised estimation and validation strategy, explaining the role of robust inference, break analysis, and out-of-sample testing.

Multiple Choice Questions with Answers

Question 1

Which condition is central to unbiased OLS estimation in the standard linear model?

A. E(ε|X)=0

B. All regressors are normally distributed

C. The dependent variable has zero variance

D. The sample contains no time dimension

Correct Answer: A. E(ε|X)=0

Question 2

What does the fixed-effects transformation primarily eliminate?

A. All time-varying shocks

B. Time-invariant unit-specific effects

C. Measurement error in the dependent variable

D. Random sampling variation

Correct Answer: B. Time-invariant unit-specific effects

Question 3

If portfolio weights are w and the covariance matrix is Σ, portfolio variance is:

A. w′Σw

B. w′Σ^{-1}w

C. Σw′w^{-1}

D. ∑w_iσ_i

Correct Answer: A. w′Σw

Question 4

What is a common consequence of heteroskedasticity for conventional OLS inference?

A. OLS coefficients are always biased

B. Conventional standard errors may be incorrect

C. The sample mean becomes zero

D. All regressors become independent

Correct Answer: B. Conventional standard errors may be incorrect

Question 5

In logistic regression, p=1/[1+exp(−x′β)] represents:

A. A conditional probability between zero and one

B. A residual variance

C. A covariance matrix

D. A guaranteed causal effect

Correct Answer: A. A conditional probability between zero and one

Question 6

Which issue most directly causes OLS inconsistency through correlation between a regressor and the error term?

A. Endogeneity

B. Multicollinearity

C. Large sample size

D. Unit conversion

Correct Answer: A. Endogeneity

Question 7

A process with a unit root is generally described as:

A. Strictly homoskedastic

B. Nonstationary

C. Normally distributed by definition

D. Free from serial dependence

Correct Answer: B. Nonstationary

Question 8

What does value-at-risk summarize?

A. The mean of every possible loss

B. A specified quantile of a loss distribution

C. The probability that a model is correct

D. The covariance between two assets

Correct Answer: B. A specified quantile of a loss distribution

Question 9

The likelihood function is maximized in order to:

A. Choose parameters that make observed data most plausible

B. Force every residual to equal zero

C. Remove all omitted variables

D. Guarantee causal identification

Correct Answer: A. Choose parameters that make observed data most plausible

Question 10

A high variance inflation factor is primarily evidence of:

A. Autocorrelation

B. Multicollinearity

C. A structural break

D. A unit root in the dependent variable

Correct Answer: B. Multicollinearity

Financial Analytics & Risk Modeling
Chapter 3 of 14