如何根据负二项模型制作调整均值的箱形图

Question

对于我使用glm.nb()中的R函数制作的模型，您如何在两个协议之间制作经过调整而非原始的方法（使用CI）的方框图。

例如（在这里容忍重量不算例子）：

library(MASS)
summary(glm.nb(formula = weight ~ 1 + Diet + Time, subset = levels(Chick) == "1", data = ChickWeight))

因此，对于Chick 1，默认情况下，饮食1与饮食2之间存在平均比较;饮食1 vs.饮食3;饮食1与饮食4相比，时间是另一个预测因素，即随着小鸡变老而调整比较方法。

如果我想从R获得那些调整后的方法，而不是使用tapply或其他原始方法：

with(data = subset(x = ChickWeight, subset = levels(Chick) == "1"),
     expr = tapply(X = weight, INDEX = Diet, FUN = mean))

我想最终得到类似的东西：

ggplot(data = subset(x = ChickWeight, subset = levels(Chick) == "1"),
       mapping = aes(x = Diet, y = weight)) + geom_boxplot()