数据分析:多组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
相关推荐
BTU_YC13 分钟前
Neo4j查询计划完全指南:读懂数据库的“执行蓝图“
数据库·neo4j
非极限码农19 分钟前
Neo4j图数据库上手指南
大数据·数据库·数据分析·neo4j
mit6.82443 分钟前
[C# starter-kit] 命令/查询职责分离CQRS | MediatR |
java·数据库·c#
咋吃都不胖lyh1 小时前
SQL-多对多关系
android·mysql·数据分析
苏打水com1 小时前
数据库进阶实战:从性能优化到分布式架构的核心突破
数据库·后端
shan~~2 小时前
linux达梦数据库操作
linux·数据库·chrome
liliangcsdn2 小时前
LLM时代基于unstructured解析非结构化pdf
linux·服务器·数据分析
在云上(oncloudai)2 小时前
深入解析 Amazon Athena:云上高效数据分析的关键引擎
数据挖掘·数据分析
武文斌773 小时前
项目学习总结:LVGL图形参数动态变化、开发板的GDB调试、sqlite3移植、MQTT协议、心跳包
linux·开发语言·网络·arm开发·数据库·嵌入式硬件·学习
CoderIsArt3 小时前
SQLite架构
数据库·sqlite