Практическое занятие №11. Трехфакторный дисперсионный анализ

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

Поскольку факторов три, модель комплекса, где строчки – градации одного фактора, а столбцы – другого не подходит. При построении комплекса по трём факторам эти факторы приходится располагать слева направо (как это принято в иерархическом комплексе).

В табл. 11.1. предлагается упрощенная общая схема трехфакторного дисперсионного комплекса при следующих заданных условиях:

1) число градаций по фактору А – 2 (А 1 и А 2);

2) число градаций по фактору В – 2 (В 1 и В 2);

3) число градаций по фактору С – 4 (С 1, С 2, С 3 и С 4);

4) число наблюдений по сочетанию градаций факторов А, В и С, то есть, nijk . – 2.

Все факторы имеют фиксированные градации. Комплекс равномерный.

После построения дисперсионного комплекса необходимо определить основные параметры, суммы и компоненты формул:

основные параметры

объем комплекса (N) = 64

число градаций по фактору А (а) = 2

число градаций по фактору В (b) = 2

число градаций по фактору С (с) = 4

число наблюдений по градации фактора А (ni ...)=32

число наблюдений по градации фактора В (n . j ..) =32

число наблюдений по градации фактора С (n .. k .) = 16

число наблюдений по сочетанию градаций АВ (nij ..) = 16

число наблюдений по сочетанию градаций АC (ni . k .) = 8

число наблюдений по сочетанию градаций ВC (n . jk .) = 8

число наблюдений по сочетанию градаций AВC (nijk .) = 4

отдельные наблюдения - xijkl

основные суммы

    сумма квадратов всех наблюдений;

     сумма всех наблюдений;

 

 


 

Таблица 11.1. Трехфакторный равномерный дисперсионный комплекс (фактор А – 2 градации; фактор В – 2 градации; фактор С – 4 градации; nijk. – 2)

 

А

В

С

xijkl

Σi…

ni…

xi…

Σ.j..

n.j..

x.j..

Σ..k.

n..k.

x..k.

Σij..

nij..

xij..

Σi.k.

ni.k.

xi.k.

Σ.jk.

n.jk.

x.jk.

Σijk.

nijk.

xijk.

Σ.…

n.…

x.…

А1

В1

С1

x1111

x1112

Σ1

n1

x1

Σ.1..

n.1..

x.1..

Σ..1.

n..1.

x..1.

Σ11..

n11..

x11..

Σ1.1.

n1.1.

x1.1.

Σ.11.

n.11.

x.11.

Σ111.

n111.

x111.

Σ.…

n.…

n.…

С2

x1121

x1122

Σ112.

n112.

x112.

С3

x1131

x1132

Σ1.2.

n1.2.

x1.2.

Σ.12.

n.12.

x.12.

Σ113.

n113.

x113.

С4

x1141

x1142

Σ114.

n114.

x114.

В2

С1

x1211

x1212

Σ..2.

n..2.

x..2.

Σ12..

n12..

x12..

Σ1.3.

n1.3.

x1.3.

Σ.13.

n.13.

x.13.

Σ121.

n121.

x121.

С2

x1221

x1222

Σ122.

n122.

x122.

С3

x1231

x1232

Σ1.4.

n1.4.

x1.4.

Σ.14.

n.14.

x.14.

Σ123.

n123.

x123.

С4

x1241

x1242

Σ124.

n124.

x124.

А2

В1

С1

x2111

x2112

Σ2

n2

x2

Σ.2..

n.2..

x.2..

Σ..3.

n..3.

x..3.

Σ13..

n13..

x13..

Σ2.1.

n2.1.

x2.1.

Σ.21.

n.21.

x.21.

Σ211.

n211.

x211.

С2

x2121

x2122

Σ212.

n212.

x212.

С3

x2131

x2132

Σ2.2.

n2.2.

x2.2.

Σ.22.

n.22.

x.22.

Σ213.

n213.

x213.

С4

x2141

x2142

Σ214.

n214.

x214.

В2

С1

x2211

x2212

Σ..4.

n..4.

x..4.

Σ14..

n14..

x14..

Σ2.3.

n2.3.

x2.3.

Σ.23.

n.23.

x.23.

Σ221.

n221.

x221.

С2

x2221

x2222

Σ222.

n222.

x222.

С3

x2231

x2232

Σ2.4.

n2.4.

x2.4.

Σ.24.

n.24.

x.24.

Σ223.

n223.

x223.

С4

x2241

x2242

Σ224.

n224.

x224.

 

 


компоненты формул для вычисления SS:

  частное от деления квадрата суммы всех наблюдений на объем комплекса

