Методи дисперсного аналізу. Множинне порівняння: процедура Тьюкі-Крамера

Усі люди від природи прагнуть знання. (Арістотель. Метафізика)

Дисперсійний аналіз

Вступний огляд

У цьому розділі ми розглянемо основні методи, припущення та термінологію дисперсійного аналізу.

Зазначимо, що у англомовної літературі дисперсійний аналіз зазвичай називається аналізом варіації. Тому, для стислості, нижче ми іноді використовуватимемо термін ANOVA (An alysis o f va riation) для звичайного дисперсійного аналізу та термін MANOVAдля багатовимірного дисперсійного аналізу У цьому розділі ми послідовно розглянемо основні ідеї дисперсійного аналізу ( ANOVA), коваріаційного аналізу ( ANCOVA), багатовимірного дисперсійного аналізу ( MANOVA) та багатовимірного коваріаційного аналізу ( MANCOVA). Після короткого обговорення переваг аналізу контрастів і апостеріорних критеріїв розглянемо припущення, на яких засновані методи дисперсійного аналізу. Ближче до кінця цього розділу пояснюються переваги багатовимірного підходу для аналізу повторних вимірів, порівняно з традиційним одновимірним підходом.

Основні ідеї

Ціль дисперсійного аналізу.Основною метою дисперсійного аналізу є дослідження значущості різницю між середніми. Глава (Глава 8) містить коротке запровадження дослідження статистичної значимості. Якщо ви просто порівнюєте середні у двох вибірках, дисперсійний аналіз дасть той самий результат, що й звичайний t- критерій для незалежних вибірок (якщо порівнюються дві незалежні групи об'єктів чи спостережень) або t- критерій для залежних вибірок (якщо порівнюються дві змінні одному й тому безлічі об'єктів чи спостережень). Якщо ви мало знайомі з цими умовами, радимо звернутися до вступного огляду глави (Глава 9).

Звідки походить назва Дисперсійний аналіз? Може здатися дивним, що порівняння середніх називається дисперсійним аналізом. Насправді це пов'язано з тим, що при дослідженні статистичної значущості відмінності між середніми ми насправді аналізуємо дисперсії.

Розбиття суми квадратів

Для вибірки обсягу n вибіркова дисперсія обчислюється як сума квадратів відхилень від середнього вибіркового, поділена на n-1 (обсяг вибірки мінус одиниця). Таким чином, при фіксованому обсязі вибірки n дисперсія є функція суми квадратів (відхилень), що позначається для стислості, SS(Від англійської Sum of Squares - Сума Квадратів). В основі дисперсійного аналізу лежить поділ (або розбиття) дисперсії на частини. Розглянемо наступний набір даних:

Середні дві групи істотно різні (2 і 6 відповідно). Сума квадратів відхилень всерединікожної групи дорівнює 2. Складаючи їх, отримуємо 4. Якщо тепер повторити ці обчислення без урахуваннягрупової приналежності, тобто якщо обчислити SSвиходячи із загального середнього цих двох вибірок, то отримаємо 28. Іншими словами, дисперсія (сума квадратів), заснована на внутрішньогруповій мінливості, призводить до набагато менших значень, ніж при обчисленні на основі загальної мінливості (щодо загальної середньої). Причина цього, очевидно, полягає у суттєвій різниці між середніми значеннями, і ця різниця між середніми і пояснює існуючу різницю між сумами квадратів. Справді, якщо використовувати для аналізу наведених даних модуль Дисперсійний аналіз, будуть отримані такі результати:

Як видно з таблиці, загальна сума квадратів SS=28 розбита у сумі квадратів, обумовлену внутрішньогруповиймінливістю ( 2+2=4 ; див. другий рядок таблиці) та суму квадратів, обумовлену різницею середніх значень. (28-(2+2)=24; див. перший рядок таблиці).

SS помилок таSS ефект.Внутрішньогрупова мінливість ( SS) зазвичай називається дисперсією помилки.Це означає, що зазвичай під час проведення експерименту вона може бути передбачена чи пояснена. З іншого боку, SS ефекту(або міжгрупову мінливість) можна пояснити різницею між середніми значеннями в групах, що вивчаються. Іншими словами, приналежність до певної групи пояснюєміжгрупову мінливість, т.к. нам відомо, що ці групи мають різні середні значення.

Перевірка важливості.Основні ідеї перевірки статистичної значимості обговорюються у розділі Елементарні поняття статистики(Глава 8). У цьому розділі пояснюються причини, через які багато критеріїв використовують ставлення поясненої і непоясненої дисперсії. Приклад такого використання є сам дисперсійний аналіз. Перевірка значущості в дисперсійному аналізі полягає в порівнянні дисперсії, обумовленої міжгруповим розкидом (названої середнім квадратом ефектуабо MSефект) та дисперсії, обумовленої внутрішньогруповим розкидом (названою середнім квадратом помилкиабо MSпомилка). Якщо вірна нульова гіпотеза (рівність середніх у двох популяціях), можна очікувати порівняно невелике різницю у вибіркових середніх через випадкової мінливості. Тому при нульовій гіпотезі внутрішньогрупова дисперсія практично співпадатиме із загальною дисперсією, підрахованою без урахування групою належності. Отримані внутрішньогрупові дисперсії можна порівняти за допомогою F- критерію, що перевіряє, чи справді відношення дисперсій значно більше 1. У розглянутому вище прикладі F- критерій показує, що різницю між середніми статистично значимо.

Основна логіка дисперсійного аналізу.Підсумовуючи, можна сказати, що метою дисперсійного аналізу є перевірка статистичної значущості різниці між середніми (для груп чи змінних). Ця перевірка проводиться з допомогою аналізу дисперсії, тобто. за допомогою розбиття загальної дисперсії (варіації) на частини, одна з яких обумовлена ​​випадковою помилкою (тобто внутрішньогруповою мінливістю), а друга пов'язана з різницею середніх значень. Остання компонент дисперсії потім використовується для аналізу статистичної значущості відмінності між середніми значеннями. Якщо ця відмінність значуща, нульова гіпотеза відкидається і приймається альтернативна гіпотеза існування різниці між середніми.

Залежні та незалежні змінні.Змінні, значення яких визначається за допомогою вимірювань у ході експерименту (наприклад, бал, набраний під час тестування), називаються залежнимизмінними. Змінні, якими можна керувати під час проведення експерименту (наприклад, методи навчання чи інші критерії, що дозволяють розділити спостереження групи) називаються факторамиабо незалежнимизмінними. Докладніше ці поняття описані у розділі Елементарні поняття статистики(Глава 8).

Багатофакторний дисперсійний аналіз

У розглянутому вище простому прикладі ви могли б відразу обчислити t-критерій для незалежних вибірок, використовуючи відповідну опцію модуля Основні статистики та таблиці.Отримані результати, звісно, ​​збігатимуться з результатами дисперсійного аналізу. Однак дисперсійний аналіз містить гнучкі та потужні технічні засоби, які можуть бути використані для більш складних досліджень.

Безліч факторів.Світ за своєю природою складний і багатовимірний. Ситуації, коли деяке явище повністю описується однією змінною, надзвичайно рідкісні. Наприклад, якщо ми намагаємося навчитися вирощувати великі помідори, слід розглядати фактори, пов'язані з генетичною структурою рослин, типом ґрунту, освітленістю, температурою тощо. Таким чином, при проведенні типового експерименту доводиться мати справу з великою кількістю факторів. Основна причина, через яку використання дисперсійного аналізу краще повторного порівняння двох вибірок при різних рівнях факторів за допомогою t- критерію, полягає в тому, що дисперсійний аналіз більш ефективнийі для малих вибірок, більш інформативний.

Управління факторами.Припустимо, що у розглянутому вище прикладі аналізу двох вибірок ми додамо ще один фактор, наприклад, Стать- Gender. Нехай кожна група складається з 3 чоловіків та 3 жінок. План цього експерименту можна подати у вигляді таблиці 2 на 2:

До проведення обчислень можна помітити, що в цьому прикладі загальна дисперсія має принаймні три джерела:

(1) випадкова помилка (внутрішньогрупова дисперсія),

(2) мінливість, пов'язана з приналежністю до експериментальної групи, та

(3) мінливість, обумовлена ​​статтю об'єктів спостереження.

(Зазначимо, що існує ще одне можливе джерело мінливості – взаємодія факторів, який ми обговоримо пізніше). Що станеться, якщо ми не включатимемо стать –genderяк фактор при проведенні аналізу та обчислимо звичайний t-Критерій? Якщо ми обчислюватимемо суми квадратів, ігноруючи стать -gender(тобто об'єднуючи об'єкти різної статі в одну групу при обчисленні внутрішньогрупової дисперсії, отримавши при цьому суму квадратів для кожної групи рівну SS=10, і загальну суму квадратів SS= 10+10 = 20), то отримаємо більше значення внутрішньогрупової дисперсії, ніж при більш точному аналізі з додатковим розбиттям на підгрупи підлозі - gender(при цьому внутрішньогрупові середні дорівнюватимуть 2, а загальна внутрішньогрупова сума квадратів дорівнює SS = 2+2+2+2 = 8). Ця різниця пов'язана з тим, що середнє значення для чоловіків - malesменше, ніж середнє значення для жінок –female, і це різницю у середніх значеннях збільшує сумарну внутригрупповую мінливість, якщо чинник статі не враховується. Управління дисперсією помилки збільшує чутливість (потужність) критерію.

На цьому прикладі видно ще одну перевагу дисперсійного аналізу порівняно зі звичайним t-Крітерієм для двох вибірок. Дисперсійний аналіз дозволяє вивчати кожен чинник, керуючи значеннями інших чинників. Це насправді і є основною причиною його більшої статистичної потужності (для отримання значних результатів потрібні менші обсяги вибірок). Тому дисперсійний аналіз навіть на невеликих вибірках дає статистично більш значущі результати, ніж простий t- критерій.

Ефекти взаємодії

Існує ще одна перевага застосування дисперсійного аналізу порівняно із звичайним t- критерієм: дисперсійний аналіз дозволяє виявити взаємодіяміж факторами і, отже, дозволяє вивчати складніші моделі. Для ілюстрації розглянемо ще один приклад.

Головні ефекти, попарні (двофакторні) взаємодії.Припустимо, що є дві групи студентів, причому психологічно студенти першої групи налаштовані виконання поставлених завдань і більш цілеспрямовані, ніж студенти другої групи, що з більш лінивих студентів. Розіб'ємо кожну групу випадково навпіл і запропонуємо одній половині в кожній групі складне завдання, а іншій - легке. Після цього виміряємо, як напружено студенти працюють над цими завданнями. Середні значення для цього (вигаданого) дослідження показані в таблиці:

Який висновок можна зробити із цих результатів? Чи можна зробити висновок, що: (1) над складним завданням студенти працюють більш напружено; (2) цілеспрямовані студенти працюють наполегливіше, ніж ліниві? Жодне з цих тверджень не відбиває сутність систематичного характеру середніх, наведених у таблиці. Аналізуючи результати, правильніше було б сказати, що над складними завданнями працюють наполегливіше лише цілеспрямовані студенти, тоді як над легкими завданнями тільки ліниві працюють наполегливіше. Тобто характер студентів та складність завдання взаємодіючиміж собою впливають на витрачається зусилля. Це приклад парної взаємодіїміж характером студентів та складністю завдання. Зазначимо, що твердження 1 та 2 описують головні ефекти.

