data(eusilcS)
data(eusilcP)
## Not run:
# ## approx. 20 seconds computation time
# inp <- specifyInput(data=eusilcS, hhid="db030", hhsize="hsize", strata="db040", weight="db090")
# inp <- simStructure(data=inp, method="direct", basicHHvars=c("age", "rb090"))
# inp <- simCategorical(inp, additional=c("pl030", "pb220a"), method="multinom",nr_cpus=1)
#
# margins <- as.data.frame(
# xtabs(rep(1, nrow(eusilcP)) ~ eusilcP$region + eusilcP$gender + eusilcP$citizenship))
# colnames(margins) <- c("db040", "rb090", "pb220a", "freq")
# inp <- addKnownMargins(inp, margins)
# str(inp)
# ## End(Not run)
Run the code above in your browser using DataLab