#r #sas #lme4
Вопрос:
меня зовут НОЗОМИ. Я биостатист в Японии.
Я впервые пытаюсь задавать вопросы на английском языке.
Я приношу извинения, если мой английский трудно понять.
Я использовал 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]