ГОСУДАРСТВЕННЫЙ СТАНДАРТ СОЮЗА ССР
ПРИКЛАДНАЯ О АТИСТИКА
ПРАВИЛА ОПРЕДЕЛЕНИЯ ОЦЕНОК И ДОВЕРИТЕЛЬНЫХ ГРАНИЦ ДЛЯ ПАРАМЕТРОВ ГАММА-РАСПРЕДЕЛЕНИЯ
ГОСТ 11.011-83
Издание официальное
ГОСУДАРСТВЕННЫЙ КОМИТЕТ СССР ПО СТАНДАРТАМ Москва
ГОСУДАРСТВЕННЫЙ СТАНДАРТ СОЮЗА ССР
ПРИКЛАДНАЯ СТАТИСТИКА
ПРАВИЛА ОПРЕДЕЛЕНИЯ ОЦЕНОК И ДОВЕРИТЕЛЬНЫХ ГРАНИЦ ДЛЯ ПАРАМЕТРОВ Г АММА-РАСП РЕД ЕЛ ЕН ИЯ
ГОСТ 11.011-83
Издание официальное
МОСКВА — 1 985
Таблица 4 |
X |
G (х) |
X |
G (х) |
X |
G (х) |
X |
G (*) |
—5,0 |
0,2116 |
—3,2 |
0,3267 |
—1,4 |
0,6407 |
0.4 |
1
1,9650 |
—4,9 |
0,2159 |
—3,1 |
0,3365 |
—1,3 |
0,6725 |
0,5 |
2,1249 |
—4, 8 |
0,2205 |
—3,0 |
0,3468 |
—1,2 |
0,7073 |
0,6 |
2,3000 |
—4,7 |
0,2252 |
-2,9 |
0,3579 |
—1,1 |
0,7444 |
0,7 |
2,4936 |
—4,6 |
0,2302 |
—2,8 |
0,3694 |
—1,0 |
0,7848 |
0,8 |
2,7075 |
—4,5 |
0,2351 |
-2,7 |
0,3817 |
—0,9 |
0,8289 |
0,9 |
2,9340 |
—4,4 |
0,2403 |
—2,6 |
0,3947 |
—0,8 |
0,8768 |
1,0 |
3,2031 |
—4,3 |
0,2458 |
-2,5 |
0,4085 |
—0,7 |
0,9291 |
Ы |
3,4905 |
—4,2 |
0,2516 |
—2,4 |
0,4231 |
—0,6 |
0,9862 |
1,2 |
3,8078 |
—4,1 |
0,2576 |
—2,3 |
0,4389 |
—0,5 |
1,0486 |
1,3 |
4,1582 |
—4,0 |
0,2639 |
—2,2 |
0,4557 |
—0,4 |
1,1168 |
1,4 |
4,5455 |
—3,9 |
0,2705 |
-2,1 |
0,4734 |
-0,3 |
1,1914 |
1,5 |
4,9728 |
—3,8 |
0,2774 |
—2,0 |
0,4926 |
—0,2 |
1,2733 |
1,6 |
5,4448 |
—3,7 |
0,2846 |
-1,9 |
0,5130 |
—0,1 |
1,3630 |
1,7 |
5,9667 |
—3,6 |
0,2922 |
-1,8 |
0,5349 |
0 |
1,4617 |
1,8 |
6,5433 |
—3,5 |
0,3002 |
— 1,7 |
0,5586 |
0,1 |
1,5701 |
1,9 |
7,1795 |
—3,4 |
0,3086 |
— 1,6 |
0,5839 |
0,2 |
1,6893 |
|
|
—3,3 |
0,3174 |
-1,5 |
0,6112 |
0,3 |
1,8205 |
|
|
|
* х + 0,556 + V (*+0,556)4-5,327 2,663 |
4.3. При —10,0^л:<—5,0 используют формулу
При х<—10,0 используют формулу
—х—0,4238
При лг<—75,0 разрешается использовать вместо (18) более простую формулу
а*=---. (19)
X
4.4. При х>1,9 следует использовать формулу
а# = е*+0,5. (20)
4.5. Для построения доверительных границ следует вычислить оценку а* (а*) среднего квадратического отклонения оценки а* параметра а, определенной в соответствии с п. 4.1, по формуле
«•(«*)=-—(21)
У пЦа*)
где /(а*) определяют по а* в соответствии с пп. 4.6—4.8,
ГОСТ 11.011-83 Стр. 9
Таблица 5 |
а |
I (а) |
а |
Па) |
а |
/(а) |
а |
/ (а) |
0,20 |
26,2674 |
1,28 |
1.1587 |
1,72 |
0,7813 |
2,80 |
0,4284 |
0,30 |
12,2454 |
1,30 |
1,1343 |
1,74 |
0,7698 |
2,90 |
0,4110 |
0,40 |
7,2754 |
1,32 |
1,1108 |
1,76 |
0,7585 |
3,00 |
0,3949 |
0,50 |
4,9348 |
1,34 |
1,0882 |
1,78 |
0,7476 |
3,10 |
0,3802 |
0,60 |
3,6361 |
1,36 |
1*0664 |
1,80 |
0,7370 |
3,20 |
0,3664 |
0,70 |
2,8340 |
1,38 |
1,0455 |
1,82 |
0,7266 |
3,30 |
0,3536 |
0,80 |
2,2995 |
1,40 |
1,0254 |
1,84 |
0,7166 |
3,40 |
0,3416- |
0,90 |
I,9226 |
1,42 |
1,0059 |
1,86 |
0,7068 |
3,50 |
0,3304- |
1,00 |
1,6449 |
1,44 |
0,9872 |
1,88 |
0,6973 |
3,60 |
0,3199 |
1,02 |
1,5981 |
1,46 |
0,9691 |
1,90 |
0,6880 |
3,70 |
0,3100 |
1,04 |
1,5537 |
1,48 |
0,9517 |
1,92 |
0,6789 |
3,80 |
0,3008 |
1,06 |
1,5115 |
1,50 |
0,9348 |
1,94 |
0,6701 |
3,90 |
0,2921 |
1,08 |
1,4715 |
1,52 |
0,9185 |
1,96 |
0,6615 |
4,00 |
0,2838 |
1,10 |
1,4333 |
1,54 |
0,9027 |
1,98 |
0,6531 |
4,10 |
0,2761 |
1,12 |
1,3970 |
1,56 |
0,8875 |
2,00 |
0,6449 |
4,20 |
0,2688 |
1,14 |
1,3623 |
1,58 |
0,8727 |
2,10 |
0,6069 |
4,30 |
0,2616 |
1,16 |
1,3292 |
1,60 |
0,8584 |
2,20 |
0,5730 |
4,40 |
0,2551 |
1,18 |
1,2976 |
1,62 |
0,8446 |
2,30 |
0,5426 |
4,50 |
0,2488 |
1,20 |
1,2674 |
1,64 |
0,8321 |
2,40 |
0,5152 |
4,60 |
0,2427 |
1,22 |
1,2385 |
1,66 |
0,8181 |
2,50 |
0,4904 |
4,70 |
0,2370 |
1,24 |
1,2107 |
1,68 |
0,8055 |
2,60 |
0,4678 |
4,80 |
0,2315 |
1,26 |
1,1842 |
1,70 |
0,7932 |
2,70 |
0,4472 |
4,90 |
0,2264 |
|
|
|
|
|
|
5,00 |
0,2213 |
|
4.6. При 0,20=^а*<:5,0 значение / (а*) определяют по а* с помощью табл. 5. Если значение а* попадает между двумя значениями а в табл. 5, то применяют дробно-линейную интерполяцию (в соответствии с разд. 9).
Примечание. Для определения I(а*) при а*, табл. 5, может быть использована рекуррентная формула
/<«)-/(«-О-_1_.
4.7. При а* <0,20 применяют формулу 1(а*)= 1,6449+
4.8. При 5,0<а* <10,0 применяют формулу
а*—0,5 12(д*—0,5)3
При а*>10,0 применяют формулу
т/„*ч 1
4.9. Верхнюю доверительную границу аВ1 соответствующую доверительной вероятности 1 — а, вычисляют по формуле
ав=а*+яав*(а*),
где а* находят в соответствии с п. 4.1, о* (а*)—по п. 4.5, квантиль и а стандартного нормального распределения, соответствующая доверительной вероятности 1 — а, приведена в табл. 3.
4.10. Нижнюю доверительную границу аЯу соответствующую доверительной вероятности 1 — а, вычисляют по формуле
(27)
где а*, «а, а* (а*) — те же, что и в п. 4.9.
Примечание. Если нижняя доверительная граница аН) вычисленная по формуле (27), отрицательна, то полагают ан = 0.
4.11. Двусторонние доверительные границы определяют в соответствии с п. 2.9.
5. ВЫБОР МЕТОДА ОЦЕНИВАНИЯ ПРИ НЕИЗВЕСТНЫХ ПАРАМЕТРАХ ФОРМЫ И МАСШТАБА И ИЗВЕСТНОМ ПАРАМЕТРЕ СДВИГА
5.1. Метод оценивания выбирают в зависимости от величины относительной или абсолютной погрешности наблюдений, объема выборки и оценки параметра формы.
5.2. Оценку а параметра формы а, используемую при выборе метода оценивания, определяют методом моментов по следующему алгоритму:
вычисляют среднее арифметическое х наблюдений хи хъ хп по формуле
вычисляют выборочную дисперсию s2 наблюдений хи х% ..., хп
по формуле
(29)
вычисляют оценку а параметра формы а 'по формуле
гост ii.aii—аз стр. и
5.3. При известной относительной погрешности наблюдений б метод оценивания выбирают следующим образом:
по оценке а параметра формы, определенной в соответствии с п. 5.2, вычисляют величину
A = a( 2а—а); (31)
по относительной погрешности наблюдений б и объему выборки п вычисляют величину
В = ; (32),
ЬУ п
если А>В, то точечные оценки и доверительные границы определяют по методу моментов (разд. 6);
если А^В, то применяют метод максимального правдоподобия (разд. 7).
5.4. При известной абсолютной погрешности А применяют правила п. 5.3 с величиной б, вычисленной по формуле
где х определяют по п. 5.2.
6. ОЦЕНИВАНИЕ ПАРАМЕТРОВ ФОРМЫ И МАСШТАБА МЕТОДОМ МОМЕНТОВ ПРИ ИЗВЕСТНОМ ПАРАМЕТРЕ СДВИГА
6.1. Точечную оценку а параметра формы а определяют в соответствии с п. 5.2.
6.2. Точечную оценку Ь (параметра масштаба b вычисляют то формуле
ь = —, (34)
а
где а: и a определяют в соответствии с п. 5.2.
6.3. Для определения доверительных границ параметра формы
следует вычислить оценку а (а) среднего квадратического отклонения оценки а 'параметра а, определенной в соответствии с п. 6.1, по формуле
Стр. 12 ГОСТ 11.011-83
6.4. Верхнюю доверительную границу ав, соответствующую доверительной вероятности 1 —> а, вычисляют по формуле
ав=а+иао(а),
где а определяют в соответствии с п. 5.2 о (а)—с п. 6.3, квантиль иа стандартного нормального распределения, соответствующая доверительной вероятности 1 — а, приведена в табл. 3.
6.5. Нижнюю доверительную границу аи, соответствующую доверительной вероятности 1 — а, вычисляют по формуле
ап=а—ила(а), (37)
где а, а (а) и иа —те же, что и в п. 6.4.
Примечание. Если нижняя доверительная граница аН) вычисленная по формуле (37), отрицательна, то полагают ап—0.
6.6. Для определения доверительных границ параметра масштаба b следует вычислить оценку а (b) среднего квадратического отклонения оценки b параметра й, определенной в соответствии с п. 6.2. Величину а (й) вычисляют по формуле
где а определяют в соответствии с п. 5.4.
6.7. Верхнюю доверительную границу йв, соответствующую доверительной вероятности 1 —■ а, вычисляют по формуле
Ьв = Ь+иав(Ь), (39)
где й определяют в соответствии с п. 6.2, а(й) —по п. 6.6, квантиль и а стандартного нормального распределения, соответствующая доверительной вероятности 1 — а, приведена в табл. 3.
6.8. Нижнюю доверительную границу йн, соответствующую доверительной вероятности 1 — а, вычисляют по формуле
Ьа=Ь—иаз (Ь),
где Ь, а (й) и иа — те же, что и в п. 6.7.
Примечание. Если нижняя доверительная граница Ьп, вычисленная по формуле (40), отрицательна, то полагают Ьл = 0.
ГОСТ 11.011-83 Стр. 13
6.9. Двусторонние доверительные границы для параметров формы а и масштаба Ь определяют в соответствии с п. 2.9.
7. ОЦЕНИВАНИЕ ПАРАМЕТРОВ ФОРМЫ И МАСШТАБА МЕТОДОМ МАКСИМАЛЬНОГО ПРАВДОПОДОБИЯ ПРИ ИЗВЕСТНОМ ПАРАМЕТРЕ СДВИГА
7.1. Точечную оценку а* параметра формы а определяют по следующему алгоритму:
вычисляют
(4.)
где х — среднее арифметическое наблюдений х, х2, .... хп (х определяют по п. 5.2);
вычисляют оценку а* по формуле
а* = Н(х), (42)
в соответствии с пп. 7.2; 7.3.
7.2. При х^0,25 значение а* определяют по х с помощью табл. 6. Если значение х, вычисленное в соответствии с п. 7.1, попадает между двумя значениями х в табл. 6, то применяют дробно-линейную интерполяцию (в соответствии с разд. 9).
Таблица 6 |
X |
я (х) |
X |
И (х) |
X |
Н (х) |
X |
Н(х) |
0,010 |
50,17 |
0,030 |
*# 41.4
16,83 |
0,050 |
10,16 |
0,070 |
7,31 |
0,011 |
45,63 |
0,031 |
16,29 |
0,051 |
9,97 |
0,071 |
7,21 |
0,012 |
41,84 |
0,032 |
15,79 |
0,052 |
9,78 |
0,072 |
7,11 |
0,013 |
38,63 |
0,033 |
15,32 |
0,053 |
9,60 |
0,073 |
7,01 |
0,014 |
35,88 |
0,034 |
14,87 |
0,054 |
9,42 |
0,074 |
6,92 |
0,015 |
33,51 |
0,035 |
14,45 |
0,055 |
9,25 |
0,075 |
6,83 |
0,016 |
31,42 |
0,036 |
14,05 |
0,056 |
9,09 |
0,076 |
6,74 |
0,017 |
29,59 |
0,037 |
13,68 |
0,057 |
8,94 |
0,077 |
6,66 |
0,018 |
27,94 |
0,038 |
13,32 |
0,058 |
8,78 |
0,078 |
6,57 |
0,019 |
26,49 |
0,039 |
12,99 |
0,059 |
8,64 |
0,079 |
6,49 |
0,020 |
25,17 |
0,040 |
12,66 |
0,060 |
8,50 |
0,080 |
6,41 |
0,021 |
23,98 |
0,041 |
12,36 |
0,061 |
8,36 |
0,081 |
6,34 |
0,022 |
22,90 |
0,042 |
12,07 |
0,062 |
8,23 |
0,082 |
6,26 |
0,023 |
21,91 |
0,043 |
11,79 |
0,063 |
8,10 |
0,083 |
6,19 |
0,024 |
21,00 |
0,044 |
11,53 |
0,064 |
7,98 |
0,084 |
6,11 |
0,025 |
20,17 |
0,045 |
11,28 |
0,065 |
7,86 |
0,085 |
6,04 |
0,026 |
19,40 |
0,046 |
11,03 |
0,066 |
7,74 |
0,086 |
5,98 |
0,027 |
18,68 |
0,047 |
10,80 |
0,067 |
7,63 |
0,087 |
5,91 |
0,028 |
18,02 |
0,048 |
10,58 |
0,068 |
7,52 |
0,088 |
5,84 |
0,029 |
17,41 1 |
1 0,049 |
10,37 1 |
1 0,069 |
7,41 |
0,089 |
5,78 |
|
Продолжение табл. 6 |
X |
И (х) |
X |
Я <*) |
X |
Я (х) |
X |
И (х) |
0,090 |
5,72 |
0,37 |
1,50 |
0,74 |
0,802 |
2,1 |
0,3242 |
0,091 |
5,66 |
0,38 |
1,46 |
0,75 |
0,792 |
2,2 |
0,3115 |
0,092 |
5,60 |
0,39 |
1,43 |
0,76 |
0,783 |
2,3 |
0,2999 |
0,093 |
5,54 |
0,40 |
1,39 |
0,77 |
0,774 |
2,4 |
0,2891 |
0,094 |
5,48 |
0,41 |
1,36 |
0,78 |
0,765 |
2,5 |
0,2791 |
0,095 |
5,42 |
0,42 |
1,33 |
0,79 |
0,757 |
2,6 |
0,2698 |
0,096 |
5,37 |
0,43 |
1,30 |
0,80 |
0,748 |
2,7 |
0,2612 |
0,097 |
5,32 |
0,44 |
1,28 |
0,81 |
0,740 |
2,8 |
0,2532 |
0,098 |
5,26 |
0,45 |
1,25 |
0,82 |
0,732 |
2,9 |
0,2456 |
0,099 |
5,21 |
0,46 |
1,23 |
0,83 |
0,725 |
3,0 |
0,2386 |
0,10 |
5,16 |
0,47 |
1,20 |
0,84 |
0,717 |
3,1 |
0,2321 |
0,11 |
4,71 |
0,48 |
1,18 |
0,85 |
0,710 |
3,2 |
0,2258 |
0,12 |
4,33 |
0,49 |
1,16 |
0,86 |
0,702 |
3,3 |
0,2197 |
0,13 |
4,01 |
0,50 |
1,14 |
0,87 |
0,695 |
3,4 |
0,2141 |
0,14 |
3,73 |
0,51 |
1,12 |
0,88 |
0,688 |
3,5 |
0,2089 |
0,15 |
3,49 |
0,52 |
1,10 |
0,89 |
0,682 |
4,0 |
0,1861 |
0,16 |
3,28 |
0,53 |
1,08 |
0,90 |
0,675 |
4,5 |
0,1680 |
0,17 |
3,10 |
0,54 |
1,06 |
0,91 |
0,668 |
5,0 |
0,1532 |
0,18 |
2,93 |
0,55 |
1,04 |
0,92 |
0,662 |
6,0 |
0,1306 |
0,19 |
2,79 |
0,56 |
1,03 |
0,93 |
0,656 |
7,0 |
0,1141 |
0,20 |
2,66 |
0,57 |
1,01 |
0,94 |
0,650 |
8,0 |
0,1013 |
0,21 |
2,54 |
0,58 |
0,996 |
0,95 |
0,644 |
9,0 |
0,0913 |
0,22 |
2,43 |
0,59 |
0,981 |
0,96 |
0,638 |
10,0 |
0,0831 |
0,23 |
2,33 |
0,60 |
0,966 |
0,97 |
0,632 |
И |
0,0781 |
0,24 |
2,24 |
0,61 |
0,952 |
0,98 |
0,627 |
12 |
0,0721 |
0,25 |
2,15 |
0,62 |
0,938 |
0,99 |
0,621 |
13 |
0,0656 |
0,26 |
2,07 |
0,63 |
0,925 |
1,00 |
0,6157 |
14 |
0,0613 |
0,27 |
2,00 |
0,64 |
0,912 |
1,1 |
0,5666 |
15 |
0,0573 |
0,28 |
1,94 |
0,65 |
0,899 |
1,2 |
0,5254 |
20 |
0,0442 |
0,29 |
1,88 |
0,66 |
0,887 |
1,3 |
0,4902 |
30 |
0,0303 |
0,30 |
1,82 |
0,67 |
0,876 |
1,4 |
0,4598 |
40 |
0,0231 |
0,31 |
1,76 |
0,68 |
0,864 |
1,5 |
0,4332 |
50 |
0,0187 |
0,32 |
1,71 |
0,69 |
0,853 |
1,6 |
0,4098 |
60 |
0,0157 |
0,33 |
1,66 |
0,70 |
0,842 |
1,7 |
0,3889 |
70 |
0,0136 |
0,34 |
1,62 |
0,71 |
0,832 |
1,8 |
0,3702 |
80 |
0,0119 |
0,35 |
1,57 |
0,72 |
0,822 |
1,9 |
0,3534 |
90 |
0,0106 |
0,36 |
1,53 |
0,73 |
0,812 |
2,0 |
0,3381 |
100 |
0,0096 |
|
7.3. При х<0,25 значение а* для х9 не содержащихся в табл. 6, вычисляют по асимптотической формуле
1+л/1+4*
а*=-1-±-. (43)
4х
При л;<0,07 допускается использовать формулу
7.4. Точечную оценку 6* параметра масштаба Ъ вычисляют по формуле
где х определяют в соответствии с п. 5.2, а а* определяют в соответствии с п. 7.1.
7.5. Для определения доверительных границ ‘параметров а и b вычисляют оценки а* (а*) и а* (6*) средних квадратических отклонений оценок а* и Ь* параметров а и Ь, определенных по пп. 7.1 и 7.4 соответственно. Способ вычисления а* (а*) и а* (&*) при 0,20 а* ^5,0 описан в п. 7.6, при а*<0,20 — в п. 7.7, при а*> >5,0 — в п. 7.8.
7.6. При 0,20^а*^5,0 величины о* (а*) и а* (й*) вычисляют по формулам
а* |
(46) |
л[я*/(я*)—1] У |
/(а*) |
(47) |
л[а*/(а*)—1] |
Значения I (а*) определяют по а* с помощью табл. 5. Интерполяцию проводят в соответствии с п. 4.6.
7.7. При а*<0,20 величины а* (а*) и о* (6*) вычисляют по формулам (46), (47), в которых I (а*) определяют по п. 4.7.
7.8. При 5,0<а*^20,0 вычисляют а* (а*) и а* (b*) по формулам:
6(я*—0,5)2
При а*>20,0 вычисляют а* (а*) и а* (6*) по формулам:
Стр. 16 ГОСТ 11.011—83
7.9. Верхнюю доверительную границу ав, соответствующую доверительной вероятности 1 — а, вычисляют по формуле
aQ—а* + tio? *(а*), (52)
где а* вычисляют в соответствии с п. 7.1, а* (а*) определяют по пп. 7.6—7.8 (в зависимости от величины а*), квантиль стандартного нормального распределения иа соответствующая доверительной вероятности 1 — а, приведена в табл. 3.
7.10. Нижнюю доверительную границу ан, соответствующую доверительной вероятности 1 — а, вычисляют по формуле
ан=а*—иае*{а*), (53)
где а*, иа , а* (а*) —те же, что и в п. 7.9.
Примечание. Если нижняя доверительная граница ан, вычисленная по формуле (53), отрицательна, то полагают ан = 0.
7.11. Верхнюю доверительную границу Ьв параметра масштаба 6, соответствующую доверительной вероятности 1 — а, вычисляют по формуле
Ьв=Ь*+иае*{Ь *), (54)
где 6* вычисляют по п. 7.4, а* (&*) — по пп. 7.6—7.8 (в зависимости от а*), квантиль иа стандартного нормального распределения приведена в табл. 3.
7.12. Нижнюю доверительную границу Ьн параметра масштаба 6, соответствующую доверительной вероятности 1 — а, вычисляют по формуле
ЬН=Ь*—tuo*(b*)> (55)
где Ь*, иа , а* (6*)v ■— те же, что ив п. 7.11.
Примечание. Если нижняя доверительная граница 6Н» вычисленная по формуле (55), отрицательна, то полагают 6Н = 0.
7.13. Двусторонние доверительные границы определяют в соответствии с п. 2.9.
8. ОЦЕНИВАНИЕ ПАРАМЕТРОВ ПРИ НЕИЗВЕСТНОМ ПАРАМЕТРЕ СДВИГА
8.1. Стандарт устанавливает правила определения точечных оценок и доверительных границ для следующих случаев: оценивание трех неизвестных параметров (пп. 8.3—8.6); оценивание параметров масштаба и сдвига при известном параметре формы (пп. 8.7—8.10);
ГОСТ 11.011—83 Стр. 17
оценивание параметров формы и сдвига при известном параметре масштаба (пп. 8.11—8.14);
оценивание параметра сдвига при известных параметрах формы и масштаба (пп. 8.15—8.18).
8.2. Для оценивания параметров применяют метод моментов и метод наилучших асимптотически нормальных оценок. Метод моментов является более простым и требует меньшего объема вычислений, чем метод наилучших асимптотически нормальных оценок. На основе оценок метода моментов вычисляют наилучшие асимптотически нормальные оценки и доверительные границы для параметров гамма-распределения. Метод наилучших асимптотически нормальных оценок при большем объеме вычислений приводит к более точным оценкам, чем метод моментов. Указанные свойства оценок используют при выборе метода оценивания. При проведении расчетов с помощью вычислительной техники предпочтительным является метод наилучших асимптотически нормальных оценок как более точный и позволяющий определить доверительные границы.
8.3. При трех неизвестных параметрах а, Ь, с оценки метода моментов вычисляют по формулам:
а=4-^-. |
(56) |
те32 |
|
(57) |
2 sa |
с = х —аЪ, |
(58) |
где среднее арифметическое наблюдений х и выборочную дисперсию s2 определяют по п. 5.2, выборочный третий центральный момент т% вычисляют по формуле
(59)
Примечание. Если оценка параметра сдвига с, вычисленная по формуле
(58), больше минимального из наблюдений хтш, то полагают с=хтin. Если отрицателен третий центральный момент т3, вычисленный по формуле (53), то неверно, что наблюдения имеют гамма-распределение.
8.4. Если а>2,5, то целесообразно определить наилучшие асимптотические нормальные оценки в соответствии о пп. 8.4.1 — 8.4.7.
© Издательство стандартов, 1985
8.4.1. По результатам наблюдений х\, x2i и оценкам метода моментов Ь, с (п. 8.3) вычисляют величины
—■ 2 In (xi — с) j — In b
8.4.2. По оценкам метода моментов а, b, определенным по п. 8.3, и величинам d\, d2, определенным по п. 8.4.1, вычисляют
k,=-W(a)+d„ (62)
k2=l—(a—\)bd2, (63)
где Ч' (а) определяют по п. 8.4.7.
8.4.3. По оценке метода моментов й (и. 8.3) вычисляют
B(a)=A{a)(aI(a)— 1), (65)
в соответствии с пп. 8.4.3.1; 8.4.3.2.
8.4.3.1. При d^5,0 вычисляют I (а) пои. 4.6, затем вычисляют А (а) и В (а) по формулам (64), (65).
8.4.3.2. При а >5,0 вычисляют А (а) и В (а) по формулам:
Л(а)=
В(а)= |
3(2а—1)з (а—1)2
) |
(66) |
3(27—1)2—4 (7—1)2 |
|
(7—1)2 [3 (27—1)2—2^1 |
(67) |
|
3(2а—1)2—4 (а—1)2 |
Примечание. Формулы (66), (67) соответствуют подстановке в формулы (64), (65) значения / (а), вычисленного по формуле (24) (с заменой а на &)* Формулу (25) для определения значения Л (а) применять запрещается, поскольку она дает относительную погрешность 33 % (при больших й).
УДК 519.2:658.562:006.354 ГОСУДАРСТВЕННЫЙ СТАНДАРТ СОЮЗА ССР
Прикладная статистика
ПРАВИЛА ОПРЕДЕЛЕНИЯ ОЦЕНОК И ДОВЕРИТЕЛЬНЫХ ГРАНИЦ ДЛЯ ПАРАМЕТРОВ ГАММА-РАСПРЕДЕЛЕНИЯ
Applied statistics. Regulations for determination of estimates and confidence limits for parameters of gamma distribution
ОКП 0021
Постановлением Государственного комитета СССР по стандартам от 27 июня 1983 г. № 2684 срок введения установлен
с 01.01.85
Настоящий стандарт устанавливает правила определения оценок и доверительных границ для параметров гамма-распределения по совокупности результатов независимых наблюдений, полученных в процессе анализов, испытаний, измерений и т. д., если исследуемые случайные величины подчиняются тамада-распределению. Проверка согласия опытного распределения этих наблюдений с теоретическим гамма-распределением — по ГОСТ 11.006-74.
1. ОСНОВНЫЕ СВОЙСТВА ГАММА-РАСПРЕДЕЛЕНИЯ
1.1. Распределение случайной величины X называется гамма-распределением, если плотность вероятности имеет вид
—!—(а*—с)й"16-аехр(——с1, если х>с,
Г(а) V ' 1 Ь У (1)
0, если х<с.
Плотность вероятности f (х; а, Ь, с) определяется тремя параметрами а, b, с, где а>0, 6>0. При этом а является .параметром формы, b — параметром масштаба, с — параметром сдвига. Множитель — является нормировочным, Г (а) —■ гамма-функция, Г(д).
оо
1(я)=1 xa~l e~xdx. (2)
о
1.2. Частные случаи гамма-раепределения >при определенных значениях параметров имеют специальные названия. При а=1
Переиздание. Февраль 1984 г.
2—812
Стр. 2 ГОСТ 11.011-83
шеем экспоненциальное распределение (ГОСТ 11.005-74). При натуральном а и с=0 гамма-распределение — распределение Эрланга. Если случайная величина X имеет гамма-распределение с параметром формы а таким, что v = 2а — целое число, 6 = 1 и с=0, то 2Х имеет распределение хи-квадрат (%2) с v степенями свободы.
1.3. Гамма-распределение имеет следующие характеристики: математическое ожидание Е (X) =а6 +с.-дисперсию D (Х) — а2 = аЬ2;
коэффициент вариации L_ — ь^ а .
Е(Х) ab+c ’
Е(Х~-Е(Х))* 2
** = VT;
Е(Х — Е(Х)У _б_ .
характеристическую функцию Eeltx=eitc( 1
2. ОБЩИЕ ПОЛОЖЕНИЯ ОПРЕДЕЛЕНИЯ ТОЧЕЧНЫХ ОЦЕНОК И ДОВЕРИТЕЛЬНЫХ ГРАНИЦ
2.1. Точечные оценки и доверительные границы для парамет
ров гамма-распределения определяются по полной выборке объема п из гамма-распределения, т. е. по совокупности результатов наблюдений хи которые рассматриваются как не
зависимые одинаково распределенные случайные величины, имеющие гамма-распределение.
Для получения оценок и доверительных границ используются метод моментов, метод максимального правдоподобия и метод наилучших асимптотически нормальных оценок (полученных из оценок метода моментов одной итерацией при решении по методу Ньютона—Рафсона системы уравнений максимального правдоподобия) .
2.2. Стандарт устанавливает правила определения точечных оценок и доверительных границ как при известном параметре сдвига с, так и при неизвестном параметре сдвига.
Если параметр сдвига известен, то, вычитая из наблюдений известное число с, получаем полную выборку объема п из гамма-раапределения с параметром сдвига с — 0. Поэтому правила настоящего стандарта устанавливаются для случая с = 0 (разд. 3—
7). Правила определения точечных оценок и доверительных границ при 'неизвестном параметре сдвига с приведены в разд. 8.
Примечание. В отдельных случаях допускается использование правил, отличных от указанных в данном стандарте, если они больше соответствуют специфике рассматриваемой области.
2.3. Стандарт устанавливает правила определения точечных оценок параметров а, b гамма-распределения для следующих случаев (с~6):
оценка параметра масштаба при известном параметре формы;
оценка параметра формы при известном параметре масштаба;
оценка обоих неизвестных параметров.
2.4. В случае двух неизвестных параметров (с = 0) стандарт предусматривает использование двух методов — метода максимального правдоподобия и метода моментов. Метод выбирают в соответствии с правилами, установленными в разд. 5. При одном неизвестном параметре (масштаба или формы) используют метод максимального правдоподобия.
2.5. Стандарт устанавливает правила определения доверительных границ для параметров формы и масштаба гамма-распределения в следующих случаях (с=0):
1) доверительные границы параметра масштаба b при известном параметре формы а;
2) доверительные границы параметра формы а при известном параметре масштаба Ь;
3) доверительные границы параметра формы а при неизвестном параметре масштаба Ь;
4) доверительные границы параметра масштаба Ъ при неизвестном параметре а.
2.6. В случаях 2), 3), 4) п. 2.5, а также при неизвестном параметре сдвига доверительные границы определяют на основе асимптотической нормальности соответствующих точечных оценок, поэтому их можно применять лишь при объеме выборки п^Ю.
2.7. Верхние доверительные границы ав, &в, соответствующие доверительной вероятности 1 — ав, вычисляют из условий:
Р(а<<2в)= 1 — ав,
P{b<bQ)^ 1—а,
2.8. Нижние доверительные границы дн, Ьи для параметров а, Ъ соответственно при доверительной вероятности 1 — ан вычисляют из условий:
Р(а> аи)=1 — а„, P(b>hn)^\— ан.
Стр. 4 ГОСТ 11.011-83
2.9. Если сп, Яв и ЬЯу Ьв определены в соотвествии -с пп. 2.7; 2.8, то
P(aH<a<aQ)= 1 — ан —ав> P(K<b<bB)^l — ан — ав,
т. е. доверительные интервалы [а1Ъ ав] и [6Н, Ьв] накрывают истинные значения параметров а, b соответственно с доверительной вероятностью 1 — ан — ав. В частности, если 'ап='<Хв=1а, т. е. односторонние -границы соответствуют доверительной вероятности 1 — а, то двусторонний доверительный интервал соответствует доверительной вероятности 1 — 2 а. Если двусторонний доверительный интервал должен соответствовать доверительной вероятности 1 — р, следует брать в качестве его концов односторонние (нижнюю и верхнюю) доверительные границы, соответствующие доверительной вероятности 1 — р/2.
2.10. Значения доверительной вероятности (1 — а) — по ГОСТ 11.001—73.
Примечание. Наиболее часто используются значения о=0,01; 0,05; 0,1.
2.11. Примеры применения правил стандарта даны в справочном приложении 1. Свойства гамма-распределения и (правила оценивания ряда его характеристик приведены в справочном приложении 2. Теоретические основы стандарта описаны в справочном приложении 3. Условные обозначения, употребляемые в стандарте, содержатся в справочном приложении 4.
3. ОЦЕНИВАНИЕ ПАРАМЕТРА МАСШТАБА ПРИ ИЗВЕСТНЫХ ПАРАМЕТРАХ
ФОРМЫ И СДВИГА
3.1. Точечную оценку 6* параметра b при известном параметре а вычисляют по формуле
(9)
3.2. Нижнюю доверительную границу Ь1Ь соответствующую доверительной вероятности 1 — а, вычисляют по формуле
^n=rt(a, па) Ь
При т = 1-И000 значения гх (а, т) берут из табл. 1 (т = па).
ГОСТ 11.011-83 Стр. 5
Таблица 1 |
т |
г! (а, т) при а |
0,001 |
0,01 |
0,025 |
0,05 |
0,1 |
0,2 |
1 |
0,14 |
0,22 |
0,27 |
1
0,33 |
1
0,43 |
0,62 |
2 |
0,22 |
0,30 |
0,36 |
0,42 |
0,51 |
0,67 |
3 |
0,27 |
0,36 |
0,42 |
0,48 |
0,57 |
0.70 |
4 |
0,31 |
0,40 |
0,46 |
0,52 |
0,60 |
0,73 |
5 |
0,34 |
0,43 |
0,49 |
0,55 |
0,62 |
0,75 |
6 |
0,37 |
0,46 |
0,51 |
0,57 |
0,65 |
0,76 |
7 |
0,39 |
0,48 |
0,54 |
0,59 |
0,66 |
0,77 |
8 |
0,41 |
0,50 |
0,56 |
0,61 |
0,68 |
0,78 |
9 |
0,43 |
0,52 |
0,57 |
0,62 |
0,69 |
0,79 |
10 |
0,44 |
0,53 |
0,59 |
0,64 |
0,70 |
0,80 |
11 |
0,46 |
0,55 |
0,60 |
0,65 |
0,71 |
0,80 |
12 |
0,47 |
0,56 |
0,61 |
0,66 |
0,72 |
0,81 |
13 |
0,48 |
0,57 |
0,62 |
0,67 |
0,73 |
0,82 |
14 |
0,49 |
0,58 |
0,63 |
0,68 |
0,74 |
0,82 |
15 |
0,50 |
0,59 |
0,64 |
0,68 |
0,74 |
0,83 |
16 |
0,51 |
0,60 |
0,65 |
0,69 |
0,75 |
0,83 |
18 |
0,53 |
0,61 |
0,66 |
0,71 |
0,76 |
0,84 |
20 |
0,55 |
0,63 |
0,67 |
0,72 |
0,77 |
0,85 |
25 |
0,58 |
0,66 |
0,70 |
0,74 |
0,79 |
0,86 |
30 |
0,60 |
0,68 |
0,72 |
0,76 |
0,80 |
0,87 |
40 |
0,64 |
0,71 |
0,75 |
0,78 |
0,83 |
0,88 |
50 |
0,67 |
0,74 |
0,77 |
0,80 |
0,84 |
0,89 |
60 |
0,70 |
0,76 |
0,79 |
0,82 |
0,86 |
0,90 |
70 |
0,71 |
0,77 |
0,80 |
0,83 |
0,87 |
0,91 |
80 |
0,73 |
0,78 |
0,81 |
0,84 |
0,87 |
0,91 |
90 |
0,74 |
0,79 |
0,82 |
0,85 |
0,88 |
0,92 |
100 |
0,75 |
0,80 |
0,83 |
0,86 |
0,88 |
0,92 |
150 |
0,79 |
0.84 |
0,86 |
0,88 |
0,90 |
0,93 |
200 |
0,81 |
0,86 |
0,87 |
0,89 |
0,92 |
0,94 |
250 |
0,83 |
0,87 |
0,89 |
0,90 |
0,92 |
0,95 |
300 |
0,84 |
0,88 |
0,90 |
0,91 |
0,93 |
0,95 |
400 |
0,86 |
0,89 |
0,91 |
0,92 |
0,94 |
0,96 |
500 |
0,88 |
0,90 |
0,92 |
0,93 |
0,94 |
0,96 |
600 |
0,89 |
0,91 |
0,92 |
0,94 |
0,95 |
0,97 |
800 |
0,90 |
0,92 |
0,93 |
0,94 |
0,96 |
0,97 |
1000 |
0,91 |
0,93 |
0,94 |
0,95 |
0,96 |
0,97 |
|
3.3. Верхнюю доверительную границу Ьъ, соответствующую доверительной вероятности 1 — а, вычисляют по формуле
Ь»=г2(<х, па) Ь*.
При т= 1 — 1000 значения Гг (а, т) приведены в табл. 2 (т=па).
3.4. Если значение т = па^ 100 не содержится в тэбл. 1, 2, применяют линейную интерполяцию (в соответствии с разд. 9).
Таблица 2 |
т |
г2 (а, т) при а |
0,001 |
0,01 |
0,025 |
0,05 |
0,1 |
0,2 |
1 |
1000 |
100 |
39,52 |
19,5 |
9,48 |
4,48 |
2 |
44,05 |
13,5 |
8,26 |
5,63 |
3,77 |
2,42 |
3 |
15,75 |
6,88 |
4,85 |
3,66 |
2,73 |
1,95 |
4 |
9,33 |
4,85 |
3,67 |
2,93 |
2,29 |
1,74 |
5 |
6,78 |
3,91 |
3,07 |
2,54 |
2,05 |
1,62 |
6 |
5,42 |
3,36 |
2,72 |
2,29 |
1,90 |
1,54 |
7 |
4,60 |
3,00 |
2,48 |
2,13 |
1,80 |
1,48 |
8 |
4,06 |
2,75 |
2,32 |
2,01 |
1,72 |
1,43 |
9 |
3,67 |
2,57 |
2,18 |
1,92 |
1,66 |
1,40 |
10 |
3,38 |
2,42 |
2,09 |
1,83 |
1,61 |
1,37 |
11 |
3,15 |
2,31 |
2,00 |
1,78 |
1,57 |
1,34 |
12 |
2,97 |
2,21 |
1,94 |
1,73 |
1,53 |
1,33 |
13 |
2,82 |
2,13 |
1,88 |
1,69 |
1,50 |
1,31 |
14 |
2,69 |
2,06 |
1,83 |
1,65 |
1,47 |
1,29 |
15 |
2,59 |
2,01 |
1,79 |
1,62 |
1,46 |
1,28 |
16 |
2,50 |
1,96 |
1,75 |
1,59 |
1,44 |
1,27 |
18 |
2,35 |
1,87 |
1,69 |
1,55 |
1,40 |
1,25 |
20 |
2,23 |
1,81 |
1,64 |
1,51 |
1,37 |
1,24 |
25 |
2,03 |
1,68 |
1,55 |
1 ,44 |
1,33 |
1,21 |
30 |
1,89 |
1,60 |
1,48 |
1,39 |
1,29 |
1,18 |
40 |
1,72 |
1,50 |
1,40 |
1,32 |
1,24 |
1,16 |
50 |
1,65 |
1,43 |
1,35 |
1,28 |
1,21 |
1,14 |
60 |
1,56 |
1,38 |
1,31 |
1,25 |
1,19 |
1,12 |
70 |
1,51 |
1,35 |
1,28 |
1,23 |
1,18 |
1,11 |
80 |
1,47 |
1,32 |
1,27 |
1,21 |
1,16 |
1,10 |
90 |
1,44 |
1,30 |
1,25 |
1,20 |
1,15 |
1,09 |
100 |
1,39 |
1,28 |
1,23 |
1,19 |
1,14 |
1,09 |
150 |
1,30 |
1,22 |
1,18 |
1,15 |
1,12 |
1,07 |
200 |
1,26 |
1,19 |
1,15 |
1,13 |
1,10 |
1,06 |
250 |
1,23 |
1,17 |
M4 |
1,11 |
1,09 |
1,06 |
300 |
1,21 |
1,15 |
1,12 |
1,10 |
1,08 |
1,05 |
400 |
1,18 |
1,13 |
1,11 |
1,09 |
1,07 |
1,04 |
500 |
1,16 |
1,11 |
1,09 |
1,08 |
1,06 |
1,04 |
600 |
1,14 |
1,10 |
1,08 |
1,07 |
1,05 |
1,04 |
800 |
1,12 |
1,09 |
1,07 |
1,06 |
1,05 |
1,03 |
1000 |
1,11 |
1,08 |
1,06 |
1,05 |
1,04 |
1,03 |
|
3.5. При т^Ю0 значения г{ («, т) и г2 (а, т) для т, не содержащегося в табл. 1, 2, вычисляют по формулам:
М«, т)=
где и а— квантиль стандартного нормального распределения с математическим ожиданием 0 и дисперсией 1, соответствующая доверительной вероятности 1 — а:
Значения иа приведены в табл. 3.
Таблица 3 |
a |
и
а |
a |
и
а |
0,001 |
3,090 |
0,025 |
1,960 |
0,0025 |
2,807 |
0,05 |
1,645 |
0,005 |
2,576 |
0,1 |
1,282 |
0,01 |
2,326 |
0,2 |
0,842 |
|
4. ОЦЕНИВАНИЕ ПАРАМЕТРА ФОРМЫ ПРИ ИЗВЕСТНЫХ ПАРАМЕТРАХ МАСШТАБА И СДВИГА |
4.1. Точечную оценку а* параметра формы а при известном значении параметра масштаба b вычисляют по формуле
a*=G
При этом сначала вычисляют величину
(16)
Затем то х находят а* в соответствии с пп. 4.2—4.4.
4.2. При —5,0^'я^1,9 значение а* определяют по х с помощью табл. 4, в которой приведены значения функции G (х).
Если значение х, найденное по формуле (16), попадает между двумя значениями х в табл. 4, то при —0,5 применяют линей
ную интерполяцию, а при х<—0,5 — дробно-линейную интерполяцию (в соответствии с разд. 9).