Структура трехфакторного дисперсионного комплекса. Трехфакторный эксперимент значительно более сложен и по организации эксперимента и по анализу данных по сравнению с двухфакторным экспериментом. Сложность анализа данных заключается в необычной структуре дисперсионного комплекса, большом числе источников вариации и вычисляемых параметров.
Поскольку факторов три, модель комплекса, где строчки – градации одного фактора, а столбцы – другого не подходит. При построении комплекса по трём факторам эти факторы приходится располагать слева направо (как это принято в иерархическом комплексе).
В табл. 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) и не забудь поделиться с друзьями:
Сейчас читают про:
|