сумма дробей, в числителях которых квадраты сумм по градациям фактора А, в знаменателях - числа наблюдений по градациям фактора А (эта и последующие формулы используются для неравномерных комплексов, в равномерных - число наблюдений по градациям фактора А константное и выносится за знак суммы), число слагаемых дробей равно числу градаций по фактору А;

сумма дробей, в числителях которых квадраты сумм по градациям фактора В, в знаменателях - числа наблюдений по градациям фактора В, число слагаемых дробей равно числу градаций по фактору В;

сумма дробей, в числителях которых квадраты сумм по градациям фактора С, в знаменателях - числа наблюдений по градациям фактора С, число слагаемых дробей равно числу градаций по фактору С;

сумма дробей, в числителях которых квадраты сумм по сочетанию градаций факторов А и В, в знаменателях - числа наблюдений по сочетанию градаций факторов А и В,число слагаемых дробей равно числу сочетаний градаций по факторам А и В, то есть произведению аb;

сумма дробей, в числителях которых квадраты сумм по сочетанию градаций факторов А и С, в знаменателях - числа наблюдений по сочетанию градаций факторов А и С, число слагаемых дробей равно числу сочетаний градаций по факторам А и C, то есть произведению аc;

сумма дробей, в числителях которых квадраты сумм по сочетанию градаций факторов В и С, в знаменателях - числа наблюдений по сочетанию градаций факторов В и С, число слагаемых дробей равно числу сочетаний градаций по факторам B и C, то есть произведению bc;

сумма дробей, в числителях которых квадраты сумм по сочетанию градаций факторов А, В и С, в знаменателях - числа наблюдений по сочетанию градаций факторов А, В и С, число слагаемых дробей равно числу сочетаний градаций по факторам A, B и C, то есть произведению аbc;

Типы варьирования переменных в трехфакторном дисперсионном комплексе. В анализируемом комплексе различают 9 типов варьирования переменных:

1) общая вариация – варьирование отдельных наблюдений (хijkl) вокруг средней по комплексу ();

2) вариация, связанная с влиянием фактора А – варьирование средних по градациям фактора А () вокруг средней по комплексу ();

3) вариация, связанная с влиянием фактора В – варьирование средних по градациям фактора В () вокруг средней по комплексу ();

4) вариация, связанная с влиянием фактора C – варьирование средних по градациям фактора C () вокруг средней по комплексу ();

5) вариация, связанная со взаимодействием факторов А и В – варьирование средних по сочетанию градаций факторов А и В () вокруг средней по комплексу () без учета вариации, связанной с влиянием факторов А и В;

6) вариация, связанная со взаимодействием факторов А и C – варьирование средних по сочетанию градаций факторов А и C () вокруг средней по комплексу () без учета вариации, связанной с влиянием факторов А и C;

7) вариация, связанная со взаимодействием факторов B и C – варьирование средних по сочетанию градаций факторов B и C () вокруг средней по комплексу () без учета вариации, связанной с влиянием факторов B и C;

8) вариация, связанная со взаимодействием факторов А, B и C – варьирование средних по сочетанию градаций факторов А, B и C () вокруг средней по комплексу () без учета вариации, связанной с влиянием факторов А, B и C;

9) случайная вариация отдельных наблюдений () вокруг средних по сочетанию градаций факторов А, В и C ().

Суммы квадратов отклонений (SS). В трехфакторном дисперсионном анализе вычисляются 9 сумм квадратов отклонений (SS). Формулы для их вычисления следующие (комплекс неравномерный, в равномерном комплексе знаменатель в виде дроби 1/ n, где n – константа, выносится перед суммой):

Проверка корректности вычислений:

Числа степеней свободы. Формулы для вычислений чисел степеней свободы (df):

Проверка корректность вычислений чисел степеней свободы:

Средние квадраты (ms):

Критерии Фишера (F):

Стандартные значения критерия Фишера для каждого фактора и взаимодействий. Стандартные значения критерия Фишера определяются из таблиц на 5% и 1% уровнях значимости при соответствующих чисел степеней свободы по фактору (или взаимодействию) и случайных отклонений.

Сравнение эмпирических значений критерия Фишера со стандартными и формулировка статистических выводов. Если какое-либо эмпирическое значение F > F 05 и F > F 01, нулевая гипотеза об отсутствии влияния данного фактора (факторов) или взаимодействия факторов отклоняется.

Если F 05 < F < F 01, то на 5% уровне значимости нулевую гипотезу следует отклонить, а на 5% - принять. В этом случае рекомендуется продолжить экспериментальные исследования по изучению влияния данного фактора на изменчивость переменного.

Если F < F 05, то нулевая гипотеза об отсутствии влияния соответствующего фактора принимается.

Дисперсии для достоверно влияющих факторов и взаимодействий (σ 2):

Доли влияния факторов и взаимодействий (pin, %):

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

