Чи покладається ANOVA на метод моментів, а не на максимальну ймовірність?


13

Я бачу в різних місцях, що ANOVA робить свою оцінку методом моментів.

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

Так:

  1. Що кваліфікує процедури ANOVA як метод моментів?

  2. Враховуючи, що ANOVA еквівалентний OLS з категоричними прогнозами, чи не це максимальна ймовірність?

  3. Якщо ці два методи якимось чином виявляться еквівалентними у спеціальному випадку звичайної ANOVA, чи існують якісь конкретні ситуації з ANOVA, коли різниця стає важливою? Неврівноважений дизайн? Повторні заходи? Змішаний (між суб'єктами + всередині предметів) дизайн?


4
У простому налаштуванні, а під простим я маю на увазі одно- та двосторонню ANOVA, ANOVA отримується від LRT при нормальних розподілах з однаковими дисперсіями, так що максимальна ймовірність використовується. Звичайно, для нормального випадку оцінки молока та мами збігаються, тому розмежування не має великого значення. Однак у більш складних налаштуваннях, а не на виведенні LRT, ми покладаємось на регресію OLS. OLS - це лише млечні вироби за звичайного розподілу і є маминим оцінником у більш загальних умовах, якщо ми нав'язуємо ортогональність із залишками, які є.
JohnK

1
@JohnK, велике спасибі за ваш коментар, але мені важко зрозуміти деякі його частини, зокрема другу частину: що таке "складніші налаштування", в яких ANOVA покладається на OLS замість LRT (і чому)? Чому OLS не є MLE в цій "більш загальній установці" - я вважав, що нормальні помилки завжди допускаються в будь-якому сценарії ANOVA? Що має робити ортогональність із залишками? Я буду дуже вдячний, якщо ви розгорнете свій коментар у відповідь.
Амеба каже, що повернеться до Моніки

2
@amoeba. Я відредагував свою відповідь, щоб вирішити ваш перший пункт про те, що ANOVA є методом оцінювання моментів. Це справедливо лише для випадкових ефектів.
Placidia

Відповіді:


12

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

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

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

Нещодавно я написав відповідь, де виникла різниця між методами MOM та ML. Питання стосувалося оцінки моделей випадкових ефектів. На даний момент традиційний підхід ANOVA повністю розділяє компанію з максимальною оцінкою ймовірності, і оцінки наслідків вже не однакові. Якщо дизайн не врівноважений, ви також не отримуєте однакових статистичних даних F.

Ще в той день, коли статистики хотіли обчислити випадкові ефекти від розщеплення або повторних проектів вимірювань, відхилення випадкових ефектів було обчислено із середніх квадратів таблиці ANOVA. Отже, якщо у вас є ділянка з дисперсією а залишкова дисперсія - , можливо, у вас з'явиться очікуване значення середнього квадрату ("очікуваний середній квадрат", EMS) , з числом розщеплень на . Ви встановите середній квадрат, рівний його очікуванню, і вирішите для σ 2 σ 2 + n σ 2 p n ^ σ 2 bσp2σ2σ2+nσp2nσb2^. ANOVA дає метод оцінки моментів для дисперсії випадкових ефектів. Зараз ми схильні вирішувати подібні проблеми за допомогою моделей зі змішаними ефектами, і компоненти дисперсії отримують за допомогою максимальної оцінки ймовірності або REML.

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

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


2
F-тест одностороннього та двостороннього ANOVA - це тест на ймовірність, і ви можете проконсультуватися з Хоггом і Крейгом, Вступ до математичної статистики, глава 9, якщо ви мені не вірите. Більше того, F-тест, що використовується в нормальній регресії, також є LRT, численні посилання на це там.
ДжонК

1
Я перевірю. Наразі абзац я видалив.
Placidia

1
+1, дуже дякую за відповідь. Це, до речі, викликало моє запитання. Є багато речей, які я не зовсім розумію в тому, що ви написали. Я буду подорожувати за Новий рік і не матиму часу думати про це, але я повернусь до цієї теми, як тільки повернусь у січні. Тим часом я хотів би ще раз закликати @JohnK опублікувати відповідь, розширюючи питання, висловлені в коментарях. Я особливо вдячний ще за деяким технічним поясненням, можливо, з конкретним прикладом. Веселого Різдва та щасливого Нового року всім вам!
Амеба каже, що поверніть Моніку

1
З Різдвом Христовим і вам. У своєму зв’язаному дописі я просто думав про оцінку випадкової дисперсії ефекту, яка є MOM у традиційному підході. Вибачте, якщо моє зауваження звучало більш загально, ніж я задумав.
Placidia

2
Через 1,5 роки я перечитав вашу відповідь і зрозумів, що це досить зрозуміло, тому я нарешті позначив її як прийняту :-) Привіт.
амеба каже, що поверніть Моніку
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.