tests/confintrep.R

library(survey)
data(api)
dclus2<-svydesign(id=~dnum+snum, fpc=~fpc1+fpc2, data=apiclus2)
rclus2<-as.svrepdesign(dclus2)

m<-svyglm(I(comp.imp=="Yes")~1, design=dclus2, family=quasibinomial)
if(anyNA(confint(m, method="likelihood"))) stop("NA in confint")
mrep<-svyglm(I(comp.imp=="Yes")~1, design=rclus2, family=quasibinomial)
if(anyNA(confint(mrep))) stop("NA in confint")
if(anyNA(confint(mrep, method="likelihood"))) stop("NA in confint")

Try the survey package in your browser

Any scripts or data that you put into this service are public.

survey documentation built on April 9, 2024, 3:01 a.m.