ggstatsplot包可以方便快捷地自動(dòng)進(jìn)行統(tǒng)計(jì)學(xué)計(jì)算并繪出圖形。本文將以R語(yǔ)言自帶數(shù)據(jù)集ToothGrowth為例，探究不同的喂養(yǎng)劑量對(duì)生長(zhǎng)速度的影響是否有差異。

suppressWarnings(library(ggstatsplot)) #加載包 
str(ToothGrowth) #查看數(shù)據(jù)結(jié)構(gòu)

 $ '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)#顯示各組均值的置信區(qū)間，并進(jìn)行多重比較。
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
  <chr>  <chr>            <dbl> <dbl>   <dbl> <dbl>   <dbl> <chr>        <chr>        
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

正態(tài)性檢驗(yàn)和方差齊性檢驗(yàn)的p值大于0.05，數(shù)據(jù)符合正態(tài)性和方差齊性。

方差分析.png

-方差分析圖中給出了方差分析的結(jié)果，F(xiàn)值為68.40 > 1， P < 0.001, 說(shuō)明不同喂養(yǎng)劑量對(duì)生長(zhǎng)的影響具有統(tǒng)計(jì)學(xué)差異。

-多重比較中，各組p值均小于0.05，各組均值的差異均有統(tǒng)計(jì)學(xué)意義，其中dose = 2，len的均值最大。