R语言用ggstatsplot包做方差分析和绘图

ggstatsplot包可以方便快捷地自动进行统计学计算并绘出图形。本文将以R语言自带数据集ToothGrowth为例，探究不同的喂养剂量对生长速度的影响是否有差异。

attach(ToothGrowth) #绑定数据集 
str(ToothGrowth) #查看数据结构

 $ 'data.frame': 60 obs. of 3 variables: $ len : num 4.2 11.5 7.3 5.8 6.4 10 11.2 11.2 5.2 7 ...
 $ supp: Factor w/ 2 levels "OJ","VC": 2 2 2 2 2 2 2 2 2 2 ...
 $ dose: num 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ...

suppressWarnings(library(ggstatsplot)) #加载包 
ggbetweenstats(ToothGrowth, dose, len, mean.ci = T, pairwise.comparisons = T)#显示各组均值的置信区间，并进行多重比较。

统计结果.png

统计结果图中，正态性检验和方差气性检验的p值大于0.05，数据符合正态性和方差齐性。

方差分析.png

方差分析图中给出了方差分析的结果，F值为68.40 > 1， P < 0.001, 说明不同喂养剂量对生长的影响具有统计学差异。
多重比较中，各组p值均小于0.05，各组均值的差异均有统计学意义，其中dose = 2，len的均值最大。