Monte Carlo, Markov Chains, Metropolis Hastings.
theta_accept = theta_accept[1:] prior_accept = prior_accept[1:] likelihood_accept = likelihood_accept[1:] Here is the first 20 θ that is accepted: theta_accept[:20] Out[6]: array([0.4182804 , 0.27898128, 0.87778638, 0.56595161, 0.70156645, 0.49273126, 0.47311259, 0.75905831, 0.40404581, 0.4929927 ,