дисперсионный анализ повторных измерений, исследование качества жизни Опции

[mamalita]

Добрый день! Кратко изложу суть исследования. Исследуем две группы один вид хирургического лечения - основная группа (n=48), другой вид - контрольная группа (n=55). В обеих группах исследуем качество жизни по шкале FACT-L до операции, через 3, 6, 9 и 12 месяцев после операции - 8 показателей на каждый период. Количество больных с течением времени сокращается 3 месяца 47 и 53 человека, 6 месяцев - 46 и 52, 9 - 45 и 50 человек, 12 - 44 и 46 человек. Программа Statistica 6. Распределение показателей качества жизни отличается от нормального. Поэтому анализ проводила непараметрический.С помощью критерия Манна-Уитни до операции все 8 показателей статистически значимо не отличаются, что собственно нам и было нужно. Попарно сравнивала между собой каждый показатель из 8 на каждый промежуток времени, чтобы доказать что качество жизни в основной группе выше на данный промежуток времени, чем в контрольной (или нет, как получалось). Ранговый ДА Фридмана и конкордации Кенделла использовала для каждой группы отдельно для оценки в динамике. р<0,001. А дальше как? сравниваем критерием Вилкоксона попарно до операции-через 3 месяца, через 3 месяца - через 6 месяцев и т.д.? Нужна ли поправка Бонферони? Основной вопрос такой: правильно ли я делала? Можно ли использовать дисперсионный анализ для повторных измерений? (я попробовала его сделать, там где различия для непараметрического метода были значимы, стали в ДА не значимы, с чем это связано?) Еще прочитала что если распределение отличается от нормального, дисперсии не равны можно воспользоваться модулем GLM - так ли это (результаты в нем получились в точности такие же как RM-ANOVA. Помогите как правильно сделать.

[nokh]

Полного непараметрического аналога RM-ANOVA нет, поэтому видится 2 пути, которыми вы собственно и идёте.
1). Попарные сравнения. На каждом сроке: группу - с группой (Манн - Уитни), для каждой группы: срок с последующим (Вилкоксон для разностей пар). Бонферрони просится. Если различия высоко статистически значимы - поправка не страшна. Если же нет - вообще всё может свести на нет. Но в любом случае в такой схеме теряется информация о различиях динамик. В дисперсионном анализе всех данных она бы вышла в виде взаимодействия факторов Группа х Срок. Нужно ли пытаться её достать - вам виднее, но скорее всего - ДА.
2). Дисперсионный анализ. Пакет Statistica не позволит обработать все ваши данные, т.к. это дисперсионный комплекс с пропусками: раз люди терялись в ходе эксперимента - в схеме с повторными измерениями появятся пропущенные ячейки. В принципе такой комплекс тоже можно обработать грамотно, но для этого нужен другой софт и знания, которыми я не обладаю (да и единства подхода в такой ситуации нет: различные схемы модифицированного расчёта сумм квадратов отклонений, импутация пропущенных значений). Если исходить из того что есть, нужно обрезать материал - я бы так и сделал. Т.е. оставлять в анализе только тех людей, данные по которым есть для всех сроков. Ваша ситуация не "аховая": не так уж много потерялось. Полагаю, что выигрыш в мощности от совместного анализа максимально возможного количества данных с лихвой окупит такую потерю.
GLM - это обобщение дисперсионного и регрессионного анализов, т.е. параметрическая техника, а значит должна быть нормальность (но не исходных данных, а остатков). Если окажется, что они распределены ненормально - не страшно: можно преобразовать исходные данные - этой теме не форуме уделялось много внимания.
Итак, в рамках такого подхода факторы будут: Группа (фиксированный), Срок (фиксированный) и Индивид внутри Группы (случайный, обозначается так: Индивид (Группа)). Т.е. имеем смешанную модель трёхфакторного дисперсионного анализа с пересекающимися и иерархическими факторами. В ходе дисперсионного анализа в модуле GLM оцениваться будут следующие эффекты:
Группа
Срок
Индивид (Группа)
Взаимодействие Группа х Срок
Взаимодействие Срок х Индивид(Группа)
В Statistica важно не забыть отметить на вкладке случайный фактор, а также нужно верно задать иерархический фактор Индивид (Группа) и взаимодействие Срок х Индивид(Группа) - описывать долго, проще "методом тыка". Попробуйте. Кое-что полезного можно найти здесь:
http://forum.disser.ru/index.php?showtopic...D1%8C&st=15 (сообщение 16).