Learn R Programming

vardpoor (version 0.6.2)

linarpt: Linearization of at-risk-of-poverty threshold

Description

Estimates the at-risk-of-poverty threshold (defined as percentage (usualy 60%) of equalised disposable income after social transfers quantile (usualy median)) and computes linearized variable for variance estimation.

Usage

linarpt(Y,
  id = NULL,
  weight = NULL,
  sort = NULL,
  Dom = NULL,
  period = NULL,
  dataset = NULL,
  percentage = 60,
  order_quant=50,
  var_name="lin_arpt")

Arguments

Y
Study variable (for example equalised disposable income after social transfers). One dimensional object convertible to one-column data.table or variable name as character, column number.
id
Optional variable for unit ID codes. One dimensional object convertible to one-column data.table or variable name as character, column number.
weight
Optional weight variable. One dimensional object convertible to one-column data.table or variable name as character, column number.
sort
Optional variable to be used as tie-breaker for sorting. One dimensional object convertible to one-column data.table or variable name as character, column number.
Dom
Optional variables used to define population domains. If supplied, linearization of at-risk-of-poverty threshold is done for each domain. An object convertible to data.table or variable names as character vector, column numbers as numeric vec
period
Optional variable for survey period. If supplied, linearization of at-risk-of-poverty threshold is done for each survey period. Object convertible to data.table or variable names as character, column numbers as numeric vector.
dataset
Optional survey data object convertable to data.table.
percentage
A numeric value in range $\left[ 0,100 \right]$ for $p$ in the formula for at-risk-of-poverty threshold computation: $$\frac{p}{100} \cdot Z_{\frac{\alpha}{100}}.$$ For example, to compute poverty threshold equal to 60% of some income quantile, $p$ shou
order_quant
A numeric value in range $\left[ 0,100 \right]$ for $\alpha$ in the formula for at-risk-of-poverty threshold computation: $$\frac{p}{100} \cdot Z_{\frac{\alpha}{100}}.$$ For example, to compute poverty threshold equal to some percentage of median income
var_name
A character specifying the name of the linearized variable.

Value

  • A list with three objects are returned:
  • quantileA data.table containing the estimated value of the quintile used for at-risk-of-poverty threshold estimation.
  • valueA data.table containing the estimated at-risk-of-poverty threshold (in percentage).
  • linA data.table containing the linearized variables of the at-risk-of-poverty threshold (in percentage).

Details

The implementation strictly follows the Eurostat definition.

References

Working group on Statistics on Income and Living Conditions (2004) Common cross-sectional EU indicators based on EU-SILC; the gender pay gap. EU-SILC 131-rev/04, Eurostat. Guillaume Osier (2009). Variance estimation for complex indicators of poverty and inequality. Journal of the European Survey Research Association, Vol.3, No.3, pp. 167-195, ISSN 1864-3361, URL http://ojs.ub.uni-konstanz.de/srm/article/view/369. Jean-Claude Deville (1999). Variance estimation for complex statistics and estimators: linearization and residual techniques. Survey Methodology, 25, 193-203, URL http://www5.statcan.gc.ca/bsolc/olc-cel/olc-cel?lang=eng&catno=12-001-X19990024882.

See Also

linarpr, incPercentile, varpoord , vardcrospoor, vardchangespoor

Examples

Run this code
data(eusilc)

dati <- data.table(IDd = 1:nrow(eusilc), eusilc)

# Full population
d1 <- linarpt(Y="eqIncome", id="IDd", weight = "rb050", Dom = NULL,
             dataset = dati, percentage = 60, order_quant=50)
d1$value

# By domains
d2 <- linarpt(Y="eqIncome", id="IDd", weight = "rb050", Dom = "db040",
              dataset = dati, percentage = 60, order_quant=50)
d2$value

Run the code above in your browser using DataLab