vardomh(Y, H, PSU, w_final, ID_household,
id = NULL, Dom = NULL, period = NULL,
N_h = NULL, fh_zero = FALSE, PSU_level = TRUE,
Z = NULL, dataset = NULL,
X = NULL, periodX = NULL, X_ID_household = NULL,
ind_gr = NULL, g = NULL, datasetX = NULL,
q = rep(1, if (is.null(datasetX))
nrow(data.frame(X)) else nrow(datasetX)),
confidence = .95, outp_lin = FALSE,
outp_res = FALSE)data.frame or variable names as character, column numbers or logical vector with only one TRUE value (length of the vector has to be the same as the column count of datasetdata.frame or variable name as character, column number or logical vector with only one TRUE value (length of the vector has to be the same as the columndata.frame or variable name as character, column number or logical vector with only one TRUE value (length of the vector has to be the same as the cdata.frame or variable name as character, column number or logical vector with only one TRUE value (length of the vector has to be the same as the column count of data.frame or variable name as character, column number or logical vector with only one TRUE value (length of the vector has to be the same as the data.frame or variable name as character, column number or logical vector with only one TRUE value (length of the vector has to be the same as data.frame or variable names as character, column numbers or logical vector (length of the vector has to be thdata.frame or variable names as character vector, column numbers or logical vector (length of the vector has todata.frame or variable names as character, column numbers or logical vector (length of the vector has to be the same as the column count of dataset).data.frame.data.frame or variable names as character, column numbers or logical vector (length of the vector has to be the same as the column count of data.frame or variable names as character, column numbers or logical vector (length of tdata.frame or variable name as character, column number or logical vector with only one TRUE value (length of the vector has to be the same as the data.frame or variable name as character, column number or logical vector with only one TRUE value (length of data.frame or variable name as character, column number or logical vector with only one TRUE value (length of the vector has to be the same as tdata.frame.data.frame or variable name as character, column number or logical vector with only one TRUE value (length of theTRUE linearized values of the ratio estimator will be printed out.TRUE estimated residuals of calibration will be printed out.data.table containing the linearized values of the ratio estimator with id and PSU.data.table containing the estimated residuals of calibration with id and PSU.data.table, which containing variables:
respondent_count - the count of respondents,
pop_size - the estimated size of population,
n_nonzero - the count of respondents, who answers are larger than zero,
estim - the estimated value,
var - the estimated variance,
se - the estimated standard error,
rse - the estimated relative standard error (coefficient of variation),
cv - the estimated relative standard error (coefficient of variation) in percentage,
absolute_margin_of_error - the estimated absolute margin of error,
relative_margin_of_error - the estimated relative margin of error,
CI_lower - the estimated confidence interval lower bound,
CI_upper - the estimated confidence interval upper bound,
var_srs_HT - the estimated variance of the HT estimator under SRS,
var_cur_HT - the estimated variance of the HT estimator under current design,
var_srs_ca - the estimated variance of the calibrated estimator under SRS,
deff_sam - the estimated design effect of sample design,
deff_est - the estimated design effect of estimator,
deff - the overall estimated design effect of sample design and estimator,
n_eff - the effective sample size.domain, lin.ratio, residual_est,
var_srs, variance_estdata(eusilc)
dataset <- data.frame(1:nrow(eusilc),eusilc)
colnames(dataset)[1] <- "IDd"
aa<-vardomh(Y="eqIncome", H="db040", PSU="db030", w_final="rb050",
ID_household="db030", id="rb030", Dom = "db040",
period = NULL, N_h = NULL, Z = NULL, dataset = dataset, X = NULL,
X_ID_household = NULL, g = NULL, datasetX = NULL,
q = rep(1, if (is.null(datasetX))
nrow(as.data.frame(H)) else nrow(datasetX)),
confidence = .95, outp_lin=TRUE, outp_res=TRUE)
dataset2 <- copy(dataset)
dataset$period=1
dataset2$period=2
dataset <-data.frame(rbind(dataset, dataset2))
aa2<-vardomh(Y="eqIncome", H="db040", PSU="db030", w_final="rb050",
ID_household="db030", id="rb030", Dom = "db040",
period = "period", N_h = NULL, Z = NULL, dataset = dataset,
X = NULL, X_ID_household = NULL,
g = NULL, datasetX = NULL,
q = rep(1, if (is.null(datasetX))
nrow(as.data.frame(H)) else nrow(datasetX)),
confidence = .95, outp_lin=TRUE, outp_res=TRUE)Run the code above in your browser using DataLab