# -*- coding: utf-8 -*-
"""
Created on Mon Jan  3 13:03:41 2022

@author: Sim
"""
import scipy.stats as st, numpy as np

drug1 = np.array([7, 8, 7, 4, 7, 7,  9, 9, 7, 4])
drug2 = np.array([6, 7, 8, 3, 5, 6, 10, 9, 5, 3])

diff = drug1-drug2
lbd, ubd = st.t.interval(0.95, len(diff)-1, loc=np.mean(diff), 
                         scale=st.sem(diff))
print(lbd, ubd)

tval, pval = st.ttest_rel(drug1, drug2, alternative='two-sided')
print(tval, pval)