双向方差分析的Levene检验

import pandas as pd
from scipy.stats import levene

# 示例数据
data = {
    'FactorA': ['A1', 'A1', 'A2', 'A2', 'A1', 'A1', 'A2', 'A2'],
    'FactorB': ['B1', 'B2', 'B1', 'B2', 'B1', 'B2', 'B1', 'B2'],
    'Value': [20, 21, 19, 22, 20, 19, 21, 23]
}

df = pd.DataFrame(data)

# 分组数据
group1 = df[(df['FactorA'] == 'A1') & (df['FactorB'] == 'B1')]['Value']
group2 = df[(df['FactorA'] == 'A1') & (df['FactorB'] == 'B2')]['Value']
group3 = df[(df['FactorA'] == 'A2') & (df['FactorB'] == 'B1')]['Value']
group4 = df[(df['FactorA'] == 'A2') & (df['FactorB'] == 'B2')]['Value']

# 执行Levene检验
stat, p_value = levene(group1, group2, group3, group4)
print('Levene统计量:', stat)
print('P值:', p_value)

# 根据P值判断方差是否齐性
if p_value < 0.05:
    print("拒绝方差齐性的零假设，方差不齐")
else:
    print("接受方差齐性的零假设，方差齐")