소수의 사이트에서 일련의 "브로큰 스틱"월간 사례 수로 구성된 데이터 세트가 있습니다. 두 가지 기술로 단일 요약 견적을 얻으려고합니다.
기술 1 : 0/1 표시기 변수가있는 Poisson GLM에 "브로큰 스틱"을 맞추고 시간 및 시간 ^ 2 변수를 사용하여 시간 추세를 제어합니다. 0/1 인디케이터 변수의 추정치와 SE는 꽤 직선적 인 위 / 아래 모멘트 방법을 사용하거나 R의 tlnise 패키지를 사용하여 "Bayesian"추정치를 얻습니다. 이것은 Peng 및 Dominici가 대기 오염 데이터를 사용하여 수행하는 것과 유사하지만 더 적은 수의 사이트 (약 12 개)를 사용합니다.
기술 2 : 시간별 추세에 대한 일부 사이트 별 제어를 포기하고 선형 혼합 모델을 사용하십시오. 특별히:
lmer(cases ~ indicator + (1+month+I(month^2) + offset(log(p)), family="poisson", data=data)
내 질문에는 이러한 추정치에서 나오는 표준 오류가 포함됩니다. 기법 1의 표준 오차는 실제로 매월 설정 한 시간이 아닌 주 단위를 사용하므로 정확도 가 더 높아야합니다.이 방법은 모멘트 방법의 경우 ~ 0.206, tlnise의 경우 ~ 0.306의 표준 오차를 갖습니다.
lmer 방법은 ~ 0.09의 표준 오류를 제공합니다. 효과 추정치가 합리적으로 가깝기 때문에 혼합 모형이 훨씬 더 효율적일 때 서로 다른 요약 추정치에서 제로화되는 것은 아닙니다.
기대하는 것이 합리적입니까? 그렇다면 왜 혼합 모델이 훨씬 더 효율적인가? 이것이 일반적인 현상입니까, 아니면이 모델의 특정 결과입니까?