Skip to Content
МетодичкиМатематическая статистикаПриложениеП7. Условное математическое ожидание

П7. Условное математическое ожидание

Условная вероятность

Определение. Условной вероятностью события AA при условии события BB (предполагается, что P(B)>0\mathsf P(B) > 0) называется

P(A∣B)=P(A∩B)/P(B).\mathsf P(A \mid B) = \mathsf P\bigl(A \cap B\bigr)/\mathsf P(B).

Формула полной вероятности. Если события A1,A2,…A_1, A_2, \dots не пересекаются и ⋃Ak=Ω\bigcup A_k = \Omega, то для любого события BB

P(B)=∑kP(B∩Ak)=∑kP(B∣Ak) P(Ak).\mathsf P(B) = \sum_k\mathsf P\bigl(B \cap A_k\bigr) = \sum_k\mathsf P(B \mid A_k)\,\mathsf P(A_k).

По кирпичикам: условная вероятность — это «пересчёт кассы» после того, как стало известно, что BB случилось. Все исходы вне BB выкидываем, оставшиеся вероятности растягиваем, чтобы в сумме снова была единица, — делим на P(B)\mathsf P(B). Формула полной вероятности — разбор по случаям: Ω\Omega разрезали на куски AkA_k, в каждом куске посчитали вероятность BB, сложили с весами кусков.

Смысл по-пацански: «какова вероятность, что пацан с района — спортсмен, если известно, что он ходит в качалку?» — условная вероятность. «Какова доля спортсменов на районе?» — взвешиваем по всем «качалка / не качалка / ушёл в армию» — формула полной вероятности.

Дискретный случай

Пусть случайные величины ξ\xi и η\eta имеют дискретные распределения.

Определение. Условное математическое ожидание ξ\xi при условии события {ω ⁣:η(ω)=yj}\{\omega\colon \eta(\omega) = y_j\} есть

M(ξ∣η=yj)=∑ixi P(ξ=xi∣η=yj).\mathsf M(\xi \mid \eta = y_j) = \sum_ix_i\,\mathsf P(\xi = x_i \mid \eta = y_j).

Определение. Условным математическим ожиданием M(ξ∣η)\mathsf M(\xi \mid \eta) называется случайная величина, принимающая на множествах вида {ω ⁣:η(ω)=yj}\{\omega\colon \eta(\omega) = y_j\} значения M(ξ∣η=yj)\mathsf M(\xi \mid \eta = y_j).

По кирпичикам: M(ξ∣η=yj)\mathsf M(\xi \mid \eta = y_j) — число (среднее ξ\xi среди тех исходов, где η=yj\eta = y_j). А M(ξ∣η)\mathsf M(\xi \mid \eta) — случайная величина: функция от η\eta, «лучший прогноз ξ\xi, если знаешь только η\eta». Её закон — значения M(ξ∣η=yj)\mathsf M(\xi \mid \eta = y_j) с вероятностями P(η=yj)\mathsf P(\eta = y_j).

Условное среднее на двух кубиках

Mξ = 72,00η = 1 (1/36)3,33η = 2 (3/36)4,80η = 3 (5/36)6,29η = 4 (7/36)7,78η = 5 (9/36)9,27η = 6 (11/36)

Бросаем два кубика: ξ — сумма очков, η — наибольшее из двух. Столбики — M(ξ | η = y), посчитанные по всем 36 равновозможным исходам. M(ξ | η) — случайная величина: она принимает значение M(ξ | η = y) с вероятностью P(η = y) (подписи под столбиками). Её среднее — свойство 1 из П7: M(M(ξ | η)) = Σ M(ξ | η = y)P(η = y) = 7,000 = Mξ = 7.

Случай с плотностью

Пусть случайные величины ξ\xi и η\eta имеют совместную плотность p(ξ,η)(x,y)p_{(\xi,\eta)}(x, y) (см. П8).

Определение. Условная плотность ξ\xi при условии η\eta есть

