atkinson.md {acid}

Atkinson Index for a Mixture of Income Distributions

Description

This function uses Monte Carlo methods to estimate the Atkinson index for a mixture of two continuous income distributions and a point mass for zero-incomes.

Usage

atkinson.md(n, epsilon = 1, dist1, dist2, theta, p0, p1, p2, 
dist.para.table, zero.approx)

Arguments

n

sample size used to estimate the Atkinson index.

epsilon

inequality aversion parameter as denoted by Atkinson (1970). The default is epsilon=1.

dist1

character string with the name of the first continuous distribution used. Must be listed in dist.para.table. Must be equivalent to the respective function of that distribution, e.g. norm for the normal distribution.

dist2

character string with the name of the second continuous distribution used. Must be listed in dist.para.table. Must be equivalent to the respective function of that distribution, e.g. norm for the normal distribution.

theta

vector with the parameters of dist1 and dist2. Order must be the same as in the functions for the distributions.

p0

scalar with probability mass for the point mass.

p1

scalar with probability mass for dist1.

p2

scalar with probability mass for dist2.

dist.para.table

a table of the same form as dist.para.t with distribution name, function name and number of parameters.

zero.approx

a scalar which replaces zero-incomes, such that the Atkinson index involving a logarithm return finite values.

Value

AIM

the selected Atkinson inequality measure.

epsilon

the inequality aversion parameter used.

y

a vector with the simulated incomes to estimate the entropy measure.

y2

a vector with the zero-replaced simulated incomes to estimate the entropy measure.

zero.replace

a logical vector indicating whether a zero has been replaced.

stat

a vector with the simulated group the observation was chosen from. 0 is the point mass, 1 dist1 and 2 dist2.

Author(s)

Alexander Sohn

References

Atkinson, A.B. (1970): On the Measurment of Inequality, in: Journal of Economic Theory, Vol. 2(3), pp. 244-263.

See Also

ineq , atkinson , atkinson.den

Examples

theta<-c(2,1,5,2)
x<- c(rgamma(50000,2,1),rgamma(50000,5,2))
para<-1

data(dist.para.t)
atkinson.md(10000,para,"gamma","gamma",theta,0,0.5,0.5,dist.para.t,zero.approx=1)$AIM
atkinson(x,1)