41、协方差分析（ANCOVA）及稳健方法详解

协方差分析（ANCOVA）及稳健方法详解

1. 更新ANCOVA模型以包含交互项

在进行ANCOVA分析时，有时需要更新原始模型以包含交互项。可以使用 update() 函数来实现这一目的。假设我们有一个原始的ANCOVA模型 viagraModel ，它已经包含了 partnerLibido 和 dose 变量。若要添加交互项，可以使用以下代码：

hoRS<-update(viagraModel, .~. + partnerLibido:dose)

这里， .~. 表示保持与之前相同的结果变量和预测变量， + partnerLibido:dose 表示添加交互项。执行上述命令创建 hoRS 对象后，可使用 Anova() 函数获取III型平方和：

Anova(hoRS, type="III")

以下是输出结果示例：
| | Sum Sq | Df | F value | Pr(>F) |
| — | — | — | — | — |
| (Intercept) | 53.542 | 1 | 21.9207 | 9.323e - 05 *** |
| partnerLibido | 17.182 | 1 | 7.034