@@ -24,10 +24,16 @@ def target(target_val, pop, factor, value):
2424 weights * factors ["APOPSNR" ][year ],
2525 weights * factors ["ARETS" ][year ],
2626 )
27- single_returns = np .where ((data ["mars" ] == 1 ) & (data ["filer" ] == 1 ), s006 , 0 )
28- joint_returns = np .where ((data ["mars" ] == 2 ) & (data ["filer" ] == 1 ), s006 , 0 )
27+ single_returns = np .where (
28+ (data ["mars" ] == 1 ) & (data ["filer" ] == 1 ), s006 , 0
29+ )
30+ joint_returns = np .where (
31+ (data ["mars" ] == 2 ) & (data ["filer" ] == 1 ), s006 , 0
32+ )
2933 hh_returns = np .where ((data ["mars" ] == 4 ) & (data ["filer" ] == 1 ), s006 , 0 )
30- returns_w_ss = np .where ((data ["e02400" ] > 0 ) & (data ["filer" ] == 1 ), s006 , 0 )
34+ returns_w_ss = np .where (
35+ (data ["e02400" ] > 0 ) & (data ["filer" ] == 1 ), s006 , 0
36+ )
3137 dep_exemptions = (
3238 np .where (data ["mars" ] == 2 , data ["XTOT" ] - 2 , data ["XTOT" ] - 1 ) * s006
3339 )
@@ -45,27 +51,39 @@ def target(target_val, pop, factor, value):
4551 # wage distribution
4652 wage1 = np .where (data ["agi" ] <= 10000 , data ["e00200" ], 0 ) * s006
4753 wage2 = (
48- np .where ((data ["agi" ] > 10000 ) & (data ["agi" ] <= 20000 ), data ["e00200" ], 0 )
54+ np .where (
55+ (data ["agi" ] > 10000 ) & (data ["agi" ] <= 20000 ), data ["e00200" ], 0
56+ )
4957 * s006
5058 )
5159 wage3 = (
52- np .where ((data ["agi" ] > 20000 ) & (data ["agi" ] <= 30000 ), data ["e00200" ], 0 )
60+ np .where (
61+ (data ["agi" ] > 20000 ) & (data ["agi" ] <= 30000 ), data ["e00200" ], 0
62+ )
5363 * s006
5464 )
5565 wage4 = (
56- np .where ((data ["agi" ] > 30000 ) & (data ["agi" ] <= 40000 ), data ["e00200" ], 0 )
66+ np .where (
67+ (data ["agi" ] > 30000 ) & (data ["agi" ] <= 40000 ), data ["e00200" ], 0
68+ )
5769 * s006
5870 )
5971 wage5 = (
60- np .where ((data ["agi" ] > 40000 ) & (data ["agi" ] <= 50000 ), data ["e00200" ], 0 )
72+ np .where (
73+ (data ["agi" ] > 40000 ) & (data ["agi" ] <= 50000 ), data ["e00200" ], 0
74+ )
6175 * s006
6276 )
6377 wage6 = (
64- np .where ((data ["agi" ] > 50000 ) & (data ["agi" ] <= 75000 ), data ["e00200" ], 0 )
78+ np .where (
79+ (data ["agi" ] > 50000 ) & (data ["agi" ] <= 75000 ), data ["e00200" ], 0
80+ )
6581 * s006
6682 )
6783 wage7 = (
68- np .where ((data ["agi" ] > 75000 ) & (data ["agi" ] <= 100_000 ), data ["e00200" ], 0 )
84+ np .where (
85+ (data ["agi" ] > 75000 ) & (data ["agi" ] <= 100_000 ), data ["e00200" ], 0
86+ )
6987 * s006
7088 )
7189 wage8 = np .where (data ["agi" ] > 100_000 , data ["e00200" ], 0 ) * s006
@@ -118,15 +136,27 @@ def target(target_val, pop, factor, value):
118136 target_name = "SS_return"
119137 rhs_vars ["returns_w_ss" ] = targets [year ][target_name ] - returns_w_ss .sum ()
120138 target_name = "Dep_return"
121- rhs_vars ["dep_exemptions" ] = targets [year ][target_name ] - dep_exemptions .sum ()
122- rhs_vars ["interest" ] = target (targets [year ]["INTS" ], apopn , aints , interest .sum ())
123- rhs_vars ["dividend" ] = target (targets [year ]["DIVS" ], apopn , adivs , dividend .sum ())
139+ rhs_vars ["dep_exemptions" ] = (
140+ targets [year ][target_name ] - dep_exemptions .sum ()
141+ )
142+ rhs_vars ["interest" ] = target (
143+ targets [year ]["INTS" ], apopn , aints , interest .sum ()
144+ )
145+ rhs_vars ["dividend" ] = target (
146+ targets [year ]["DIVS" ], apopn , adivs , dividend .sum ()
147+ )
124148 rhs_vars ["biz_income" ] = target (
125149 targets [year ]["SCHCI" ], apopn , aschci , biz_income .sum ()
126150 )
127- rhs_vars ["biz_loss" ] = target (targets [year ]["SCHCL" ], apopn , aschcl , biz_loss .sum ())
128- rhs_vars ["cap_gain" ] = target (targets [year ]["CGNS" ], apopn , acgns , cap_gain .sum ())
129- rhs_vars ["pension" ] = target (targets [year ]["Pension" ], apopn , atxpy , pension .sum ())
151+ rhs_vars ["biz_loss" ] = target (
152+ targets [year ]["SCHCL" ], apopn , aschcl , biz_loss .sum ()
153+ )
154+ rhs_vars ["cap_gain" ] = target (
155+ targets [year ]["CGNS" ], apopn , acgns , cap_gain .sum ()
156+ )
157+ rhs_vars ["pension" ] = target (
158+ targets [year ]["Pension" ], apopn , atxpy , pension .sum ()
159+ )
130160 rhs_vars ["sch_e_income" ] = target (
131161 targets [year ]["SCHEI" ], apopn , aschei , sch_e_income .sum ()
132162 )
@@ -136,15 +166,33 @@ def target(target_val, pop, factor, value):
136166 rhs_vars ["ss_income" ] = target (
137167 targets [year ]["SS" ], apopsnr , asocsec , ss_income .sum ()
138168 )
139- rhs_vars ["ucomp" ] = target (targets [year ]["UCOMP" ], apopn , aucomp , ucomp .sum ())
140- rhs_vars ["wage1" ] = target (targets [year ]["wage1" ], apopn , awage , wage1 .sum ())
141- rhs_vars ["wage2" ] = target (targets [year ]["wage2" ], apopn , awage , wage2 .sum ())
142- rhs_vars ["wage3" ] = target (targets [year ]["wage3" ], apopn , awage , wage3 .sum ())
143- rhs_vars ["wage4" ] = target (targets [year ]["wage4" ], apopn , awage , wage4 .sum ())
144- rhs_vars ["wage5" ] = target (targets [year ]["wage5" ], apopn , awage , wage5 .sum ())
145- rhs_vars ["wage6" ] = target (targets [year ]["wage6" ], apopn , awage , wage6 .sum ())
146- rhs_vars ["wage7" ] = target (targets [year ]["wage7" ], apopn , awage , wage7 .sum ())
147- rhs_vars ["wage8" ] = target (targets [year ]["wage8" ], apopn , awage , wage8 .sum ())
169+ rhs_vars ["ucomp" ] = target (
170+ targets [year ]["UCOMP" ], apopn , aucomp , ucomp .sum ()
171+ )
172+ rhs_vars ["wage1" ] = target (
173+ targets [year ]["wage1" ], apopn , awage , wage1 .sum ()
174+ )
175+ rhs_vars ["wage2" ] = target (
176+ targets [year ]["wage2" ], apopn , awage , wage2 .sum ()
177+ )
178+ rhs_vars ["wage3" ] = target (
179+ targets [year ]["wage3" ], apopn , awage , wage3 .sum ()
180+ )
181+ rhs_vars ["wage4" ] = target (
182+ targets [year ]["wage4" ], apopn , awage , wage4 .sum ()
183+ )
184+ rhs_vars ["wage5" ] = target (
185+ targets [year ]["wage5" ], apopn , awage , wage5 .sum ()
186+ )
187+ rhs_vars ["wage6" ] = target (
188+ targets [year ]["wage6" ], apopn , awage , wage6 .sum ()
189+ )
190+ rhs_vars ["wage7" ] = target (
191+ targets [year ]["wage7" ], apopn , awage , wage7 .sum ()
192+ )
193+ rhs_vars ["wage8" ] = target (
194+ targets [year ]["wage8" ], apopn , awage , wage8 .sum ()
195+ )
148196
149197 model_vars = [
150198 "single_returns" ,
0 commit comments