documenting run_simulation.py

Python/CURRENT_DEBUGGED_VERSION/run_simulation.py

@@ -1,6 +1,6 @@
 '''
 ------------------------------------------------------------------------
-Last updated 6/19/2015
+Last updated 7/16/2015
 
 This will run the steady state solver as well as time path iteration.
 
@@ -46,7 +46,6 @@
 J            = number of different ability groups
 T            = number of time periods until steady state is reached
 lambdas  = desired percentiles of ability groups
-scal         = scalar multiplier used in SS files to make the initial value work
 starting_age = age of first members of cohort
 ending age   = age of the last members of cohort
 E            = number of cohorts before S=1
@@ -62,9 +61,6 @@
                period
 g_y_annual   = annual growth rate of technology
 g_y          = growth rate of technology for one cohort
-slow_work    = time at which chi_n starts increasing from 1
-chi_n_multiplier = scalar which is increased to force the labor
-               distribution to 0
 maxiter   = Maximum number of iterations that TPI will undergo
 mindist_SS   = Cut-off distance between iterations for SS
 mindist_TPI   = Cut-off distance between iterations for TPI

Python/OUTPUT/SSinit/ss_init_tpi_vars.pkl

@@ -0,0 +1,43 @@
+(dp1
+S'nssmat_init'
+p2
+cnumpy.core.multiarray
+_reconstruct
+p3
+(cnumpy
+ndarray
+p4
+(I0
+tS'b'
+tRp5
+(I1
+(L10L
+L2L
+tcnumpy
+dtype
+p6
+(S'f8'
+I0
+I1
+tRp7
+(I3
+S'<'
+NNNI-1
+I-1
+I0
+tbI00
+S"\xb2|\xd2{\x98\xb0\xce?IURIA\xff\x88?\x11\xbdj?\x1e\xf7\xca?b(\x8fm\x155k?\xf4\x8diR\xd9c\xc9?;\xff\xf9in\xa9i?\x1c\xa7v\xe7\x1dX\xc8?\xb8H\xf43\xd7\x024?\xe2\xd1_\x82\xb7t\xc7?\xda\xdd\xc7\x8a\x7f\xa4_?FfP\xce9l\xc6?\xee\x1e\xea\xda\x1d\xd1T?W\xdf\xc2Bp*\xc6?\x8a3\n|\xb5\xadV?\x1dI3B\x08t\xc6?\xb7|A3\xbaHA?\x7f\xeb\x17Y\x97,\xc7?\xe5]r\xc0iQ'?\x16+\xad~$0\xc8?\xc4Ev\xef\xc3\x99!?"
+tbsS'bssmat_init'
+p8
+g3
+(g4
+(I0
+tS'b'
+tRp9
+(I1
+(L10L
+L2L
+tg7
+I00
+S'v\xcf\x18\xdc\x83\xf8\xc1?8\x1c\x10m\t\x1c\xf8?C\xb4\x05jw\xa5\xc6?\xb1\xe5\x02I\x82\xec\x07@\xf9\xcc__\xa8\xea\xc7?\x89\x81\xdcH1y|@W\xd4\xfd\xb3\x1d\xdc\xc7?\xb1\xdf\xce3\xae\x81{\xc0\x0e\x95j\x9e\xe4\xc6\xc6?n\xd1\x04\xfd\x12\xc7\x1e@\x1c\xa6\xb6\x1d[\x05\xcb?\x83\xe9\xaf\x15\xdb\xf3"@\xfeF\x1d\xc4`3\xd0?\xcf8oPX\xe0&@\x98\\\x0b;\xe1\xad\xd4?wIk\x1du*+@l\xddU\xa4h\xbc\xdb?]\x0c\xf9m\x9f\xe8-@@\xefW\xb5,\xd8\xe2?\xe2\x12\xd3&\x97\xbc\xb1@'
+tbs.

Python/OUTPUT/SSinit/ss_init_vars.pkl

@@ -0,0 +1,201 @@
+(dp1
+S'BQss'
+p2
+cnumpy.core.multiarray
+_reconstruct
+p3
+(cnumpy
+ndarray
+p4
+(I0
+tS'b'
+tRp5
+(I1
+(L2L
+tcnumpy
+dtype
+p6
+(S'f8'
+I0
+I1
+tRp7
+(I3
+S'<'
+NNNI-1
+I-1
+I0
+tbI00
+S'X\xaf\x99\x9a\x0eO\x9c?bP\xa3\x83\x04bj@'
+tbsS'chi_b'
+p8
+g3
+(g4
+(I0
+tS'b'
+tRp9
+(I1
+(L10L
+L2L
+tg7
+I00
+S'\x00\x00\x00\x00\x00\x00\xf0?\x00\x00\x00\x00\x00j\xf8@\x00\x00\x00\x00\x00\x00\xf0?\x00\x00\x00\x00\x00j\xf8@\x00\x00\x00\x00\x00\x00\xf0?\x00\x00\x00\x00\x00j\xf8@\x00\x00\x00\x00\x00\x00\xf0?\x00\x00\x00\x00\x00j\xf8@\x00\x00\x00\x00\x00\x00\xf0?\x00\x00\x00\x00\x00j\xf8@\x00\x00\x00\x00\x00\x00\xf0?\x00\x00\x00\x00\x00j\xf8@\x00\x00\x00\x00\x00\x00\xf0?\x00\x00\x00\x00\x00j\xf8@\x00\x00\x00\x00\x00\x00\xf0?\x00\x00\x00\x00\x00j\xf8@\x00\x00\x00\x00\x00\x00\xf0?\x00\x00\x00\x00\x00j\xf8@\x00\x00\x00\x00\x00\x00\xf0?\x00\x00\x00\x00\x00j\xf8@'
+tbsS'cssmat'
+p10
+g3
+(g4
+(I0
+tS'b'
+tRp11
+(I1
+(L10L
+L2L
+tg7
+I00
+S'\x99W\x0e\x8b\xf3;\xe7?\x85Py;\x93Nz@\xf2\\\xf5\x98\x90\xe4\xe7?\xe7\x1b\xb3\x91\xc6Jz@H\xd0b\t\x1a\xa7\xe7?\x00c>\xc0\x82$F\xc0X\xac%\x01\x15I\xe7?e[\xbe\xcd\x1ct\x94@x\xda\x9e\x88\xbe\xf7\xe6?\xdc\x1f\xcc\xe2zR\x11\xc0:\xb4p=a\xdf\xe6?g\xa5\xa2\xeb\x03>z@\xc5\xd528\xf8U\xe6?\x0b\xba`\x15\xd38z@\xd27E\xe9\xf7\x81\xe5?\xb4\x891i\n3z@@\x98\xbc`v~\xe4?)\x0b\x85"\xf3<z@\xda\xbf\xc3\xb6\x17k\xe3?\x07}\xf9~Y\x92\xb0\xc0'
+tbsS'Kss'
+p12
+cnumpy.core.multiarray
+scalar
+p13
+(g7
+S'\xe6@L}|?l@'
+tRp14
+sS'bssmat_s'
+p15
+g3
+(g4
+(I0
+tS'b'
+tRp16
+(I1
+(L10L
+L2L
+tg7
+I00
+S'\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00v\xcf\x18\xdc\x83\xf8\xc1?8\x1c\x10m\t\x1c\xf8?C\xb4\x05jw\xa5\xc6?\xb1\xe5\x02I\x82\xec\x07@\xf9\xcc__\xa8\xea\xc7?\x89\x81\xdcH1y|@W\xd4\xfd\xb3\x1d\xdc\xc7?\xb1\xdf\xce3\xae\x81{\xc0\x0e\x95j\x9e\xe4\xc6\xc6?n\xd1\x04\xfd\x12\xc7\x1e@\x1c\xa6\xb6\x1d[\x05\xcb?\x83\xe9\xaf\x15\xdb\xf3"@\xfeF\x1d\xc4`3\xd0?\xcf8oPX\xe0&@\x98\\\x0b;\xe1\xad\xd4?wIk\x1du*+@l\xddU\xa4h\xbc\xdb?]\x0c\xf9m\x9f\xe8-@'
+tbsS'chi_n'
+p17
+g3
+(g4
+(I0
+tS'b'
+tRp18
+(I1
+(L10L
+tg7
+I00
+S'\x00\x00\x00\x00\x00\x00\x14@\x00\x00\x00\x00\x00\x00\x18@\x00\x00\x00\x00\x00\x00\x1c@\x00\x00\x00\x00\x00\x00 @\x00\x00\x00\x00\x00\x00"@\x00\x00\x00\x00\x00\x00$@\x00\x00\x00\x00\x00\x00&@\x00\x00\x00\x00\x00\x00(@\x00\x00\x00\x00\x00\x00*@\x00\x00\x00\x00\x00\x00,@'
+tbsS'bssmat_splus1'
+p19
+g3
+(g4
+(I0
+tS'b'
+tRp20
+(I1
+(L10L
+L2L
+tg7
+I00
+S'v\xcf\x18\xdc\x83\xf8\xc1?8\x1c\x10m\t\x1c\xf8?C\xb4\x05jw\xa5\xc6?\xb1\xe5\x02I\x82\xec\x07@\xf9\xcc__\xa8\xea\xc7?\x89\x81\xdcH1y|@W\xd4\xfd\xb3\x1d\xdc\xc7?\xb1\xdf\xce3\xae\x81{\xc0\x0e\x95j\x9e\xe4\xc6\xc6?n\xd1\x04\xfd\x12\xc7\x1e@\x1c\xa6\xb6\x1d[\x05\xcb?\x83\xe9\xaf\x15\xdb\xf3"@\xfeF\x1d\xc4`3\xd0?\xcf8oPX\xe0&@\x98\\\x0b;\xe1\xad\xd4?wIk\x1du*+@l\xddU\xa4h\xbc\xdb?]\x0c\xf9m\x9f\xe8-@@\xefW\xb5,\xd8\xe2?\xe2\x12\xd3&\x97\xbc\xb1@'
+tbsS'euler_savings'
+p21
+g3
+(g4
+(I0
+tS'b'
+tRp22
+(I1
+(L10L
+L2L
+tg7
+I00
+S'\x00\x00\x00\x00\x00\xd8O=\x1e\xcf\xd6\x98L\xbcG\xc0\x00\x00\x00\x00\x00\xf7R\xbd[O3P0\x10\x1a\xc0\x00\x00\x00\x00\x00\xa0\x08=](b\xacAR\xec\xbe\x00\x00\x00\x00\x00@5=\x95G\x81\x8f@\xa2\x8a?\x00\x00\x00\x00\x00\xa0;\xbd\x97\xed3\xf1G\x1b\xe0\xbf\x00\x00\x00\x00\x00t1=\xa7\xd2\x07JFy\xd2\xbf\x00\x00\x00\x00\x00.0=B\xcaC\xcaV#\xc7\xbf\x00\x00\x00\x00\x00\xf0\x1f=\x99\x9cW\x83\xa9w\xbe\xbf\x00\x00\x00\x00\x00\x98\x00=\xe0\x88-\xfb\x91(\xb9\xbf\x00\x00\x00\x00\x00\x80J\xbdG\xf3m\xd6*a\xb0\xbe'
+tbsS'T_Hss'
+p23
+g13
+(g7
+S':l\x92E\xebJ\xc1?'
+tRp24
+sS'Yss'
+p25
+g13
+(g7
+S'\xbd\t\xcb\xfd\xdd\xa7\xe3?'
+tRp26
+sS'wss'
+p27
+g13
+(g7
+S'\xff\x85\x7fF\xb551@'
+tRp28
+sS'euler_labor_leisure'
+p29
+g3
+(g4
+(I0
+tS'b'
+tRp30
+(I1
+(L10L
+L2L
+tg7
+I00
+S'\x00\x00\x00\x00\x00\xa0B=\xbf%\xc8Q\x15\xd4\x8e\xbf\x00\x00\x00\x00\x00\x80\\\xbd\xa9\x15\x1f\xa7\x9djT\xbf\x00\x00\x00\x00\x00\xa0I=\xb3g\xb1\x8f\xe9\x9cX\xbf\x00\x00\x00\x00\x00\xa0B=;\xcd\ry\x1d\xca\xe9\xbe\x00\x00\x00\x00\x00\x803\xbd \xfa z\x08{\xf2?\x00\x00\x00\x00\x00@J=?b\td(\xe52\xbf\x00\x00\x00\x00\x00@3=\xc0\x10\xab\x14\xd2\xc98\xbf\x00\x00\x00\x00\x00\x00-\xbdKcf\xdb\xab\xb0\r\xbf\x00\x00\x00\x00\x00@H\xbd\xc5\x92`\xc6Z\x8f\xda\xbe\x00\x00\x00\x00\x00\x00O\xbdB\xc0\xf8\xad9\xc2\xd0\xbe'
+tbsS'Lss'
+p31
+g13
+(g7
+S'\xb1\xfc\xb0\x00| \x9a?'
+tRp32
+sS'nssmat'
+p33
+g3
+(g4
+(I0
+tS'b'
+tRp34
+(I1
+(L10L
+L2L
+tg7
+I00
+S"\xb2|\xd2{\x98\xb0\xce?IURIA\xff\x88?\x11\xbdj?\x1e\xf7\xca?b(\x8fm\x155k?\xf4\x8diR\xd9c\xc9?;\xff\xf9in\xa9i?\x1c\xa7v\xe7\x1dX\xc8?\xb8H\xf43\xd7\x024?\xe2\xd1_\x82\xb7t\xc7?\xda\xdd\xc7\x8a\x7f\xa4_?FfP\xce9l\xc6?\xee\x1e\xea\xda\x1d\xd1T?W\xdf\xc2Bp*\xc6?\x8a3\n|\xb5\xadV?\x1dI3B\x08t\xc6?\xb7|A3\xbaHA?\x7f\xeb\x17Y\x97,\xc7?\xe5]r\xc0iQ'?\x16+\xad~$0\xc8?\xc4Ev\xef\xc3\x99!?"
+tbsS'theta'
+p35
+g3
+(g4
+(I0
+tS'b'
+tRp36
+(I1
+(L2L
+tg7
+I00
+S'\xbd\x18\x8a>>\xccl?yN\xbe\x00e\xe1a?'
+tbsS'factor_ss'
+p37
+g13
+(g7
+S'\x12\r\xea0\x1cD\x1eA'
+tRp38
+sS'bssmat'
+p39
+g3
+(g4
+(I0
+tS'b'
+tRp40
+(I1
+(L9L
+L2L
+tg7
+I00
+S'v\xcf\x18\xdc\x83\xf8\xc1?8\x1c\x10m\t\x1c\xf8?C\xb4\x05jw\xa5\xc6?\xb1\xe5\x02I\x82\xec\x07@\xf9\xcc__\xa8\xea\xc7?\x89\x81\xdcH1y|@W\xd4\xfd\xb3\x1d\xdc\xc7?\xb1\xdf\xce3\xae\x81{\xc0\x0e\x95j\x9e\xe4\xc6\xc6?n\xd1\x04\xfd\x12\xc7\x1e@\x1c\xa6\xb6\x1d[\x05\xcb?\x83\xe9\xaf\x15\xdb\xf3"@\xfeF\x1d\xc4`3\xd0?\xcf8oPX\xe0&@\x98\\\x0b;\xe1\xad\xd4?wIk\x1du*+@l\xddU\xa4h\xbc\xdb?]\x0c\xf9m\x9f\xe8-@'
+tbsS'rss'
+p41
+g13
+(g7
+S'\x924o?*\x11\xa9\xbf'
+tRp42
+s.

Python/OUTPUT/Saved_moments/labor_data_moments.pkl

@@ -0,0 +1,29 @@
+(dp1
+S'labor_dist_data'
+p2
+cnumpy.core.multiarray
+_reconstruct
+p3
+(cnumpy
+ndarray
+p4
+(I0
+tS'b'
+tRp5
+(I1
+(L80L
+tcnumpy
+dtype
+p6
+(S'f8'
+I0
+I1
+tRp7
+(I3
+S'<'
+NNNI-1
+I-1
+I0
+tbI00
+S"\xdeM\xd6r\x9c\xa2\xd4?\xc1\x8aV^\x8cC\xd5?!\xbc\x89\xea9i\xd6?\x8cT\xb3\xaa\xfds\xd7?\x8e'`+\x15X\xd8?\x8a\x87\xc8\xf8\x89\xca\xd8?o\xad\xfd\x99\x8a \xd9?\xee(\xb3{\xbd\x1e\xd9?\xc5\x853\x8e\xfbF\xd9?\xa1}\xae\xc4vg\xd9?6\xd9\xcd\x1aJ\\\xd9?\x86B(\x19ek\xd9?E\xd4\xd7\xfd\xccn\xd9?\x0f=s\xcc\xb4h\xd9?\x9c\x1d\x9d\xaa\x1d\x93\xd9?\xe8\xfb\xe5\xef\xcd\xa4\xd9?\xaf\xe9a\xe7\x1a\xb1\xd9?\xd3\xed\xbc$\x92\xbd\xd9?\xcf\xd2\xcd\x01\xd3\xc9\xd9?\x8e\xddU\x8c\xdb\xba\xd9?y\x96<\x90m\xd3\xd9?\xd0\xa1\x8f\xa7\xae\xde\xd9?\x89\xca\xf7\xc3\xb9\xd9\xd9?\x9d,\x0c\x8b0\xf3\xd9?\xab\x98\x95\x03\xd8\xec\xd9?;\xa84\xc0)\xf2\xd9?8\x83Q\x7f\x0c&\xda?\xbb\xba\xe2\x06h\xff\xd9?\xf9\xb6\xf7\x7f\x0c\t\xda?\xd07\x11\x89P\xeb\xd9?\x99I\x8c{\xff\xbe\xd9?Q\xd8E\x13\x9f\xa1\xd9?|\x918\r\xbcg\xd9?Sx\x05TH9\xd9?\xd1\x9b`\xd0\xa5\xec\xd8?\xa4\xfe\xc2\x15\x00\xa1\xd8?t4\x8a\xfd\xac\x1f\xd8?\xe3\xf5\x8b\xa5\xc0\xe7\xd7?M\xd7\xd1\xcb\xd2\x83\xd7?R\xb9\x8d0\x82\xf3\xd6?#\xb9\xb6$\xd6V\xd6?o\x91\xfe3\x0b\xe8\xd4?\x03\xb3\xf5\xa0\x87\xa5\xd3?\xaf\x01=\xdc5\xe4\xd2?\x0f\x05\xf3\xc3D\xa7\xd1?\xb0\xad\x9eeCz\xd0?\nX\xe0\xf4P\xc4\xcf?\x1e\xae\xc6\xf5\xb0\xf6\xce?\xfe\x11p\x9aJ\xb5\xcd?X\xbf\xfd\xc83\t\xcd?\x1e\xed\x13IC%\xcd?B3\xab\xab\x91V\xcc?\xb89\x0ct\x07%\xcc?n\x13\xa5V\xebb\xcc?]\xfa\x00\xd7\xc1)\xcc?\x0ex7FR\xbc\xcb?\n\r\xf1\xee\x8c\xd9\xcb?\n\r\xf1\xee\x8c\xd9\xcb?\x1278i:\xaa\xcb?\x1aa\x7f\xe3\xe7z\xcb?!\x8b\xc6]\x95K\xcb?)\xb5\r\xd8B\x1c\xcb?1\xdfTR\xf0\xec\xca?8\t\x9c\xcc\x9d\xbd\xca?@3\xe3FK\x8e\xca?H]*\xc1\xf8^\xca?P\x87q;\xa6/\xca?X\xb1\xb8\xb5S\x00\xca?`\xdb\xff/\x01\xd1\xc9?g\x05G\xaa\xae\xa1\xc9?o/\x8e$\\r\xc9?wY\xd5\x9e\tC\xc9?\x7f\x83\x1c\x19\xb7\x13\xc9?\x86\xadc\x93d\xe4\xc8?\x8e\xd7\xaa\r\x12\xb5\xc8?\x96\x01\xf2\x87\xbf\x85\xc8?\x9e+9\x02mV\xc8?\xa6U\x80|\x1a'\xc8?\xae\x7f\xc7\xf6\xc7\xf7\xc7?\xb5\xa9\x0equ\xc8\xc7?"
+tbs.