Skip to Content

Таблицы

В конце книги — шесть таблиц: T1–T5 и T7 (таблицы с номером 6 в книге нет). Здесь они не сфотографированы, а посчитаны: каждая клетка вычисляется в браузере, поэтому годятся любые уровни и степени свободы, а не только напечатанные. Числа совпадают с книжными в пределах последнего знака.

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

По кирпичикам: любая таблица здесь отвечает на один из двух вопросов.

  • Прямой: дано значение статистики — какова вероятность хвоста? (T2.)
  • Обратный: дан уровень pp — какое значение отсекает хвост? Это квантиль ypy_p: P(ξ≤yp)=p\mathsf P(\xi \le y_p) = p. (T3, T4, T5, T7.)

Критерий уровня α\alpha отвергает гипотезу, если статистика вылезла за квантиль уровня 1−α1 - \alpha (односторонний) или 1−α/21 - \alpha/2 (двусторонний).

T1. Случайные числа

Фрагменты таблиц случайных чисел: равномерные двузначные числа 00–99 и нормальные N(0,1)\mathcal N(0, 1) с двумя знаками. Ими в книге «вручную» моделируются выборки и случайные разбиения (гл. 2, гл. 24). Откуда такие таблицы берутся и чем их заменили датчики — гл. 2 и 3.

Таблица 1. Случайные числа

1100973253376520135863467354876
2375420480564894742962480524037
3084226895319645093032320902560
4990190252909376707153831131165
5128079997080157361476403236653
6660657471734072768503669736170
7310601080545571824063530342614
8852697760202051656926866574818
9635733213505325470489055357548
10737964575303529647783580834282
11985201776714905686072210940558
12118050543139808277325072568248
13834529963406288980831374670078
14886854020086507584013676667951
15995946734887517649699182608928
16654811767417468509505804776974
17801243563517727080154531822374
18743509981777402772144323600210
19699162680366252291483693687203
20098932050514225685144642756788
−1,75−0,33−1,260,321,530,35−0,96−0,060,42−1,08
−0,290,091,70−1,09−0,44−0,290,25−0,54−1,380,32
−0,930,130,630,901,41−0,88−0,100,230,130,37
−0,45−0,240,071,031,73−0,06−1,49−0,08−2,36−0,99
0,51−0,880,49−1,30−0,270,76−0,360,19−1,080,53

T1: фрагменты таблиц случайных чисел — равномерно распределённые двузначные числа 00–99 и нормально распределённые N(0, 1) с двумя знаками. По ним в книге моделируются выборки (например, в примере 2 гл. 24 разбиение на подгруппы начинается с чисел 10, 09, 73, 25 первой строки и 37, 54 второй). Кнопка «Новая таблица» заполняет ту же сетку датчиком createRng — так делают сейчас вместо печатных таблиц (гл. 2).

T2. Нормальный закон

Таблица хвоста 1−Φ(x)1 - \Phi(x) стандартного нормального закона. Строка — целые и десятые xx, столбец — сотые. Для двустороннего уровня хвост удваивается: P(∣Z∣≥x)=2(1−Φ(x))\mathsf P(|Z| \ge x) = 2(1 - \Phi(x)).

Ходовые квантили, которые полезно помнить наизусть:

pp0,900,950,9750,990,9950,999
xpx_p1,2821,6451,9602,3262,5763,090

Таблица 2. Хвост нормального закона

1 − Φ(1,96) = 0,0250; двусторонний 0,0500
x0,000,010,020,030,040,050,060,070,080,09
0,00,5000,4960,4920,4880,4840,4800,4760,4720,4680,464
0,10,4600,4560,4520,4480,4440,4400,4360,4330,4290,425
0,20,4210,4170,4130,4090,4050,4010,3970,3940,3900,386
0,30,3820,3780,3740,3710,3670,3630,3590,3560,3520,348
0,40,3450,3410,3370,3340,3300,3260,3230,3190,3160,312
0,50,3090,3050,3020,2980,2950,2910,2880,2840,2810,278
0,60,2740,2710,2680,2640,2610,2580,2550,2510,2480,245
0,70,2420,2390,2360,2330,2300,2270,2240,2210,2180,215
0,80,2120,2090,2060,2030,2000,1980,1950,1920,1890,187
0,90,1840,1810,1790,1760,1740,1710,1690,1660,1640,161
1,00,1590,1560,1540,1520,1490,1470,1450,1420,1400,138
1,10,1360,1330,1310,1290,1270,1250,1230,1210,1190,117
1,20,1150,1130,1110,1090,1070,1060,1040,1020,1000,099
1,30,0970,0950,0930,0920,0900,0890,0870,0850,0840,082
1,40,0810,0790,0780,0760,0750,0740,0720,0710,0690,068
1,50,0670,0660,0640,0630,0620,0610,0590,0580,0570,056
1,60,0550,0540,0530,0520,0510,0490,0480,0470,0460,046
1,70,0450,0440,0430,0420,0410,0400,0390,0380,0380,037
1,80,0360,0350,0340,0340,0330,0320,0310,0310,0300,029
1,90,0290,0280,0270,0270,0260,0260,0250,0240,0240,023
2,00,0230,0220,0220,0210,0210,0200,0200,0190,0190,018
2,10,0180,0170,0170,0170,0160,0160,0150,0150,0150,014
2,20,0140,0140,0130,0130,0130,0120,0120,0120,0110,011
2,30,0110,0100,0100,0100,0100,0090,0090,0090,0090,008
2,40,0080,0080,0080,0080,0070,0070,0070,0070,0070,006
2,50,0060,0060,0060,0060,0060,0050,0050,0050,0050,005
2,60,0050,0050,0040,0040,0040,0040,0040,0040,0040,004
2,70,0030,0030,0030,0030,0030,0030,0030,0030,0030,003
2,80,0030,0020,0020,0020,0020,0020,0020,0020,0020,002
2,90,0020,0020,0020,0020,0020,0020,0020,0010,0010,001
3,00,0010,0010,0010,0010,0010,0010,0010,0010,0010,001