Стандартные значения критерия Тьюки (Q05). Критерии Тьюки определяются из таблицы по числу сравниваемых средних и числу степеней свободы случайной вариации.

Наименьшая существенная разность (НСР05):

НСР05 (А)=

НСР05 (В)=

НСР05 (С)=

НСР05 (АВ)=

НСР05 (АС)=

НСР05 (ВС)=

НСР05 (АВС)=

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

Оценка достоверности различий между средними по градациям факторов и взимодействий. Строятся матрицы разностей между средними по градациям факторов и взаимодействий. Затем разности сравниваются по абсолютной величине со значениями НСР05.

Результативная таблица (табл. 11.2.):

Таблица 11.2. Результаты трехфакторного дисперсионного анализа данных об изменчивости…:

Источник вариации

SS

df

ms

σ 2

F

F 05

F 01

pin %

НСР05

Общая

SSy

dfy

σy 2

pyin %

фактор А

SSA

dfA

msA

σA 2

FA

F 05

F 01

pAin %

НСР A 05

фактор В

SSB

dfB

msB

σB 2

FB

F 05

F 01

pBin %

НСР B 05

фактор С

SSC

dfC

msC

σC 2

FC

F 05

F 01

pCin %

НСР C 05

взаимодействие АВ

SSAB

dfAB

msAB

σAB 2

FAB

F 05

F 01

pABin %

НСР AB 05

взаимодействие АС

SSAC

dfAC

msAC

σAC 2

FAC

F 05

F 01

pACin %

НСР AC 05

взаимодействие ВС

SSBC

dfBC

msBC

σBC 2

FBC

F 05

F 01

pBCin %

НСР BC 05

взаимодействие АВС

SSABC

dfABC

msABC

σABC 2

FABC

F 05

F 01

pABCin %

НСР ABC 05

Случайная

SSz

df z

msz

σe 2

pein %

Выводы по трехфакторному дисперсионному анализу.

Пример 1. Изучали изменчивость урожайности томата в защищенном грунте (кг/м2) при малообъемной технологии в зависимости от трех факторов: массы субстрата (фактор А), интенсивности орошения (фактор В) и дозы удобрений (N,P,K) – фактор С:

Масса субстрата (фактор А) Интенсивность орошения (фактор В) Доза удобрений: N,P,K (фактор С)

Урожайность (кг/м2)

А1

В1

С1 17 18 21 20
С2 19 21 22 22
С3 20 22 23 23
С4 20 23 24 25

В2

С1 19 22 24 23
С2 21 22 25 24
С3 21 24 26 25
С4 23 23 25 25

А2

В1

С1 24 26 29 28
С2 27 28 32 30
С3 26 30 31 32
С4 30 29 33 32

В2

С1 29 30 32 33
С2 31 32 35 34
С3 34 35 38 37
С4 36 38 39 36

Какое влияние оказывают перечисленные факторы и их взаимодействие на урожайность томата в защищенном грунте?

Решение:

1. Копируем содержание таблицы в Excel

2. Строим трехфакторный дисперсионный комплекс

3. Определим основные параметры, суммы и компоненты формул:

объем комплекса (N) = 64

число градаций по фактору А (а) = 2

число градаций по фактору В (b) = 2

число градаций по фактору С (с) = 4

число наблюдений по градации фактора А (ni...)=32

число наблюдений по градации фактора В (n.j..) =32

число наблюдений по градации фактора С (n..k.) = 16

число наблюдений по сочетанию градаций АВ (nij..) = 16

число наблюдений по сочетанию градаций АC (ni.k.) = 8

число наблюдений по сочетанию градаций ВC (n.jk.) = 8

число наблюдений по сочетанию градаций AВC (nijk.) = 4

Σxijkl2 =48698

4. Суммы квадратов отклонений отдельных наблюдений от средних (SS):


Трехфакторный дисперсионный комплекс изменчивости урожайности томата (кг/м2) зависимости от массы субстрата (фактор А), интенсивности орошения (фактор В) и дозы удобрений (фактор С)

 

А

 В

 С

кг/м2

Σi…

ni…

xi…

Σ.j..

n.j..

x.j..

Σ..k.

n..k.

x..k.

Σij..

nij..

xij..

Σi.k.

ni.k.

xi.k.

Σ.jk.

n.jk.

x.jk.

Σijk.

nijk.

xijk.

Σ.…

n.…

x.…

А1

В1

С1

17

18

21

20

712 А1

32

22,3

807 В1

32

25,2

395 С1

16

24,69

340 А1В1

16

21,3

164 А1С1

8

20,5

183 В1С1

8

22,88

76

4

19,0

1728

64

27,0

С2

19

21

22

22

84

4

21,0

С3

20

22

23

23

176 А1С2

8

22,0

201 В1С2

8

25,13

88

4

22,0

С4

20

23

