vhat2: Estimated variance of a total estimation, with strata

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

Estimated variance of a total estimation, with strata

Description

Estimated variance of a total estimation, with strata

Usage

vhat2(data, weights, mh, Nh, strata = NULL)

Arguments

data	un data frame con il sottoinsieme dei dati per cui si vuole stimare la varianza della stima di un totale. Da notare la differenza rispetto all'argomento data della funzione vhat_strata, dove data rappresenta tutto il dataset, mentre qui rappresenta solo il sottoinsieme dei dati la funzione indicatrice vale uno.
weights	stringa con il nome della colonna dei dati che contiene i pesi individuali
mh	nome della colonna dei dati che contiene il numero di osservazioni dello strato. Deve essere un vettore numerico con lo stesso valore ripetuto il numero di righe dei dati.
Nh	nome della colonna dei dati che contiene la stima degli effettivi dello strato. Come mh, deve essere lo stesso valore ripetuto n volte, dove n è il numero di righe dei dati.
strata	character string with the name of the column that contains the strata

Examples

str(d)

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
str(d)
vhat2(data = d[d$xcat == "cat2", ], weights = "w", mh = "mh", Nh = "Nh", strata = "s")

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