如何在R coefplot::multiplot中按系数而不是模型进行着色

multiplot <- function(..., coef_colors = NULL) {
  plots <- list(...)
  num_plots <- length(plots)
  
  if (is.null(coef_colors)) {
    coef_colors <- rainbow(num_plots)
  }
  
  par(mfrow = c(ceiling(num_plots/2), 2))
  
  for (i in 1:num_plots) {
    plot <- plots[[i]]
    coef_color <- coef_colors[i]
    
    coefplot::coefplot(plot, col = coef_color)
  }
}

# 假设有两个模型 model1 和 model2
model1 <- lm(Sepal.Length ~ Sepal.Width + Petal.Length, data = iris)
model2 <- lm(Sepal.Length ~ Sepal.Width + Petal.Width, data = iris)

# 定义颜色向量
coef_colors <- c("red", "blue")

# 使用修改后的 multiplot 函数绘制系数图
multiplot(model1, model2, coef_colors = coef_colors)