如何用R对两个回归的斜率进行Welch t检验？

# 加载必要的库
library(tidyverse)

# 示例数据
set.seed(123)
data <- data.frame(
  x = rnorm(100),
  group = sample(c("A", "B"), 100, replace = TRUE),
  y = rnorm(100)
)

# 拟合两个回归模型
model_a <- lm(y ~ x, data = data[data$group == "A", ])
model_b <- lm(y ~ x, data = data[data$group == "B", ])

# 提取斜率
slope_a <- coef(model_a)[2]
slope_b <- coef(model_b)[2]

# 创建一个数据框来存储斜率
slopes <- data.frame(
  group = c("A", "B"),
  slope = c(slope_a, slope_b)
)

# 执行Welch t检验
t_test_result <- t.test(slope ~ group, data = slopes, var.equal = FALSE)

# 打印t检验结果
print(t_test_result)