Tutorials › CAS Exam 9 - Risk Management for Actuaries › Estimating and Smoothing Transition Matrices

CAS Exam 9 · Credit Risk Modeling and Management · Tutorial 724 of 995

Estimating and Smoothing Transition Matrices

Estimate transition probabilities from cohort counts, smooth sparse cells without violating structural constraints, and choose sampling weights that match the question being measured.

Advanced 12 min read

What You'll Learn

  • Estimate each transition-matrix row from origin-state cohort counts using the multinomial maximum-likelihood method.
  • Handle zero-count cells correctly and distinguish unobserved transitions from structurally impossible ones.
  • Apply Dirichlet prior smoothing and interpret its pseudo-counts and assumptions.
  • Compare raw and smoothed estimates for sparse transition data.
  • Diagnose how exposure-weighted and equally weighted cohorts produce different transition and default estimates.
  • Identify denominator, follow-up, and repeated-observation errors in transition studies.

Why Estimation Choices Matter

A transition matrix summarizes how obligors move between credit states over a stated interval. Its entries may look like simple proportions, but those proportions depend on which obligors enter the study, how their starting and ending ratings are recorded, and how much weight each observation receives. Sparse data add another challenge: a transition not observed in a finite sample receives an estimated probability of zero under the raw frequency method, even when the transition is plausible.

The objective is not to make every estimated cell positive. It is to estimate the transition process consistently with the data and the intended population, while making uncertainty and modeling choices explicit. The methods here concern estimation of a one-period matrix from cohort observations.

Definition: Let the credit states be \(1,\ldots,K\). For a specified observation horizon, \(N_{ij}\) is the number of obligors observed in state \(i\) at the start and state \(j\) at the end. The origin-state cohort size is \(N_i=\sum_{j=1}^K N_{ij}\). The transition probability \(p_{ij}\) is the probability that an obligor starting in state \(i\) is in state \(j\) at the end of the horizon. A transition matrix has entries \(p_{ij}\geq 0\), with \(\sum_{j=1}^K p_{ij}=1\) for every origin state \(i\).

Rows represent starting states and columns represent ending states. All counts in a row must use the same defined horizon and eligible origin cohort. For example, if the horizon is one year, a rating at the start and a rating one year later define a transition observation. A name with no valid ending observation cannot simply be counted as a nondefault: the treatment of withdrawals, mergers, missing ratings, and incomplete follow-up must be specified.

A default state is often modeled as absorbing, meaning an obligor already in default remains in that state. This is a structural restriction, not merely a transition that happened to be absent in the sample. If default is absorbing, the default row is fixed at probability one for default-to-default, and smoothing must not assign probability to transitions out of default.

Raw Cohort Estimates

Suppose each obligor contributes one transition observation to a particular origin-state cohort. Conditional on starting in state \(i\), the destination counts have a multinomial model with probabilities \(p_{i1},\ldots,p_{iK}\). The maximum-likelihood estimate for that row is the observed destination proportion.

$$ \widehat{p}_{ij}=\frac{N_{ij}}{N_i}, \qquad N_i>0. $$

Each estimated row sums to one because \(\sum_j N_{ij}=N_i\). A row with \(N_i=0\) has no empirical estimate: its probabilities cannot be obtained by dividing by zero. It requires additional data, a justified external assumption, or a separately specified prior or model.

Derivation: For one origin state, the likelihood as a function of its destination probabilities is proportional to \(\prod_j p_{ij}^{N_{ij}}\), subject to nonnegative probabilities summing to one. Any destination with a positive count must have positive probability for the likelihood to be positive. Let \(J_+=\{j:N_{ij}>0\}\). Assigning positive probability to a zero-count destination leaves less probability for destinations with positive counts and cannot improve the likelihood, so the maximum sets \(p_{ij}=0\) for \(j\notin J_+\). On the remaining destinations, maximize \(\sum_{j\in J_+}N_{ij}\log p_{ij}\) subject to \(\sum_{j\in J_+}p_{ij}=1\). The Lagrange equations give \(N_{ij}/p_{ij}=\lambda\), so \(p_{ij}=N_{ij}/\lambda\). Summing over \(J_+\) yields \(\lambda=\sum_{j\in J_+}N_{ij}=N_i\). Thus the maximum is \(N_{ij}/N_i\) for positive-count destinations and zero for zero-count destinations. This boundary treatment is necessary because \(\log(p_{ij})\) is not defined at zero.

Worked Example: Estimating a Three-State Matrix

A hypothetical lender follows each obligor for one year from the start-of-year rating. Its states are A, BBB, and D, where D denotes default and is assumed absorbing. The observed transition counts are:

Starting stateEnd AEnd BBBEnd DRow total
A186024
BBB212620
D001010

Method. Apply the multinomial maximum-likelihood estimate separately to each origin-state row. The A row total is \(18+6+0=24\), so its estimates are \(18/24=0.7500\), \(6/24=0.2500\), and \(0/24=0.0000\). The BBB row total is \(2+12+6=20\), giving \(2/20=0.1000\), \(12/20=0.6000\), and \(6/20=0.3000\).

