数据分析:多组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
相关推荐
paopaokaka_luck15 分钟前
基于SpringBoot+Vue的酒类仓储管理系统
数据库·vue.js·spring boot·后端·小程序
薛晓刚1 小时前
哪个领域数据库最难替换?
数据库
过期的秋刀鱼!1 小时前
用“做饭”理解数据分析流程(Excel三件套实战)
数据挖掘·数据分析·excel·powerbi·数据分析入门
芷栀夏2 小时前
基于Anything LLM的本地知识库系统远程访问实现路径
数据库·人工智能
软件2052 小时前
【redis使用场景——缓存——数据淘汰策略】
数据库·redis·缓存
ChinaRainbowSea2 小时前
9-2 MySQL 分析查询语句:EXPLAIN(详细说明)
java·数据库·后端·sql·mysql
时序数据说2 小时前
Java类加载机制及关于时序数据库IoTDB排查
java·大数据·数据库·物联网·时序数据库·iotdb
deeper_wind2 小时前
MySQL数据库基础(小白的“升级打怪”成长之路)
linux·数据库·mysql
加勒比海涛2 小时前
Spring Cloud Gateway 实战:从网关搭建到过滤器与跨域解决方案
数据库·redis·缓存
belldeep2 小时前
java:如何用 JDBC 连接 TDSQL 数据库
java·数据库·jdbc·tdsql