Односторонний анализ случайных эффектов в SAS: PROC GLM или MIXED?

Я пытаюсь провести простой односторонний дисперсионный анализ со случайными эффектами в SAS. Я хочу знать, значительно ли дисперсия населения отличается от нуля или нет.

На сайте idre UCLA указано, что используется PROC MIXED. следующим образом:

proc mixed data = in.hsb12 covtest noclprint;
   class school;
   model mathach = / solution;
   random intercept / subject = school;
run;

Это имеет смысл для меня, учитывая мой предыдущий опыт использования PROC MIXED.

Однако в тексте Мюррея Логана «Биостатистический план и анализ с использованием R» он говорит, что для однофакторного дисперсионного анализа фиксированные и случайные эффекты не различаются, и проводит (в R) «стандартный» однофакторный дисперсионный анализ, хотя он проверяет дисперсия, а не средства. Я обнаружил, что в SAS его процедура R эквивалентна использованию любого из следующих действий:

  1. PROC ANOVA
  2. PROC GLM (то же, что и ANOVA, но с GLM вместо ANOVA)
  3. PROC GLM с оператором RANDOM

Значения p для трех вышеупомянутых моделей одинаковы, но отличаются от модели PROC MIXED, используемой Калифорнийским университетом в Лос-Анджелесе. Для моих данных разница составляет p=0,2508 и p=0,3138. Хотя выводы в этом случае не меняются, меня не очень устраивает эта разница.

Может ли кто-нибудь дать совет, какой из них более подходит, а также почему существует эта разница?


sas anova random-effects
person Meg    schedule 05.10.2014    source источник
comment
Какие параметры вы используете в вызове PROC ANOVA?   -  person Kostya    schedule 16.10.2014
comment
Нет, только PROC ANOVA, за которым следуют операторы CLASS и MODEL.   -  person Meg    schedule 16.10.2014

Ответы (1)


arrow_upward
1
arrow_downward

Для вашей модели разница между PROC ANOVA и PROC MIXED связана только с числовым шумом (REML-оценка PROC MIXED). Однако значения p, упомянутые в вашем вопросе, соответствуют разным тестам. Чтобы получить значение F, используя выходные данные COVTEST в PROC MIXED, вам необходимо пересчитать MS_groups с учетом неравных размеров выборки (либо вручную, как описано на стр. 231 http://bio.classes)..ucsc.edu/bio286/MIcksBookPDFs/QK08.PDF или просто используя PROC MIXED с та же спецификация фиксированной модели, что и в PROC ANOVA). Этот документ (http://isites.harvard.edu/fs/docs/icb.topic1140782.files/S98.pdf) содержит несколько примеров использования PROC MIXED в дополнение к руководству по SAS.

person Kostya    schedule 20.10.2014
comment
Очень полезно, спасибо. Мне нужно будет найти время, чтобы прочитать эти ссылки, чтобы понять, почему PROC MIXED не может обрабатывать несбалансированные данные, а PROC ANOVA и GLM в порядке. - person Meg; 23.11.2014