置信系数1-α=0.95
ltx=[221,244,243,288,233,220,210,258,245]
lty=[268,213,188,189,217,207]
def compute(ltx):
le=len(ltx)
sum=0
S=0
for i in range(0,le):
sum+=ltx[i]
avg=sum/le
print('avg:',avg)
for i in range(0,le):
S+=(ltx[i]-avg)*(ltx[i]-avg)
s=S/(le-1)
print('S²:',s)
return s,avg
S1,x=compute(ltx)
S2,y=compute(lty)
n1=len(ltx)
n2=len(lty)
v=n1+n2-2
Sw=(((n1-1)*S1+(n2-1)*S2)/v)**0.5
print("自由度v:",v)
print("Sw:",Sw)
#查表t₀.₀₂₅[v]=2.1604
N=(1/n1+1/n2)**0.5
print("N:",N)
se1=x-y-2.1604*Sw*N
se2=x-y+2.1604*Sw*N
print('μ₁-μ₂的置信区间为',(se1,se2))點擊查看更多內容
為 TA 點贊
評論
評論
共同學習，寫下你的評論
評論加載中...
作者其他優質文章
正在加載中
感謝您的支持，我會繼續努力的～
掃碼打賞，你說多少就多少
贊賞金額會直接到老師賬戶
支付方式
打開微信掃一掃，即可進行掃碼打賞哦