Взаємодія вищих порядків.У той час, як пояснити попарні взаємодії ще порівняно легко, взаємодії вищих порядків пояснити значно складніше. Уявімо, що в аналізований вище приклад введено ще один фактор стать -Genderі ми отримали наступну таблицю середніх значень:

Які висновки можна зробити з отриманих результатів? Графіки середніх дозволяють легко інтерпретувати складні ефекти. Модуль дисперсійного аналізу дозволяє будувати ці графіки практично одним клацанням мишки.

Зображення на графіках внизу являє собою трифакторну взаємодію, що вивчається.

Дивлячись на графіки, можна сказати, що у жінок існує взаємодія між характером та складністю тесту: цілеспрямовані жінки працюють над важким завданням більш напружено, ніж над легким. У чоловіків ця ж взаємодія має зворотний характер. Видно, що опис взаємодії між факторами стає більш заплутаним.

Загальний спосіб опису взаємодій.У випадку взаємодія між чинниками описується як зміни одного ефекту під впливом іншого. У розглянутому вище прикладі двофакторну взаємодію можна описати як зміну головного ефекту фактора, що характеризує складність завдання під впливом фактора, що описує характер студента. Для взаємодії трьох факторів із попереднього параграфа можна сказати, що взаємодія двох факторів (складності завдання та характеру студента) змінюється під впливом статі – Gender. Якщо вивчається взаємодія чотирьох чинників, можна сказати, взаємодія трьох чинників, змінюється під впливом четвертого чинника, тобто. існують різні типи взаємодій різних рівнях четвертого чинника. Виявилося, що в багатьох областях взаємодія п'яти чи навіть більшої кількості факторів не є чимось незвичайним.

Складні плани

Міжгрупові та внутрішньогрупові плани (плани з повторними вимірами)

При порівнянні двох різних груп зазвичай використовується t- критерій для незалежних вибірок (з модуля Основні статистики та таблиці). Коли порівнюються дві змінні на тому самому безлічі об'єктів (спостережень), використовується t-Критерій для залежних вибірок. Для дисперсійного аналізу також важливо залежні чи ні вибірки. Якщо є повторні вимірювання тих самих змінних (за різних умов або в різний час) для тих самих об'єктів, то говорять про наявність фактора повторних вимірів(називається також внутрішньогруповим фактором,оскільки з оцінки його значимості обчислюється внутригрупповая сума квадратів). Якщо порівнюються різні групи об'єктів (наприклад, чоловіки та жінки, три штами бактерій тощо), то різниця між групами описується міжгруповий фактор.Способи обчислення критеріїв значущості для двох описаних типів факторів різні, але їх загальна логіка та інтерпретації збігається.

Між- та внутрішньогрупові плани.У багатьох випадках експеримент вимагає включення до плану і міжгрупового фактора, і фактора повторних вимірів. Наприклад, вимірюються математичні навички студентів жіночої та чоловічої статі (де стать -Gender-міжгруповий фактор) на початку та наприкінці семестру. Два виміри навичок кожного студента утворюють внутрішньогруповий фактор (фактор повторних вимірів). Інтерпретація головних ефектів та взаємодій для міжгрупових факторів та факторів повторних вимірювань збігається, і обидва типи факторів можуть, очевидно, взаємодіяти між собою (наприклад, жінки набувають навичок протягом семестру, а чоловіки їх втрачають).

Неповні (гніздові) плани

У багатьох випадках можна знехтувати ефектом взаємодії. Це відбувається або коли відомо, що у популяції ефект взаємодії відсутній, або коли здійснення повного факторногоплану неможливо. Наприклад, вивчається вплив чотирьох добавок до палива на витрату пального. Вибираються чотири автомобілі та чотири водії. Повний факторнийЕксперимент вимагає, щоб кожна комбінація: добавка, водій, автомобіль – з'явилися хоча б один раз. Для цього потрібно не менше 4 x 4 x 4 = 64 груп випробувань, що потребує надто великих часових витрат. Крім того, навряд чи існує взаємодія між водієм та добавкою до палива. Зважаючи на це, можна використовувати план Латинські квадрати,в якому міститься лише 16 груп випробувань (чотири добавки позначаються буквами A, B, C та D):

Латинські квадрати описані в більшості книг з планування експериментів (наприклад, Hays, 1988; Lindman, 1974; Milliken and Johnson, 1984; Winer, 1962), і тут вони не будуть детально обговорюватися. Зазначимо, що латинські квадрати це неnолніплани, у яких беруть участь в повному обсязі комбінації рівнів чинників. Наприклад, водій 1 керує автомобілем 1 тільки з добавкою А водій 3 керує автомобілем 1 тільки з добавкою С. Рівні фактора добавок ( A, B, C і D) вкладені в комірки таблиці автомобіль x водій –як яйця в гнізда. Це мнемонічне правило корисне для розуміння природи гніздових чи вкладенихпланів. Модуль Дисперсійний аналізнадає прості методи аналіз планів такого типу.

Коваріаційний аналіз

Основна ідея

В розділі Основні ідеїкоротко обговорювалася ідея управління факторами та те, яким чином включення адитивних факторів дозволяє зменшувати суму квадратів помилок та збільшувати статистичну потужність плану. Все це може бути поширене і на змінні з безперервним безліччю значень. Коли такі безперервні змінні включаються в план як фактори, вони називаються підступами.

Фіксовані коваріати

Припустимо, що порівнюються математичні навички двох груп студентів, які навчалися за двома різними підручниками. Припустимо, що є дані про коефіцієнт інтелекту (IQ) для кожного студента. Можна припустити, що коефіцієнт інтелекту пов'язаний з математичними навичками та використовувати цю інформацію. Для кожної з двох груп студентів можна визначити коефіцієнт кореляції між IQ і математичними навичками. Використовуючи цей коефіцієнт кореляції, можна виділити частку дисперсії в групах, що пояснюється впливом IQ і нез'ясовну частку дисперсії (див. також Елементарні поняття статистики(глава 8) та Основні статистики та таблиці(Глава 9)). Частка дисперсії, що залишилася, використовується при проведенні аналізу як дисперсія помилки. Якщо є кореляція між IQ та математичними навичками, то можна суттєво зменшити дисперсії помилки SS/(n-1) .

Вплив коваріат наF- критерій. F-критерій оцінює статистичну значущість відмінності середніх значень у групах, при цьому обчислюється відношення міжгрупової дисперсії ( MSефект) до дисперсії помилок ( MSerror) . Якщо MSerrorзменшується, наприклад, при врахуванні фактора IQ, значення Fзбільшується.

Безліч коваріат.Міркування, використані вище для однієї кваріати (IQ), легко поширюються на кілька коваріат. Наприклад, крім IQ, можна включити вимір мотивації, просторового мислення тощо. Замість звичайного коефіцієнта кореляції у своїй використовується множинний коефіцієнт кореляції.

Коли значенняF -Критерію зменшується.Іноді введення коваріату в план експерименту зменшує значення F-критерія . Зазвичай це вказує на те, що коваріати корелюються не тільки із залежною змінною (наприклад, математичними навичками), а й з факторами (наприклад, із різними підручниками). Припустимо, що IQ вимірюється наприкінці семестру, після майже річного навчання двох груп студентів за двома різними підручниками. Хоча студенти розбивалися на групи випадковим чином, може виявитися, що відмінність підручників настільки велика, що і IQ та математичні навички у різних групах сильно відрізнятимуться. У цьому випадку, каварыати не лише зменшують дисперсію помилок, а й міжгрупову дисперсію. Іншими словами, після контролю за різницею IQ у різних групах, різниця в математичних навичках вже буде несуттєвою. Можна сказати інакше. Після “виключення” впливу IQ, ненавмисно виключається вплив підручника на розвиток математичних навичок.

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

Взаємодія між коваріатами та факторами.Також як досліджується взаємодія між факторами, можна досліджувати взаємодію між коваріатами та між групами факторів. Припустимо, що один із підручників особливо підходить для розумних студентів. Другий підручник для розумних студентів з'їде, а для менш розумних студентів цей же підручник важкий. В результаті є позитивна кореляція між IQ і результатом навчання в першій групі (розумніші студенти, краще результат) і нульова або невелика негативна кореляція в другій групі (чим розумніший студент, тим менш ймовірне придбання математичних навичок з другого підручника). У деяких дослідженнях ця ситуація обговорюється як приклад порушення припущень коварійного аналізу. Однак так як у модулі Дисперсійний аналіз використовуються найзагальніші способи коваріаційного аналізу, можна, зокрема, оцінити статистичну значущість взаємодії між факторами та коваріатами.

Змінні коваріати

У той час, як фіксовані кваріати обговорюються в підручниках досить часто, змінні кваріати згадуються набагато рідше. Зазвичай, під час проведення експериментів з повторними вимірами, нас цікавлять розбіжності у вимірах тих самих величин у різні моменти часу. Зокрема, нас цікавить значимість цих відмінностей. Якщо одночасно з вимірюваннями залежних змінних проводиться вимірювання коваріат, можна обчислити кореляцію між коваріатою та залежною змінною.

Наприклад, можна вивчати інтерес до математики та математичні навички на початку та в кінці семестру. Цікаво було б перевірити, чи корельовані між собою зміни на користь математики зі зміною математичних навичок.

Модуль Дисперсійний аналізв STATISTICAавтоматично оцінює статистичну значущість зміни коваріат у тих планах, де це можливо.

Багатомірні плани: багатовимірний дисперсійний та коварійний аналіз

Міжгрупові плани

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

Наприклад, проводиться дослідження двох різних підручників. При цьому вивчаються успіхи студентів у вивченні фізики та математики. У цьому випадку є дві залежні змінні і потрібно з'ясувати, як впливають на них одночасно два різні підручники. Для цього можна скористатися багатовимірним дисперсійним аналізом (MANOVA). Замість одновимірного Fкритерію, використовується багатовимірний Fкритерій (l-критерій Вілкса), заснований на порівнянні матриці коваріаційної матриці помилок і міжгрупової матриці коваріаційної.

Якщо залежні змінні корелированы між собою, це кореляція повинна враховуватися при обчисленні критерію значимості. Очевидно, якщо один і той самий вимір повторюється двічі, то нічого нового отримати при цьому не можна. Якщо до наявного виміру додається корельований з ним вимір, то виходить деяка нова інформація, але при цьому нова змінна містить надмірну інформацію, яка відображається в коваріації між змінними.

Інтерпретація результатів.Якщо загальний багатовимірний критерій значимий, можна зробити висновок, що відповідний ефект (наприклад, тип підручника) значимий. Однак постають такі питання. Чи впливає тип підручника на покращення лише математичних навичок, лише фізичних навичок, або одночасно на покращення тих та інших навичок. Насправді, після отримання значущого багатовимірного критерію для окремого головного ефекту або взаємодії досліджується одномірний Fкритерій. Іншими словами, окремо досліджуються залежні змінні, які роблять внесок у значущість багатовимірного критерію.

Плани з повторними вимірами

Якщо вимірюються математичні та фізичні навички студентів на початку семестру та наприкінці, то це і є повторні виміри. Вивчення критерію значущості у таких планах це логічний розвиток одновимірного випадку. Зауважимо, що методи багатовимірного дисперсійного аналізу зазвичай також використовуються для дослідження значущості одновимірних факторів повторних вимірів, що мають більш як два рівні. Відповідні застосування будуть розглянуті пізніше у цій частині.

