数学建模习题9.3

`import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns
from scipy import stats

file_path = '9.3.xlsx'
data = pd.read_excel(file_path, header=0)
labs = ['Lab1', 'Lab2', 'Lab3', 'Lab4', 'Lab5', 'Lab6', 'Lab7']
assert all(lab in data.columns for lab in labs), "Excel文件中的列名与实验室名称不匹配"

plt.figure(figsize=(12, 8))
sns.boxplot(data=data)
plt.xticks(rotation=45)
plt.title('Boxplot of Chlorpheniramine Content Measurements by Laboratory')
plt.xlabel('Laboratory')
plt.ylabel('Chlorpheniramine Content (mg)')
plt.show()

进行方差分析

f_stat, p_value = stats.f_oneway(data['Lab1'], data['Lab2'], data['Lab3'], data['Lab4'], data['Lab5'], data['Lab6'], data['Lab7'])
print(f"F-statistic: {f_stat}")
print(f"P-value: {p_value}") # 修正括号

alpha = 0.05
if p_value < alpha:
print("在显著水平α=0.05下，各实验室测量的扑尔敏有效含量的均值有显著差异。")
else:
print("在显著水平α=0.05下，各实验室测量的扑尔敏有效含量的均值没有显著差异。")
print("学号：05")`