pξ∣η(x,y)={p(ξ,η)(x,y)/pη(y),если pη(y)>0,0,если pη(y)=0.p_{\xi|\eta}(x, y) = \begin{cases}p_{(\xi,\eta)}(x, y)/p_\eta(y), & \text{если } p_\eta(y) > 0,\\ 0, & \text{если } p_\eta(y) = 0.\end{cases}

Определение. Условное математическое ожидание ξ\xi при условии события {ω ⁣:η(ω)=y}\{\omega\colon \eta(\omega) = y\} задаётся формулой

M(ξ∣η=y)=∫−∞∞x pξ∣η(x,y) dx.\mathsf M(\xi \mid \eta = y) = \int_{-\infty}^\infty x\,p_{\xi|\eta}(x, y)\,dx.

Определение. Условное математическое ожидание M(ξ∣η)\mathsf M(\xi \mid \eta) — случайная величина, принимающая на множествах {ω ⁣:η(ω)=y}\{\omega\colon \eta(\omega) = y\} значения M(ξ∣η=y)\mathsf M(\xi \mid \eta = y).

(Определение M(ξ∣η)\mathsf M(\xi \mid \eta) в общем случае — через теорему Радона — Никодима — есть в курсах теории вероятностей, на которые ссылается Лагутин.)

По кирпичикам: у непрерывной η\eta событие {η=y}\{\eta = y\} имеет вероятность ноль, и делить на него по первой формуле нельзя. Выход — делить плотности: «вертикальный срез» совместной плотности на высоте yy, нормированный так, чтобы площадь под ним была 1. В виджете это делается «руками»: облако режется на узкие полосы по η\eta, в каждой берётся среднее ξ\xi. Для двумерного нормального закона с корреляцией ρ\rho (и единичными дисперсиями) получается прямая M(ξ∣η=y)=ρy\mathsf M(\xi \mid \eta = y) = \rho y — линия регрессии (гл. 21, 22).

Условное математическое ожидание: средние по полосам

ηξ

1500 точек (η, ξ) двумерного нормального закона с ρ = 0,60; по горизонтали η, по вертикали ξ. Облако режется на вертикальные полосы, в каждой полосе ξ усредняется — красные точки. Это «эмпирическое» M(ξ | η ∈ полоса). Синяя прямая — точное M(ξ | η = y) = ρy (П7, для нормального закона). Чем уже полосы, тем ближе к определению M(ξ | η = y) — и тем меньше точек в полосе, тем больше дрожат средние.

Смысл по-пацански: хочешь угадать вес пацана, зная его рост, — бери всех пацанов того же роста и усредняй их вес. Для нормального закона эти средние ложатся на прямую — вот откуда в статистике столько прямых линий.

Свойства условного математического ожидания

  1. M(M(ξ∣η))=Mξ\mathsf M\bigl(\mathsf M(\xi \mid \eta)\bigr) = \mathsf M\xi.
  2. Если ξ\xi не зависит от η\eta, то M(ξ∣η)=Mξ\mathsf M(\xi \mid \eta) = \mathsf M\xi (п. н.).

Для любой борелевской функции φ(x)\varphi(x) имеют место равенства

  1. M(φ(η)ξ∣η)=φ(η) M(ξ∣η)\mathsf M(\varphi(\eta)\xi \mid \eta) = \varphi(\eta)\,\mathsf M(\xi \mid \eta) (п. н.),
  2. M[M(ξ∣η)∣φ(η)]=M(ξ∣φ(η))\mathsf M\bigl[\mathsf M(\xi \mid \eta) \mid \varphi(\eta)\bigr] = \mathsf M(\xi \mid \varphi(\eta)) (п. н.).

По кирпичикам:

  • Свойство 1 — «формула полного математического ожидания», то же, что формула полной вероятности, но для средних: усреднил по кускам — усреднил по всему.
  • Свойство 3 — если знаешь η\eta, то φ(η)\varphi(\eta) для тебя константа, её можно вынести.
  • Свойство 4 — «башенное»: прогноз по более бедной информации φ(η)\varphi(\eta) от прогноза по богатой η\eta — это просто прогноз по бедной.

На свойствах 1 и 4 стоит теорема Рао — Блэкуэлла (гл. 10): усреднение оценки по достаточной статистике не меняет среднего, а дисперсию только уменьшает.

Смысл по-пацански: условное среднее — это лучший расклад, который можно сделать, зная часть информации. Свойство 1 говорит: если потом усреднить свои прогнозы по всем возможным подсказкам, в среднем ты не ошибся. Свойство 4: если тебе подсказку огрубили, твой прогноз — это огрублённый прогноз умника, который подсказку знал целиком.

Проверка

Бросают два кубика. Найди M(сумма∣максимум=3)\mathsf M(\text{сумма} \mid \text{максимум} = 3).

Загрузка редактора…

Проверь себя

Тест

Что такое M(ξ | η) — число или случайная величина?

Тест

Если ξ не зависит от η, чему равно M(ξ | η)?
Обновлено