Підсумовування значень змінних та багатовимірний дисперсійний аналіз

Навіть досвідчені користувачі одновимірного та багатовимірного дисперсійного аналізу часто утрудняються, отримуючи різні результати при застосуванні багатовимірного дисперсійного аналізу, наприклад, для трьох змінних, і при застосуванні одновимірного дисперсійного аналізу до суми цих трьох змінних, як до однієї змінної.

Ідея підсумовуваннязмінних у тому, кожна змінна містить у собі деяку істинну змінну, що й досліджується, і навіть випадкову помилку виміру. Тому при усередненні значень змінних помилка вимірювання буде ближче до 0 для всіх вимірювань і усереднене значень буде більш надійним. Насправді, в цьому випадку застосування дисперсійного аналізу до суми змінних є розумним і є потужним методом. Однак якщо залежні змінні за своєю природою багатовимірні, підсумовування значень змінних є недоречним.

Наприклад, нехай залежні змінні складаються з чотирьох показників успіху у суспільстві. Кожен показник характеризує незалежну сторону людської діяльності (наприклад, професійний успіх, успішність у бізнесі, сімейний добробут тощо). Додавання цих змінних подібне до додавання яблука і апельсина. Сума цих змінних не буде відповідним одновимірним показником. Тому з такими даними потрібно поводитися як з багатовимірними показниками багатовимірному дисперсійному аналізі.

Аналіз контрастів та апостеріорні критерії

Чому порівнюються окремі множини середніх?

Зазвичай гіпотези щодо експериментальних даних формулюються непросто у термінах основних ефектів чи взаємодій. Прикладом може бути така гіпотеза: деякий підручник підвищує математичні навички лише в студентів чоловічої статі, тоді як інший підручник приблизно однаково ефективний обох статей, проте менш ефективний чоловікам. Можна передбачити, що ефективність підручника взаємодіє зі статтю студента. Однак цей прогноз стосується також природивзаємодії. Очікується значне різницю між статями для учнів з однієї книжці і майже залежні від статі результати для які у інших книжці. Такий тип гіпотез зазвичай досліджується за допомогою аналізу контрастів.

Аналіз контрастів

Якщо говорити коротко, аналіз контрастів дозволяє оцінювати статистичну значимість деяких лінійних комбінацій ефектів складного плану. Аналіз контрастів є головним і обов'язковим елементом будь-якого складного плану дисперсійного аналізу. Модуль Дисперсійний аналізмає досить різноманітні можливості аналізу контрастів, які дозволяють виділяти та аналізувати будь-які типи порівнянь середніх.

Апостеріорніпорівняння

Іноді внаслідок обробки експерименту виявляється несподіваний ефект. Хоча у більшості випадків творчий дослідник зможе пояснити будь-який результат, це не дає можливостей для подальшого аналізу та отримання оцінок для прогнозу. Ця проблема є однією з тих, для яких використовуються апостеріорні критерії, тобто критерії, які не використовують апріорнігіпотези. Для ілюстрації розглянемо такий експеримент. Припустимо, що у 100 картках записані числа від 1 до 10. Опустивши всі ці картки в шапку, ми випадково вибираємо 20 разів по 5 карток, і обчислюємо кожної вибірки середнє значення (середнє чисел, записаних на картки). Чи можна очікувати, що знайдуться дві вибірки, у яких середні значення значно відрізняються? Це дуже правдоподібно! Вибираючи дві вибірки з максимальним і мінімальним середнім, можна отримати різницю середніх, що відрізняється від різниці середніх, наприклад, перших двох вибірок. Цю різницю можна дослідити, наприклад, за допомогою аналізу контрастів. Якщо не вдаватися в деталі, то існує кілька так званих апостеріорнихкритеріїв, які засновані в точності на першому сценарії (взяття екстремальних середніх із 20 вибірок), тобто ці критерії засновані на виборі найбільш відмінних середніх для порівняння всіх середніх значень у плані. Ці критерії застосовуються для того, щоб суто випадково не отримати штучний ефект, наприклад, виявити значну різницю між середніми, коли його немає. Модуль Дисперсійний аналізпропонує широкий вибір таких критеріїв. Коли в експерименті, пов'язаному з кількома групами, трапляються несподівані результати, то використовуються апостеріорніпроцедури на дослідження статистичної значимості отриманих результатів.

Сума квадратів типу I, II, III та IV

Багатомірна регресія та дисперсійний аналіз

Існує тісний взаємозв'язок між методом багатовимірної регресії та дисперсійним аналізом (аналізом варіацій). І в тому, і в іншому методі досліджується лінійна модель. Якщо говорити коротко, то практично всі плани експерименту можна досліджувати за допомогою багатовимірної регресії. Розглянемо наступний простий міжгруповий 2 x 2 план.

Стовпці А та В містять коди, що характеризують рівні факторів А та В, стовпець АxВ містить добуток двох стовпців А та В. Ми можемо аналізувати ці дані за допомогою багатовимірної регресії. Змінна DVвизначається як залежна змінна, змінні від Aдо AxBяк незалежні змінні. Дослідження значущості для коефіцієнтів регресії співпадатиме з обчисленнями у дисперсійному аналізі значимості головних ефектів факторів Aі Bта ефекту взаємодії AxB.

Незбалансовані та збалансовані плани

При обчисленні кореляційної матриці для всіх змінних, наприклад для даних, зображених вище, можна помітити, що головні ефекти факторів Aі Bта ефект взаємодії AxBнекорельовані. Цю властивість ефектів називають також ортогональністю. Говорять, що ефекти Aі B - ортогональніабо незалежніодин від одного. Якщо всі ефекти в плані ортогональні один одному, як у наведеному вище прикладі, то кажуть, що план збалансований.

Збалансовані плани мають “хорошу властивість”. Обчислення під час аналізу таких планів дуже прості. Усі обчислення зводяться до обчислення кореляції між ефектами та залежними змінними. Так як ефекти ортогональні, приватні кореляції (як у повній багатовимірноїрегресії) не обчислюються. Однак у реальному житті плани не завжди збалансовані.

Розглянемо реальні дані з нерівним числом спостережень у осередках.

Якщо закодувати ці дані як вище та обчислити кореляційну матрицю для всіх змінних, то виявиться, що фактори плану корелюються один з одним. Фактори в плані тепер не ортогональні і такі плани називаються незбалансованими.Зауважимо, що в прикладі, що розглядається, кореляція між факторами повністю пов'язана з відмінністю частот 1 і -1 в стовпцях матриці даних. Іншими словами, плани експериментів з нерівними обсягами осередків (точніше, непропорційними обсягами) будуть незбалансованими, це означає, що головні ефекти та взаємодії змішуватимуться. І тут для обчислення статистичної значущості ефектів необхідно повністю обчислювати багатовимірну регресію. Тут є кілька стратегій.

Сума квадратів типу I, II, III та IV

Сума квадратів типуIіIII. Для вивчення значущості кожного фактора в багатовимірній моделі можна обчислювати приватну кореляцію кожного фактора за умови, що всі інші фактори вже враховані в моделі. Можна також вводити фактори в модель покроковим способом, фіксуючи всі фактори, що вже введені в модель і ігноруючи всі інші фактори. Взагалі, в цьому і полягає різниця між типом IIIі типомIсуми квадратів (ця термінологія була введена в SAS, див. наприклад, SAS, 1982; детальне обговорення можна також знайти в Searle, 1987, стор. 461; Woodward, Bonett, and Brecht, 1990, стор. 216; або Milliken and Johnson, 1984, стор 138).

Сума квадратів типуІІ.Наступна “проміжна” стратегія формування моделі полягає: у контролі всіх основних ефектів щодо значимості окремого головного ефекту; у контролі всіх основних ефектів та всіх попарних взаємодій, коли досліджується значимість окремої попарної взаємодії; у контролі всіх основних ефектів всіх попарних взаємодій та всіх взаємодій трьох факторів; щодо окремого взаємодії трьох чинників тощо. Суми квадратів для ефектів, що обчислюються таким способом, називаються типомIIсуми квадратів. Отже, типIIсуми квадратів контролює всі ефекти того ж порядку та нижче, ігноруючи всі ефекти вищого порядку.

Сума квадратів типуIV. Нарешті, для деяких спеціальних планів із пропущеними осередками (неповними планами) можна обчислювати, так звані, типу IVсуми квадратів. Цей метод обговорюватиметься пізніше у зв'язку з неповними планами (планами з пропущеними осередками).

Інтерпретація гіпотези про суму квадратів типу I, II та III

Суму квадратів типуIIIнайлегше інтерпретувати. Нагадаємо, що суми квадратів типуIIIдосліджують ефекти після контролю всіх інших ефектів. Наприклад, після знаходження статистично значущого типуIIIефекту для фактора Aу модулі Дисперсійний аналіз, можна сказати, що існує єдиний значущий ефект фактора Aпісля введення всіх інших ефектів (факторів) і відповідно інтерпретувати цей ефект. Ймовірно, у 99% усіх додатків дисперсійного аналізу саме цей тип критерію цікавить дослідника. Цей тип суми квадратів зазвичай обчислюється в модулі Дисперсійний аналізза замовчуванням, незалежно від того вибрано опцію Регресійний підхідчи ні (стандартні підходи прийняті в модулі Дисперсійний аналізобговорюються нижче).

Значні ефекти, отримані за допомогою сум квадратів типуабо типуIIсуми квадратів інтерпретувати не так просто. Найкраще їх інтерпретувати в контексті покрокової багатовимірної регресії. Якщо при використанні суми квадратів типуIголовний ефект фактора В виявився значимим (після включення в модель фактора А, але перед додаванням взаємодії між А і В), можна зробити висновок, що існує значний головний ефект фактора В, за умови, що немає взаємодії між факторами А і В. (Якщо при використання критерію типуIII, фактор В також виявився значним, то можна зробити висновок, що існує значний головний ефект фактора B, після введення в модель всіх інших факторів та їх взаємодій).

У термінах маргінальних середніх гіпотези типуIі типуIIзазвичай немає простої інтерпретації. У таких випадках кажуть, що не можна інтерпретувати значущість ефектів, розглядаючи лише маргінальні середні. Швидше представлені pзначень середніх мають відношення до складної гіпотези, яка комбінує середні та обсяг вибірки. Наприклад, типIIгіпотези для фактора А в простому прикладі плану 2 x 2, що раніше розглядаються (див. Woodward, Bonett, and Brecht, 1990, стор 219):

nij- Число спостережень в осередку

uij- Середнє значення в осередку

n. j- маргінальне середнє

Якщо не вдаватися в деталі (детальніше див. Milliken and Johnson, 1984, глава 10), то ясно, що це не прості гіпотези і в більшості випадків жодна з них не має особливого інтересу у дослідника. Однак існують випадки, коли гіпотези типуIможуть бути цікавими.

Обчислювальний підхід у модулі, що приймається за умовчанням Дисперсійний аналіз

