2015-07-16

R로 Linear Mixed-effects Model 통계 분석

x가 y에 영향을 미치는 지 알고 싶어서 실험을 합니다.

여러 subjects로 실험을 하는데, 그런데 subject마다 y의 평균값이 높기도 하고 낮기도 하고 하는 식으로 다를 수 있습니다.

하지만 어쨋거나 x가 y에 linear하게 영향을 미치는지 확인하고 싶습니다.

이럴때, Linear Mixed-effects Model로 분석을 합니다.

이때 x가 fixed effect이고 subjectID는 random effect입니다.

R에서 하는 법:
  1. R에 lme4 package를 사용합니다: library(lme4)라고 입력합니다.
  2. lmer함수를 사용합니다: lmer(y~x+ (1|subjectID), data=dtableName)
단순 linear model이아니고 generalized linear model로 분석하고 싶으면,
lmer대신 glmer을 쓰면 되고요. 이때는 family를 세팅해줘야합니다.