「医学统计学之一」方差齐性检验：python-scipy实现

△ 内容：

1 概念和举例子。

2 python相关知识，scipy的统计学包stats使用。

△ 方差齐性检验：

1 对两个或者以上的研究总体的总体平均数的差异进行显著性检验。

2通俗地说，就是两组或多组之间如果进行其他数据的比较，那么前提是这两组或者多组的方差是齐性的，什么意思，就是它们之间没有统计学上的的差异，这样子后续的研究才有意义。

△ 举例子1：口服阿奇霉素冲剂治疗轻型支气管肺炎的疗效观察

1 随机分为实验组和对照组，收集资料，分布如下。

2 那么就要分析这两组是不是具有统计学上的差异，最好是统计学上无明显差异，这样后续的研究才有意义。

3 代码

from scipy.stats import levene# 数据，一行一行数据# 实验组x=[35,41,5.1,121,31]# 对照组y=[33,46,5.3,118,30]# Levene检验，适用于多组方差的比较，且对正态性没要求out= levene(x, y)print(out)# 结果#LeveneResult(statistic=0.0006253614373741172, pvalue=0.9806617265486998)# P约等于0.98，＞0.05，无差异，具有方差齐性

△ scipy包的stats的levene：

1 还可以用于多组的方差齐性的比较。

2 举例2：三组xyz的方差齐性的比较，代码如下：

# 同时熟悉python和levene的2种的使用方法from scipy.stats import levene# 模拟数据，x、y、z代表不同组的样本。x=[11,22,33,44]y=[12,23,34,45]z=[10,20,30,40]# 方法一stat, p = levene(x, y, z)print(stat, p)# 结果#0.03508771929824561 0.9656521217710577# 方法二out= levene(x, y, z)print(out)# 结果#LeveneResult(statistic=0.03508771929824561, pvalue=0.9656521217710577)# 以上结果均提示P≈0.966＞0.05，统计学上无差异。

△ 赋诗一首（巩固一下原创性）：

医生科研要努力，

统计数据需分析，

小编带你巧学习，

服务病人才专一。

===========================

未完待续。。。。。。

===========================