data(eusilc)
dati <- data.table(IDd=1:nrow(eusilc),eusilc)
# Full population
dd <- linqsr(Y="eqIncome", id="IDd", weight="rb050",
Dom=NULL, dataset= dati, alpha=20)
dd$value
# By domains
dd <- linqsr(Y="eqIncome", id="IDd", weight="rb050",
Dom="db040", dataset= dati, alpha=20)
dd$value
Run the code above in your browser using DataLab