在进行ANCOVA时，确保至少有两个组别以及一个与预测变量相关的协变量变量。这是为了比较不同组别之间的差异是否受到协变量的影响。因此，需要进行对比。

示例代码：

创建一个数据集

library(dplyr) library(ggplot2) set.seed(123) df <- data.frame(group = sample(c("A", "B"), 50, replace = TRUE), x = rnorm(50, 10, 2), y = rnorm(50, 20, 3))

运行ANCOVA

model <- lm(y ~ group + x, data = df) summary(model)

对组A和组B进行对比

contrast <- c(-1, 1) contrast.matrix <- makeContrasts(contrast = contrast, levels = levels(df$group)) fit <- contrasts.fit(model, contrast.matrix) fit_summary <- summary(fit) fit_summary

可视化效果

ggcontrast(fit, contrast = "contrast", axis.title = "Group Contrast")