Simulating counterfactuals