Теоретический материал   Reading materials  

Web-версия учебного курса "Основы математической статистики"

Раздел 3.7.  Дисперсионный анализ. Проверка гипотезы об отсутствии влияния факторов на параметр (продолжение)

На следующем примере покажем, как при небольшом числе опытов можно проследить влияние сразу нескольких факторов. Поставим вопрос: влияет ли на светоотдачу лампы (параметр) конструктивное исполнение горелки (первый фактор) и тепловой режим работы горелки (второй фактор), задаваемый наличием и конструктивным исполнением внешнего стеклянного баллона, окружающего горелку? (Как видно из этого примера, факторы не обязательно имеют численное выражение, важно, чтобы они произвольно устанавливались или контролировались экспериментатором). Испытаны семь горелок мощностью 250 Вт. Эти горелки одинаковы по наполнению ртутью и излучающими добавками, но несколько отличаются по межэлектродным расстояниям, диаметрам и точности центровки электродов. Каждая горелка испытана в трех тепловых режимах: режим 1 - открытая горелка, режим 2 - в стандартном заводском баллоне, режим 3 - в баллоне специальной конструкции, обеспечивающей более равномерное распределение температуры вдоль стенок горелки. Результаты испытаний представлены в таблице 3.6. Числа в таблице - это ток регистрирующего фотоэлемента в микроамперах, т.е. величина, пропорциональная светоотдаче.

Введем, как и раньше обозначения:

- среднее по строке (в нашем случае среднее из трех значений светоотдачи по всем режимам для j-той горелки);

- среднее по столбцу, в нашем случае среднее значение светоотдачи по всем горелкам для i-того режима.

- генеральное среднее всех чисел в таблице

Таблица 3.6 Светоотдача образцов ламп (в условных единицах)

Горелка    j

i      Режим  (v уровней)

r

j

1

2

3

у

1

6.9

7.8

9.5

8.07

р

2

10.1

10.5

13.0

11.2

о

3

9.7

10.5

11.9

10.7

в

4

9.1

10.9

11.6

10.5

н

5

8.5

9.0

10.1

9.2

е

6

9.0

10.0

12.2

10.4

й

7

9.2

9.8

10.4

9.8

8.9

9.78

11.24

9.98

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

Q1 - характеризует разброс между строками, т.е. влияние первого фактора (у нас - конструкции ); Q2 - характеризует разброс между столбцами, т.е. влияние второго фактора (теплового режима). В нашем случае:

Q1 = 3[(8.07 - 9.98)2 + (11.2 - 9.98)2 + ... + (10.4 - 9.98)2] = 20.37 (семь слагаемых),

Q2 = 7[(8.9 - 9.98)2 + (9.78 - 9.98)2 + (11.24 - 9.98)2 = 19.17.

Теперь вычисляем сумму квадратов Q, которая характеризует отклонение всех значений таблицы от генерального среднего (см. 3.12 ):

Q = (6.9 - 9.98)2 + (7.8 - 9.98)2 + (9.8 - 9.98)2 = 42.02 (21 слагаемое).

Можно показать, что, как и в случае однофакторного анализа, разность
Q
3 = Q - Q1 - Q2                                                                                              (3.22)

характеризует отклонение от среднего не из-за влияния факторов, а из-за случайных погрешностей эксперимента (к ним следует отнести и влияние факторов, не контролируемых в данном эксперименте).

Вычисляем Q3: Q3 = 42.02 - 20.37 - 19.16 = 2.49.

Если гипотеза об отсутствии влияния факторов верна, то, как и в предыдущем случае, отношение всех Q к общей дисперсии распределено по закону Χ2, причем число степеней свободы для Q / σ2 очевидно равно rv-1, для Q1 / σ2  -  r-1, для Q2 / σ2 - ν-1, тогда на "долю" Q3 / σ2 остается rν-1-(r-1)-(ν-1) = (r-1)(ν-1) степеней свободы.

По этим значениям находим соответствующие оценки дисперсий:

Если факторы не влияют на результат то все это- оценки одной и той же дисперсии. Если влияют, то в S12 отразилось влияние первого фактора, в S22 - второго. Проверка влияния факторов ведется по критерию Фишера.

, m=1, 2. Если Fmэ>Fq, то m-тый фактор влияет на результат.

Выбрав α = 0.05 имеем:

Для первого фактора при к1 = 6, к2 = 12, Fq = 2.996, F1э = 3.395 / 0.207 = 16.38. Гипотеза о том, что первый фактор (конструкция) не влияет на светоотдачу, отвергается.

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

Для второго фактора при к1 = 2, к2 = 12, Fq = 3.88, F2э = 9.58 / 0.207 = 46.2. Гипотеза о том, что второй фактор (тепловой режим ) не влияет на светоотдачу также отвергается.

Следовательно, тепловой режим заметно влияет на светоотдачу.

Как видно, результат статистической обработки даже небольшой серии измерений приводит к важным практическим результатам, показывает необходимость дальнейшей работы по оптимизации разрабатываемого изделия и позволяет наметить пути достижения стабильных параметров.

Следующий пример иллюстрирует трехфакторный дисперсионный анализ.

Предположим, мы хотим выяснить влияние различных компонент наполнения на выход излучения металлогалоидных ламп. Изготовим 16 образцов горелок, в которых варьируются: 1) дозировка NaJ (первый фактор);

