getSampleVarTotalSTSI:

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

Usage

getSampleVarTotalSTSI(sampleDesign, y)

Arguments

sampleDesign
y

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 (sampleDesign, y) 
{
    getStratumVarTotal <- function(sample) {
        N <- sum(sample$weights)
        n <- nrow(sample)
        return(N^2 * (1 - n/N) * var(sample$y)/n)
    }
    sample <- data.frame(y = y, weights = sampleDesign$weights)
    strataVar <- sapply(split(sample, sampleDesign$strata), getStratumVarTotal)
    sampleVar <- sum(strataVar)
    return(sampleVar)
  }

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