《概率论与数理统计》小鸡增肥例子
巴特利特检验
import scipy.stats as sp
import numpy as np
from matplotlib import pyplot as plt
#chicks
a1=np.array([73,9,60,1,2,12,9,28])+1000
a2=np.array([107,92,-10,109,90,74,122,1])+1000
a3=np.array([93,29,80,21,22,32,29,48])+1000
sp.bartlett(a1, a2, a3)
out
BartlettResult(statistic=3.386318395372458, pvalue=0.18393751017215176)