За замовчуванням, якщо не зазначено опцію Регресійний підхід, модуль Дисперсійний аналізвикористовує модель середніх по осередках. Для цієї моделі характерно, що суми квадратів для різних ефектів обчислюються для лінійних комбінацій середніх значень осередків. У повному факторному експерименті це призводить до сум квадратів, які збігаються з сумами квадратів, які раніше обговорювали як тип III. Однак у опції Сплановані порівняння(у вікні Результати дисперсійного аналізу), користувач може перевіряти гіпотезу щодо будь-якої лінійної комбінації зважених або незважених середніх по осередках. Таким чином, користувач може перевіряти не лише гіпотези типуIIIале гіпотези будь-якого типу (включаючи типIV). Цей загальний підхід є особливо корисним, коли досліджуються плани з пропущеними осередками (так звані неповні плани).

Для повних факторних планів цей підхід корисно також використовувати у випадках, коли хочуть аналізувати зважені маргінальні середні. Наприклад, припустимо, що в аналізованому раніше простому 2 x 2 плані, потрібно порівняти виважені (за рівнями фактора B) маргінальні середні для фактора А. Це буває корисним, коли розподіл спостережень по осередках не готувалося експериментатором, а будувалося випадково, і ця випадковість відображається у розподілі числа спостережень за рівнями фактора B у сукупності.

Наприклад, є фактор – вік вдів. Можлива вибірка респондентів розбита на дві групи: молодше 40 років та старше 40 (фактор В). Другий чинник (фактор А) у плані - отримували чи ні соціальну підтримку вдови у певному агентстві (при цьому одні вдови були обрані випадково, інші служили як контроль). У цьому випадку розподіл удів за віком у вибірці відображає дійсний розподіл вдів за віком у сукупності. Оцінка ефективності групи соціальної підтримки вдів з всім вікомбуде відповідати виважене середнє для двох вікових груп (з вагами, що відповідають числу спостережень у групі).

Сплановані порівняння

Зауважимо, що сума запроваджених коефіцієнтів контрастів не обов'язково дорівнює 0 (нулю). Натомість програма автоматично вносити поправки, щоб відповідні гіпотези не змішувалися із загальним середнім.

Для ілюстрації цього повернемося знову до простого 2 x 2 плану, розглянутого раніше. Нагадаємо, що числа спостережень у осередках цього незбалансованого плану -1, 2, 3, і 1. Припустимо, що ми хочемо порівняти зважені середні маргінальні для фактора А (зважені з частотою рівнів фактора В). Можна ввести коефіцієнти розмаїття:

Зауважимо, що ці коефіцієнти не дають у сумі 0. Програма встановлюватиме коефіцієнти так, що в сумі вони даватиму 0, і при цьому зберігатимуться їх відносні значення, тобто:

1/3 2/3 -3/4 -1/4

Ці контрасти порівнюватимуть зважені середні для фактора А.

Гіпотези про головне середнє.Гіпотеза, у тому, що ні зважене головне середнє дорівнює 0 може досліджуватися з допомогою коефіцієнтів:

Гіпотеза про те, що зважене головне середнє 0 перевіряється за допомогою:

У жодному разі програма не здійснює коригування коефіцієнтів контрастів.

Аналіз планів із пропущеними осередками (неповні плани)

Факторні плани, що містять порожні осередки (обробка комбінацій осередків, у яких немає спостережень), називаються неповними. У таких планах деякі фактори зазвичай не ортогональні і деякі взаємодії не можуть бути обчислені. Взагалі немає кращого методу аналізу таких планів.

Регресійний підхід

У деяких старих програмах, які ґрунтуються на аналізі планів дисперсійного аналізу за допомогою багатовимірної регресії, фактори в неповних планах за замовчуванням задаються звичайним чином (начебто план повний). Потім проводиться багатовимірний регресійний аналіз цих фіктивно закодованих чинників. На жаль, цей метод призводить до результатів, які дуже важко, або навіть неможливо, інтерпретувати, оскільки неясно, як кожен ефект бере участь у лінійній комбінації середніх значень. Розглянемо наступний приклад.

Якщо виконуватиметься багатовимірна регресія виду Залежна змінна = Константа + Фактор A + Фактор B, то гіпотеза про значимість факторів A і B в термінах лінійних середніх комбінацій виглядає так:

Фактор A: Осередок A1, B1 = Осередок A2, B1

Фактор B: Осередок A1, B1 = Осередок A1, B2

Цей випадок простий. У складніших планах неможливо фактично визначити, що точно досліджуватиметься.

Середні осередки, підхід дисперсійного аналізу , гіпотези типу IV

Підхід, який рекомендується в літературі і який здається кращим – дослідження осмислених (з точки зору дослідницьких завдань) апріорнихгіпотез про середні, що спостерігаються в осередках плану. Докладне обговорення цього підходу можна знайти в Dodge (1985), Heiberger (1989), Milliken and Johnson (1984), Searle (1987) або Woodward, Bonett, and Brecht (1990). Суми квадратів, асоційовані з гіпотезами про лінійну комбінацію середніх у неповних планах, що досліджують оцінки частини ефектів, називаються також сумами квадратів IV.

Автоматична генерація гіпотез типуIV. Коли багатофакторні плани мають складний характер пропущених осередків, бажано визначити ортогональні (незалежні) гіпотези, дослідження яких еквівалентне дослідженню основних ефектів чи взаємодій. Були розвинені алгоритмічні (обчислювальні) стратегії (засновані на псевдозворотній матриці плану) для генерування ваги для таких порівнянь. На жаль, остаточні гіпотези визначаються не єдиним чином. Звичайно, вони залежать від порядку, в якому були визначені ефекти і рідко допускають просту інтерпретацію. Тому рекомендується уважно вивчити характер пропущених осередків, потім формулювати гіпотези типуIV, які найбільш змістовно відповідають цілям дослідження. Потім дослідити ці гіпотези, використовуючи опцію Сплановані порівнянняу вікні Результати. Найлегший шлях задати порівняння у цьому випадку – вимагати введення вектора контрастів для всіх факторів разому вікні Сплановані порівняння.Після виклику діалогового вікна Сплановані порівняннябудуть показані всі групи поточного плану та позначені ті, що пропущені.

Пропущені осередки та перевірка специфічного ефекту

Існує кілька типів планів, у яких розташування пропущених осередків невипадково, але ретельно сплановано, що дозволяє проводити простий аналіз головних ефектів не торкаючись інших ефектів. Наприклад, коли необхідна кількість комірок у плані недоступна, часто використовуються плани. Латинські квадратидля оцінювання основних ефектів кількох чинників із великою кількістю рівнів. Наприклад, 4 x 4 x 4 x 4 факторний план потребує 256 осередків. У той же час можна використовувати Греко-латинський квадратдля оцінки головних ефектів, маючи лише 16 осередків у плані (глава Планування експерименту, том IV містить детальний опис таких планів). Неповні плани, в яких головні ефекти (і деякі взаємодії) можуть бути оцінені за допомогою простих лінійних середніх комбінацій, називаються збалансованими неповними планами.

У збалансованих планах стандартний (за замовчуванням) метод генерування контрастів (ваг) для головних ефектів і взаємодій буде проводити аналіз таблиці дисперсій, в якій суми квадратів для відповідних ефектів не змішуються один з одним. Опція Специфічний ефектвікна Результатибуде генерувати пропущені контрасти, записуючи нуль у пропущені комірки плану. Відразу після того, як буде запрошено опцію Специфічний ефектдля користувача, який вивчає деяку гіпотезу, з'являється таблиця результатів із фактичними вагами. Зауважимо, що у збалансованому плані, суми квадратів відповідних ефектів обчислюються тільки, якщо ці ефекти ортогональні (незалежні) всім іншим головним ефектам та взаємодіям. В іншому випадку потрібно скористатися опцією Сплановані порівняннявивчення змістовних порівнянь між середніми.

Пропущені осередки та об'єднані ефекти/члени помилки

Якщо опція Регресійний підхіду стартовій панелі модуля Дисперсійний аналізне вибрано, то при обчисленні суми квадратів для ефектів використовуватиметься модель середніх за комірками (установка за замовчуванням). Якщо план не збалансований, то при поєднанні неортогональних ефектів (див. вище обговорення опції Пропущені осередки та специфічний ефект) можна отримати суму квадратів, що складається з неортогональних (або перекриваються) компонентів. Отримані при цьому результати зазвичай не інтерпретовані. Тому треба бути дуже обережним під час виборів та реалізації складних неповних експериментальних планів.

Існує багато книг із детальним обговоренням планів різного типу. (Dodge, 1985; Heiberger, 1989; Lindman, 1974; Milliken and Johnson, 1984; Searle, 1987; Woodward and Bonett, 1990), але така інформація лежить поза межами цього підручника. Проте пізніше в цьому розділі буде продемонстровано аналіз різного типу планів.

Припущення та ефекти порушення припущень

Відхилення від припущення щодо нормальності розподілів

Припустимо, що залежна змінна виміряна у числовій шкалі. Припустимо також, що залежна змінна має нормальний розподіл усередині кожної групи. Дисперсійний аналізмістить широкий набір графіків та статистик для обґрунтування цього припущення.

Ефекти порушення.Взагалі Fкритерій дуже стійкий до відхилення від нормальності (докладні результати див. у роботі Lindman, 1974). Якщо ексцес більший за 0, то значення статистики Fможе стати дуже маленьким. Нульова гіпотеза у своїй приймається, хоча може бути й неправильна. Ситуація змінюється на протилежну, коли ексцес менший за 0. Асиметрія розподілу зазвичай незначно впливає на Fстатистику. Якщо кількість спостережень у осередку досить велика, то відхилення від нормальності не має особливого значення в силу центральної граничної теореми, відповідно до якої, розподіл середнього значення близький до нормального, незалежно від початкового розподілу. Детальне обговорення стійкості FСтатистики можна знайти в Box and Anderson (1955), або Lindman (1974).

Однорідність дисперсії

Припущення.Передбачається, що дисперсії різних груп плану однакові. Це припущення називається припущенням про однорідності дисперсії.Згадаймо, що на початку цього розділу, описуючи обчислення суми квадратів помилок, ми робили підсумовування всередині кожної групи. Якщо дисперсії у двох групах відрізняються один від одного, то додавання їх не дуже природне і не дає оцінки загальної внутрішньогрупової дисперсії (оскільки в цьому випадку загальної дисперсії взагалі не існує). Модуль Дисперсійний аналіз -ANOVA/MANOVAмістить великий набір статистичних критеріїв виявлення відхилення припущень однорідності дисперсії.

Ефекти порушення.Ліндман (Lindman 1974, стор 33) показує, що Fкритерій цілком стійкий щодо порушення припущень однорідності дисперсії ( неоднорідністьдисперсії, див. також Box, 1954a, 1954b; Hsu, 1938).

Спеціальний випадок: кореленість середніх та дисперсій.Бувають випадки, коли Fстатистика може ввести в оману.Це буває, коли в осередках плану середні значення корелюються з дисперсією. Модуль Дисперсійний аналіздозволяє будувати діаграми розсіювання дисперсії або стандартного відхилення щодо середніх виявлення такої кореляції. Причина, через яку така кореляція небезпечна, полягає в наступному. Уявімо, що є 8 осередків у плані, 7 з яких мають майже однакове середнє, а в одному осередку середнє набагато більше за інших. Тоді Fкритерій може виявити статистично значущий ефект. Але припустимо, що у осередку з великим середнім значенням і дисперсія значно більше інших, тобто. середнє значення і дисперсія в осередках залежні (що більше середнє, то більше вписувалося дисперсія). І тут велике середнє значення ненадійно, оскільки може бути викликано великий дисперсією даних. Однак Fстатистика, заснована на об'єднаноюдисперсії всередині осередків, фіксуватиме велике середнє, хоча критерії, засновані на дисперсії у кожному осередку, в повному обсязі відмінності середніх вважатимуть значимими.

