数据分析:多组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
相关推荐
l1t1 小时前
利用DeepSeek实现服务器客户端模式的DuckDB原型
服务器·c语言·数据库·人工智能·postgresql·协议·duckdb
大数据CLUB5 小时前
基于spark的澳洲光伏发电站选址预测
大数据·hadoop·分布式·数据分析·spark·数据开发
MarkHard1236 小时前
如何利用redis使用一个滑动窗口限流
数据库·redis·缓存
island13147 小时前
【Redis#10】渐进式遍历 | 数据库管理 | redis_cli | RES
数据库·redis·bootstrap
心想事成的幸运大王7 小时前
Redis的过期策略
数据库·redis·缓存
倔强的石头_7 小时前
CentOS 上安装KingbaseES(ISO包)详细教程
数据库
2401_897930068 小时前
使用Docker轻松部署Neo4j图数据库
数据库·docker·neo4j
诗句藏于尽头8 小时前
Django模型与数据库表映射的两种方式
数据库·python·django
智数研析社8 小时前
9120 部 TMDb 高分电影数据集 | 7 列全维度指标 (评分 / 热度 / 剧情)+API 权威源 | 电影趋势分析 / 推荐系统 / NLP 建模用
大数据·人工智能·python·深度学习·数据分析·数据集·数据清洗
寻星探路9 小时前
数据库造神计划第六天---增删改查(CRUD)(2)
java·大数据·数据库