Example #1
0
def income_prob_density_func(gender, race, age):
	(inc_med, inc_mean) = income_table.get_income(gender, race, age)
	k = (.2*(inc_med/inc_mean) + .8)/(3-3*(inc_med/inc_mean))
	theta = inc_mean/k
#	func = lambda x: x**(k-1)*exp(-x/theta)/(gamma(x)*theta**k)
	func = lambda x: x**(k-1)*exp(-x/theta)/(Gamma(x)*theta**k)
	return func
Example #2
0
def income_cumulative_prob_inverse(gender, race, age):
	(inc_med, inc_mean) = income_table.get_income(gender, race, age)
	k = (.2*(inc_med/inc_mean) + .8)/(3-3*(inc_med/inc_mean))
	theta = inc_mean/k
	func = lambda prob: gammaincinv(k, prob)*theta
	return func
Example #3
0
def income_cumulative_prob_func(gender, race, age):
	(inc_med, inc_mean) = income_table.get_income(gender, race, age)
	k = (.2*(inc_med/inc_mean) + .8)/(3-3*(inc_med/inc_mean))
	theta = inc_mean/k
	func = lambda x: gammainc(k, x/theta)
	return func