具有系数均值和置信区间的条形图

Question

我想创建一个类似于所讨论的here的情节。

但是，不同之处在于我有28个以上的回归，因此，我想创建一个类似的图形，该图形显示每个变量及其置信区间的平均系数。需要说明的是：我不希望每个变量有28个系数。我想要每个变量1个系数（平均值）及其置信区间。

我总共有28个因变量和21个解释变量。

我该怎么做？

MRE

ami_data <- read.table("http://static.lib.virginia.edu/statlab/materials/data/ami_data.DAT")
names(ami_data) <- c("TOT","AMI","GEN","AMT","PR","DIAP","QRS")

mlm1 <- lm(cbind(TOT, AMI) ~ GEN + AMT + PR + DIAP + QRS, data = ami_data)

在上面的回归中，有两个因变量（TOT和AMI；在我的情况下，我有28个）。