Monte Carlo, Markov Chains, Metropolis Hastings.
np.std(theta_accept) print('std of theta, ', sd) df = 10 se = sd / np.sqrt(df) print ('standard error of theta, ', se) t = np.abs((mean-0.5) / se) print('t = ', t) p = (1-scipy.stats.t.cdf(t, 10))*2 print('p value from t test is, ', p) # output mean of theta, 0.4967416562316055 std of theta,