2017-12-22 26 views
0

는, 나는 다음과 같은 명령을 사용하면, 내 종속 변수 (thkbins) 내 두 개의 예측 변수 (cc & tv 사이에 가능한 각 조합에 대한 logits을 얻을 수 있다는 사실을 알고) :R은 + 다단계 로지스틱 회귀 + 그룹 고정 효과 STATA에서

melogit thkbins cc#tv || school:, 

enter image description here

R 비슷한 출력을 생성하는 방법이 있나요? 나는 glmer 명령을 lme4 패키지에서 사용하고 있으며 상호 작용 용어로 출력을 얻을 수는 있지만 Stata에서 정확히 생성 할 수있는 것은 아닙니다.

model1 <- glmer(thkbin ~ cc + tv + cc*tv + (1|school), 
      data=thkdata, family = binomial, nAGQ = 7) 
summary(model1) 

enter image description here

답변

1

나는 패키지 ordinal (튜토리얼 here)에서 clmm을 사용 :

model<-clmm(DepVar~IndVar+(1|WithinVar),data=df) 

도움이 되었기를 바랍니다.