我正在嘗試使用具有各自置信區間的負二項式隨機效應模型繪制估計均值，其中資料和擬合定義如下：
d.AD <- data.frame(treatment = gl(12,12),
treatment2 = gl(4,1,36),
counts = rpois(144, 4.03),
treatment3 = gl(3,4,36),
ID = gl(12,12))
d.AD$treatment2 = as.integer(d.AD$treatment2)
d.AD$treatment3 = as.factor(d.AD$treatment3)
d.AD$ID = as.factor(d.AD$ID)
library(ggplot2)
library(lme4)
library(emmeans)
mode1A <- glmer.nb(counts~treatment3 poly(treatment2, 2) (1|ID),
data=d.AD, verbose=F)
使用所使用的置信區間進行平均的函式是emmeans包的emmeans。但是，看到圖形的形狀是倒置的，也就是說，我希望 x 軸成為我的 y 軸，反之亦然。
emm1 = emmeans(mode1A,~"treatment3", type="response")
emm1df = as.data.frame(emm1)
plot(emm1)
geom_point(size=2.2)
labs(x = 'Response', y = 'Treatment')
scale_shape_manual(NULL, values = 18)
theme(legend.title = element_text(size = 25,color = "black"),
legend.text = element_text(size = 25,color = "black"),
axis.title = element_text(size = 25,color = "black"),
axis.text.x = element_text(color = "black", hjust=1),
axis.text.y = element_text(color = "black", hjust=1),
axis.text = element_text(size = 25,color = "black"),
plot.title = element_blank(),
strip.text.x = element_text(size = 15,color = "black"))
我想要一些接近這個例子的東西。
uj5u.com熱心網友回復：
您可以執行以下操作：
emm1df |>
ggplot(aes(x = treatment3,
y = response,
label = round(response,2)
)
)
geom_point(size = 2)
geom_errorbar(aes(ymin = asymp.LCL,
ymax = asymp.UCL),
width = 0
)
geom_label(nudge_y = (emm1df$asymp.UCL-emm1df$response)*1.2)
ylim(c(3,5))
labs(x = "Treatment",
y = "Response")
uj5u.com熱心網友回復：
嘗試
plot(emm1, horizontal = FALSE)
見`? plot.emmGrid
轉載請註明出處，本文鏈接：https://www.uj5u.com/gongcheng/487256.html
上一篇：如何繪制隨時間變化的比較直方圖？