getSampleStrata:

In BERENZ/RISQ: What the package does (short line)

Usage

getSampleStrata(sampleWeights, nMaxStrata = 20)

Arguments

sampleWeights
nMaxStrata

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (sampleWeights, nMaxStrata = 20) 
{
    weights <- sort(unique(sampleWeights))
    indices <- seq(along = sampleWeights)
    strata <- factor(seq(along = weights))
    if (length(weights) <= nMaxStrata) {
        values <- merge(data.frame(weight = sampleWeights, index = indices), 
            data.frame(weight = weights, stratum = strata), all.x = TRUE)
        sampleStrata <- values[order(values$index), "stratum"]
    }
    else sampleStrata <- factor(rep(1, length(sampleWeights)))
    return(sampleStrata)
  }

BERENZ/RISQ documentation built on May 5, 2019, 10:25 a.m.