|
1 | | -import json |
2 | 1 | import numpy as np |
3 | 2 | import pandas as pd |
4 | 3 | import os |
@@ -1093,19 +1092,22 @@ def model_fit_table( |
1093 | 1092 | model_val = tpi_output["r"][t] |
1094 | 1093 | elif target_desc == r"Capital share of output": |
1095 | 1094 | model_val = ( |
1096 | | - 1 - tpi_output["r"][t] * tpi_output["K"][t] / tpi_output["Y"][t] |
| 1095 | + 1 |
| 1096 | + - tpi_output["r"][t] * tpi_output["K"][t] / tpi_output["Y"][t] |
1097 | 1097 | ) |
1098 | 1098 | elif target_desc == r"Labor share of output": |
1099 | | - model_val = tpi_output["w"][t] * tpi_output["L"][t] / tpi_output["Y"][t] |
1100 | | - # Fiscal moments |
1101 | | - elif target_desc == r"Revenue to GDP ratio $(T/Y)$": |
1102 | 1099 | model_val = ( |
1103 | | - tpi_output["total_tax_revenue"][t] / tpi_output["Y"][t] |
| 1100 | + tpi_output["w"][t] * tpi_output["L"][t] / tpi_output["Y"][t] |
1104 | 1101 | ) |
| 1102 | + # Fiscal moments |
| 1103 | + elif target_desc == r"Revenue to GDP ratio $(T/Y)$": |
| 1104 | + model_val = tpi_output["total_tax_revenue"][t] / tpi_output["Y"][t] |
1105 | 1105 | elif target_desc == r"Gov't consumption to GDP ratio $(G/Y)$": |
1106 | 1106 | model_val = tpi_output["G"][t] / tpi_output["Y"][t] |
1107 | 1107 | elif target_desc == r"Pension outlays to GDP ratio $(Pension/Y)$": |
1108 | | - model_val = tpi_output["agg_pension_outlays"][t] / tpi_output["Y"][t] |
| 1108 | + model_val = ( |
| 1109 | + tpi_output["agg_pension_outlays"][t] / tpi_output["Y"][t] |
| 1110 | + ) |
1109 | 1111 | elif target_desc == r"Infrastructure spending to GDP ratio $(I_g/Y)$": |
1110 | 1112 | model_val = tpi_output["I_g"][t] / tpi_output["Y"][t] |
1111 | 1113 | elif target_desc == r"Debt to GDP ratio $(D/Y)$": |
@@ -1141,7 +1143,9 @@ def model_fit_table( |
1141 | 1143 | # Demographic moments |
1142 | 1144 | elif target_desc == r"Fraction 65+": |
1143 | 1145 | model_val = ( |
1144 | | - params.omega[t, -35:].sum() # NOTE: not flexible with S, E changes |
| 1146 | + params.omega[ |
| 1147 | + t, -35: |
| 1148 | + ].sum() # NOTE: not flexible with S, E changes |
1145 | 1149 | / params.omega[t].sum() |
1146 | 1150 | ) |
1147 | 1151 | elif target_desc == r"Pop growth rate": |
|
0 commit comments