
Биометрия методичка
.pdf
Таблица 6.2 – Результаты расчетов коэффициентов асимметрии и эксцесса
Сорта яровой пшеницы |
Коэффициент КА |
Коэффициент КЕ |
Задворская |
0,438 |
0,792 |
Янтарная |
0,084 |
0,804 |
Новобелицкая |
1,433 |
1,502 |
Солнечная |
0,429 |
0,063 |
Далее переходим к ОДА и применяем команды Данные /
Анализ данных / Однофакторный дисперсионный анализ.
Появляется диалоговое окно, показанное на Рис. 6.2.
В окошко «Входной интервал» вводим весь дисперсионный комплекс (с номерами групп), для этого обводим его курсором. Группировку ставим «по строкам», т. е. по группам поросят. Флажок в окошке «Метки в первом столбце» вводит в итоговую таблицу номера групп. Максимальное требуемое значение уровня значимости (0,05) вводим в окошко «Альфа». Далее указываем ячейку для выходного интервала (например, G3) и нажимаем кнопку «ОК».
Рис. 6.2 – Диалоговое окно инструмента «Однофакторный дисперсионный анализ»
41

Результаты применения инструмента «Однофакторный дисперсионныйанализ» представленынаРис. 6.3.
Вэтойтаблицеприводятся:
средние значения признака («Среднее») и его вариансы («Дисперсия») погруппам;
«SS» – суммыквадратов;
«df» – числостепенейсвободы;
«MS» – средние квадраты (значения SS, которые разделены насоответствующиеdf);
«F» – фактическое значение F-критерия (отношение первой ячейкистолбца«MS» ковторой).
Критическое значение F-критерия приведено на уровне
значимости p < 0,05, которое было указано в диалоговом окне.
Рис. 6.3 – Результаты применения инструмента «Однофакторный дисперсионный анализ»
Из данных Рис. 6.3 видно, что влияние фактора (сортовые особенности) на урожайность яровой пшеницы оказалось значимым на уровне p < 7,78148 · 10–18 (это столбик «Р-значение» на Рис. 6.3). Показателем силы влияния фактора является отношение первой ячейки столбика SS к послед-
ней (Итого): ηх2 = 2563,6616 / 2603,5776 = 0,985. Это свиде-
тельствует о том, что влияние фактора (сортовые особенности) на урожайность яровой пшеницы очень существенно.
42
Задания Задание 1. Изучали влияние улучшенной схемы кормле-
ния различных пород свиней мясо-сального направления. Исследования проводили в шестикратной повторности. Результаты исследований – масса свиней 100-дневного возраста – приведены в Таблице 6.4.
Необходимо установить достоверность и силу влияния фактора – улучшенной схемы кормления – на массу свиней.
Таблица 6.4 – Масса свиней различных пород 100-дневного возраста
Повторность |
|
Масса свиней, кг |
|
||
Муромская |
Миргородская |
Ливенская |
Кемеровская |
||
|
|||||
1 |
51,6 |
49,2 |
51,4 |
55,0 |
|
2 |
50,5 |
48,5 |
49,7 |
51,4 |
|
3 |
53,3 |
46,8 |
48,6 |
49,7 |
|
4 |
52,4 |
47,4 |
51,2 |
50,6 |
|
5 |
54,0 |
45,8 |
50,5 |
52,9 |
|
6 |
51,9 |
46,2 |
51,2 |
54,8 |
Задание 2. Изучали отличия сортов голубики высокорослой по массе ягоды. Исследования проводили в шестикратной повторности. Результаты исследований приведены вТаблице 6.5.
Необходимо установить достоверность отличия сортов и силу влияния фактора.
Таблица 6.5 – Отличия сортов голубики высокорослой по массе ягоды
Повторность |
|
|
Масса ягоды, г |
|
|
|
Патриот |
Бонус |
Дюк |
Река |
Элизабет |
Чандлер |
|
1 |
1,9 |
2,4 |
2,5 |
1,8 |
2,4 |
3,0 |
2 |
1,6 |
2,5 |
1,9 |
1,6 |
2,0 |
2,4 |
3 |
1,8 |
2,8 |
2,1 |
1,8 |
2,1 |
2,3 |
4 |
1,7 |
3,0 |
2,0 |
1,7 |
1,9 |
2,6 |
5 |
1,8 |
2,9 |
2,4 |
1,6 |
2,3 |
2,8 |
6 |
1,7 |
3,1 |
2,3 |
1,8 |
2,2 |
2,7 |
43
Лабораторная работа № 7 Многофакторный дисперсионный анализ
Двухфакторный дисперсионный анализ
Двухфакторный дисперсионный анализ используют для исследования влияния на результативный признак одновременно двух факторов. При этом выбор факторов ограничен требованием полной независимости их между собой. Факторы могут иметь как количественную, так и качественную градацию.
Пример. Изучали урожайность двух сортов крыжовника и их устойчивость против вредоносного действия крыжовникового пилильщика. Получили следующие результаты (Таблица 7.1). Необходимо провести двухфакторный дисперсионный анализ данных с помощью пакета «Анализ данных» программы Ms Excel.
Таблица 7.1 – Урожайность двух сортов крыжовника и их устойчивость против вредоносного действия крыжовникового пилильщика
Урожай с отдельных кустов крыжовника, кг |
Сорта крыжовника (А) |
||
английский желтый |
малахитовый |
||
|
(А1) |
(А2) |
|
|
5,4 |
4,6 |
|
|
6,5 |
4,3 |
|
Не поврежденных пилильщиком (В1) |
4,2 |
3,9 |
|
4,9 |
4,3 |
||
|
|||
|
6,1 |
3,5 |
|
|
5,2 |
3,6 |
|
|
5,2 |
3,9 |
|
|
4,3 |
3,1 |
|
Поврежденных пилильщиком (В2) |
3,5 |
3,2 |
|
4,5 |
3,8 |
||
|
|||
|
4,2 |
2,6 |
|
|
5,2 |
4,2 |
Целью исследований было установление достоверности влияния на урожай крыжовника каждого из двух факторов в отдельности, а также степень влияния их совместного действия. Эти исследования были организованы в виде двухфакторного равномерного дисперсионного комплекса с шестикратной повторяемостью, поскольку объемы выборок были неизменными (n = 6). Поэтому следует применить двухфакторный дисперсионный анализ с повторениями.
44

