> dd Group y 1 A 4.439524 2 A 4.769823 3 A 6.558708 4 A 5.070508 5 A 5.129288 6 A 6.715065 7 A 5.460916 8 A 3.734939 9 A 4.313147 10 A 4.554338 11 B 16.224082 12 B 15.359814 13 B 15.400771 14 B 15.110683 15 B 14.444159 16 B 16.786913 17 B 15.497850 18 B 13.033383 19 B 15.701356 20 B 14.527209
> dd Group y 1 A 4.439524 2 A 4.769823 3 A 6.558708 4 A 5.070508 5 A 5.129288 6 A 6.715065 7 A 5.460916 8 A 3.734939 9 A 4.313147 10 A 4.554338 11 B 16.224082 12 B 15.359814 13 B 15.400771 14 B 15.110683 15 B 14.444159 16 B 16.786913 17 B 15.497850 18 B 13.033383 19 B 15.701356 20 B 14.527209 21 C 13.932176 22 C 14.782025 23 C 13.973996 24 C 14.271109 25 C 14.374961 26 C 13.313307 27 C 15.837787 28 C 15.153373 29 C 13.861863 30 C 16.253815
2.1 箱线图+散点图
p = ggboxplot(dd,x = "Group",y = "y",color = "Group",add = "jitter") p
> dd Group1 Group2 y 1 A X 4.439524 2 A X 4.769823 3 A X 6.558708 4 A X 5.070508 5 A X 5.129288 6 A X 6.715065 7 A X 5.460916 8 A X 3.734939 9 A X 4.313147 10 A X 4.554338 11 B X 9.224082 12 B X 8.359814 13 B X 8.400771 14 B X 8.110683 15 B X 7.444159 16 B X 9.786913 17 B X 8.497850 18 B X 6.033383 19 B X 8.701356 20 B X 7.527209 21 C X 5.932176 22 C X 6.782025 23 C X 5.973996 24 C X 6.271109 25 C X 6.374961 26 C X 5.313307 27 C X 7.837787 28 C X 7.153373 29 C X 5.861863 30 C X 8.253815 31 A Y 15.426464 32 A Y 14.704929 33 A Y 15.895126 34 A Y 15.878133 35 A Y 15.821581 36 A Y 15.688640 37 A Y 15.553918 38 A Y 14.938088 39 A Y 14.694037 40 A Y 14.619529 41 B Y 17.305293 42 B Y 17.792083 43 B Y 16.734604 44 B Y 20.168956 45 B Y 19.207962 46 B Y 16.876891 47 B Y 17.597115 48 B Y 17.533345 49 B Y 18.779965 50 B Y 17.916631 51 C Y 17.253319 52 C Y 16.971453 53 C Y 16.957130 54 C Y 18.368602 55 C Y 16.774229 56 C Y 18.516471 57 C Y 15.451247 58 C Y 17.584614 59 C Y 17.123854 60 C Y 17.215942
3.1 绘制分组箱线图
p = ggboxplot(dd,x = "Group1",y="y",color = "Group2", add = "jitter") p
3.2 增加P值
p + stat_compare_means(aes(group = Group2),method = "t.test")
3.3 修改为显著性结果
p + stat_compare_means(aes(group = Group2),method = "t.test",label = "p.signif")
3.4 将分组数据分开绘制
p = ggboxplot(dd,x = "Group2",y="y",color = "Group1", add = "jitter",facet.by = "Group1") p
3.5 分组显示统计检验
p + stat_compare_means(method = "t.test")
3.6 分组显示显著性结果
p + stat_compare_means(method = "t.test",label = "p.signif",label.y = 17)
4. 单因素直方图绘制
直方图+标准误,之前用ggplot2需要很长的代码,这里有更好的方案。
4.1 直方图+标准误
p = ggbarplot(dd,x = "Group1",y = "y",add = "mean_se",color = "Group1") p
4.2 直方图+标准误+显著性
p + stat_compare_means(method = "anova",,label.y = 15)+ stat_compare_means(comparisons = my_comparisons)
5. 单因素折线图绘制
5.1 折线图+标准误
p = ggline(dd,x = "Group1",y = "y",add = "mean_se") p
5.2 折线图+标准误+显著性
p + stat_compare_means(method = "anova",,label.y = 15)+ stat_compare_means(comparisons = my_comparisons)
6. 二因素直方图绘制
6.1 直方图+标准误
p = ggbarplot(dd,x = "Group1",y = "y",add = "mean_se",color = "Group2", position = position_dodge(0.8)) p
6.2 直方图+标准误+显著性
p + stat_compare_means(aes(group=Group2), label = "p.signif")
7. 二因素折线图绘制
7.1 折线图+标准误
p = ggline(dd,x = "Group1",y = "y",add = "mean_se",color = "Group2", position = position_dodge(0.8)) p
7.2 折线图+标准误+显著性
p + stat_compare_means(aes(group=Group2), label = "p.signif")