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

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

suppressWarnings(library(ggstatsplot)) #加载包 
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 ...

ggbetweenstats(ToothGrowth, dose, len, mean.ci = T, pairwise.comparisons = T)#显示各组均值的置信区间，并进行多重比较。
Note: 95% CI for effect size estimate was computed with 100 bootstrap samples.
# A tibble: 3 x 9
  group1 group2 mean.difference    se t.value    df p.value significance p.value.label
                                         
1 0.5    1                 9.13 0.997    6.48  38.0       0 ***          p <= 0.001   
2 0.5    2                15.5  0.929   11.8   36.9       0 ***          p <= 0.001   
3 1      2                 6.36 0.918    4.9   37.1       0 ***          p <= 0.001   
Note: Shapiro-Wilk Normality Test for len : p-value = 0.109
Note: Bartlett's test for homogeneity of variances for factor dose: p-value = 0.717

正态性检验和方差齐性检验的p值大于0.05，数据符合正态性和方差齐性。

方差分析.png

-方差分析图中给出了方差分析的结果，F值为68.40 > 1， P < 0.001, 说明不同喂养剂量对生长的影响具有统计学差异。

-多重比较中，各组p值均小于0.05，各组均值的差异均有统计学意义，其中dose = 2，len的均值最大。