Решение. Открываем новую книгу (или новый рабочий лист) в ППП Excel и копируем Таблицу 7.1. Затем применя-
ем команды Данные / Анализ данных / Двухфакторный дисперсионный анализ с повторениями. Появляется диало-
говое окно, которое показано на Рис. 7.1.
Рис. 7.1 – Диалоговое окно инструмента «Двухфакторный дисперсионный анализ с повторениями»
В окошко «Входной интервал» вводим весь двухфакторный дисперсионный комплекс с надписями, обводя его курсором. В окошко «Число строк для выборки» вводим объем каждой выборки – 6. Максимально допустимое значение уровня значимости (0,05) вводим в окошко «Альфа». Далее указываем ячейку выходного интервала (например, G11) и нажимаем кнопку «ОК». Результаты применения инструмента «Двухфакторный дисперсионный анализ с повторениями» приведены на Рис. 7.2. Разбиение Рис. 7.2 на секторы «B1», «B2», «Итого» отражает влияние факторов A и B по отдельности, а также влияние их совместного действия – AB.
Источники вариации:
«Выборка» – фактор B;
«Столбцы» – фактор A;
«Взаимодействие» – A и B совместно;
«Внутри» – случайная вариация признака внутри каждой выборки.
45

Рис. 7.2 – Результаты применения инструмента «Двухфакторный дисперсионный анализ с повторениями»
Фактическое значение F-критерия (столбик «F») определяется из столбика «МS» (средние квадраты) как отношение значения соответствующей ячейки столбика «MS» к значению последней «Внутри». «Р-значение» – это значимость влияния факторов в соответствии с F-критерием. Критические значения F-критерия приведены на уровне значимости p < 0,05, которое указано в диалоговом окне.
Как видно из Рис. 7.2, на результативный признак достоверно влияет каждый из факторов A и B, но отсутствует влияние совместного действия факторов – АВ (поскольку
F (0,402414) меньше F-критическое (4,3512435)).
Силу влияния факторов можно определить из 1-го столбика раздела «Дисперсионный анализ» «SS» (сумма квадратов), разделив соответствующее значение SS на суммарное
(20,0783333).
46

