BiometriaKnuga
.pdfРозділ 3.6 ЕкологіяРегресійний аналіз в біологічних і екологічних дослідженнях |
|
125 |
|
||
|
|
|
6.8. Контрольні питання
1.В чому полягає регресійний аналіз даних спостережень біологічних об'єктів?
2.Які існують види кривих регресій і які їх особливості?
3.Як характеризується прямолінійна регресія, яка формула її апроксимації?
4.Як характеризується криволінійна регресія?
5.Як визначається помилка репрезентативності прямолінійної регресії?
6.Як визначається параболічна регресія, які основні її види?
7.В чому полягає гіперболічна регресія?
8.Які основні властивості коефіцієнту регресії?
9.Як визначається критерій достовірності різниці між коефіцієнтами регресії?
10.Які особливості застосування регресійного аналізу при дослідженнях біологічних об'єктів?
126
Розділ 7.
Дисперсійний
аналіз
Біометрія
7.1. Дисперсія як критерій порівняльної оцінки варіаційних рядів
Дисперсія (від лат. dispersio – розсіювання) є мірою відхилення показників варіаційного ряду від середнього арифметичного значення цього ряду (x). Як вже визначалось раніше, біологічні об'єкти ніколи не виступають як незалежні представники флори чи фауни – вони завжди пов'язані з оточуючим середовищем, складною системою стосунків, вертикальних та горизонтальних зв'язків. На кожний біологічний об'єкт здійснюють сумарний вплив багато факторів як біотичних, так і абіотичних. Цей вплив проявляється у формуванні його індивідуальних особливостей, що пов'язано з дією механізмів адаптації.
Сума індивідуальних ознак біологічних об'єктів формує особливості їх угрупувань. При проведенні біологічних та екологічних досліджень часто необхідно з'ясувати: яким чином той чи інший фактор вплинув на підсилення або послаблення тієї чи іншої характеристики. Це питання набуває особливо важливого значення у зв'язку з необхідністю встановлювати силу впливу окремих антропогенних дій в ситуації величезного антропогенного навантаження на природне середовище.
Завдання ускладнюється тим, що живий організм реагує у своєму розвитку одночасно на зміну декількох факторів, внаслідок чого визначення характеру впливу одного з них є досить важким, а іноді і невирішуваним. Так, наприклад, вивчаючи вплив технології доглядів за молодими рослинами
Розділ 3.7 ДисперсійнийЕкологія аналіз |
127 |
|
деревних порід, дослідник мусить враховувати вплив факторів – зміна аерації корененасиченого шару грунту, зміну температури, зміну вологості і т.д. Вивчаючи вплив різних доз або різного складу добрив на підвищення урожайності пшениці, жита або іншої сільськогосподарської культури, термінів і способів внесення їх у грунт, також необхідно враховувати вплив відповідного комплексу факторів, частину з яких дослідник може регулювати за власним бажанням, а прояв іншої частини від нього не залежить. Вивчаючи вплив нового виду домішок до корму худоби на показники її продуктивності, також необхідно враховувати одночасну дію інших умов оточуючого середовища. Вирішення поставлених завдань можливо лише порівнянням результатів досліду, в яких даний фактор виступає у відповідній градації. А оскільки всі фактори, що впливають на результат, виступають в комплексі, то застосування знайомого нам методичного прийому парного порівняння рядів розподілу стає невиправдано важким, а іноді, неможливим.
Для подальшого викладення приймемо, що причини зміни характеристик біологічних об'єктів звуться факторами, а ті характеристики біологічних об'єктів, які змінюються під їх впливом, звуться результативними ознаками.
Вивчаючи дію окремого фактору F1, треба мати на увазі, що для цього необхідно відокремлювати зміни характеристик, покликані дією цього фактору – результативні ознаки (F1) від результативних ознак (F2), (F3), ..., (Fn). Цій меті і відповідає дисперсійний аналіз.
Його сутність полягає в наступному. Під час дослідження одержуємо декілька вибіркових сукупностей, в яких фактор (-и), який нас цікавить (назвемо його фокусним), діє з різною силою. Обрахувавши біометричні показники вибіркових сукупностей, ми знаємо, що дія фокусного фактора, як і дія інших факторів, вплинула на отримані показники, але поки що не можемо сказати в якій мірі.
Часто можна зустріти такий розподіл факторів в дисперсійному аналізі, коли вони позначаються як регульовані і нерегульовані. В цьому випадку має місце проведення так званого умовного експерименту, коли ще на стадії планування дослідження передбачається: скільки градацій фактору слід ввести, регулюючи його силу (кількість доглядів за деревостанами, дозу добрив і т.ін.). Але аналіз може проводитись і стосовно неантропогенних факторів – вміст хімічних елементів в грунті, рівень грунтових вод та ін. В цьому випадку дослідник відбирає випадкові об'єкти з сукупностей із
128 |
|
Біометрія |
|
||
|
|
заздалегідь відомою силою впливу нерегульованого фактору. Тому в цьому посібнику замість терміну "регульований фактор" використається термін "фокусний фактор", що позначає фактор, для якого досліджується статистична значимість ступеня впливу на формування результативної ознаки. Але в будь-якому випадку результати дисперсійного аналізу можуть бути достовірними лише за умов випадково зроблених незалежних спостережень у відповідній їх кількості.
В 1925 році англійський вчений Р.А.Фішер запропонував для вирішення цієї проблеми метод, який ми сьогодні і називаємо "дисперсійним аналізом". Користуючись ним, ми передбачаємо, що дисперсія результативної ознаки всього варіаційного комплексу (який ми далі будемо звати дисперсійним) покликана впливом всієї суми факторів на біологічні об'єкти – складові вибіркових сукупностей і може бути розкладена на факторну дисперсію, покликану впливом фокусного фактору, і залишкову, покликану впливом всіх інших факторів.
Зіставляючи факторну дисперсію з остаточною і порівнюючи результат із теоретичним значенням критерію F-Фішера для прийнятого рівня значимості (частіше – 5 %, рідше 1 % або 0,1 %) та ступенів вільності факторної та остаточної дисперсії, яке знаходять по спеціальній таблиці, розробленій Фішером, ми дізнаємось – чи є статистично значимою різниця між цими видами дисперсій, тобто чи достовірна різниця між впливом фокусного фактору і всіх інших.
Види дисперсій, які фігурують в дисперсійному аналізі, отримаються наступним чином.
Загальна дисперсія (σзаг2), покликана впливом на складові дисперсійного комплексу всіх факторів, вже відома нам з розділу 7 дисперсія, основана на сумі квадратів відхилень окремих варіант (xi) від загальної середньої (x) – Σ(xi
–x)2 (позначається Dy). Число степенів вільності для неї дорівнює N – 1, де N
–загальна кількість варіант дисперсійного комплексу (позначається Ky).
Якщо фокусний фактор значимо впливає на результативну ознаку, то групи показників будуть відрізнятись одна від одної тим сильніше, чим сильніше вплив фокусного фактору. Це знайде свій прояв в різниці групових середніх (xi) вибіркових сукупностей і наприкінці – в степені розсіяння цих групових середніх навколо загальної середньої арифметичної. Ця степінь розсіяння позначається як Σ(xi – x)2niDx, де ni–кількість ознак по градаціях.
Розділ 3.7 ДисперсійнийЕкологія аналіз |
129 |
|
Відхилення зводяться у квадрат заради запобігання взаємного знищення від'ємних і додатних результатів. Відношення факторіальної суми квадратів відхилень групових середніх від загальної середньої до відповідної кількості ступенів вільності і є факторною (міжгруповою) дисперсією (σфакт2). Число ступенів вільності для факторної дисперсії дорівнює а – 1, де a – кількість градацій фокусного фактору (позначається Kx).
Те, що в середині градацій фокусного фактора має місце відхилення окремих групових варіант (xi) від групових середніх (xi), свідчить про вплив на результативну ознаку і інших факторів, окрім фокусного. Сума квадратів цих відхилень по всіх градаціях комплексу (Dz = Σ(Σ(xi – x)2)) поділена на відповідну кількість ступенів вільності (Kz = Σa(ni – 1), де а – кількість градацій фокусного фактору, а ni – кількість показників в окремій градації) є остаточною (випадковою) дисперсією (σост2).
Слід відмітити, що існує математична залежність між сумами квадратів відхилень, яка виражається як
або
Dy = Dx + Dz
ітака ж сама залежність між числом ступенів вільності:
∑(x − x)2 = ∑(xi − x)2 ni + ∑(∑(xi − xi )2 )
K y = Kx + Kz або |
(N −1) = (a −1) + ∑a(ni −1). |
Ці рівняння можна використати як для перевірки точності визначених сум квадратів відхилень і ступенів вільності, так і для знаходження третього члена рівняння по двох відомих.
Дисперсійний аналіз є одним з найдосконаліших біометричних методів. За його допомогою можна аналізувати вплив не лише одного фактору, а двох,
130 |
|
Біометрія |
|
||
|
|
трьох і таке інше. В такому випадку факторна дисперсія розглядається як сума дисперсій, покликаних дією фокусних факторів F1, F2 , F3, ... , Fn. В цьому посібнику ми обмежимось розглядом лише дисперсійного аналізу однофакторних комплексів. Дисперсійний аналіз дозволяє вивчати і неоднорідні за біологічною природою сукупності об'єктів (різні популяції, особини різного полу, віку і т.ін.). Це обумовлює його високу цінність для біологічних досліджень. Особлива привабливість дисперсійного аналізу в тому, що навіть при великій кількості спостережень він дозволяє обійтися без громіздких розрахунків.
Дисперсійний аналіз може застосовуватись також і для встановлення однорідності вибіркових сукупностей (дисперсії, яких повинні бути однаковими). У разі встановлення однорідності сукупності можуть бути об'єднані в одну. Більша кількість спостережень в такому разі дозволить робити більш грунтовні статистичні висновки, підвищувати точність обрахування результатів експерименту.
130 |
|
Біометрія |
|
7.2. Критерії дисперсійного аналізу
При проведенні дисперсійного аналізу одночасно обробляються дані декількох вибіркових сукупностей, які складають дисперсійний комплекс. Всі дані при неавтоматизованому обрахуванні оформлюються в спеціальну робочу таблицю. Один з вимірів таблиці відповідає градаціям фокусного фактору (наприклад, доза ліків: "10 мг", "20 мг", "30 мг" ... ), інший вимір таблиці відповідає повторностям виміру результативних ознак. Варіанти таблиць (при дотриманні наведеного правила) можуть бути різними (що можна побачити в наведених далі прикладах).
Дисперсійні комплекси можуть бути ортогональними (такими, в яких значення результативних ознак рівномірно групуються за градаціями комплексів – порівну або пропорційно) та неортогональними (в яких вони групуються нерівномірно).
Оцінка достовірності різниці результативних ознак (або впливу фокусного фактору) в дисперсійному аналізі проводиться за критерієм F- Фішера. Цей критерій являє собою відношення фактичної дисперсії (σфакт2) до залишкової дисперсії (σост2). Таким чином фактичний критерій Fфакт = σ
факт2/σост2. Після визначення Fфакт його значення порівнюється з табличним
Розділ 3.7 ДисперсійнийЕкологія аналіз |
|
131 |
|
||
|
|
|
значенням Fтабл, які представлені за відповідними ступенями вільності (додатки 3, 4).
Нульова гіпотеза полягає в припущенні відсутності впливу організованого фактора на результативну ознаку. Вона відкидається, якщо Fфакт буде рівним або більше Fтабл. В даному випадку приймається, що вплив організованого фактора на результативну ознаку не є випадковим, тобто стверджується наявність достовірної різниці генеральних параметрів за градаціями дисперсійного комплексу.
Якщо значення Fфакт буде менше значення Fтабл, то це означає, що нульова гіпотеза підтверджується, тобто та різниця, яка зафіксована між груповими середніми, є статистично недостовірною. По встановленні достовірності статистично значимого впливу фокусного фактору на результативну ознаку, якщо необхідно встановити, яка з градацій фактору впливає на неї найбільш значимо, слід попарно порівняти середні.
Далі наводиться алгоритм дисперсійного аналізу, який відображає логічну сутність послідовності його проведення (табл. 8).
Для визначення дисперсій Dy, Dx, Dz застосовуються такі робочі формули:
|
= ∑x2 |
|
(∑x)2 |
; Dx = ∑ |
(∑x )2 |
(∑x)2 |
||||
Dy |
− |
|
|
i |
− |
|
; |
|||
|
N |
|
|
|||||||
|
|
|
|
|
|
n |
N |
|||
|
|
|
|
|
|
|
i |
|
|
|
|
= ∑x2 |
− ∑ |
(∑x )2 |
|
|
|
|
|||
Dz |
|
i |
. |
|
|
|
|
|||
|
|
|
|
|
|
|||||
|
|
|
|
n |
|
|
|
|
||
|
|
|
|
|
i |
|
|
|
|
|
Тут х – варіанти дисперсійного комплексу, хі – варіанти, що входять в склад окремих градацій або груп дисперсійного комплексу; x – загальна середня арифметична; xі – групові середні арифметичні; nі – кількість варіант по градаціях фокусного фактору.
Розділ 3.7 ДисперсійнийЕкологія аналіз |
|
131 |
|
||
|
|
|
7.3. Приклади дисперсійного аналізу
Далі наведено приклад аналізу рівномірного однофакторного комплексу. Такий вид аналізу є досить часто поширеним в біометрії. Досліджувався вплив кількості доглядів (розпушень грунту) на ріст 2-х річних сіянців сосни звичайної після висаджування їх в грунт. Дослід організовано на однорідній території з виділенням дослідних площ розміром 0,5 га. Сукупність рослин на кожній дослідній площі є відповідною генеральною сукупністю. Вибіркові
132 |
|
|
|
Біометрія |
|
|
|
||
|
||||
Таблиця 8. Алгоритм дисперсійного аналізу однофакторних комплексів |
||||
|
|
|
|
|
номер |
|
Зміст організаційної або математичної дії |
Результат дії |
|
дії |
|
|
|
|
1 |
Групування вибіркових матеріалів у |
Таблиця дисперсійного |
||
|
комбінаційну таблицю дисперсійного |
комплексу |
||
2 |
Визначення значень: середнього арифметичного |
Середні значення (x) і (xi) |
||
|
всього комплексу (x) і групових середніх за |
|
|
|
|
градаціями організованого фактора (xi) |
|
|
|
3 |
Визначення загальної суми квадратів відхилень |
Значення Dy |
||
|
(Dy), тобто суми квадратів відхилень варіант від |
|
|
|
|
загальної середньої: Dy=(x– x)2 |
|
|
|
4 |
Визначення міжгрупової суми квадратів |
ЗначенняDx |
||
|
відхилень, яка дорівнює сумі квадратів |
|
|
|
|
відхилень групових середніх від загальної |
|
|
|
|
середньої з урахуванням статистичної ваги (ni) |
|
|
|
|
групових середніх: |
|
|
|
|
• |
у випадку рівних чисел варіант в |
|
|
|
|
градаціях комплексу – Dx = ni Σ (xi – x)2; |
|
|
|
• |
у випадку різної кількості варіант в |
|
|
|
|
градаціях комплексу – |
|
|
5 |
Визначення внутрішньогрупової суми квадратів, |
ЗначенняDz |
||
|
тобто суми квадратів відхилень групових |
|
|
|
|
варіант від групових середніх: |
|
|
|
6 |
Визначення дисперсій (середніх квадратів |
Значення дисперсій: |
||
|
відхилень): |
σзаг2; σфакт2; σост2 |
||
|
• |
загальна: σзаг2= Dy / (N-1); |
|
|
|
• |
факторна: σфакт2= Dx / (a-1), |
|
|
|
|
де а – кількість груп; |
|
|
|
• |
остаточна: σост2= Dz / (N-a) |
|
|
7 |
Визначення фактичного значення критерію |
Значення фактичного |
||
|
|
|
критерію F |
|
8 |
Порівняння фактичного значення критерію F з |
Висновки про відхилення |
||
|
його табличним (стандартним) значенням для |
або визнання нульової |
||
|
відповідного рівня значимості (р) і даних чисел |
гіпотези |
||
|
|
|
|
|