Такий характер даних (велике середнє і велика дисперсія) - часто зустрічається, коли є спостереження, що різко виділяються. Одне або два різко виділяються спостережень сильно зміщують середнє значення і дуже збільшують дисперсію.

Однорідність дисперсії та підступів

Припущення.У багатовимірних планах, з багатовимірними залежними вимірами, також застосовуються припущення однорідності дисперсії, описані раніше. Однак так як існують багатовимірні залежні змінні, то потрібно так само, щоб їх взаємні кореляції (коваріації) були однорідними по всіх осередках плану. Модуль Дисперсійний аналізпропонує різні способи перевірки цих припущень.

Ефекти порушення. Багатовимірний аналог F- критерію - λ-критерій Вілкса. Не так багато відомо про стійкість (робастність) λ-критерію Вілкса щодо порушення зазначених вище припущень. Тим не менш, так як інтерпретація результатів модуля Дисперсійний аналізґрунтується зазвичай на значущості одновимірних ефектів (після встановлення значущості загального критерію), обговорення робастності стосується переважно одномірного дисперсійного аналізу. Тому має бути уважно досліджено значущість одновимірних ефектів.

Спеціальний випадок: підступний аналіз.Особливо серйозні порушення однорідності дисперсії/коваріацій можуть відбуватися, коли до плану включаються коваріати. Зокрема, якщо кореляція між коваріатами і залежними вимірами різна в різних осередках плану, може бути неправильне тлумачення результатів. Слід пам'ятати, що в коваріаційному аналізі, по суті, проводиться регресійний аналіз усередині кожного осередку для того, щоб виділити ту частину дисперсії, яка відповідає коваріату. Припущення про однорідність дисперсії/коваріації передбачає, що цей регресійний аналіз проводиться за наступного обмеження: всі регресійні рівняння (нахили) всім осередків однакові. Якщо це не передбачається, то можуть виникнути великі помилки. Модуль Дисперсійний аналізмає кілька спеціальних критеріїв для перевірки цього припущення. Можна порадити використовувати ці критерії, щоб переконатися, що регресійні рівняння для різних осередків приблизно однакові.

Сферичність та складна симетрія: причини використання багатовимірного підходу до повторних вимірів у дисперсійному аналізі

У планах, що містять фактори повторних вимірювань з більш ніж двома рівнями, застосування одновимірного дисперсійного аналізу потребує додаткових припущень: припущення складної симетрії та припущення сферичності. Ці припущення рідко виконуються (див. нижче). Тому в останні роки багатовимірний дисперсійний аналіз завоював популярність у таких планах (обидва підходи поєднані в модулі Дисперсійний аналіз).

Припущення про складну симетріюПрипущення складної симетрії у тому, що дисперсії (загальні внутригрупповые) і ковариации (по групам) щодо різних повторних вимірів однорідні (однакові). Це достатня умова для того, щоб одномірний критерій F для повторних вимірювань був обґрунтованим (тобто видані F-значення в середньому відповідали F-розподілу). Проте в даному випадку ця умова не є необхідною.

Припущення про сферичність.Припущення про сферичність є необхідною та достатньою умовою того, щоб F-критерій був обґрунтованим. Воно у тому, що у груп всі спостереження незалежні і однаково розподілені. Природа цих припущень, а також вплив їх порушень зазвичай не дуже добре описані в книгах дисперсійного аналізу - ця буде описана в наступних параграфах. Там буде показано, що результати одновимірного підходу можуть відрізнятися від результатів багатовимірного підходу, і буде пояснено, що це означає.

Необхідність незалежності гіпотез.Загальний спосіб аналізу даних у дисперсійному аналізі – це припасування моделі. Якщо щодо моделі, що відповідає даним, є деякі апріорнігіпотези, то дисперсія розбивається для перевірки цих гіпотез (критерії основних ефектів, взаємодій). З погляду обчислень, цей підхід генерує кілька контрастів (множина порівнянь середніх у плані). Однак якщо контрасти не незалежні один від одного, розбиття дисперсій стає беззмістовним. Наприклад, якщо два контрасти Aі Bтотожні і виділяється відповідна їм частина з дисперсії, то та сама частина виділяється двічі. Наприклад, безглуздо і безглуздо виділяти дві гіпотези: "середнє в осередку 1 вище середнього в осередку 2" і "середнє в осередку 1 вище середнього в осередку 2". Отже, гіпотези мають бути незалежні або ортогональні.

Незалежні гіпотези при повторних вимірах.Загальний алгоритм, реалізований у модулі Дисперсійний аналіз, намагатиметься для кожного ефекту генерувати незалежні (ортогональні) контрасти. Для фактора повторних вимірювань ці контрасти задають безліч гіпотез щодо різницьміж рівнями аналізованого фактора. Однак якщо ці різниці корелюються всередині груп, то результуючі контрасти не є більш незалежними. Наприклад, у навчанні, де учні вимірюються тричі за один семестр, може статися, що зміни між 1 і 2 виміром негативно корелюють зі зміною між 2 та 3 вимірами суб'єктів. Ті, хто більшу частину матеріалу освоїв між 1 і 2 вимірами, освоюють меншу частину протягом того часу, який пройшов між 2 і 3 виміром. Насправді, для більшості випадків, де дисперсійний аналіз використовуються при повторних вимірах, можна припустити, що зміни за рівнями корелюються суб'єктами. Однак коли це трапляється, припущення про складну симетрію та припущення про сферичність не виконуються і незалежні контрасти не можуть бути обчислені.

Вплив порушень та способи їх виправлення.Коли припущення про складну симетрію або сферичність не виконуються, дисперсійний аналіз може видати помилкові результати. До того, як були розроблені багатовимірні процедури, було запропоновано кілька припущень для компенсації порушень цих припущень. (див., наприклад, роботи Greenhouse & Geisser, 1959 та Huynh & Feldt, 1970). Ці методи досі широко використовуються (тому вони представлені в модулі Дисперсійний аналіз).

Підхід багатовимірного дисперсійного аналізу до повторних вимірів.Загалом проблеми складної симетрії та сферичності відносяться до того факту, що безліч контрастів, включених у дослідження ефектів факторів повторних вимірів (з числом рівнів більшим, ніж 2) не незалежні один від одного. Однак їм не обов'язково бути незалежними, якщо використовується багатовимірнийкритерій одночасної перевірки статистичного значимості двох чи більше контрастів чинника повторних вимірів. Це є причиною того, що методи багатовимірного дисперсійного аналізу стали частіше використовуватися для перевірки значущості факторів одновимірних повторних вимірів з більш ніж 2 рівнями. Цей підхід широко поширений, тому що він, у загальному випадку, не вимагає припущення про складну симетрію та припущення про сферичність.

Випадки, в яких підхід багатовимірного дисперсійного аналізу не може бути використаний.Існують приклади (плани), коли підхід багатовимірного дисперсійного аналізу може бути застосований. Зазвичай це випадки, коли є невелика кількість суб'єктів у плані та багато рівнів у факторі повторних вимірів. Тоді для проведення багатовимірного аналізу може бути замало спостережень. Наприклад, якщо є 12 суб'єктів, p = 4 фактора повторних вимірювань, і кожен фактор має k = 3 рівнів. Тоді взаємодія 4-х факторів "витрачатиме" (k-1) P = 2 4 = 16 степенів свободи. Проте є лише 12 суб'єктів, отже, у цьому прикладі багатовимірний тест може бути проведено. Модуль Дисперсійний аналізсамостійно виявить ці спостереження та обчислить лише одномірні критерії.

Відмінності в одновимірних та багатовимірних результатах.Якщо дослідження включає велику кількість повторних вимірювань, можуть виникнути випадки, коли одновимірний підхід дисперсійного аналізу до повторних вимірювань дає результати, які сильно відрізняються від тих, які були отримані при багатовимірному підході. Це означає, що різниці між рівнями відповідних повторних вимірів корелюються суб'єктами. Іноді цей факт представляє певний самостійний інтерес.

Багатомірний дисперсійний аналіз та структурне моделювання рівнянь

В останні роки моделювання структурних рівнянь стало популярним як альтернатива багатовимірному аналізу дисперсії (див. наприклад, Bagozzi and Yi, 1989; Bagozzi, Yi, and Singh, 1991; Cole, Maxwell, Arvey, and Salas, 1993). Цей підхід дозволяє перевіряти гіпотези не тільки про середні в різних групах, але так само і про кореляційні матриці залежних змінних. Наприклад, можна послабити припущення про однорідність дисперсії та підступів і явно включити в модель для кожної групи дисперсії та підступності помилки. Модуль STATISTICAМоделювання структурними рівняннями (SEPATH) (див. Том III) дозволяє проводити такий аналіз.

Дисперсійний аналіз

Курсова робота з дисципліни: «Системний аналіз»

Виконавець студент гр. 99 ІСЕ-2 Жбанов В.В.

Оренбурзький державний університет

Факультет інформаційних технологій

Кафедра прикладної інформатики

м. Оренбург-2003

Вступ

Мета роботи: ознайомитись із таким статистичним методом, як дисперсійний аналіз.

Дисперсійний аналіз (від латинського Dispersio – розсіювання) – статистичний метод, що дозволяє аналізувати вплив різних факторів на зміну, що досліджується. Метод був розроблений біологом Р. Фішером у 1925 році та застосовувався спочатку для оцінки експериментів у рослинництві. Надалі з'ясувалась загальнонаукова значущість дисперсійного аналізу для експериментів у психології, педагогіці, медицині та ін.

Метою дисперсійного аналізу є перевірка значущості різницю між середніми з допомогою порівняння дисперсій. Дисперсію вимірюваної ознаки розкладають на незалежні доданки, кожне з яких характеризує вплив того чи іншого фактора або їхньої взаємодії. Подальше порівняння таких доданків дозволяє оцінити значущість кожного фактора, що вивчається, а також їх комбінації /1/.

При істинності нульової гіпотези (про рівність середніх у кількох групах спостережень, вибраних з генеральної сукупності), оцінка дисперсії, пов'язаної із внутрішньогруповою мінливістю, має бути близька до оцінки міжгрупової дисперсії.

Під час проведення дослідження ринку нерідко постає питання про сумісність результатів. Наприклад, проводячи опитування щодо споживання будь-якого товару в різних регіонах країни, необхідно зробити висновки, наскільки дані опитування відрізняються або не відрізняються один від одного. Порівнювати окремі показники не має сенсу і тому процедура порівняння та подальшої оцінки проводиться за деякими усередненими значеннями та відхиленнями від цієї усередненої оцінки. Вивчається варіація ознаки. За міру варіації може бути прийнята дисперсія. Дисперсія σ 2 – міра варіації, яка визначається як середня з відхилень ознаки, зведених у квадрат.

Насправді часто виникають завдання більш загального характеру – завдання перевірки суттєвості відмінностей середніх вибіркових кількох сукупностей. Наприклад, потрібно оцінити вплив різної сировини на якість виробленої продукції, вирішити завдання вплив кількості добрив на врожайність с/г продукції.

