R
!!!rstatix
https://www.wenyanet.com/opensource/ja/6005734135624a73b01c238b.html
https://github.com/kassambara/rstatix


{{pre
install.packages("rstatix")
library(rstatix)
}}

!get_summary_stats(omega.dat2)

!mshapiro_test(omega.dat2$MDD)

!pairwise_wilcox_test
{{pre
pairwise_wilcox_test(omega.dat2, MHD~year)

	.y.	group1 group2	n1	n2	statistic	p	p.adj	p.adj.signif

1	MHD	2	3	699	589	144470.5	2.36e-20	4.72e-20	****
2	MHD	2	ns	699	462	70514.5	1.28e-59	3.84e-59	****
3	MHD	3	ns	589	462	91477.5	6.77e-20	6.77e-20	****

}}

!wilcox_effsize
{{pre
wilcox_effsize(omega.dat2, MHD~year)

	.y.	group1	group2	effsize	n1	n2	magnitude

1	MHD	2	3	0.2575914	699	589	small
2	MHD	2	ns	0.4779203	699	462	moderate
3	MHD	3	ns	0.2816649	589	462	small

}}
!!ggpubr
{{pre
library(ggpubr)
ggboxplot(omega.dat2, x = "year", y = "MHD")+
  stat_pvalue_manual(
    pairwise.test.result, label = "p.adj", 
    y.position = c(6, 7, 8)
}}
{{ref_image MHDyear.jpg}}