Множественные сравнения линейных регрессий |
Здравствуйте, гость ( Вход | Регистрация )
Множественные сравнения линейных регрессий |
27.02.2011 - 16:26
Сообщение
#1
|
|
Группа: Пользователи Сообщений: 244 Регистрация: 28.08.2009 Пользователь №: 6286 |
Например, надо сравнить три линейных регрессии. Выборки неравных объемов. Омнибусный тест показывает наличие различий. Как в таких случаях поступают? Существуют ли техники применения апостериорных критериев для сравнения линейных регрессий?
|
|
28.02.2011 - 15:53
Сообщение
#2
|
|
Группа: Пользователи Сообщений: 244 Регистрация: 28.08.2009 Пользователь №: 6286 |
Может быть через General Linear Models это можно как-то решить?
|
|
1.03.2011 - 05:25
Сообщение
#3
|
|
Группа: Пользователи Сообщений: 1202 Регистрация: 13.01.2008 Из: Челябинск Пользователь №: 4704 |
А почему не хотите простой поправкой Бонферрони обойтись? Выборок здесь очень мало (3), поэтому она не будет слишком консервативной. Вообще мне эта поправка не нравится, но в данном случае использовал бы её.
|
|
1.03.2011 - 09:46
Сообщение
#4
|
|
Группа: Пользователи Сообщений: 244 Регистрация: 28.08.2009 Пользователь №: 6286 |
А почему не хотите простой поправкой Бонферрони обойтись? Выборок здесь очень мало (3), поэтому она не будет слишком консервативной. Вообще мне эта поправка не нравится, но в данном случае использовал бы её. Я думал об этом. В смысле, как по Гланцу: сравнивать попарно, но считать критическим уровень значимости не 0,05, а 0,05/3=0,017. Мне бы именно как-то так и надо ("правильно" уменьшить уровень значимости). Поскольку еще есть такая же задача, только с четырьмя регрессиями и гетероскедастичными остатками, и я там думаю с помощью dummy и весовой переменной сравнивать. Тогда введение поправки при попарных сравнениях было бы выходом. Просто, поскольку Гланц пишет, что все апостериорные критерии, как и Бонферрони, являются модификациями критерия Стьюдента, я подумал, что может можно примерно также уменьшать, но только не по Бонферони, а по Fisher LSD или Unequal N HSD . Nokh, а если опираться на Бонферрони, не подскажете какой-нибудь авторитетный источник (посильнее Гланца, можно на англ.), где бы было сказано про максимально возможное число сравнений (Гланц пишет про 8, на форуме читал про 5, а где-то попадалась статья "Кошмар Бонферрони", там вообще запугали); чтобы сослаться можно было. Сообщение отредактировал Pinus - 1.03.2011 - 09:48 |
|
1.03.2011 - 22:09
Сообщение
#5
|
|
Группа: Пользователи Сообщений: 377 Регистрация: 18.08.2008 Из: Москва Златоглавая Пользователь №: 5224 |
максимально возможное число сравнений (Гланц пишет про 8, на форуме читал про 5, а где-то попадалась статья "Кошмар Бонферрони", там вообще запугали) Попарные сравнения с поправкой Бонферони сильно снижают мощность критерия с увеличением количества сравнений и не дают найти различия там, где они на самом деле есть. Если у вас нет обобщенного критерия типа дисперсионного анализа, сравнивающего все группы одновременно, то вы можете использовать Бонферони, расплачиваясь за кажущуюся простоту потерей мощности. Для обоснования применения поправки Бонферони рекомендую задать минимальную мощность критерия, на которую вы еще согласны, и рассчитать для этой мощности максимальное количество попарных сравнений. Для расчетов можно использовать метод Монте-Карло (случайных испытаний), если не получится вывести точную формулу.Просто включи мозги => http://doctorstat.narod.ru
|
|
18.03.2011 - 16:13
Сообщение
#6
|
|
Группа: Пользователи Сообщений: 244 Регистрация: 28.08.2009 Пользователь №: 6286 |
Пишут
http://www.statsoft.ru/home/portal/applica...comparisons.htm что для смягчения поправки Бонферрони нужно использовать обобщенную оценку внутригрупповой дисперсии. Число степеней свободы при таком подходе для критерия Стьюдента равно f = m*(n - 1), где n - объем групп, а для групп разного объема число степеней свободы будет равно суммарной численности всех групп N минус количество групп m (что в случае m>2 превышает обычное число степеней свободы для критерия Стьюдента, равное суммарной численности двух непосредственно сравниваемых групп минус 2). Правильно ли я понимаю, что кроме того, что при введении поправки мы делаем более высоким уровень значимости, следует еще и p-уровень определять с другим числом степеней свободы (f = m*(n - 1))? |
|