|
Так, |
|
|
для фактора В: |
ηВ |
= 0,160; |
|
для фактора А: |
ηА |
= 0,418; |
|
|
для AВ: |
ηAВ |
= 0,008; |
|
для случайных: |
ηсл |
= 0,412. |
Из рисунка видно, что наибольшее влияние на результативный признак оказывает действие фактора А – влияние сортовых особенностей крыжовника.
Задания
В Таблицу 7.2 занесены данные о влиянии породных свойств и качества пчелиных маток на яйценоскость потомства. Необходимо провести двухфакторный дисперсионный анализ данных и сделать соответствующие выводы.
Таблица 7.2 – Яйценоскость пчел
Породы пчел |
|
А1 |
|
|
А2 |
|
|
А3 |
|
Матки |
В1 |
В2 |
В3 |
В1 |
В2 |
В3 |
В1 |
В2 |
В3 |
|
14 |
14 |
15 |
14 |
20 |
16 |
10 |
16 |
12 |
|
10 |
14 |
17 |
17 |
18 |
17 |
12 |
12 |
15 |
Яйценоскость |
15 |
15 |
16 |
20 |
18 |
20 |
8 |
10 |
8 |
дочерей, шт. |
14 |
17 |
13 |
16 |
16 |
17 |
13 |
16 |
12 |
|
16 |
15 |
14 |
15 |
20 |
18 |
18 |
10 |
16 |
|
12 |
16 |
15 |
20 |
19 |
20 |
17 |
13 |
18 |
Получены данные (Таблица 7.3) о массе трёх пород поросят. Необходимо провести двухфакторный дисперсионный анализ данных и сделать соответствующие выводы о влиянии породных свойств.
Таблица 7.3 – Масса разных пород поросят
Породы свиней |
Кемеровская |
|
Муромская |
Миргородская |
||||||
Свиноматки |
С1 |
С2 |
С3 |
С1 |
С2 |
С3 |
С1 |
С2 |
С3 |
|
|
23 |
20 |
25 |
24 |
20 |
24 |
25 |
20 |
17 |
|
|
28 |
28 |
29 |
21 |
19 |
21 |
22 |
17 |
20 |
|
Масса поросят, кг |
17 |
23 |
27 |
16 |
19 |
23 |
24 |
19 |
19 |
|
21 |
22 |
19 |
17 |
18 |
27 |
27 |
14 |
21 |
||
|
||||||||||
|
18 |
21 |
22 |
24 |
20 |
24 |
21 |
17 |
19 |
|
|
18 |
22 |
23 |
19 |
22 |
18 |
23 |
19 |
23 |
47

Лабораторная работа № 8 Корреляционный анализ
Если исследуется зависимость результативного признака Y от нескольких признаков (факторов) X1, X2, ..., Xn, то на первом этапе сравнивают между собой выборки из этих признаков (разумеется, одинакового объема).
Для оценки силы линейной связи между признаками в ППП MS Excel применяют инструмент «Корреляция». Он позволяет рассчитать выборочные коэффициенты парной корреляции между результирующим признаком Y и каждым из признаков X1, X2, ..., Xn, а также между всеми парами признаков X1, X2, ..., Xn – корреляционную матрицу.
Пример 1. Нужно исследовать корреляционную зависимость между жирностью молока Y и надоем X за лактационный период для коров чёрно-пёстрой породы.
Соответствующие данные приведены в Таблице 8.1.
Таблица 8.1 – Надой и жирность молока коров чёрно-пёстрой породы за лактационный период
Номер животного |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
Надой X, ц |
40,5 |
38,1 |
41,3 |
38,5 |
41,1 |
39,3 |
39,5 |
40,3 |
40,5 |
38,7 |
39,5 |
38,7 |
Жирность Y, % |
3,5 |
3,9 |
3,4 |
3,7 |
3,5 |
3,8 |
3,7 |
3,5 |
3,3 |
3,8 |
3,5 |
3,7 |
Решение. Открываем новую книгу (или новый рабочий лист) в ППП Excel и копируем Таблицу 8.1. Потом применя-
ем команды Данные / Анализ данных / Корреляция. Появля-
ется диалоговое окно, которое показано на Рис. 8.1.
Вокошко «Входной интервал» вводим данные таблицы
снадписями, обводя их курсором. При этом группирование делаем «По строкам».
Чтобы надписи попали в итоговую таблицу, выставляем флажок в окошке «Метки в первой строке». Указываем ячейку выходного интервала (например, В36) и нажимаем кнопку
«ОК».
Результаты приведены на Рис. 8.2.
48

