用python做方差分析_python做单因素方差分析

方差分析的主要功能就是验证两组样本，或者两组以上的样本均值是否有显著性差异，即均值是否一样。

这里有两个大点需要注意：①方差分析的原假设是：样本不存在显著性差异(即，均值完全相等)；②两样本数据无交互作用(即，样本数据独立)这一点在双因素方差分析中判断两因素是否独立时用。

原理：

方差分析的原理就一个方程：SST=SS组间+SSR组内  (全部平方和=组间平方和+组内平方和)

说明：方差分析本质上对总变异的解释。

组间平方和=每一组的均值减去样本均值

组内平方和=个体减去每组平方和

方差分析看的最终结果看的统计量是：F统计量、R2。

其中：g为组别个数，n为每个组内数据长度。

python实现：

from scipy import stats

from statsmodels.formula.api import ols

from statsmodels.stats.anova import anova_lm

from statsmodels.stats.multicomp import pairwise_tukeyhsd

import warnings

warnings.filterwarnings("ignore")

import itertools

df2=pd.DataFrame()

df2['group']=list(itertools.repeat(-1.,9))+ list(ite