2) давление Ar (второй фактор); 3) дозировка TlJ (третий фактор). Выполним эксперимент по плану, представленному в таблице 3.7. Условные обозначения уровней факторов здесь примут конкретные значения. Для первого фактора (NaJ): а- соответствует 5 мг; b - соответствует 20 мг; с- соответствует 30 мг; d- соответствует 40 мг. Для второго фактора ( давление Ar): А - 10 мм рт.ст.; В - 15 мм рт.ст.; С - 20 мм рт.ст.; D - 25 мм рт.ст. Для третьего фактора (TlJ): α - 2 мг; β - 4 мг; γ - 6 мг; δ - 8 мг. Результат измерения светоотдачи (в относительных единицах) приведен в той же таблице

Таблица 3.7. Светоотдача ламп в условных единицах

Фактор 1

j

Фактор 2

i  n-уровней

n

j

A

B

C

D

у

р

a

α

14

β

17

γ

22

δ

39

23

о

в

b

δ

45

α

11

β

15

γ

18

22.25

н

е

c

γ

19

δ

41

α

12

β

14

21.5

й

d

β

13

γ

23

δ

36

α

10

20.5

22.7

23

21.2

20/25

21.8

Результаты надо понимать, например, так: горелка с 5 мг NaJ, 5 мм рт.ст. Ar и 8 мг TlJ дала 45 ед. светоотдачи. Как и раньше вычислим средние по строкам и по столбцам и занесем эти данные в таблицу. Кроме того, вычислим средние по уровням третьего фактора
(p = α, β, γ, δ).

Находим также генеральное среднее всех значений в таблице и поместим его в правом нижнем углу таблицы.

Оценка влияния факторов опять выполняется с помощью сравнения соответствующих оценок дисперсий, которые вычисляются из соответствующих сумм квадратов:



, характеризующих влияние первого (Q1), второго (Q2) и третьего (Q3) факторов соответственно и "остаточной суммы квадратов".

Q4 = Q - Q1 - Q2 - Q3,                                                                                     (3.23)

где Q - сумма квадратов отклонений всех значений в таблице от генерального среднего.

Число степеней свободы для каждой оценки, характеризующей влияние факторов равно n-1, для Q оно, очевидно, равно n2-1, поэтому для "остаточной" оценки остается n2-1-3(n-1) степеней свободы. Проверка гипотезы об отсутствии влияния фактора с номером m (m = 1, 2, 3) проводится по критерию Фишера, экспериментальное значение которого определяется отношением

, где (m = 1, 2, 3), .

Критическое значение Fq для всех факторов при выбранном уровне значимости одинаково и берется при числе степеней свободы к1= n-1, к2= n2-1-3(n-1).

В нашем примере:

Q1 = 4[(23 - 21.8)2 + (22.25 - 21.8)2 + (21.5 - 21.8)2 + (20.5 - 21.3)2] = 13,

Q2 = 4[(22.7 - 21.8)2 + (23 - 21.8)2 + (21.2 - 21.8)2 + (20.25 - 21.8)2] = 17,

Q3 = 4[(11.7 - 21.8)2 + (14.7 - 21.8)2 + (20.5 - 21.8)2 + (40.25 - 21.8)2] = 1968,

Q = (14 - 21.8)2 + (17 - 21.8)2 + ... + (10 - 21.8)2 = 2047 (16 слагаемых).

Q4 = 2047 - 13 - 17 - 1968 = 49.

к1 = 3, к2 = 6, соответственно:

S12 = Q1 / 3 = 4.33, S22 = Q2 / 3 = 5.66, S32 = Q3 / 3 = 656, S42 = Q4 / 6 = 8.16.

F1э = 0.53, F2э = 0.69, F3э = 80.4; Fq =  4.77 (при к1 = 3, к2 = 6, α = 0.05).

Делаем вывод, что факторы первый и второй не влияют на светоотдачу, т.е. дозировки NaJ и аргона при вариации их в изученных пределах не сказываются на светоотдаче (это, разумеется, не исключает возможности заметить влияние этих факторов, варьируя их в более широких пределах). В то же время дозировка TlJ при изменении ее в диапазоне от 2 до 8 мг заметно влияет на светоотдачу. Следовательно, в дальнейших экспериментах интервал варьирования количества NaJ и буферного газа следует расширить, если необходимо найти, когда они станут влиять на параметр, а характер влияния TlJ на светоотдачу требуется проследить более детально.