T2: y = 1 − Φ(x) — вероятность верхнего «хвоста» стандартного нормального закона. Строка — целые и десятые x, столбец — сотые. Двусторонний уровень вдвое больше: P(|Z| ≥ x) = 2(1 − Φ(x)). Обратная задача — квантиль x_p = Φ⁻¹(p): x₀,₉₅ = 1,645, x₀,₉₇₅ = 1,960, x₀,₉₉₅ = 2,576. Ползунок даёт хвост для любого x, в том числе далеко за краем книжной таблицы.

Смысл по-пацански: 1,96 — это «два сигма» для тех, кто любит точность. Вылезла нормированная статистика за 1,96 по модулю — на уровне 5 % гипотеза идёт лесом.

Проверка

Найди 1−Φ(1,96)1 - \Phi(1{,}96).

T3. Распределение хи-квадрат

Квантили ypy_p закона χk2\chi^2_k с kk степенями свободы: P(χk2≤yp)=p\mathsf P(\chi^2_k \le y_p) = p. Нужны в критерии хи-квадрат (гл. 18), в доверительных интервалах для дисперсии (гл. 11), в критериях Фридмана и Краскела — Уоллиса при больших выборках (гл. 16, 17).

При больших kk хорошо работает приближение Уилсона — Хилферти: yp≈k(1−29k+xp29k)3y_p \approx k\bigl(1 - \frac2{9k} + x_p\sqrt{\frac2{9k}}\bigr)^3, где xpx_p — нормальный квантиль из T2.

Таблица 3. Квантили χ² с k степенями свободы

уровень p:
k = 10, p = 0,950: y_p = 18,307
k0,9000,9500,9750,9900,9950,999
12,713,845,026,637,8810,83
24,615,997,389,2110,6013,82
36,257,819,3511,3412,8416,27
47,789,4911,1413,2814,8618,47
59,2411,0712,8315,0916,7520,52
610,6412,5914,4516,8118,5522,46
712,0214,0716,0118,4820,2824,32
813,3615,5117,5320,0921,9526,12
914,6816,9219,0221,6723,5927,88
1015,9918,3120,4823,2125,1929,59
1218,5521,0323,3426,2228,3032,91
1522,3125,0027,4930,5832,8037,70
2028,4131,4134,1737,5740,0045,31
2534,3837,6540,6544,3146,9352,62
3040,2643,7746,9850,8953,6759,70
4051,855,859,363,766,873,4
6074,479,183,388,492,099,6
120140,2146,6152,2159,0163,6173,6

p-квантиль y_p: P(ξ ≤ y_p) = p. При больших k удобно приближение Уилсона — Хилферти: y_p ≈ k(1 − 2/(9k) + x_p√(2/(9k)))³. В отличие от книжной таблицы, здесь годятся любые p и k: квантиль считается обращением функции распределения.

Проверка

Найди квантиль χ102\chi^2_{10} уровня 0,95 (два знака).

T4. Распределение Стьюдента

Квантили tt-распределения с kk степенями свободы. Закон симметричен, поэтому y1−p=−ypy_{1-p} = -y_p. Нужны в доверительных интервалах для среднего нормальной выборки (гл. 11), в двухвыборочном критерии Стьюдента (гл. 14), в регрессии (гл. 21). При k→∞k \to \infty квантили сходятся к нормальным из T2; уже при k=30k = 30 разница с 1,960 — меньше 0,09.

Таблица 4. Квантили распределения Стьюдента