Іноді дисперсійний аналіз застосовується, щоб встановити однорідність кількох сукупностей (дисперсії цих сукупностей однакові за припущенням; якщо дисперсійний аналіз покаже, як і математичні очікування однакові, то цьому сенсі сукупності однорідні). Однорідні ж сукупності можна об'єднати в одну і тим самим отримати про неї повнішу інформацію, отже, і більш надійні висновки /2/.

1 Дисперсійний аналіз

1.1 Основні поняття дисперсійного аналізу

У процесі спостереження за об'єктом, що досліджується, якісні фактори довільно або заданим чином змінюються. Конкретна реалізація фактора (наприклад, певний температурний режим, обране обладнання чи матеріал) називається рівнем чинника чи способом обробки. Модель дисперсійного аналізу з фіксованими рівнями факторів називають моделлю I, модель із випадковими факторами – моделлю II. Завдяки варіюванню фактора можна дослідити його вплив на величину відгуку. Нині загальна теорія дисперсійного аналізу розроблено моделей I.

Залежно від кількості факторів, що визначають варіацію результативної ознаки, дисперсійний аналіз поділяють на однофакторний та багатофакторний.

Основними схемами організації вихідних даних із двома та більше факторами є:

Перехресна класифікація, притаманна моделей I, у яких кожен рівень одного чинника поєднується при плануванні експерименту з кожною градацією іншого фактора;

Ієрархічна (гніздова) класифікація, характерна моделі II, у якій кожному випадковому, навмання обраному значенню одного чинника відповідає своє підмножина значень другого чинника.

Якщо одночасно досліджується залежність відгуку якісних і кількісних чинників, тобто. факторів змішаної природи, то використовується підступний аналіз /3/.

Таким чином, дані моделі відрізняються між собою способом вибору рівнів фактора, що, очевидно, насамперед впливає на можливість узагальнення отриманих експериментальних результатів. Для дисперсійного аналізу однофакторних експериментів відмінність цих двох моделей не така істотна, проте в багатофакторному дисперсійному аналізі воно може виявитися дуже важливим.

При проведенні дисперсійного аналізу повинні виконуватися такі статистичні припущення: незалежно від рівня фактора величини відгуку мають нормальний (Гаусовський) закон розподілу та однакову дисперсію. Така рівність дисперсій називається гомогенністю. Таким чином, зміна способу обробки позначається лише на положенні випадкової величини відгуку, що характеризується середнім або медіаною значенням. Тому всі спостереження відгуку належать сімейному зрушенню нормальних розподілів.

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

При невідомому законі розподілу величин відгуку використовують непараметричні (найчастіше рангові) методи аналізу.

В основі дисперсійного аналізу лежить поділ дисперсії на частини чи компоненти. Варіацію, обумовлену впливом фактора, покладеного в основу угруповання, характеризує міжгрупова дисперсія 2 . Вона є мірою варіації приватних середніх за групами

де k – число груп;

n j - число одиниць у j-ій групі;

Варіацію, обумовлену впливом інших факторів, характеризує у кожній групі внутрішньогрупова дисперсія j 2 .

Між загальною дисперсією σ 0 2 внутрішньогруповою дисперсією σ 2 і міжгруповою дисперсією

Внутрішньогрупова дисперсія пояснює вплив неврахованих при групуванні факторів, а міжгрупова дисперсія пояснює вплив факторів угруповання на середнє значення групи /2/.

1.2 Однофакторний дисперсійний аналіз

Однофакторна дисперсійна модель має вигляд:

x ij = μ + F j + ε ij , (1)

де х ij - значення досліджуваної зміною, отриманої на i-му рівні фактора (i = 1,2, ..., т) c j-м порядковим номером (j = 1,2, ..., n);

Fi – ефект, обумовлений впливом i-го рівня фактора;

ε ij – випадкова компонента, чи обурення, викликане впливом неконтрольованих чинників, тобто. варіацією зміною всередині окремого рівня.

Основні передумови дисперсійного аналізу:

Математичне очікування обурення ε ij дорівнює нулю будь-яких i, тобто.

M(ε ij) = 0; (2)

Обурення ε ij взаємно незалежні;

Дисперсія змінної x ij (або обурення ε ij) постійна для

будь-яких i, j, тобто.

D(ε ij) = σ 2; (3)

Змінна x ij (або обурення ε ij) має нормальний закон

розподілу N(0; 2).

Вплив рівнів фактора може бути фіксованим або систематичним (модель I), так і випадковим (модель II).

Нехай, наприклад, необхідно з'ясувати, чи є суттєві різницю між партіями виробів за деяким показником якості, тобто. перевірити вплив на якість одного фактора – партії виробів. Якщо включити у дослідження всі партії сировини, вплив рівня такого чинника систематичне (модель I), а отримані висновки застосовні лише до окремих партій, які залучалися щодо дослідження. Якщо ж включити лише відібрану випадково частину партій, вплив фактора випадкове (модель II). У багатофакторних комплексах можлива змішана модель III, у якій одні чинники мають випадкові рівні, інші – фіксовані.

Численні порівняння.

Досить часто виникає завдання порівняння двох груп (вибірок), а кількох – так звана завдання множинних порівнянь . Наприклад, різних вікових, професійних, соціальних верств населення, чи впливу різних доз препарату, методів діагностики тощо.

При дослідженні нового лікарського препарату нас може цікавити не лише питання достовірності впливу препарату на величину деякого параметра Х (наприклад, систолічного тиску або рівня холестерину), але й значущий вплив різних доз препарату. Тобто нам треба вивчити вплив фактора Дози(F 1), а вибірки (групи) складені із значень параметра Х для різних доз називаються – рівні фактора дози . Також нас може цікавити фактор Часу (F 2) , коли застосовують препарат (ранок, день вечір або дрібніша градація – це рівні фактора Часу ) та вплив фактора Вік пацієнта (F 3) і т.д.

У цьому випадку рекомендується на початковому етапі провести факторний дисперсійнийаналіз, який дозволяє відповісти на запитання, чи важлива хоча б одна відмінність між групами, що порівнюються.Нульова гіпотеза про рівність всіх середніх (медіан) перевіряється за параметричним критерієм Фішера або непараметричний аналог – критерієм Фаркела-Уоліса.

Здавалося б навіщо потрібний дисперсійний аналіз, якщо існує такий прекрасний і зрозумілий статистичний критерій, як т-КритерійСтьюдента?

Головне обмеження т-критеріяперед дисперсійним аналізомполягає в тому, що перший призначений для парних порівнянь, тобто ситуації, коли у нас є тільки дві вибіркиі він потребує поправки намножинні порівняння, якщо у нас більше двох груп (вибірок). По-друге, уявимо, якщо у нас 6 груп і ми шукаємо статистично значущі різницю між ними, скільки попарних порівняньу такому разі потрібно зробити?

Проблема численних порівнянь.

Якщо на одному і тому жНабір даних виконується багаторазова перевірка гіпотез, то при перевірці кожної статистичної гіпотези закладається можливість помилки першого роду (тобто відхилення вірної нульової гіпотези). Чим більше гіпотез ми перевіряємо на тих самих даних, тим більша ймовірність припуститися як мінімум однієї такої помилки. Це явище називають ефектом множинних порівнянь(multiple testing) ...

Нехай маємо три вибірки А, У і З. Порівняння за критерієм Стьюдента, Рівень значимості – a=0,05.

При порівнянні груп A і В ризик помилитися з ймовірністю 5%. хоча б в одномуіз цих трьох порівнянь складе:

Це набагато вище за 0,05.

Очевидно, що подальше збільшення числа гіпотез, що перевіряються, неминуче супроводжуватиметься і зростанням у кожному окремому тесті помилки 2 роду (зниження потужності критерію).

У такому разі необхідно використовувати критерій, який призначений для ситуацій, коли порівнюється велика кількість груп і який нам дасть єдину відповідь на всі групи, що вивчаються. дисперсійний аналіз.

Дисперсійний аналіз.

У зарубіжній літературі дисперсійний аналіз часто позначається як ANOVA , Що перекладається як аналіз варіативності (Analysis of Variance).

Дисперсійний аналіз було розроблено англійським математиком – статистиком Р. Фішером у 1918 р. для обробки результатів агрономічних дослідів щодо виявлення умов отримання максимального врожаю різних сортів сільськогосподарських культур.

Дисперсійний аналіз – це статистичний метод,призначений для оцінки впливу різних факторів на результат експерименту, а також подальшого планування аналогічних експериментів. Тобто в результаті проведення експерименту ми отримуємо кілька вибірок для кожного фактора (рівні фактора), які ми повинні порівнювати.

Умови застосування дисперсійного аналізу ANOVA

Перед тим як приступити до застосування дисперсійного аналізу,який призначений для мінімізації ризику неправильної оцінки помилки 1 роду у разі множинних порівнянь необхідно переконатися у дотриманні низки умов:

1. Безперервнийтип даних, дискретні дані менш бажані.

2. Незалежніміж собою вибірки.

3. Нормальний розподілознаки у статистичних сукупностях, у тому числі вилучені вибірки. ( Дисперсійний аналіз- Це параметричний критерій)

4. Рівність дисперсійознаки, що вивчається в статистичних сукупностях, з яких вилучені вибірки.

5. Незалежні спостереженняу кожній із вибірок.

Якщо дані не підкоряються нормальному розподілу, то при аналізі можна використовувати два способи: застосування різних арифметичних перетворень до досягнення нормальності розподілу і далі вже застосовувати дисперсійний аналіз, або використовувати критерій Краскела-Уолліса (Kruskal-Wallis H-test), Іноді його також називають непараметричним дисперсійним аналізом.

Метод дисперсійного аналізу базується на припущенні про те, що якщо на об'єкт (групу випробуваних) впливає кілька незалежних факторіві їх вплив складається, то загальну дисперсію значень ознаки, що характеризує об'єкт (групу випробуваних) можна розкласти на суму дисперсій , що виникають в результаті впливу кожного окремого фактор А , а також обумовлених випадковими впливами (залишкова дисперсія). Порівняння дисперсій зумовлених впливом різних факторів із випадковою (залишковою) дисперсією дозволяє оцінити значущість вкладу кожного з факторів, тобто оцінити достовірністьцих впливів.

Іноді дисперсійний аналіз застосовується, щоб встановити однорідність кількох сукупностей. А саме: якщо дисперсії цих сукупностей однакові (за припущенням) і якщо математичні очікування однакові, то можна вважати, що сукупності однорідні. Однорідні сукупності можна об'єднати в одну і тим самим спростити дослідження та отримати повнішу інформацію, більш надійні висновки.

За допомогою дисперсійного аналізу дослідники з дитячої лікарні у місті Пітсбурзі (США) отримали перші докази того, що при ГРВІ підвищується рівень гістаміну, що проявляється у сечі на 2-5 добу вірусної інфекції, спричиненої вірусом А. (Гістамін – основний ендогенний (внутрішній) медіатор (фактор) алергічних реакцій).

за числу факторіввплив яких досліджується, розрізняють однофакторний і багатофакторний дисперсійний аналіз.

Як було вже зазначено, дисперсійний метод тісно пов'язаний зі статистичними угрупованнями і передбачає, що сукупність, що вивчається, підрозділена на групи за факторними ознаками, вплив яких повинен бути вивчений.

На основі дисперсійного аналізу проводиться:

1. оцінка достовірності відмінностей у групових середніх за однією факторною ознакою або декількома;

2. оцінка достовірності взаємодій факторів;

3. оцінка окремих відмінностей між парами середніх.

