vhat3: Estimated variance of a total estimation, with strata (vector...

In gibonet/vhatbfs: Implementation of Confidence Intervals for Swiss Structural Survey Data

Estimated variance of a total estimation, with strata (vector version)

Description

Estimated variance of a total estimation, with strata (vector version)

Usage

vhat3(x, mh, Nh, strata)

Arguments

x	un vettore numerico con gli elementi di uno strato (di solito i pesi di campionamento)
mh	vettore con il numero di osservazioni dello strato (può essere un vettore di lunghezza uno oppure un vettore più lungo ma con un valore unico ripetuto).
Nh	vettore con la stima degli effettivi dello strato (può essere un vettore di lunghezza uno oppure un vettore più lungo ma con un valore unico ripetuto). In pratica è la somma dei pesi di campionamento dello strato.
strata	vettore che contiene le informazioni relative agli strati

Examples

tmp <- n_strata(d, weights = "w", strata = "s")
tmp$s <- 1L:3L
d <- merge(d, tmp, by = "s", all.x = TRUE)

# Adesso d ha le colonne mh e Nh
head(d)

k <- with(d, xcat == "cat2")
vhat3(d$w[k], d$mh[k], d$Nh[k], d$s[k])

gibonet/vhatbfs documentation built on Nov. 24, 2024, 3:25 p.m.