Plotting counterfactuals is almost working, except that random shocks
differ between original simulation and new simulation