Description
Calibrates and Reweights Units in Samples.
Description
Provides user-friendly tools for calibration in survey sampling. The package is production-oriented, and its interface is inspired by the famous popular macro 'Calmar' for SAS, so that 'Calmar' users can quickly get used to 'icarus'. In addition to calibration (with linear, raking and logit methods), 'icarus' features functions for calibration on tight bounds and penalized calibration.
README.md
Icarus
Icarus (Icarus Calibrates And Reweights Units in Samples) is an R package providing useful functions for calibration and reweighting estimators in survey sampling. The former name of this package was gaston.
Cite
To cite Icarus in publications use: Rebecq, Antoine (2017). Icarus: an R package for calibration in survey sampling. R package version 0.3.0.
Install
You can use the following instruction to install icarus (from CRAN):
install.packages("icarus")
However, if you wish to install the latest version of icarus, you can use devtools and install directly from this github repo:
install.packages("devtools")
library(devtools)
install_github("haroine/icarus")
Short example
In this example, we perform calibration (with the "raking" method) on the test dataset data_employees included in icarus:
library(icarus)
N <- 300 ## Population size
## Compute the Horvitz-Thompson estimator (returns 1.666667)
weightedMean(data_employees$movies, data_employees$weight, N)
## Add calibration margins
mar1 <- c("category",3,80,90,60)
mar2 <- c("sex",2,140,90,0)
mar3 <- c("department",2,100,130,0)
mar4 <- c("salary", 0, 470000,0,0)
margins <- rbind(mar1, mar2, mar3, mar4)
## Compute calibration weights
wCal <- calibration(data=data_employees, marginMatrix=margins, colWeights="weight"
, method="raking", description=FALSE)
## Value of the calibrated estimator: 2.471917
weightedMean(data_employees$movies, wCal, N)