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", "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
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) svrepratio(~alive, ~arrests, design=scd2jkn) svrepratio(~alive, ~arrests, design=scd2jknf) data(api) dclus1<-svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc) rclus1<-as.svrepdesign(dclus1) svrepmean(~api00, rclus1) svreptotal(~enroll, rclus1)