В основі застосування дисперсійного аналізу лежить закон розкладання дисперсій (варіацій) ознаки на складові.

Загальна варіація D про результативну ознаку при угрупуванні може бути розкладена на наступні складові:

1. на міжгрупову D м пов'язану з групувальною ознакою;

2. на залишкову(внутрішньогрупову) D B , не пов'язану з групувальною ознакою.

Співвідношення між цими показниками виражається так:

D про = D м + D ст. (1.30)

Розглянемо застосування дисперсійного аналізу з прикладу.

Припустимо, потрібно довести, чи терміни посіву впливають на врожайність пшениці. Вихідні дослідні дані дисперсійного аналізу представлені в табл. 8.

Таблиця 8

У цьому прикладі N = 32, K = 4, l = 8.

Визначимо загальну сумарну варіацію врожайності, яка є сумою квадратів відхилень індивідуальних значень ознаки від загальної середньої:

де N - Число одиниць сукупності; Y i - Індивідуальні значення врожайності; Y o - загальна середня врожайність по всій сукупності.

Для визначення міжгрупової сумарної варіації, що визначає варіацію результативної ознаки за рахунок фактора, що вивчається, необхідно знати середні значення результативної ознаки по кожній групі. Ця сумарна варіація дорівнює сумі квадратів відхилень групових середніх величин від загальної середньої величини ознаки, виваженої на число одиниць сукупності в кожній із груп:

Внутрішньогрупова сумарна варіація дорівнює сумі квадратів відхилень індивідуальних значень ознаки від групових середніх за кожною групою, що підсумовується за всіма групами сукупності.

Вплив фактора на результативну ознаку проявляється у співвідношенні між D м і D в: чим сильніший вплив фактора на величину досліджуваної ознаки, тим більше D м і менше D в.

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

Обсяг варіації вже встановлено, тепер необхідно визначити кількість ступенів волі варіації. Число ступенів свободи - Це число незалежних відхилень індивідуальних значень ознаки від його середнього значення. Загальна кількість ступенів свободи, що відповідає загальній сумі квадратів відхилень у дисперсійному аналізі, розкладається за складовими варіації. Так, загальній сумі квадратів відхилень D о відповідає число ступенів свободи варіації, що дорівнює N – 1 = 31. Груповій варіації D м відповідає число ступенів свободи варіації, що дорівнює K – 1 = 3. K=28.

Тепер, знаючи суми квадратів відхилень та кількість ступенів свободи, можна визначити дисперсії для кожної складової. Позначимо ці дисперсії: d м – групові та d у – внутрішньогрупові.

Після обчислення цих дисперсій приступимо до встановлення значущості впливу фактора на результативну ознаку. Для цього знаходимо відношення: d M / d B = F ф,

Величина F ф, звана критерієм Фішера , Порівнюється з табличним, F табл. Як було зазначено, якщо F ф > F табл, вплив чинника на результативний ознака доведено. Якщо F ф< F табл то можно утверждать, что различие между дисперсиями находится в пределах возможных случайных колебаний и, следовательно, не доказывает с достаточной вероятностью влияние изучаемого фактора.

Теоретична величина пов'язана з ймовірністю, і в таблиці її значення наводиться за певного рівня ймовірності судження. У додатку є таблиця, що дозволяє встановити можливу величину F при ймовірності судження, що найчастіше використовується: рівень ймовірності «нульової гіпотези» – 0,05. Замість ймовірностей «нульової гіпотези» таблицю можна назвати таблицею для ймовірності 0,95 суттєвості впливу фактора. Підвищення рівня ймовірності вимагає порівняння вищого значення F табл.

Величина F табл залежить також від числа ступенів свободи двох дисперсій, що порівнюються. Якщо число ступенів свободи прагне нескінченності, то F табл прагне одиниці.

Таблиця значень F табл побудована наступним чином: у стовпцях таблиці вказано ступеня свободи варіації для більшої дисперсії, а в рядках – ступеня свободи для меншої (внутрішньогрупової) дисперсії. Величина F знаходиться на перетині стовпця та рядки відповідних ступенів свободи варіації.

Так, у прикладі F ф = 21,3/3,8 = 5,6. Табличне значення F табл для ймовірності 0,95 і ступенів свободи, відповідно рівних 3 і 28, F табл = 2,95.

Значення F ф отримане досвіді, перевищує теоретичне значення навіть ймовірності 0,99. Отже, досвід із ймовірністю понад 0,99 доводить вплив досліджуваного чинника на врожайність, т. е. досвід вважатимуться надійним, доведеним, отже, терміни посіву істотно впливають на врожайність пшениці. Оптимальним терміном посіву слід вважати період з 10 по 15 травня, тому що саме при цьому строку посіву отримано найкращі результати врожайності.

Нами розглянуто методику дисперсійного аналізу при групуванні за однією ознакою та випадковим розподілом повторностей усередині групи. Однак часто буває так, що дослідна ділянка має якісь відмінності в родючості ґрунту і т. д. Тому може виникнути така ситуація, що більша кількість ділянок одного з варіантів потрапить на кращу частину, і її показники будуть завищені, а іншого варіанта – на найгіршу частину, і результати в цьому випадку, природно, будуть гіршими, тобто занижені.

Щоб виключити варіювання, яке викликається причинами, що не належать до досвіду, треба з внутрішньогрупової (залишкової) дисперсії вичленувати дисперсію, розраховану за повторностями (блоками).

Загальна сума квадратів відхилень поділяється в цьому випадку вже на 3 складові:

D про = D м + D повт + D зуп. (1.33)

Для нашого прикладу сума квадратів відхилень, викликана повторностями, дорівнюватиме:

Отже, що випадкова сума квадратів відхилень дорівнюватиме:

D ост = D в - D повт; D зост = 106 - 44 = 62.

Для залишкової дисперсії число ступенів свободи дорівнюватиме 28 – 7 = 21. Результати дисперсійного аналізу представлені у табл. 9.

Таблиця 9

Оскільки фактичні значення F-критерію для ймовірності 0,95 перевищують табличні, вплив строків посіву і повторностей на врожайність пшениці слід вважати істотним. Розглянутий спосіб побудови досвіду, коли ділянка попередньо ділиться на блоки з відносно вирівняними умовами, а варіанти, що перевіряються, розподіляються всередині блоку у випадковому порядку, називається способом рендомізованих блоків.

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

Двофакторний дисперсійний аналіз відрізняється від двох однофакторних тим, що він може відповісти на такі питання:

1. 1який вплив обох факторів разом?

2. яка роль поєднання цих факторів?

Розглянемо дисперсійний аналіз досвіду, у якому слід виявити вплив як термінів посіву, а й сортів на врожайність пшениці (табл. 10).

Таблиця 10. Дані досвіду щодо впливу строків посіву та сортів на врожайність пшениці

- Це сума квадратів відхилень індивідуальних значень від загальної середньої.

Варіація щодо спільного впливу строків посіву та сорту

– це сума квадратів відхилень середніх по підгрупах від загальної середньої, зважених число повторностей, т. е. на 4.

Обчислення варіації щодо впливу лише строків посіву:

Залишкова варіація визначається як різниця між загальною варіацією та варіацією щодо спільного впливу досліджуваних факторів:

D ост = D про - D пс = 170 - 96 = 74.

Усі розрахунки можна оформити як таблиці (табл. 11).

Таблиця 11. Результати дисперсійного аналізу

Результати дисперсійного аналізу показують, що вплив факторів, що вивчаються, тобто термінів посіву і сорту, на врожайність пшениці суттєво, оскільки F-критерії фактичні по кожному з факторів значно перевищують табличні, знайдені для відповідних ступенів свободи, і при цьому з досить високою ймовірністю (р = 0,99). Вплив поєднання чинників у разі відсутня, оскільки чинники незалежні друг від друга.

Аналіз впливу трьох факторів на результат ведеться за таким же принципом, що і для двох факторів, тільки в цьому випадку три дисперсії за факторами і чотири дисперсії по поєднанню факторів. Зі збільшенням числа факторів різко збільшується обсяг розрахункових робіт і, крім того, важко оформляти вихідну інформацію в комбінаційну таблицю. Тому навряд чи доцільно вивчати вплив багатьох факторів на результат із використанням дисперсійного аналізу; краще взяти менше їх число, але вибрати найбільш істотні чинники з погляду економічного аналізу.

Нерідко досліднику доводиться мати справу з так званими непропорційними дисперсійними комплексами, тобто такими, в яких не дотримується пропорційності чисельностей варіантів.

У таких комплексах варіація сумарної дії факторів не дорівнює сумі варіації за факторами та варіації поєднання факторів. Вона відрізняється на величину, яка залежить від ступеня зв'язків між окремими факторами, що виникають унаслідок порушення пропорційності.

І тут виникають труднощі щодо ступеня впливу кожного чинника, оскільки сума приватних впливів не дорівнює сумарному впливу.

p align="justify"> Одним із способів приведення непропорційного комплексу до єдиної структури є спосіб його заміни пропорційним комплексом, в якому частоти усереднені по групах. Коли така заміна зроблена, завдання вирішується за принципами пропорційних комплексів.

У практичній діяльності лікарів при проведенні медико-біологічних, соціологічних та експериментальних досліджень виникає потреба встановити вплив факторів на результати вивчення стану здоров'я населення, при оцінці професійної діяльності, ефективності нововведень.