The D row is constrained by the absorbing-default assumption: its estimated probabilities are \(0/10=0\), \(0/10=0\), and \(10/10=1\). Therefore the matrix estimate, in A, BBB, D order for rows and columns, is:

$$ \widehat{P}= \begin{pmatrix} 0.7500&0.2500&0.0000\\ 0.1000&0.6000&0.3000\\ 0.0000&0.0000&1.0000 \end{pmatrix}. $$

Check. The row sums are \(0.7500+0.2500+0.0000=1.0000\), \(0.1000+0.6000+0.3000=1.0000\), and \(0+0+1=1\). The estimated one-year default proportion from BBB is \(6/20=0.3000\).

Common mistake avoided. A zero observed count in the A-to-D cell gives an unsmoothed estimate of \(0/24=0.0000\), not a conclusion that A-to-D default is impossible. The estimate records the sample frequency; it does not establish that the underlying transition probability is exactly zero.

Dirichlet Smoothing for Sparse Rows

When a row is sparse, an analyst may use a Dirichlet prior to stabilize its estimates. For a row with \(K_i\) allowed destinations, assign a common prior parameter \(\alpha>0\) to each allowed probability. Conditional on the data, the posterior parameters are \(N_{ij}+\alpha\). The posterior mean is a smoothed estimate:

$$ \widetilde{p}_{ij} =\frac{N_{ij}+\alpha}{N_i+K_i\alpha}, \qquad j\text{ an allowed destination}. $$

The denominator includes \(K_i\alpha\), not \(K\alpha\) automatically: it is the number of destinations permitted by the model for that row. If a transition is structurally impossible, it is excluded from the allowed set and receives probability zero. For a fixed absorbing default row with only D allowed, the formula gives \((N_{DD}+\alpha)/(N_{DD}+\alpha)=1\).

The formula is the posterior mean under the stated Dirichlet prior, not a universal correction. It can also be interpreted as adding \(\alpha\) prior pseudo-counts to each allowed destination. A larger \(\alpha\) pulls the estimate more strongly toward the prior mean, which is uniform when the same \(\alpha\) is used for every allowed destination. This does not mean the pseudo-counts are additional observed obligors.

A posterior mean should not be confused with a maximum a posteriori estimate. The posterior mean above is valid for every \(\alpha>0\); the mode has a different formula and may lie on the boundary. When reporting smoothing, state the method, the prior parameters, and which destinations are allowed.

Worked Example: Smoothing a Sparse Rating Row

A hypothetical portfolio has 12 one-year observations from A: 9 remain A, 3 move to BBB, and none move to D. Assume all three destinations are possible from A. Use a symmetric Dirichlet prior with \(\alpha=0.5\), and report the posterior-mean estimates to four decimal places.

Method. There are \(K_i=3\) allowed destinations and \(N_i=9+3+0=12\). The smoothed denominator is \(12+3(0.5)=13.5\). The three numerators are \(9+0.5=9.5\), \(3+0.5=3.5\), and \(0+0.5=0.5\). Thus:

$$ \widetilde{p}_{A,A}=\frac{9.5}{13.5}=0.7037,\qquad \widetilde{p}_{A,BBB}=\frac{3.5}{13.5}=0.2593,\qquad \widetilde{p}_{A,D}=\frac{0.5}{13.5}=0.0370. $$

Before rounding, the fractions sum to \((9.5+3.5+0.5)/13.5=13.5/13.5=1\). The displayed rounded estimates also sum to \(0.7037+0.2593+0.0370=1.0000\). The raw maximum-likelihood default estimate was \(0/12=0\); smoothing gives a positive estimate because default was specified as an allowed transition and the prior assigns it positive probability.

Common mistake avoided. Adding \(0.5\) to the three numerators but leaving the denominator at 12 would produce a row summing to \(13.5/12\), not one. Smoothing adds the total pseudo-count \(3(0.5)=1.5\) to the denominator as well.

The reason the posterior mean has this form is that the multinomial likelihood contributes a factor \(p_{ij}^{N_{ij}}\), while the Dirichlet prior contributes a factor proportional to \(p_{ij}^{\alpha-1}\). Their product has powers \(N_{ij}+\alpha-1\), giving a Dirichlet posterior with parameters \(N_{ij}+\alpha\). A Dirichlet distribution with parameters \(\beta_j\) has mean \(\beta_j/\sum_k\beta_k\), which yields the formula above. The prior and the multinomial model are assumptions; smoothing does not remedy a poorly defined cohort or a misspecified set of possible destinations.

Sampling Choices Change the Estimated Matrix

The simple row proportions estimate the transition distribution for the observations that entered that row. Whether this is the desired target depends on the sampling design. Pooling all obligors gives larger cohorts more influence. Giving each cohort equal weight answers a different question: the average transition distribution across cohorts, regardless of their sizes. Neither weighting is automatically correct.

The distinction is important when cohort sizes and transition experience vary. To estimate transitions for a population represented by all observed obligors, pooling counts is generally appropriate if observations are representative and each obligor is counted once in the defined period. If the target is an equally weighted average across years or separately sampled groups, compute group-specific transition proportions and then average them using the intended weights.

