Блог

меня зовут НОЗОМИ. Я биостатист в Японии.
Я впервые пытаюсь задавать вопросы на английском языке.
Я приношу извинения, если мой английский трудно понять.

Я использовал lmer() в пакете lme4 для анализа модели смешанных эффектов.
Стандартные ошибки коэффициентов каждого уровня случайного эффекта отличались от результатов, рассчитанных SAS.
В качестве примечания, точечные оценки были близки между R и SAS.

Я загружу сценарии и выходные данные R и SAS, поэтому, пожалуйста, скажите мне, как описать сценарий в R, чтобы он соответствовал стандартной ошибке в SAS.

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

[Сценарий R]

 model <- lmer(y ~ f   (1|r),data=data)
coef_random <- tidy(model, effects="ran_vals", conf.int=T)
 

[Вывод R]

[Сценарий SAS]

 proc mixed data=data ;
 class r ;
 model y = f / cl ;
 random r / cl ;
run ;
 

[Вывод SAS]