Статистическая обработка результатов измерений. Практикум. Учебное пособие
.pdf
двух нормальных распределений с известными дисперсиями(известными из предшествующего опыта или найденными теоретически). Подобная задача возникает, например, при сравнении качества изделий, изготовленных на разных установках.
Пусть генеральные совокупности X и Y распределены нормально, по независимым выборкам, объемы которых соответственно равны n и m, извлеченным из этих совокупностей, определены выборочные средние x и y .
Нулевая гипотеза, которую требуется проверить по выборочным средним при заданном уровне значимостиa, запишется следующим образом:
H0 : M (X )= M (Y ),
т. е. генеральные средние (математические ожидания) рассматриваемых совокупностей равны между собой.
Если окажется, что нулевая гипотеза справедлива, т.е. генеральные средние одинаковы, то различие выборочных средних незначимо и объясняется случайными причинами. Например, если физические величины А и В имеют одинаковые истинные размеры, а средние арифметические x и y результатов измерений
этих величин различны, то это различие незначимое. Если нулевая гипотеза отвергнута, т.е. генеральные средние неодинаковы, то различие выборочных средних значимо и не может быть объяснено случайными причинами, а объясняется тем, что сами генеральные средние (математические ожидания) различны. Например, если среднее арифметическое x результатов измерений физической величины А значимо отличается от среднего арифметического y результатов измерений физической величиныВ, то
это означает, что истинные размеры (математические ожидания) этих величин различны.
В качестве критерия проверки нулевой гипотезы принимают случайную величину
|
|
|
|
|
|
|
|
Z = |
X - Y |
|
|
||||
D(X ) D(Y ). |
|||||||
|
|
+ |
|
|
|
||
|
|
|
m |
||||
|
n |
||||||
31
Если нулевая гипотеза выполняется, то критерий Z – нормированная случайная величина, т.е. ее математическое ожидание М(Z)=0 и стандартное отклонение s(Z)=1.
Для проверки гипотезы о различии между средними(математическим ожиданиями) двух нормальных распределений с известными дисперсиями используется режим Двухвыборочный z- тест для средних. В диалоговом окне данного режима(рис. 7) задаются следующие параметры.
Рис. 7
1.Интервал переменной 1 – вводится ссылка на ячейки, содержащие результаты наблюдений величины X. Диапазон данных должен состоять из одного столбца или из одной строки.
2.Интервал переменой 2 – вводится ссылка на ячейки, содержащие результаты наблюдений величины Y. Диапазон данных должен состоять из одного столбца или из одной строки.
3.Гипотетическая средняя разность– вводится число,
равное предполагаемой разности средних(математических ожиданий) изучаемых генеральных совокупностей. Значение 0 указывает на то, что проверяется гипотеза Н0: ax=ay.
4.Дисперсия переменной 1 (известная) – вводится извест-
ное значение дисперсии генеральной совокупности величины X.
5.Дисперсия переменной 2 (известная) – вводится извест-
ное значение дисперсии генеральной совокупности величины Y.
32
6.Флажок Метки – устанавливается в активное состояние, если первая строка (столбец) во входном диапазоне содержит заголовки. Если заголовки отсутствуют, флажок следует деактивизировать, в этом случае будут автоматически созданы стандартные названия для данных выходного диапазона.
7.Альфа – вводится уровень значимости a, равный вероятности возникновения ошибки первого рода.
8.Переключатель Выходной интервал / Новый рабочий лист / Новая рабочая книга. В положении Выходной интервал
активизируется поле, в которое необходимо ввести ссылку на левую верхнюю ячейку выходного диапазона. Размер выходного диапазона будет определен автоматически и на экране появится сообщение в случае возможного наложения выходного диапазона на исходные данные. В положении Новый рабочий лист открывается новый лист, в который начиная с ячейки А1 вставляются результаты анализа. Если необходимо задать имя открываемого рабочего листа, его вводят в поле, расположенное напротив соответствующего положения переключателя. В положении Новая рабочая книга открывается новая книга, на первом листе которой начиная с ячейки А1 вставляются результаты анализа.
Пример расчета
Выборочные данные о диаметре валиков(мм), обработанных на автоматах 1 и 2, приведены в табл. 9, сформированной на рабочем листе Microsoft Excel. По выборке объема n=14 найден средний размер диаметра валиков x =124,49 мм, обработанных на автомате №1 (ячейка D38 содержит формулу =СРЗНАЧ(D24:D37)). По выборке объемаm=9 найден средний размер диаметра валиков y =124,62 мм, обработанных на автома-
те № 2 (ячейка Е38 содержит формулу =СРЗНАЧ(Е24:Е32)).
|
|
|
Таблица 9 |
|
|
|
|
|
C |
D |
E |
23 |
Номер |
Автомат № 1 |
Автомат № 2 |
24 |
1 |
124,55 |
126,73 |
25 |
2 |
123,08 |
124,42 |
26 |
3 |
125,37 |
124,92 |
33
|
|
|
Продолжение табл. 9 |
|
|
|
|
|
C |
D |
E |
27 |
4 |
126,91 |
124,70 |
28 |
5 |
126,80 |
123,45 |
29 |
6 |
127,60 |
123,29 |
30 |
7 |
121,72 |
125,76 |
31 |
8 |
124,65 |
122,45 |
32 |
9 |
126,64 |
125,90 |
33 |
10 |
123,37 |
|
34 |
11 |
123,96 |
|
35 |
12 |
122,46 |
|
36 |
13 |
122,23 |
|
37 |
14 |
123,53 |
|
38 |
Среднее |
124,49 |
124,62 |
39 |
Дисперсия |
3,66 |
1,92 |
|
|
|
|
|
Кроме того, предварительным анализом установлено, что |
||||
размер диаметра валиков, обработанных на каждом автомате, |
|||||
имеет |
нормальный |
закон |
распределения |
с |
дисперсие |
sx2=3,66 мм2 для автомата |
1 и sy2=1,92 мм2 для автомата |
2. |
|||
В ячейках D39 и Е39 содержатся соответственно следующие |
|||||
формулы: =ДИСП(D24:D37) и =ДИСП(Е24:Е32). |
|
|
|||
Можно ли на уровне значимостиa = 0,05 объяснить различие выборочных средних случайными причинами? Или, иными
словами, при уровне значимостиa = 0,05 требуется проверить
гипотезу Н0: ax = ay.
Для решения задачи используем режим работыДвухвыборочный z-тест для средних. Устанавливаем параметры в одноименном диалоговом окне (см. рис. 6).
1.Интервал переменной 1 – $D23:$D37.
2.Интервал переменной 2 – $Е23:$Е32.
3.Гипотетическая средняя разность – 0.
4.Дисперсия переменной 1 (известная) – 3,66.
5.Дисперсия переменной 2 (известная) – 1,92.
6.Метки – устанавливаем флажок в активное состояние.
7.Альфа – 0,05.
8.Выходной интервал – $C$41.
34
Рассчитанные |
в |
данном |
режиме |
показатели |
представлены |
||||
в табл. 10. |
|
|
|
|
|
Таблица 10 |
|||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
D |
|
|
E |
|
41 |
|
Двухвыборочный z-тест для средних |
|
|
|
|
|||
42 |
|
|
|
|
|
|
|
|
|
43 |
|
|
|
|
Автомат № 1 |
Автомат № 2 |
|
||
44 |
|
|
Среднее |
|
124,49 |
|
124,62 |
|
|
|
|
|
|
|
|||||
45 |
|
Известная дисперсия |
3,66 |
|
|
1,92 |
|
||
|
|
|
|
|
|
|
|
||
|
|
Наблюдения |
14 |
|
|
9 |
|
||
46 |
|
|
|
|
|||||
47 |
|
Гипотетическая разность средних |
0 |
|
|
|
|
||
|
|
|
z |
|
-0,19 |
|
|
|
|
48 |
|
|
|
|
|
|
|
||
49 |
|
P(Z £ z) одностороннее |
0,423826596 |
|
|
|
|||
|
|
|
|
||||||
50 |
|
z критическое одностороннее |
1,64 |
|
|
|
|
||
|
|
|
|
|
|||||
51 |
|
P(Z £ z) двухстороннее |
0,847653193 |
|
|
|
|||
|
|
|
|
|
|
|
|
||
52 |
|
z критическое двухстороннее |
1,96 |
|
|
|
|
||
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
Нулевая гипотеза имеет вид Н0: ax = ay, следовательно, альтернативная ей гипотеза – Н1: ax ¹ ay. В этом случае стоим двустороннюю критическую область (рис. 8), исходя из требования, чтобы вероятность попадания критерия в эту область в предположении справедливости нулевой гипотезы была равна принятому уровню значимости a.
Наибольшая мощность критерия(вероятность попадания критерия в критическую область при справедливости конкурирующей гипотезы) достигается тогда, когда «левая» и «правая» критические точки выбраны так, что вероятность попадания критерия в каждый из двух интервалов критической области равна
a 2 : P(Z < zлев.кр )=a 2 , |
P(Z > zправ.кр )= a 2 . |
|
a/2 |
0 |
a/2 |
|
||
-zкр |
zкр |
Z |
|
Рис. 8
Поскольку z – нормированная нормальная случайная величина, критические точки симметричны относительно 0.
35
В результате расчета получено: расчетное значение критерия zр = -0,19; критическая область (-¥; -1,96) È (1,96; +¥) (рис. 9).
a/2 |
-0,19 |
0 |
a/2 |
|
|
Z -1,96 1,96
Рис. 9
Так как расчетное значение критерия не попадает в критическую область, нет оснований отвергнуть нулевую гипотезу, следовательно, различие выборочных средних обусловлено случайными причинами.
Проверка гипотезы о равенстве дисперсий двух нормальных распределений
Дисперсия является важной характеристикой совокупности результатов измерений и характеризует точность выполнения технологических процессов, точность средств измерения и .тп. Пусть даны генеральные совокупности X и Y, распределенные по нормальному закону, и из них извлечены независимые выборки, объемы которых соответственно равныn и m. В математической статистике доказывается, что если гипотеза H0 : s2X = s2Y выполняется, то величина
F= σ2X
σ2Y
имеет F-распределение с k = n-1 и l = m-1 числом степеней свободы. Величину F называют дисперсионным отношением Фишера (или просто статистикой Фишера) и используют в качестве критерии для проверки данной нулевой гипотезы. При этом в качест-
ве альтернативной рассматривается гипотезаН1 : s2X < s2Y, если s2X < s2Y, или гипотеза Н1 s2X > s2Y, если s2X > s2Y.
Для проверки нулевой гипотезы используется режим рабо-
ты Двухвыборочный F-тест для дисперсий. Диалоговое окно режима представлено на рис. 10. Параметры, которые здесь зада-
36
ются, аналогичны режиму работы Двухвыборочный z-тест для средних, только отсутствуют поляДисперсия переменной 1
(известная), Дисперсия переменной 2 (известная) и Гипотетическая средняя разность.
Рис. 10
Пример расчета
Результаты измерения массовой доли влаги сахарного -пе ченья, %, изготовленного из различных партий сырья, представлены в табл. 11, сформированы на листе Microsoft Excel.
|
|
|
Таблица 11 |
|
|
|
|
|
C |
D |
E |
2 |
Номер |
Партия сырья 1 |
Партия сырья 2 |
3 |
1 |
8,6 |
8,7 |
4 |
2 |
8,6 |
9,4 |
5 |
3 |
9,5 |
9,7 |
6 |
4 |
9,0 |
8,6 |
7 |
5 |
8,5 |
9,0 |
8 |
6 |
9,5 |
9,5 |
9 |
7 |
9,6 |
8,1 |
10 |
8 |
8,5 |
9,4 |
11 |
9 |
8,8 |
8,7 |
12 |
10 |
|
8,9 |
13 |
11 |
|
9,5 |
14 |
12 |
|
9,4 |
15 |
13 |
|
9,9 |
|
|
|
|
37
Рассчитаем выборочные дисперсии этих совокупностей с использованием функции ДИСП, получим sX2 = 0,21 %; sY2 = 0,27 %. Проверим, можно ли на уровне значимости a = 0,05
считать статистически незначимым отличие между оценками дисперсий двух совокупностей.
Для решения задачи используем режим работыДвухвыборочный F-тест для дисперсий. Значения параметров в соответствующем диалоговом окне показаны на рис. 11, а рассчитанные данные в табл. 12.
Рис. 11
|
|
|
|
Таблица 12 |
|
|
|
|
|
|
C |
D |
|
E |
20 |
Двухвыборочный |
F-тест для дисперсии |
|
|
21 |
|
|
|
|
|
|
Партия сырья1 |
|
Партия сырья2 |
22 |
|
|||
23 |
Среднее |
8,96 |
|
9,14 |
24 |
Дисперсия |
0,21 |
|
0,27 |
|
Наблюдения |
9 |
|
13 |
25 |
|
|||
|
df |
8 |
|
12 |
26 |
|
|||
27 |
F |
0,80 |
|
|
|
P(F<=f) одностороннее |
0,39 |
|
|
28 |
|
|
||
29 |
F критическое одностороннее |
0,30 |
|
|
Получено расчетное значение F-критерия Fр = 0,80, а критическая область образуется левосторонним интервалом (0; 0,30).
38
Так как Fр не попадает в критическую область, гипотезу о равенстве дисперсий не отвергаем.
Порядок выполнения работы
1. В соответствии с требованиями ГОСТ Р15810-2014 массовая доля влаги заварных пряников должна находиться в переделах от 8,5 до 16,0 %. Результаты измерений данного показателя представлены в табл. 13.
Таблица 13
Номер |
|
|
Влажность изделий, % |
|
|
|||
изде- |
Вари- |
Вари- |
Вари- |
Вари- |
Вари- |
Вари- |
Вари- |
Вари- |
лия |
ант |
ант |
ант |
ант |
ант |
ант |
ант |
ант |
|
1, 9 |
2, 10 |
3, 11 |
4, 12 |
5, 13 |
6, 14 |
7, 15 |
8, 16 |
1 |
11,9 |
10,8 |
12,6 |
13,8 |
13,7 |
14,4 |
9,7 |
12,0 |
2 |
13,6 |
11,0 |
11,5 |
10,3 |
10,1 |
11,1 |
11,4 |
9,8 |
3 |
11,6 |
11,8 |
12,5 |
11,9 |
11,9 |
11,9 |
13,9 |
12,2 |
4 |
12,1 |
11,7 |
14,7 |
13,3 |
15,2 |
11,5 |
14,3 |
10,4 |
5 |
12,9 |
13,4 |
14,6 |
12,2 |
11,7 |
13,1 |
11,8 |
13,2 |
6 |
10,6 |
11,3 |
10,5 |
11,9 |
12,3 |
12,3 |
11,9 |
14,9 |
7 |
10,2 |
11,4 |
9,2 |
14,0 |
10,8 |
11,5 |
13,2 |
12,9 |
8 |
13,3 |
13,0 |
10,7 |
11,0 |
13,1 |
12,7 |
11,2 |
12,0 |
9 |
12,5 |
13,0 |
12,5 |
11,2 |
14,6 |
12,9 |
12,4 |
13,3 |
10 |
13,3 |
11,5 |
11,2 |
13,6 |
10,9 |
10,4 |
13,2 |
13,1 |
11 |
14,9 |
14,0 |
13,9 |
12,4 |
12,3 |
12,8 |
12,3 |
11,0 |
12 |
10,2 |
13,3 |
12,8 |
13,0 |
12,6 |
11,1 |
13,8 |
11,9 |
13 |
11,3 |
11,3 |
11,8 |
11,8 |
11,7 |
13,3 |
12,9 |
11,6 |
14 |
13,9 |
10,2 |
13,0 |
12,2 |
12,4 |
11,5 |
11,6 |
13,3 |
15 |
13,3 |
12,8 |
13,1 |
14,3 |
12,7 |
13,0 |
14,5 |
11,9 |
16 |
13,5 |
12,5 |
13,7 |
12,1 |
11,4 |
13,6 |
11,6 |
12,9 |
17 |
13,0 |
11,9 |
11,8 |
10,7 |
14,7 |
11,6 |
12,4 |
12,0 |
18 |
15,7 |
13,8 |
13,4 |
13,9 |
12,5 |
13,0 |
12,0 |
10,8 |
19 |
13,8 |
12,0 |
10,7 |
13,2 |
13,2 |
12,8 |
12,8 |
11,5 |
20 |
13,2 |
12,9 |
13,6 |
9,0 |
12,9 |
14,1 |
10,2 |
12,4 |
21 |
13,9 |
14,3 |
12,6 |
12,5 |
14,5 |
12,2 |
11,5 |
13,4 |
22 |
12,3 |
13,5 |
12,8 |
14,0 |
12,2 |
11,1 |
10,5 |
9,2 |
23 |
12,4 |
12,7 |
12,3 |
12,5 |
12,6 |
14,1 |
12,4 |
11,0 |
24 |
10,6 |
12,5 |
13,4 |
10,5 |
11,6 |
12,9 |
12,0 |
11,9 |
25 |
12,1 |
11,8 |
11,3 |
13,1 |
12,4 |
13,8 |
12,6 |
12,4 |
26 |
12,6 |
12,6 |
12,3 |
10,8 |
12,6 |
12,5 |
13,2 |
11,4 |
27 |
13,0 |
12,6 |
12,4 |
13,3 |
14,5 |
10,7 |
10,8 |
13,0 |
28 |
13,7 |
11,2 |
9,7 |
11,2 |
12,5 |
12,8 |
11,1 |
11,3 |
39
Окончание табл. 13
Номер |
Вари- |
Вари- |
Вари- |
Вари- |
Вари- |
Вари- |
Вари- |
Вари- |
изде- |
ант |
ант |
ант |
ант |
ант |
ант |
ант |
ант |
лия |
1, 9 |
2, 10 |
3, 11 |
4, 12 |
5, 13 |
6, 14 |
7, 15 |
8, 16 |
29 |
12,9 |
12,4 |
11,7 |
12,4 |
11,2 |
11,2 |
13,6 |
14,3 |
30 |
11,5 |
12,1 |
12,4 |
14,5 |
13,0 |
14,3 |
12,5 |
14,2 |
31 |
11,5 |
11,1 |
13,7 |
13,4 |
10,2 |
13,8 |
12,1 |
10,9 |
32 |
10,0 |
11,3 |
12,9 |
13,4 |
12,2 |
12,2 |
12,6 |
11,9 |
33 |
13,4 |
11,7 |
15,1 |
13,2 |
13,0 |
11,1 |
13,3 |
14,0 |
34 |
11,4 |
11,9 |
10,7 |
13,9 |
13,6 |
12,5 |
9,7 |
12,9 |
35 |
12,9 |
13,8 |
11,2 |
13,7 |
11,6 |
11,4 |
10,2 |
13,9 |
36 |
11,1 |
14,2 |
10,9 |
11,8 |
10,7 |
14,2 |
12,3 |
12,7 |
37 |
13,5 |
14,6 |
12,2 |
11,1 |
11,6 |
13,3 |
13,3 |
12,9 |
38 |
10,5 |
10,6 |
12,1 |
11,5 |
11,5 |
13,4 |
11,3 |
12,2 |
39 |
12,0 |
14,1 |
10,7 |
11,2 |
12,1 |
11,0 |
10,6 |
9,7 |
40 |
13,0 |
11,9 |
12,8 |
9,8 |
13,7 |
10,0 |
11,7 |
11,1 |
41 |
12,2 |
12,5 |
14,4 |
11,7 |
11,3 |
11,4 |
13,0 |
13,8 |
42 |
10,6 |
13,5 |
14,2 |
11,5 |
12,4 |
11,3 |
13,1 |
12,6 |
43 |
12,4 |
11,1 |
11,9 |
13,2 |
12,3 |
13,5 |
11,1 |
13,2 |
44 |
13,8 |
10,9 |
12,8 |
10,2 |
12,4 |
13,5 |
13,5 |
12,4 |
45 |
11,4 |
13,0 |
12,9 |
13,2 |
11,2 |
13,0 |
14,2 |
12,6 |
46 |
11,8 |
12,6 |
14,5 |
9,4 |
12,3 |
12,2 |
13,5 |
13,8 |
47 |
14,3 |
10,1 |
13,3 |
12,7 |
13,2 |
12,9 |
12,1 |
14,3 |
48 |
12,8 |
12,3 |
12,3 |
12,2 |
10,3 |
12,8 |
12,5 |
13,1 |
49 |
12,1 |
13,3 |
13,9 |
14,6 |
10,7 |
13,0 |
12,0 |
11,4 |
50 |
13,2 |
15,1 |
13,9 |
12,5 |
11,7 |
12,0 |
11,0 |
12,4 |
2.Открыть новую рабочую книгу в табличном процессоре Microsoft Excel. На первом листе рабочей книги заполнить результаты измерений.
3.Рассортировать полученную последовательность чисел
(кнопка АЯ
на панели инструментов) и разбить ее на несколько
интервалов (8-9 интервалов), записать исходные данные для расчета в виде таблицы (аналогично данным в табл. 7).
4.Проверить гипотезу о нормальном законе распределения последовательности случайных чисел.
5.Массовая доля фруктового сырья в пастильных изделиях нормируется в стандарте ГОСТ6441-2014. В зефире и пастиле допускается не менее 11 % фруктового сырья. Результаты изме-
40