уровень p:
k = 10, p = 0,950: y_p = 1,812
k0,9000,9500,9750,9900,9950,999
13,086,3112,7131,8263,66318,31
21,892,924,306,969,9222,33
31,642,353,184,545,8410,21
41,532,132,783,754,607,17
51,482,022,573,364,035,89
61,441,942,453,143,715,21
71,411,892,363,003,504,79
81,401,862,312,903,364,50
91,381,832,262,823,254,30
101,371,812,232,763,174,14
121,361,782,182,683,053,93
151,341,752,132,602,953,73
201,331,722,092,532,853,55
251,321,712,062,492,793,45
301,311,702,042,462,753,39
401,301,682,022,422,703,31
601,301,672,002,392,663,23
1201,291,661,982,362,623,16

p-квантиль y_p: P(ξ ≤ y_p) = p. Распределение Стьюдента симметрично: y₁₋ₚ = −y_p; двусторонний интервал уровня 0,95 берёт y₀,₉₇₅. При k → ∞ квантили стремятся к нормальным (1,645; 1,960; 2,326; 2,576). В отличие от книжной таблицы, здесь годятся любые p и k: квантиль считается обращением функции распределения.

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

T5. Распределение Фишера — Снедекора

Квантили FF-распределения с k1k_1 степенями свободы числителя и k2k_2 — знаменателя. Нижние квантили получаются из верхних: v1−p(k1,k2)=1/vp(k2,k1)v_{1-p}(k_1, k_2) = 1/v_p(k_2, k_1). Нужны в сравнении дисперсий (гл. 24), в дисперсионном анализе (гл. 16, 17) и в общей линейной гипотезе (гл. 21).

Таблица 5. Квантили распределения Фишера — Снедекора

уровень p:
k₁ = 10, k₂ = 20, p = 0,950: y_p = 2,348
k₁0,9000,9500,9750,9900,9950,999
12,974,355,878,109,9414,82
22,593,494,465,856,999,95
32,383,103,864,945,828,10
42,252,873,514,435,177,10
52,162,713,294,104,766,46
62,092,603,133,874,476,02
72,042,513,013,704,265,69
82,002,452,913,564,095,44
91,962,392,843,463,965,24
101,942,352,773,373,855,08
121,892,282,683,233,684,82
151,842,202,573,093,504,56
201,792,122,462,943,324,29
251,762,072,402,843,204,12
301,742,042,352,783,124,00
401,711,992,292,693,023,86
601,681,952,222,612,923,70
1201,641,902,162,522,813,54

p-квантиль y_p: P(ξ ≤ y_p) = p. Нижние квантили — из соотношения v₁₋ₚ(k₁, k₂) = 1/v_p(k₂, k₁) (как в книге). Таблица снизу — для текущего k₂. В отличие от книжной таблицы, здесь годятся любые p и k: квантиль считается обращением функции распределения.

T7. Коэффициент ранговой корреляции Спирмена

Односторонние критические значения rpr_p: P(ρS≥rp)=p\mathsf P(\rho_S \ge r_p) = p при гипотезе независимости, kk — размер выборки (гл. 20). Для k≤9k \le 9 виджет пересчитывает таблицу точным перебором всех k!k! перестановок и подсвечивает расхождения с книгой: их ровно два, и оба — пустые клетки, где уровень на самом деле достижим.

Таблица 7. Критические значения коэффициента Спирмена

p0,1000,0500,0250,0100,0050,001
книга, k = 80,5240,6430,7380,8330,8810,952
перебор0,5240,6430,7380,8330,8810,952

T7: односторонние критические значения r_p коэффициента Спирмена: P(ρ_S ≥ r_p) = p при независимости, k — размер выборки. Двусторонний уровень — 2p (закон ρ_S симметричен). Пустая клетка — уровень недостижим при таком k. Для k ≤ 9 снизу — проверка точным перебором всех k! перестановок (функция из гл. 20): совпадает с книгой, только в двух клетках книга пустая, хотя уровень достижим — при k = 4, p = 0,05 и при k = 5, p = 0,01 (P(ρ_S = 1) = 1/24 ≈ 0,042 и 1/120 ≈ 0,008); они выделены цветом. В тексте §7 гл. 20 эта таблица названа T6 — в книге она под номером 7, а таблицы 6 нет.

Замечание. В тексте §7 гл. 20 книги эта таблица названа «T6», но в приложении она идёт под номером 7, а таблицы 6 нет вовсе.

Таблица в одну строчку

Всё, что выше, считается без всяких библиотек: функция распределения — интеграл от плотности, квантиль — обращение функции распределения делением отрезка пополам.

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

Проверь себя

Тест

Двусторонний критерий Стьюдента уровня 0,05 при k = 10 степенях свободы. Какой квантиль брать из T4?

Тест

Как получить нижний квантиль Фишера v₀,₀₅(5, 10), если в таблице только верхние?

Тест

Почему в T7 при k = 4 в книге пустая клетка для p = 0,05, хотя виджет показывает число?
Обновлено