R レポートなどに載せるグラフに検定結果（p値）のアスタリスク付け方

Question

前提・実現したいこと

R初心者でデータの解析をしています。下の画像のように箱ひげ図の検定結果のバーとアスタリスクを表示させたいです。

https://kimbio.info/r_t-test_asterisk_ggplot2/
上のURLのサイトを参考にやっているのですが、２群の場合はうまくできるのですが、３群の場合どのようにやればいいのかわかりません。ご教示いただけると幸いです。

該当のソースコード

library(ggplot2)
d <- read.csv(file='sample.csv')

plt <- ggplot(d, aes(y=height,x=classes)) +
  geom_boxplot(width=0.4)
plt

library("ggsignif")
plt2 <- plt +   geom_signif(comparisons = list(c("classA", "classB")),
              test = "t.test",
              na.rm = FALSE,
              map_signif_level = TRUE) + theme_minimal()
plt2

試したこと

geom_signif()のcomparisonsにlist(c("classA","classB"),c("classA","classB"))を渡したのですが２つのバーが同じ高さになってしまって重なってしまいました。

Answer 1

これが参考になるのではないでしょうか。

Rのggplot2でグラフを描いたときに有意差の記号、*とかそういうのをグラフ上に書こうとしたときの方法について - 備忘ログ

Comments

[Hiro10]

ありがとうございます！まさにこれを探してました！