数据分析:多组ANOVA和后置检验

介绍

单因素多个levels数据在做方差分析后,需要进一步做如t检验等置后检验,判断两两组间差异。

加载R包和准备数据

{r} 复制代码
# Load required R packages
library(tidyverse)
library(rstatix)
library(ggpubr)
  
data("PlantGrowth")
set.seed(1234)
PlantGrowth %>% sample_n_by(group, size = 1)

查看统计量

{r} 复制代码
PlantGrowth %>%
  group_by(group) %>%
  get_summary_stats(weight, type = "mean_sd")

使用ANOVA test比较多组的平均值

{r} 复制代码
res.aov <- PlantGrowth %>% anova_test(weight ~ group)
res.aov

组间t检验

{r} 复制代码
pwc <- PlantGrowth %>%
  pairwise_t_test(weight ~ group, p.adjust.method = "bonferroni")
pwc

p值的坐标轴参数

{r} 复制代码
pwc <- pwc %>% add_xy_position(x = "group")
pwc

可视化

{r} 复制代码
ggboxplot(PlantGrowth, x = "group", y = "weight") +
  stat_pvalue_manual(pwc, label = "p.adj", tip.length = 0, step.increase = 0.1) +
  labs(
    subtitle = get_test_label(res.aov, detailed = TRUE),
    caption = get_pwc_label(pwc))

Reference

  1. HOW TO PERFORM T-TEST FOR MULTIPLE GROUPS IN R
相关推荐
Zzz 小生1 小时前
Claude Code学习笔记(四)-助你快速搭建首个Python项目
大数据·数据库·elasticsearch
nongcunqq4 小时前
abap 操作 excel
java·数据库·excel
rain bye bye5 小时前
calibre LVS 跑不起来 就将setup 的LVS Option connect下的 connect all nets by name 打开。
服务器·数据库·lvs
阿里云大数据AI技术6 小时前
云栖实录|MaxCompute全新升级:AI时代的原生数据仓库
大数据·数据库·云原生
不剪发的Tony老师7 小时前
Valentina Studio:一款跨平台的数据库管理工具
数据库·sql
weixin_307779137 小时前
在 Microsoft Azure 上部署 ClickHouse 数据仓库:托管服务与自行部署的全面指南
开发语言·数据库·数据仓库·云计算·azure
六元七角八分7 小时前
pom.xml
xml·数据库
虚行7 小时前
Mysql 数据同步中间件 对比
数据库·mysql·中间件
奥尔特星云大使7 小时前
mysql读写分离中间件Atlas安装部署及使用
数据库·mysql·中间件·读写分离·atlas
牛马baby8 小时前
【mysql】in 用到索引了吗?
数据库·mysql·in