add code to simulate from conditional rhos in SIW model