关于置信水平,求区间的简便算法

 Ps~N(p,pq/n)

import scipy.stats as ss
import math
var1=0.25*0.75/50
std1=math.sqrt(var1)
intt=ss.norm.interval(0.99,0.25,std1)

　

 

----------------------------------------T分布------------------------------------------

 

t=ss.t.isf((1-0.95)/2,9)

 

 

 

 

 python t分布

 interval=stats.t.interval(a,b,mean,tsem)
 t分布的置信区间
 a:置信水平
 b:检验量的自由度
 mean:样本均值
 tsem:均值标准差 等同于math.sqrt(var/len(x))    np.var(x1,ddof=1)

 

from scipy import stats
import math
import numpy as np
x=[10.1,10,9.8,10.5,9.7,10.1,9.9,10.2,10.3,9.9]
x1=np.array(x)
mean=np.mean(x1)
tsem=stats.tsem(x1,ddof=1)#均值标准差
interval=stats.t.interval(0.95,len(x)-1,mean,tsem)#求总体均值置信区间

 

 

 

来源：《深入浅出统计学》