24

25

92

4

23,0

В2

С1

19

22

24

23

425  С2

16

26,56

372 А1В2

16

23,3

184 А1С3

8

23,0

207 В1С3

8

25,88

88

4

22,0

С2

21

22

25

24

92

4

23,0

С3

21

24

26

25

188 А1С4

8

23,5

216 В1С4

8

27

96

4

24,0

С4

23

23

25

25

96

4

24,0

А2

В1

С1

24

26

29

28

1016 А2

32

31,8

921 В2

32

28,8

447 С3

16

27,94

467 А2В1

16

29,2

231 А2С1

8

28,9

212 В2С1

8

26,5

107

4

26,8

С2

27

28

32

30

117

4

29,3

С3

26

30

31

32

249 А2С2

8

31,1

224 В2С2

8

28

119

4

29,8

С4

30

29

33

32

124

4

31,0

В2

С1

29

30

32

33

461 С4

16

28,81

549 А2В2

16

34,3

263 А2С3

8

32,9

240 В2С3

8

30

124

4

31,0

С2

31

32

35

34

132

4

33,0

С3

34

35

38

37

273 А2С4

8

34,1

245 В2С4

8

30,63

144

4

36,0

С4

36

38

39

36

149

4

37,3

 


Проверим корректность вычислений:

5. Вычислим числа степеней свободы (df):

Проверим корректность вычислений:

6. Вычислим средние квадраты (ms):

7. Вычислим критерии Фишера (F):

8. Определим стандартные значения критерия Фишера для каждого фактора (стат. табл. IV-V):

для фактора А: F05 (dfa =1; dfz =48)=4,04; F01 =7,21;

для фактора В: F05 (dfb =1; dfz =48)=4,04; F01 =7,21;

для фактора С: F05 (dfс =3; dfz =48)=2,79; F01 =4,21;

для взаимодействия факторов АВ: F05 (dfab =1; dfz =48)=4,04; F01 =7,21;

для взаимодействия факторов AC: F05 (dfac =3; dfz =48)=2,79; F01 =4,21;

для взаимодействия факторов ВC: F05 (dfbc =3; dfz =48)=2,79; F01 =4,21;

для взаимодействия факторов AВC: F05 (dfabc =3 dfz =48)= 2,79; F01 =4,21;

9. Сравним эмпирические значения критерия Фишера со стандартными:

Fa =396,07> F01 =7,21;

Fb =55,70> F01 =7,21;

Fc =14,19> F01 =4,21;

Fab =10,71> F01 =7,21;

Fac =1,03< F05 =2,79;

Fbc =0,29< F05 =2,79;

Fabc =1,02< F05 =2,79;

10. Сформулируем статистические выводы:

по фактору А: нулевая гипотеза отклоняется на 1% уровне значимости (фактор влияет);

по фактору В: нулевая гипотеза отклоняется на 1% уровне значимости (фактор влияет);

по фактору С: нулевая гипотеза отклоняется на 1% уровне значимости (фактор влияет);

по взаимодействию факторов АВ: нулевая гипотеза отклоняется на 1% уровне значимости (взаимодействие АВ влияет);

по взаимодействию факторов АС: нулевая гипотеза принимается на 5% уровне значимости (взаимодействие АС не влияет);

по взаимодействию факторов ВС: нулевая гипотеза принимается на 5% уровне значимости (взаимодействие ВС не влияет);

по взаимодействию факторов АВС: нулевая гипотеза принимается на 5% уровне значимости (взаимодействие АВС не влияет).

11. Вычислим дисперсии (σ2):

 (взаимодействие АС не влияет)

 (взаимодействие ВС не влияет)

 (взаимодействие АВС не влияет)

12. Вычислим доли влияния факторов (pin%):

13. Вычислим ошибки групповых средних по факторам (mx), ошибки вычисляются только для тех средних, по которым доказано достоверное влияние фактора:

14. Определяем критерий Тьюки (Q05) по стат. табл.VI:

для фактора А: Q05 (а =2; dfZ =48)=2,85

для фактора В: Q05 (b =2; dfZ =48)=2,85

для фактора С: Q05 (с=4; dfZ =48)=3,75

для взаимодействия факторов АВ: Q05 (ab =4; dfZ =48)=3,75

15. Вычисляем НСР05:

для фактора А: НСР05=

для фактора В: НСР05=

для фактора C: НСР05=

для взаимодействия факторов АВ: НСР05=

16. Сводим все данные в результативную таблицу:

Результаты трехфакторного дисперсионного анализа данных об изменчивости урожайности томата в защищенном грунте (кг/м2) в зависимости от массы субстрата (фактор А), интенсивности орошения (фактор В) и дозы удобрений (фактор С):

Источник вариации

SS

df

ms


Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями:  



double arrow