Worked Example: Exposure-Weighted Versus Equal-Cohort Default Rates

A hypothetical study observes two starting-state cohorts over the same one-year horizon. In cohort 1, 20 obligors produce 12 same-state transitions, 6 nondefault downgrades, and 2 defaults. In cohort 2, 80 obligors produce 40 same-state transitions, 24 nondefault downgrades, and 16 defaults. Compare the pooled transition row with the equal-weight average of the two cohort rows.

Method: pool counts for an obligor-weighted estimate. The combined row total is \(20+80=100\). The counts by destination are \(12+40=52\) same-state, \(6+24=30\) nondefault downgrade, and \(2+16=18\) default. The pooled row is therefore \((52/100,\;30/100,\;18/100)=(0.5200,\;0.3000,\;0.1800)\). Its default estimate is \(18/100=0.1800\).

Method: average cohort rows equally. Cohort 1 has proportions \((12/20,\;6/20,\;2/20)=(0.6000,\;0.3000,\;0.1000)\). Cohort 2 has proportions \((40/80,\;24/80,\;16/80)=(0.5000,\;0.3000,\;0.2000)\). The equal-weight averages are \(((0.6000+0.5000)/2,\;(0.3000+0.3000)/2,\;(0.1000+0.2000)/2)=(0.5500,\;0.3000,\;0.1500)\). The equal-cohort default estimate is \(0.1500\).

Check. Each result is a valid row: the pooled proportions sum to \(0.5200+0.3000+0.1800=1\), and the equal-weight averages sum to \(0.5500+0.3000+0.1500=1\). The estimates differ because the larger second cohort has a higher default rate and receives 80% of the weight in the pooled estimate, but only 50% in the equal-cohort estimate.

Common mistake avoided. Averaging the two default counts, \((2+16)/2=9\), and dividing by the pooled 100 obligors mixes incompatible weights. For an obligor-weighted target, the correct calculation is \(18/100=18\%\); for an equally weighted cohort target, it is \((2/20+16/80)/2=15\%\). State which target is intended before choosing the calculation.

Data and Model Checks

Before estimating or smoothing, verify how the observations were assembled. In particular, check that each row denominator represents obligors that actually began the period in the corresponding state. A denominator based on end-of-period counts answers a different question. Also confirm that all cohorts use the same rating definitions, horizon, and default definition, or document any changes and handle them deliberately.

  • One observation per intended unit. If an obligor appears in multiple annual cohorts, pooling all issuer-year records weights frequent observations more heavily than a cross-sectional estimate. Repeated records may also be statistically dependent; they are not automatically equivalent to independent obligors.
  • Complete follow-up. Excluding names with missing end ratings can bias estimates if missingness is associated with credit quality. Treating every missing record as a nondefault is also unjustified unless the data definition supports it.
  • Consistent treatment of defaults and exits. Specify whether default is absorbing and how exits, mergers, or other terminal events map to states. Do not smooth into transitions prohibited by the model.
  • Appropriate weights. Cohort weights should reflect the estimand: obligor-weighted, exposure-weighted, or group-weighted estimates are distinct quantities. An exposure-weighted matrix requires exposure weights, not just raw name counts.
  • Sparse cells. Report raw counts as well as smoothed probabilities. A smoothed positive cell may be driven substantially by the prior when the observed count is small.
Key takeaway: Raw transition estimates are row-specific observed proportions. Smoothing can stabilize sparse estimates, but it adds an explicit prior assumption and must respect structural restrictions. Cohort selection and weighting determine the population or average that the resulting matrix describes.

Check Your Understanding

For each problem, identify the denominator and method before calculating.

  1. A one-year cohort begins in state X. The observed destinations are 14 in X, 5 in Y, and 1 in default. Find the raw maximum-likelihood row, in X, Y, default order. (Answer: \((0.7000,\;0.2500,\;0.0500)\))
  2. A row has counts \((7,0,0)\) across three allowed destinations. Use a symmetric Dirichlet prior with \(\alpha=1\) and report the posterior-mean row. (Answer: \((0.8000,\;0.1000,\;0.1000)\))
  3. Two equally weighted cohorts have 40 and 60 obligors, with 4 and 18 defaults, respectively. Calculate the pooled obligor-weighted default rate and the equally weighted cohort default rate. (Answer: pooled \(22/100=22.00\%\); equal-cohort \((10\%+30\%)/2=20.00\%\))
  4. A row has counts \((0,3,9)\) across three allowed destinations. Smooth with \(\alpha=0.5\), and report probabilities to four decimal places. (Answer: \((0.0370,\;0.2593,\;0.7037)\))
  5. Default is absorbing, so the only allowed destination from D is D. There are 6 observed D-to-D records. Under symmetric Dirichlet smoothing with \(\alpha=0.5\) over allowed destinations, what is the estimated D row? (Answer: \((1.0000)\) for the sole allowed destination D; probabilities to all other states remain structurally zero.)
Clinical Biostats