scau华南农业大学 高等统计学实验

scau华南农业大学 高等统计学实验四

7.10

import math
from scipy.stats import t
a = [6.5, 6.6, 6.7, 6.8, 7.1, 7.3, 7.4, 7.7, 7.7, 7.7]    # 方式1
b = [4.2, 5.4, 5.8, 6.2, 6.7, 7.7, 7.7, 8.5, 9.3, 10.0]    # 方式2

n = 10
A = 0.05
averageA = 0
averageB = 0
sumA = 0
sumB = 0

for i in range(10):
    sumA += a[i]

averageA = sumA / n
temp = 0

for i in range(10):
    temp += (a[i] - averageA)*(a[i] - averageA)

sa2 = temp / (n - 1)

print("第一种排队方式等待时间方差：", sa2)
print()


# 查表
xleft = t.ppf(1-A/2, df=n-1)

left = averageA-(xleft*sa2/math.sqrt(n))
right = averageA+(xleft*sa2/math.sqrt(n))

compareA = right - left

print("第一种排队方式等待时间标准差的置信区间为：", left, "~", right)
print("----------------------------------------------------------------")
#构建第二种排队方式等待时间标准差的95%的置信区间
for i in range(10):
    sumB += b[i]

averageB = sumB / n
temp = 0

for i in range(10):
    temp += (b[i]-averageB)*(b[i]-averageB)

sb2 = temp / (n-1)
print("第二种排队方式等待时间方差：",sb2)
print()


left = averageB-(xleft*sb2/math.sqrt(n))
right = averageB+(xleft*sb2/math.sqrt(n))
compareB = right - left

print("第二种排队方式等待时间标准差的置信区间为：",left,"~",right)
print("----------------------------------------------------------------")

#比较
if compareA > compareB:
    print("第二种排队方式更好，因为它的离散程度小于第一种排队方式。")
elif compareB > compareA:
    print("第一种排队方式更好，因为它的离散程度小于第二种排队方式。")

7.13

import math
from scipy.stats import t
a = [78,63,72,89,91,49,68,76,85,55]
b = [71,44,61,84,74,51,55,60,77,39]
n = 10
c = [a[i]-b[i] for i in range(10)]


sum = 0
for i in range(n):
    sum += c[i]
d = sum / n


temp = 0
for i in range(n):
    temp += (c[i]-d)*(c[i]-d)
sd = math.sqrt(temp/(n-1))


check = t.ppf(0.95,9)
Float = check*sd/math.sqrt(n)
print(" 两种方法平均自信心得分之差的标准差：", sd)
print(" 两种方法平均自信心得分之差的95%的置信区间为:", d - Float, "~", d + Float)