in the code: add: import statistics modify : sns.pointplot(x=x, y=sample_means, **linestyle='None'**, ax=ax) modify: ci = t.interval(**confidence=.95**, df=len(simdata)-1, loc=np.mean(simdata), scale=sem(simdata))