diff --git a/HARK/ConsumptionSaving/ConsIndShockModel.py b/HARK/ConsumptionSaving/ConsIndShockModel.py index 887fc4a3d..43c76ff37 100644 --- a/HARK/ConsumptionSaving/ConsIndShockModel.py +++ b/HARK/ConsumptionSaving/ConsIndShockModel.py @@ -1732,6 +1732,13 @@ def sim_death(self): DiePrb = DiePrb_by_t_cycle[ self.t_cycle - 1 ] # Time has already advanced, so look back one + + # In finite-horizon problems the previous line gives newborns the + # survival probability of the last non-terminal period. This is okay, + # however, since they will be instantly replaced by new newborns if + # they die. + # See: https://github.com/econ-ark/HARK/pull/981 + DeathShks = Uniform(seed=self.RNG.randint(0, 2 ** 31 - 1)).draw( N=self.AgentCount )