So I\'m using pymc3 to calculate the Geweke and Gelmen-Rubin diagnostic for an MCMC chain like:
chain = sampler(data,samples=5000,mu_init=9,sigma=sigma,propos