Існує ряд статистичних методів, що дозволяють визначити силу, напрям, закономірності впливу факторів на результат у генеральній або вибірковій сукупності (розрахунок критерію I, кореляційний аналіз, регресія, Χ 2 - (критерій згоди Пірсона та ін.). Дисперсійний аналіз був розроблений та запропонований англійською вченим, математиком та генетиком Рональдом Фішером у 20-х роках XX століття.

Дисперсійний аналіз найчастіше використовують у науково-практичних дослідженнях громадського здоров'я та охорони здоров'я для вивчення впливу одного чи кількох факторів на результативну ознаку. Він заснований на принципі "відображення різноманітностей значень факторного(их) на розмаїтті значень результативної ознаки" та встановлює силу впливу фактора(ів) у вибіркових сукупностях.

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

Дисперсійний аналіз- це статистичний метод оцінки зв'язку між факторними та результативними ознаками у різних групах, відібраний випадковим чином, заснований на визначенні відмінностей (різноманіття) значень ознак. У основі дисперсійного аналізу лежить аналіз відхилень всіх одиниць досліджуваної сукупності середнього арифметичного. Як міра відхилень береться дисперсія (В)-середній квадрат відхилень. Відхилення, що викликаються впливом факторної ознаки (фактора), порівнюються з величиною відхилень, що викликаються випадковими обставинами. Якщо відхилення, викликані факторним ознакою, істотніші, ніж випадкові відхилення, вважається, що чинник істотно впливає на результативний ознака.

Для того щоб обчислити дисперсію значення відхилень кожної варіанти (кожного зареєстрованого числового значення ознаки) від середнього арифметичного зводять у квадрат. Тим самим позбавляються негативних знаків. Потім ці відхилення (різниці) підсумовують і поділяють число спостережень, тобто. усереднюють відхилення. Таким чином, набувають значення дисперсій.

Важливим методичним значенням застосування дисперсійного аналізу є правильне формування вибірки. Залежно від поставленої мети та завдань вибіркові групи можуть формуватися випадковим чином незалежно одна від одної (контрольна та експериментальна групи для вивчення деякого показника, наприклад, вплив високого артеріального тиску на розвиток інсульту). Такі вибірки називаються незалежними.

Нерідко результати впливу факторів досліджуються в однієї і тієї ж вибіркової групи (наприклад, в тих самих пацієнтів) до і після впливу (лікування, профілактика, реабілітаційні заходи), такі вибірки називаються залежними.

Дисперсійний аналіз, у якому перевіряється вплив одного чинника, називається однофакторним (одномірний аналіз). При вивченні впливу більш як одного фактора використовують багатофакторний дисперсійний аналіз (багатомірний аналіз).

Факторні ознаки - це ознаки, що впливають на досліджуване явище.
Результативні ознаки - це ознаки, які змінюються під впливом факторних ознак.

Для проведення дисперсійного аналізу можуть використовуватися як якісні (стаття, професія), так і кількісні ознаки (число ін'єкцій, хворих у палаті, число ліжко-днів).

Методи дисперсійного аналізу:

1. Метод Фішеру (Fisher) - критерій F (значення F див. у додатку N 1);
   Метод застосовується в однофакторному дисперсійному аналізі, коли сукупна дисперсія всіх значень, що спостерігаються, розкладається на дисперсію всередині окремих груп і дисперсію між групами.
2. Метод "загальної лінійної моделі".
   У його основі лежить кореляційний чи регресійний аналіз, застосовуваний багатофакторному аналізі.

Зазвичай у медико-біологічних дослідженнях використовуються лише однофакторні, максимум двофакторні дисперсійні комплекси. Багатофакторні комплекси можна досліджувати, послідовно аналізуючи одно-або двофакторні комплекси, що виділяються з усієї спостережуваної сукупності.

Умови застосування дисперсійного аналізу:

1. Завданням дослідження є визначення сили впливу одного (до 3) факторів на результат або визначення сили спільного впливу різних факторів (стаття та вік, фізична активність та харчування тощо).
2. Досліджувані чинники мають бути незалежні (непов'язані) між собою. Наприклад, не можна вивчати спільний вплив стажу роботи та віку, зростання та ваги дітей тощо. на захворюваність населення.
3. Підбір груп на дослідження проводиться рандомизированно (випадковий відбір). Організація дисперсійного комплексу із дотриманням принципу випадковості відбору варіантів називається рандомізацією (перев. з англ. – random), тобто. вибрані навмання.
4. Можна застосовувати як кількісні, і якісні (атрибутивні) ознаки.

При проведенні однофакторного дисперсійного аналізу рекомендується (необхідна умова застосування):

1. Нормальність розподілу аналізованих груп чи відповідність вибіркових груп генеральним сукупностям з нормальним розподілом.
2. Незалежність (не пов'язаність) розподілу спостережень у групах.
3. Наявність частоти (повторність) спостережень.

Нормальність розподілу визначається кривою Гаусса (Де Мавура), яку можна описати функцією у = f(х), так як вона належить до законів розподілу, що використовуються для наближеного опису явищ, які носять випадковий, імовірнісний характер. Предмет медико-біологічних досліджень - явища імовірнісного характеру, нормальний розподіл у таких дослідженнях трапляється досить часто.

Принцип застосування методу дисперсійного аналізу

Спочатку формулюється нульова гіпотеза, тобто передбачається, що досліджувані фактори не впливають на значення результативної ознаки та отримані відмінності випадкові.

Потім визначаємо, яка можливість отримати спостерігаються (чи сильніші) відмінності за умови справедливості нульової гіпотези.

Якщо ця ймовірність мала*, ми відкидаємо нульову гіпотезу і укладаємо, що результати дослідження статистично значущі. Це ще не означає, що доведено дію саме досліджуваних факторів (це питання, перш за все, планування дослідження), але все ж таки малоймовірно, що результат обумовлений випадковістю.
__________________________________
* Максимальну прийнятну можливість відкинути правильну нульову гіпотезу називають рівнем значущості і позначають α = 0,05.

При виконанні всіх умов застосування дисперсійного аналізу розкладання загальної дисперсії математично виглядає наступним чином:

D ощ. = D факт + D зуп. ,

D ощ. - загальна дисперсія значень, що спостерігаються (варіант), характеризується розкидом варіант від загального середнього. Вимірює варіацію ознаки у всій сукупності під впливом всіх факторів, що зумовили цю варіацію. Загальна різноманітність складається з міжгрупового та внутрішньогрупового;

D факт - факторна (міжгрупова) дисперсія, що характеризується відмінністю середніх у кожній групі та залежить від впливу досліджуваного фактора, за яким диференціюється кожна група. Наприклад, у групах різних за етіологічним фактором клінічного перебігу пневмонії середній рівень проведеного ліжко-дня неоднаковий – спостерігається міжгрупова різноманітність.

D зуп. - Залишкова (внутрішньогрупова) дисперсія, яка характеризує розсіювання варіант всередині груп. Відбиває випадкову варіацію, тобто. частина варіації, що відбувається під впливом неуточнених факторів і не залежить від ознаки - фактора, покладеного в основу угруповання. Варіація ознаки, що вивчається, залежить від сили впливу якихось неврахованих випадкових факторів, як від організованих (заданих дослідником), так і від випадкових (невідомих) факторів.

Тому загальна варіація (дисперсія) складається з варіації, викликаної організованими (заданими) чинниками, званими факторіальної варіацією і неорганізованими чинниками, тобто. залишковою варіацією (випадковою, невідомою).

Класичний дисперсійний аналіз проводиться за такими етапами:

1. Побудова дисперсійного комплексу.
2. Обчислення середніх квадратів відхилень.
3. Обчислення дисперсії.
4. Порівняння факторної та залишкової дисперсій.
5. Оцінка результатів з допомогою теоретичних значень розподілу Фішера-Снедекора (додаток N 1).

АЛГОРИТМ ПРОВЕДЕННЯ ДИСПЕРСІЙНОГО АНАЛІЗУ ЗА СПРОЩЕНИМ ВАРІАНТОМ

Алгоритм проведення дисперсійного аналізу за спрощеним способом дозволяє отримати ті ж результати, але розрахунки виконуються значно простіше:

І етап. Побудова дисперсійного комплексу

Побудова дисперсійного комплексу означає побудову таблиці, де було б чітко розмежовані чинники, результативний ознака і добір спостережень (хворих) у групу.

Однофакторний комплекс складається з кількох градацій одного фактора (А). Градації - це вибірки різних генеральних сукупностей (А1, А2, АЗ).

Двофакторний комплекс – складається з кількох градацій двох факторів у комбінації між собою. Етіологічні чинники захворюваністю на пневмонію ті ж (А1, А2, АЗ) у поєднанні з різними формами клінічного перебігу пневмонії (Н1 - гостре, Н2 - хронічне).

ІІ етап. Обчислення загальної середньої (М обш)

Обчислення суми варіант за кожною градацією факторів: Vj = V 1 + V 2 + V 3

Обчислення загальної суми варіант (Σ V заг) за всіма градаціями факторної ознаки: Σ V заг = Σ Vj 1 + Σ Vj 2 + Σ Vj 3

Обчислення середньої групової (М грн) факторної ознаки: М грн. = Vj / N,
де N - сума числа спостережень за всіма градаціями факторної I ознаки (?n за групами).

ІІІ етап. Розрахунок дисперсій:

За дотримання всіх умов застосування дисперсійного аналізу математична формула виглядає так:

D ощ. = D факт + D зуп.

D ощ. - загальна дисперсія, що характеризується розкидом варіант (спостережуваних значень) від загального середнього;
D факт. - факторна (міжгрупова) дисперсія, що характеризує розкид групових середніх від загального середнього;
D зуп. - залишкова (внутрішньогрупова) дисперсія, що характеризує розсіювання варіант усередині груп.

1. Обчислення факторіальної дисперсії (D факт.): D факт. = Σ h - H
2. Обчислення h проводиться за такою формулою: h = (Σ Vj) / N
3. Обчислення Н проводиться за такою формулою: H = (Σ V) 2 / N
4. Обчислення залишкової дисперсії: D зуп. = (Σ V) 2 - Σ h
5. Обчислення загальної дисперсії: D ощ. = (Σ V) 2 - Σ H

ІV етап. Розрахунок основного показника сили впливу фактора, що вивчаєтьсяПоказник сили впливу (η ​​2) факторної ознаки на результат визначається часткою факторіальної дисперсії (D факт.) у загальній дисперсії (D oбщ.), η 2 (ця) - показує яку частку займає вплив фактора, що вивчається, серед усіх інших факторів і визначається за формулою :

V етап. Визначення достовірності результатів дослідження методом Фішера проводять за такою формулою:

Порівняння критерію Фішера (F) із стандартним (табличним) F проводять за графами таблиці з урахуванням ступенів свободи:

v 1 = n - 1
v 2 = N - 1

По горизонталі визначають v 1 по вертикалі - v 2 на їх перетині визначають табличне значення F, де верхнє табличне значення р ≥ 0,05, а нижнє відповідає р > 0,01, і порівнюють з обчисленим критерієм F. Якщо значення обчисленого критерію F одно або більше табличного, результати достовірні і Н 0 не відкидається.

Умова задачі:

На підприємстві Н. підвищився рівень травматизму, у зв'язку з чим лікар провів дослідження окремих факторів, серед яких вивчався стаж роботи працюючих у цехах. Вибірки зроблено на підприємстві Н. з 4 цехів із близькими умовами та характером праці. Рівні травматизму розраховані на 100 працюючих за минулий рік.

При дослідженні фактора робочого стажу отримано такі дані:

На підставі даних проведеного дослідження було висунуто нульову гіпотезу (Н0) про вплив стажу роботи на рівень травматизму працівників підприємства А.А.

Завдання
Підтвердіть або спростуйте нульову гіпотезу методом однофакторного дисперсійного аналізу:

1. визначте силу впливу;
2. оцінити достовірність впливу фактора.

Етапи застосування дисперсійного аналізу
для визначення впливу фактора (стажу роботи) на результат (рівень травматизму)

Висновок.У вибірковому комплексі виявлено, що сила впливу стажу роботи на рівень травматизму становить 80% від загального числа інших факторів. Для всіх цехів заводу можна з ймовірністю 99,7% (13,3 > 8,7) стверджувати, що стаж роботи впливає рівень травматизму.

Таким чином, нульова гіпотеза (Н0) не відкидається і вплив стажу роботи на рівень травматизму в цехах заводу А вважається доведеним.

Значення F (критерій Фішера) стандартного при р ≥ 0,05 (верхнє значення) при р ≥ 0,01 (нижнє значення)

1. Власов В.В. Епідеміологія. - М: ГЕОТАР-МЕД, 2004. 464 с.
2. Архіпова ГЛ., Лаврова І.Г., Трошина І.М. Деякі сучасні методи статистичного аналізу у медицині. - М.: Метропостач, 1971. - 75 с.
3. Зайцев В.М., Ліфляндський В.Г., Марінкін В.І. Прикладна статистика. – СПб.: ТОВ "Видавництво ФОЛІАНТ", 2003. – 432 с.
4. Платонов А.Є. Статистичний аналіз у медицині та біології: завдання, термінологія, логіка, комп'ютерні методи. – М.: Видавництво РАМН, 2000. – 52 с.
5. Плохінський Н.А. Біометрія. - Видавництво Сибірського відділення АН СРСР Новосибірськ. – 1961. – 364 с.