as.svrepdesign {survey} | R Documentation

## Convert a survey design to use replicate weights

### Description

Creates a replicate-weights survey design object from a traditional strata/cluster survey design object. `JK1` and `JKn` are jackknife methods, `BRR` is Balanced Repeated Replicates and `Fay` is Fay's modification of this, `bootstrap` is Canty and Davison's bootstrap and `subbootstrap` is Rao and Wu's (n-1) bootstrap.

### Usage

    as.svrepdesign(design, type=c("auto", "JK1", "JKn", "BRR", "bootstrap","subbootstrap","Fay"),
    fay.rho = 0, ..., compress=TRUE)

### Arguments

`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 `jk1weights`, `jknweights`, `brrweights`, or `bootweights` .
`compress` | Use a compressed representation of the replicate weights matrix.

### Value

Object of class `svyrep.design`.

### See Also

`brrweights`, `svydesign`, `svrepdesign`, `bootweights`, `subbootweights`

### Examples

    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")

    ## with user-supplied hadamard matrix
    scd2brr1 <- as.svrepdesign(scdnofpc, type="BRR", hadamard.matrix=paley(11))

    svyratio(~alive, ~arrests, design=scd2brr)
    svyratio(~alive, ~arrests, design=scd2brr1)
    svyratio(~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)

[Package survey version 3.19 Index]