Рис. 8.1 – Диалоговое окно инструмента «Корреляция»
По результатам вычислений, представленным на Рис. 8.2 видно, что корреляция высокая и отрицательная, т. е. с увеличением надоя Х жирность молока Y существенно уменьшается.
Рис. 8.2 – Результаты расчета коэффициента корреляции
Задания
Задание 1. Изучали зависимость между массой матерей (mм) лисицимассойноворожденныхдетенышей(mд) (Таблица8.2).
Таблица 8.2 – Масса матерей лисиц (кг) и их новорожденных детенышей (кг)
mм |
10,1 |
11,2 |
12,4 |
10,3 |
10,4 |
10,1 |
11,3 |
10,2 |
13,4 |
12,3 |
14,5 |
9,8 |
12,2 |
11,8 |
13,2 |
10,8 |
10,1 |
14,6 |
13,4 |
12,8 |
mд |
0,68 |
0,9 |
1,2 |
0,7 |
0,7 |
0,65 |
0,8 |
0,6 |
0,7 |
1,1 |
1,2 |
0,5 |
0,9 |
0,7 |
1,1 |
0,7 |
0,6 |
1,1 |
1,2 |
0,9 |
Задание 2. Изучали зависимость между длиной цветоносных побегов (X) и длиной листовой пластинки (Y) земляники садовой сорта Неженка (Таблица 8.3).
49

Таблица 8.3 – Длина цветоносных побегов (X) и длина листовой пластинки (Y) земляники садовой (сорт Неженка)
X |
7,2 |
8,7 |
9,8 |
9,6 |
7,4 |
7,5 |
9,5 |
10,1 |
7,3 |
8,9 |
10,4 |
9,8 |
8,7 |
9,8 |
7,5 |
8,8 |
9,1 |
8,6 |
8,4 |
10,1 10,4 |
|
Y |
3,1 |
3,8 |
4,1 |
4,2 |
3,3 |
3,2 |
3,7 |
4,4 |
3,2 |
3,6 |
4,2 |
3,9 |
3,5 |
3,5 |
3,5 |
3,4 |
3,3 |
3,1 |
3,2 |
3,4 |
4,2 |
Задание 3. Изучали зависимость массы колоса пшеницы
(Y) сорта Мироновская от ширины (X1) и длины (X2) листовой пластинки (Таблица 8.4).
Таблица 8.4 – Масса колоса (Y), ширина (X1) и длина (X2) листовой пластинки пшеницы (сорт Мироновская)
Y, г |
62,5 |
80,2 |
74,1 |
72,3 |
68,2 |
74,6 |
79,5 |
81,3 |
64,4 |
68,7 |
78,2 |
67,5 |
63,4 |
78,5 |
80,6 |
69,4 |
79,3 |
77,8 |
62,2 |
72,6 |
79,8 |
X1, мм |
10 |
18 |
16 |
14 |
13 |
15 |
16 |
17 |
12 |
13 |
14 |
12 |
11 |
13 |
14 |
11 |
16 |
15 |
11 |
13 |
14 |
X2, см |
7,3 |
12,6 |
11,4 |
12,1 |
9,2 |
10,8 |
11,9 |
12,5 |
8,4 |
9,2 |
11,4 |
12,1 |
8,3 |
11,6 |
12,2 |
8,8 |
11,2 |
10,8 |
7,8 |
8,6 |
8,8 |
50