是否有 R function 來執行 Ancova 事后測試以檢查回歸斜率的同質性

Question

在調整了一些線性模型后，我首先要測試回歸斜率的同質性。 第二步，這是我的疑問，我想使用事后測試來兩個兩個比較斜率。

這是一個從https://www.datanovia.com/en/lessons/ancova-in-r/修改的示例

獲取數據

data("anxiety", package = "datarium")
anxiety <- anxiety[,c("id","group","t1","t3")]
names(anxiety)[c(3,4)] <- c("pretest","posttest")

plot 回歸線

ggscatter(anxiety,x="pretest",y="posttest",color="group",add="reg.line")+
    stat_regline_equation(aes(label=paste(..eq.label.., ..rr.label.., sep = "~~~~"),color = group))

檢查回歸斜率的同質性

anova_test(anxiety,posttest~group*pretest)

在這里，我們可以看到統計上不顯着的 p 值 4.15e-01

事后測試 emmeans_test 執行成對比較以確定哪些組不同。 不過，我想采用多重比較程序來確定哪些B （斜率）與其他哪些不同。

這個有 function 嗎？ 提前致謝。

Answer 1

在閱讀和搜索更多之后，我准備了一個我試圖做的分析的例子。 我希望它是有用的。 一個重要的來源是https://cran.r-project.org/web/packages/emmeans/vignettes/interactions.html

## packages
library(ggpubr)
library(rstatix)
library(emmeans)
library(data.table)

## prepare the example data
rm(list = ls())
set.seed(321)

a1 <- 0
b1 <- 1
a2 <- 1
b2 <- 1.7

x <- c(1:10)
y1 <- (a1+b1*x)+rnorm(10,0,.6)
y2 <- (a1+b1*x)+rnorm(10,0,.6)
y3 <- (a2+b1*x)+rnorm(10,0,.6)
y4 <- (a2+b2*x)+rnorm(10,0,.6)

dat <- data.frame(x=rep(x,4),y=c(y1,y2,y3,y4),group=rep(c("A","B","C","D"),each=10))

## regression and coefficients
lm.dat <- lm(y~x*group,dat)
summary(lm.dat)

## coeficients, confidence intervals and R2 by group
as.data.table(dat)[,as.list(coef(lm(y~x))),by=group]
as.data.table(dat)[,as.list(confint(lm(y~x))),by=group]
as.data.table(dat)[,list(r2=summary(lm(y~x))$r.squared),by=group]

## plots
emmip(lm.dat,group~x,cov.reduce=range)
ggscatter(dat,x="x",y="y",color="group",add="reg.line")+
    stat_regline_equation(aes(label=paste(..eq.label.., ..rr.label.., sep = "~~~~"),color = group))

## anova
anova(lm.dat)
anova_test(dat,y~x*group)

## interactions with covariates
## slopes for each group and pairwise comparisons
emtrends(lm.dat,pairwise~group,var="x")