数据分析:多组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
相关推荐
KIDAKN几秒前
Redis 分布式锁
数据库·redis·分布式
程序新视界8 分钟前
如何为MySQL中的JSON字段设置索引
数据库·mysql
Ultipa17 分钟前
查询语言的进化:SQL之后,为什么是GQL?数据世界正在改变
数据库·sql·图数据库·gql
LB211220 分钟前
SQL隐式链接显式连接
大数据·数据库·sql
隔壁阿布都1 小时前
spring boot + mybatis 使用线程池异步修改数据库数据
数据库·spring boot·mybatis
没有梦想的咸鱼185-1037-16635 小时前
【高分论文密码】大尺度空间模拟预测与数字制图
信息可视化·数据分析·r语言
MAGICIAN...8 小时前
【Redis】--持久化机制
数据库·redis·缓存
我真的是大笨蛋9 小时前
JVM调优总结
java·jvm·数据库·redis·缓存·性能优化·系统架构
步步为营DotNet10 小时前
5-2EFCore性能优化
数据库·性能优化·.net
2501_9200470311 小时前
Redis-集群
数据库·redis·bootstrap