使用anova 得到不同药物剂量下对安全性指标（a1,a2,a3,a4,a5）的影响代码R语言

在R语言中，ANOVA（Analysis of Variance，方差分析）通常用于比较两个或更多组之间数值变量的平均值是否有显著差异。如果你想研究不同药物剂量对安全性指标的影响，你可以使用`aov()`函数来进行单因素或多因素ANOVA。假设我们有一个数据框`df`，其中包含两列：`dose`（药物剂量）和`safety`（安全指标），下面是一个简单的例子：

# 假设你的数据结构如下：
# df <- data.frame(dose = c("low", "medium", "high", ..., "low", "medium", "high"), # 分为几个级别
#                   a1 = c(数值, 数值, 数值, ..., 数值, 数值, 数值), # 安全指标1
#                   a2 = c(数值, 数值, 数值, ..., 数值, 数值, 数值), # 同上
#                   a3 = c(数值, 数值, 数值, ..., 数值, 数值, 数值), # 同上
#                   a4 = c(数值, 数值, 数值, ..., 数值, 数值, 数值), # 同上
#                   a5 = c(数值, 数值, 数值, ..., 数值, 数值, 数值)) # 同上

# 对所有安全指标进行一次方差分析
model <- aov(safety ~ dose, data = df)

# 查看结果
summary(model)

这将返回一个包含F统计量、P值等信息的结果，可以帮你判断不同药物剂量是否显著影响了各个安全指标。