R语言双因素方差分析

R语言双因素方差分析

双因素方差分析（Two-way ANOVA）是一种用于探索两个自变量对因变量的影响的统计方法。在R语言中，我们可以使用多种包来执行双因素方差分析，如car、stats和afex等。本文将介绍如何使用这些包来进行双因素方差分析，并提供相应的源代码。

首先，我们需要安装并加载所需的R包。在执行下面的代码之前，请确保已经安装了这些包：

install.packages("car")
install.packages("stats")
install.packages("afex")

library(car)
library(stats)
library(afex)

接下来，我们准备一个样本数据集来进行双因素方差分析。假设我们正在研究两种不同的肥料对植物生长的影响，同时考虑了两个因素：肥料类型（fertilizer）和施肥时间（time）。我们随机地选择了16个样本，每个样本都测量了植物的生长高度（height）。下面是一个示例数据集：

# 创建数据框
data <- data.frame(
  fertilizer = rep(c("A", "B"), each = 8),
  time = factor(rep(c("T1", "T2"), each = 4, times = 2)),
  height = c(10, 12, 14, 16, 9, 11, 13, 15, 8, 10, 12, 14, 7, 9, 11, 13)
)

# 查看数据框
print(data)

在执