이것은 나의 첫 번째 게시물입니다. 이 커뮤니티에 진심으로 감사드립니다.
0으로 잘린 종단 반응 (응답 변수 = 0이 될 확률)과 평균! = 분산을 분석하려고하므로 포아송보다 음의 이항 분포가 선택되었습니다.
내가 배제한 기능 / 명령 :
아르 자형
- R의 gee () 함수는 제로 잘림이나 음의 이항 분포를 고려하지 않습니다 (MASS 패키지가로드 된 경우에도)
- R의 glm.nb ()는 다른 상관 관계 구조를 허용하지 않습니다.
- VGAM 패키지의 vglm ()은 posnegbinomial 제품군을 사용할 수 있지만 Stata의 ztnb 명령 (아래 참조)과 동일한 문제가 있습니다.
스타 타
- 데이터가 종단이 아닌 경우 Stata 패키지 ztnb를 사용하여 분석을 실행할 수 있지만 해당 명령은 내 관측치가 독립적이라고 가정합니다.
또한 다양한 방법 론적 / 철학적 이유로 GLMM을 배제했습니다.
지금은 Stata의 xtgee 명령에 동의했습니다. xtgee 사용의 추가 이점은 qic 값 (qic 명령 사용)을 계산하여 응답 변수에 가장 적합한 상관 관계 구조를 결정할 수도 있다는 것입니다.
R 또는 Stata에 1) nbinomial 제품군, 2) GEE 및 3) 잘림을 고려할 수있는 패키지 / 명령이 있으면 알고 싶지 않습니다.
당신이 가질 수있는 아이디어를 대단히 감사하겠습니다. 감사합니다.
케이스