# -*- coding: utf-8 -*-
"""
Created on Wed Jul  7 11:46:24 2021

@author: Sim
"""

import scipy.stats as st
import numpy as np
import matplotlib.pyplot as plt

def t_ci_mean(x, alpha):
    ave = np.mean(x)
    df = len(x) - 1
    lbd = ave - st.t.ppf(1-alpha/2, df)*np.std(x, ddof=1)/(len(x)**.5)
    ubd = ave + st.t.ppf(1-alpha/2, df)*np.std(x, ddof=1)/(len(x)**.5)
    return([ave, lbd, ubd])

plt.plot( [-1,100],[-3,3], visible=False)
plt.plot((-1,100), (0,0), 'k-')
for i in range(100):
  xx = st.norm.rvs(size=5)
  ave, lbd, ubd = t_ci_mean(xx, 0.05)
  print(lbd, ubd)
  if lbd*ubd > 0: 
    plt.plot([i,i], [lbd, ubd],  'r-', linewidth=1)
  else:
    plt.plot([i,i], [lbd, ubd],  'b-', linewidth=1)