as.svrepdesign {survey} | R Documentation |
Creates a replicate-weights survey design object from a traditional strata/cluster survey design object.
as.svrepdesign(design, type=c("auto", "JK1", "JKn", "BRR", "bootstrap","Fay"), fay.rho = 0, ..., compress=TRUE)
design |
Object of class survey.design |
type |
Type of replicate weights. "auto" uses JKn for
stratified, JK1 for unstratified designs |
fay.rho |
Tuning parameter for Fay's variance method |
... |
Other arguments to brrweights , controlling
the handling of strata that do not have exactly 2 clusters. |
compress |
Use a compressed representation of the replicate weights matrix. |
Object of class svrep.design
.
brrweights
, svydesign
,
svrepdesign
, bootweights
data(scd) scddes<-svydesign(data=scd, prob=~1, id=~ambulance, strata=~ESA, nest=TRUE, fpc=rep(5,6)) scdnofpc<-svydesign(data=scd, prob=~1, id=~ambulance, strata=~ESA, nest=TRUE) # convert to BRR replicate weights scd2brr <- as.svrepdesign(scdnofpc, type="BRR") scd2fay <- as.svrepdesign(scdnofpc, type="Fay",fay.rho=0.3) # convert to JKn weights scd2jkn <- as.svrepdesign(scdnofpc, type="JKn") # convert to JKn weights with finite population correction scd2jknf <- as.svrepdesign(scddes, type="JKn") svrepratio(~alive, ~arrests, design=scd2brr) svrepratio(~alive, ~arrests, design=scd2fay) svyratio(~alive, ~arrests, design=scd2jkn) svyratio(~alive, ~arrests, design=scd2jknf) data(api) ## one-stage cluster sample dclus1<-svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc) ## convert to JK1 jackknife rclus1<-as.svrepdesign(dclus1) ## convert to bootstrap bclus1<-as.svrepdesign(dclus1,type="bootstrap", replicates=100) svymean(~api00, dclus1) svytotal(~enroll, dclus1) svymean(~api00, rclus1) svytotal(~enroll, rclus1) svymean(~api00, bclus1) svytotal(~enroll, bclus1)