python 方差分析 样本量太大_python单因素方差分析实例

试验设计

探究施肥是否会对促进植株生成(植株生长以树高作为指标来衡量)。试验为： - 对照组：清水 - 实验组： 某肥料四个浓度梯度，分别是A,B,C,D，施肥一段时间之后测量树高(要控制其他变量保持一致，比如施肥之前的树高要基本保持一致，生长势基本保持一致等等)

做方差分析的时候数据需要满足正态分布；方差齐性等。正常拿到数据后需要对数据是否符合正态分布和组间方差是否一致做检验。如何来做以上两个检验今天先忽略掉，在默认拿到的数据符合条件后直接在做单因素方差分析。

模拟生成数据

使用numpy模块模拟生成5组，每组100个正态分布数据 正态分布函数参数依次是均值，标准差、数据的个数

import numpy as np

df = {'ctl':list(np.random.normal(10,5,100)),

'treat1':list(np.random.normal(15,5,100)),\

'treat2':list(np.random.normal(20,5,100)),\

'treat3':list(np.random.normal(30,5,100)),\

'treat4':list(np.random.normal(31,5,100))}

#组合成数据框

import pandas as pd

df = pd.DataFrame(df)

df.head()

ctl treat1 treat2 treat3 treat4

0 9.614605 15.719777 17.068697 23.842793 32.206690

1 7.617131