GEE (Generalized Estimating Equation)를 사용하는 계수의 신뢰 간격

Question

geepack을 사용하는 일반화 된 추정 방정식을 사용하여 선형 회귀 모델을 실행 중입니다. confint(fit) 명령이 여기에서 작동하지 않는 것 같습니다. 예를 들면 :

f2 <- geeglm(FEV1 ~ Age, data = Hospdata, family=gaussian, id=HHID) 
summary(f2) 
confint(f2) 

내가 confint(f2)을 실행에 다음과 같은 오류 메시지가 얻을 :

> confint(f2) 
Waiting for profiling to be done... 
Error in `[.data.frame`(summ$coefficients, , "Std. Error", drop = FALSE) : undefined columns selected 

를 여기에서 신뢰 구간을 찾을 수있는 방법이 있나요?

Answer 1

뭔가 :

confint.geeglm <- function(object, parm, level = 0.95, ...) { 
    cc <- coef(summary(object)) 
    mult <- qnorm((1+level)/2) 
    citab <- with(as.data.frame(cc), 
        cbind(lwr=Estimate-mult*Std.err, 
         upr=Estimate+mult*Std.err)) 
    rownames(citab) <- rownames(cc) 
    citab[parm,] 
} 

confint(gee1) 

:

library(geepack) 
data(dietox) 
dietox$Cu  <- as.factor(dietox$Cu) 
mf1 <- formula(Weight~Cu*poly(Time,3)) 
gee1 <- geeglm(mf1, data=dietox, id=Pig, 
       family=poisson("identity"),corstr="ar1") 
cc <- coef(summary(gee1)) 
citab <- with(as.data.frame(cc), 
    cbind(lwr=Estimate-1.96*Std.err, 
      upr=Estimate+1.96*Std.err)) 
rownames(citab) <- rownames(cc) 

의 편의를 위해,이를 캡슐화하는 confint 방법을 쓸 수있다

Answer 2

confint는 통계 패키지입니다. geeglm은 geepack 패키지에 있습니다.

확실한 간격을 모델 출력에 저장할 위치를 결정해야합니다.

str(f2)을 사용하거나 summary(f2)에서 유도하십시오.

또한 모델 개체를 통해 f2$ 및 탭을 자동 완성하도록 살펴보십시오.

또한, 설명서 - link을 살펴보십시오. 내가 실행 한 예제 모델이 CI를 생성하지 않았기 때문에 스스로를 구성해야 할 수도 있습니다. 매개 변수 표준 오류에서 롤백해야 할 수도 있습니다. 이 같은