Уровень значимости ggplot2 со звездочкой с скорректированными значениями p

Question

Я абсолютный новичок в R, и мне было интересно, могу ли я получить некоторую помощь в визуализации данных с помощью пакета ggplot2.

Я пытаюсь пометить уровень значимости с помощью скорректированного значения p (с поправкой Бонферрони) вместо обычного p-значения со звездочкой.

Коды, которые я пробовал: ...

stat_compare_means(aes(label=..p.adj.., group=sample_type), label = "p.signif", method = "t.test", hide.ns = TRUE)

stat_compare_means(aes(group=sample_type), p.adjust.method = "bonferroni", label = "p.signif", method = "t.test", hide.ns = TRUE)

Answer 1

См. Пакет EnvStats для этого. Пример , предоставленный авторами (стр. 1180):

library(ggplot2)
library(EnvStats)
p <- ggplot(mtcars, aes(x = factor(cyl), y = mpg, color =     
factor(cyl))) + theme(legend.position = "none")

dev.new()
p + geom_point() +
stat_n_text() + stat_mean_sd_text() +
stat_test_text() +
labs(x = "Number of Cylinders", y = "Miles per Gallon")