如何在R中进行方差分析？

Question

我有一个有三列的dataframe : data$input、data$output和data$columns。输入和输出都是连续数字，类别是离散字符。我知道ANCOVA是一种在控制类别效应的同时，分析输入对产出的影响的方法。然而，我没有找到确切的命令来做这个在线。

下面是一个关于数据外观的小例子：

   input output category

1    0.4   0.55        A

2    0.5   0.66        A

3    0.6   0.57        A

4    0.3   0.23        B

5    0.4   0.53        B

6    0.7   0.75        B

7    1.1   1.31        C

8    0.9   1.01        C

9    0.8   0.58        C

10   0.5   0.34        C

我应该输入哪些命令来执行ANCOVA？非常感谢!

Answer 1

您可以使用函数lm，它运行线性回归(最后，ANOVA和ANCOVA只是线性模型的受限版本)。

mod <- lm(输出~输入+类别，data=data)

可以使用summary函数查看输出。

摘要(国防部)

如果您确实需要输出'ANOVA样式‘，那么您可以在创建的模型上应用函数anova (我个人认为lm输出信息更丰富，但可能有特定的限制)。

anova(国防部)