トップ 差分 一覧 Farm ソース 検索 ヘルプ PDF RSS ログイン

rstatix

*disclaimer
231312

R

rstatix

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


install.packages("rstatix")
library(rstatix)

get_summary_stats(omega.dat2)


mshapiro_test(omega.dat2$MDD)


pairwise_wilcox_test

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

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

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