Session 5.5.1 – Multinomial Probability Distribution

What is it? The multinomial distribution generalizes the binomial distribution to k categories. We perform n independent trials, and in each trial the outcome falls into exactly one category \(1,2,\dots,k\) with fixed probabilities \(p_1,\dots,p_k\), where \(\sum_{i=1}^k p_i=1\). We want the probability that the counts are \(X_1=x_1,\dots,X_k=x_k\) with \(\sum_{i=1}^k x_i=n\).
\[ P(X_1=x_1,\dots,X_k=x_k) = \frac{n!}{x_1!\,x_2!\cdots x_k!}\; p_1^{x_1} p_2^{x_2}\cdots p_k^{x_k}. \]

Example (4 classes): Digital bits labeled E, G, F, P

# Label Probability \(p_i\) Count \(x_i\) \(E[X_i]=n p_i\) \(\mathrm{Var}(X_i)=n p_i(1-p_i)\)
Sums →
Checks: probs sum to 1 → ; counts sum to n →

Key properties

Excel steps (copy/paste)
  1. Put labels in A2:A5 (E,G,F,P), probabilities in B2:B5, counts in C2:C5, and n in C6 (or use =SUM(C2:C5)).
  2. Check sums: B6: =SUM(B2:B5) (should be 1), C6: =SUM(C2:C5) (should equal n)
  3. Multinomial coefficient (ways) in D2: =MULTINOMIAL(C2,C3,C4,C5)
  4. Product of \(p^x\) in E2: =PRODUCT(B2^C2, B3^C3, B4^C4, B5^C5) Log-safe: =EXP(SUM(C2*LN(B2), C3*LN(B3), C4*LN(B4), C5*LN(B5)))
  5. Final probability in F2: =D2*E2
  6. E[Xi] in D3:D6: =C$6*B2 then drag down
  7. Var(Xi) in E3:E6: =C$6*B2*(1-B2) then drag down
Tip: If B6 ≠ 1, scale probs by =B2/$B$6 etc. If C6 ≠ n, change one count or compute “last x = n − (sum others)”.