-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathinitdata.jl
254 lines (244 loc) · 9.42 KB
/
initdata.jl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
#===variable descriptions===#
# parameters and counters
# alpha # capital-share parameter in production function
# delta # depreciation rate
# sigL # elasticity of hours-disutility
# fag # first economically active age group
# rho # discount rate households
# sigma # elasticity of intertemporal substitution
# tt # time index
# age-dependent variables
# parlv0 # multiplicative shift parameter hours-disutillity
# parlv1 # additive shift parameter hours-disutillity
# time-dependent variables
# A # aggregate assets
# Cons # aggregate consumption
# CG # aggregate public consumption
# DG # public debt
# Div # dividends
# Exp # government expenditure
# Inv # investment
# K # capital stock
# LD # aggregate labor demand
# LS # aggregate labor supply
# N # aggregate population
# NB # number of newborns
# Nc # number of children
# Nr # number of retirees
# Nw # number of workers
# P # pension expenditure
# PB # primary balance
# Rev # government revenue
# TaxF # payroll tax revenue
# TFP # TFP stock
# V # value of representative firm
# Y # output (GDP)
# ab # aggregate accidental bequests
# eda # excess demand asset market
# edab # (excess demand) accidental bequests
# edg # (excess demand) government budget
# ediv # (excess demand) intervivo transfers
# edl # excess demand labor market
# edw # (excess demand) Walras' Law
# edy # excess demand goods market
# iv # aggregate received intervivo transfers
# pc # price of consumption
# qTob # Tobin's q
# r # real interest rate
# rag # retirement age group
# tauC # consumption tax rate
# tauF # payroll tax rate
# tauW # average wage tax rate
# taul # lump-sum tax rate
# tauprof # corporate tax rate
# uck # user cost of capital
# w # wage rate
# age-dependent and time-dependent variables (with v and z suffix)
# A # assets
# Cons # consumption
# HH_nonconv # counter for non-converging households
# N # population mass
# Sav # end-of-period savings
# ab # accidental bequest
# dis_tot # disutility of labor
# cG # public consumption
# ell # hours worked
# gam # conditional survival probability
# iv # intervivo transfers
# lambda # shadow price of assets
# notret # not retired indicator
# p # pension income
# pc # price of consumption
# r # real interest rate
# tauC # consumption tax rate
# tauW # wage tax rate
# taul # lump-sum tax rate
# theta # productivity (age-specific)
# w # wage rate
# y # per-period labor and pension income
#=====initialization=====#
const ncoh = tend + (nag-1)
# parameters and counters
alpha = 0.0
delta = 0.0
sigL = 0.0
fag = 0.0
rho = 0.0
sigma = 0.0
tt = 0.0
# age-dependent variables
parlv0 = zeroscol(nag)
parlv1 = zeroscol(nag)
# time-dependent variables
A = zerosrow(tend)
A0 = 0.0
Cons = zerosrow(tend)
Cons0 = 0.0
CG = zerosrow(tend)
CG0 = 0.0
DG = zerosrow(tend)
DG0 = 0.0
Div = zerosrow(tend)
Div0 = 0.0
Exp = zerosrow(tend)
Exp0 = 0.0
Inv = zerosrow(tend)
Inv0 = 0.0
K = zerosrow(tend)
K0 = 0.0
LD = zerosrow(tend)
LD0 = 0.0
LS = zerosrow(tend)
LS0 = 0.0
N = zerosrow(tend)
N0 = 0.0
NB = zerosrow(tend)
NB0 = 0.0
Nc = zerosrow(tend)
Nc0 = 0.0
Nr = zerosrow(tend)
Nr0 = 0.0
Nw = zerosrow(tend)
Nw0 = 0.0
P = zerosrow(tend)
P0 = 0.0
PB = zerosrow(tend)
PB0 = 0.0
Rev = zerosrow(tend)
Rev0 = 0.0
TaxF = zerosrow(tend)
TaxF0 = 0.0
TFP = zerosrow(tend)
TFP0 = 0.0
V = zerosrow(tend)
V0 = 0.0
Y = zerosrow(tend)
Y0 = 0.0
ab = zerosrow(tend)
ab0 = 0.0
eda = zerosrow(tend)
eda0 = 0.0
edab = zerosrow(tend)
edab0 = 0.0
edg = zerosrow(tend)
edg0 = 0.0
ediv = zerosrow(tend)
ediv0 = 0.0
edl = zerosrow(tend)
edl0 = 0.0
edw = zerosrow(tend)
edw0 = 0.0
edy = zerosrow(tend)
edy0 = 0.0
iv = zerosrow(tend)
iv0 = 0.0
pc = zerosrow(tend)
pc0 = 0.0
qTob = zerosrow(tend)
qTob0 = 0.0
r = zerosrow(tend)
r0 = 0.0
rag = zerosrow(tend)
rag0 = 0.0
tauC = zerosrow(tend)
tauC0 = 0.0
tauF = zerosrow(tend)
tauF0 = 0.0
tauW = zerosrow(tend)
tauW0 = 0.0
taul = zerosrow(tend)
taul0 = 0.0
tauprof = zerosrow(tend)
tauprof0 = 0.0
uck = zerosrow(tend)
uck0 = 0.0
w = zerosrow(tend)
w0 = 0.0
# age-dependent and time-dependent variables
Av = zerosmat(nag,tend)
Az = zerosmat(nag,ncoh)
Av0 = zeroscol(nag)
Consv = zerosmat(nag,tend)
Consz = zerosmat(nag,ncoh)
Consv0 = zeroscol(nag)
HH_nonconvv = zerosmat(nag,tend)
HH_nonconvz = zerosmat(nag,ncoh)
HH_nonconvv0 = zeroscol(nag)
Nv = zerosmat(nag,tend)
Nz = zerosmat(nag,ncoh)
Nv0 = zeroscol(nag)
Savv = zerosmat(nag,tend)
Savz = zerosmat(nag,ncoh)
Savv0 = zeroscol(nag)
abv = zerosmat(nag,tend)
abz = zerosmat(nag,ncoh)
abv0 = zeroscol(nag)
dis_totv = zerosmat(nag,tend)
dis_totz = zerosmat(nag,ncoh)
dis_totv0 = zeroscol(nag)
cGv = zerosmat(nag,tend)
cGz = zerosmat(nag,ncoh)
cGv0 = zeroscol(nag)
ellv = zerosmat(nag,tend)
ellz = zerosmat(nag,ncoh)
ellv0 = zeroscol(nag)
gamv = zerosmat(nag,tend)
gamz = zerosmat(nag,ncoh)
gamv0 = zeroscol(nag)
ivv = zerosmat(nag,tend)
ivz = zerosmat(nag,ncoh)
ivv0 = zeroscol(nag)
lambdav = zerosmat(nag,tend)
lambdaz = zerosmat(nag,ncoh)
lambdav0 = zeroscol(nag)
notretv = zerosmat(nag,tend)
notretz = zerosmat(nag,ncoh)
notretv0 = zeroscol(nag)
pv = zerosmat(nag,tend)
pz = zerosmat(nag,ncoh)
pv0 = zeroscol(nag)
pcv = zerosmat(nag,tend)
pcz = zerosmat(nag,ncoh)
pcv0 = zeroscol(nag)
rv = zerosmat(nag,tend)
rz = zerosmat(nag,ncoh)
rv0 = zeroscol(nag)
tauCv = zerosmat(nag,tend)
tauCz = zerosmat(nag,ncoh)
tauCv0 = zeroscol(nag)
tauWv = zerosmat(nag,tend)
tauWz = zerosmat(nag,ncoh)
tauWv0 = zeroscol(nag)
taulv = zerosmat(nag,tend)
taulz = zerosmat(nag,ncoh)
taulv0 = zeroscol(nag)
thetav = zerosmat(nag,tend)
thetaz = zerosmat(nag,ncoh)
thetav0 = zeroscol(nag)
wv = zerosmat(nag,tend)
wz = zerosmat(nag,ncoh)
wv0 = zeroscol(nag)
yv = zerosmat(nag,tend)
yz = zerosmat(nag,ncoh)
yv0 = zeroscol(nag)
;