如何从MuMIn model.avg()摘要中绘图

library(MuMIn)
library(ggplot2)

# 创建模型列表
model1 <- lm(y ~ x1, data = data)
model2 <- lm(y ~ x1 + x2, data = data)
model3 <- lm(y ~ x1 + x2 + x3, data = data)
models <- list(model1, model2, model3)

# 计算模型平均效应
avg <- model.avg(models)

# 提取需要绘制的数据
estimates <- avg$coefficients[, "Estimate"]
se <- avg$coefficients[, "SE"]

# 创建绘图
data <- data.frame(Model = c("Model 1", "Model 2", "Model 3"),
                   Estimate = estimates,
                   SE = se)
p <- ggplot(data, aes(x = Model, y = Estimate)) +
  geom_bar(stat = "identity", fill = "blue", width = 0.5) +
  geom_errorbar(aes(ymin = Estimate - SE, ymax = Estimate + SE),
                width = 0.2, color = "black") +
  labs(x = "Model", y = "Estimate") +
  ggtitle("Model Average Effects") +
  theme_minimal()

# 显示绘图
print(p)