双因素方差分析_R海拾遗双因素重复测量方差分析

重复测量方差分析

sunqi

2020/7/26

概述

双因素的重复测量资料方差分析

代码

数据获得

library(tidyverse)
library(ggpubr)
library(rstatix)
rm(list=ls())
set.seed(123)
data("selfesteem2", package = "datarium")
# 抽样
selfesteem2 %>% sample_n_by(treatment, size = 1)

## # A tibble: 2 x 5
##   id    treatment    t1    t2    t3
##   
## 1 3     ctr          93    92    89
## 2 3     Diet         91    91    92

# 数据含有5个变量，其中三个时间点，一个为治疗方式，一个为id
# 个案为12个，每个人进行3次测量，2种治疗

# 对数据进行长转宽
#将id和时间转换为因子
selfesteem2 %
  gather(key = "time", value = "score", t1, t2, t3) %>%
  convert_as_factor(id, time)
# 检查数据
set.seed(123)
selfesteem2 %>% sample_n_by(treatment, time, size = 1)

## # A tibble: 6 x 4
##   id    treatment time  score
##   
## 1 3     ctr       t1       93
## 2 3     ctr       t2       92
## 3 10    ctr       t3       81
## 4 2     Diet      t1      100
## 5 6     Diet      t2       75
## 6 11    Diet      t3       91

# 描述数据
selfesteem2 %>%
  group_by(treatment, time) %>%
  get_summary_stats(score, type = "mean_sd")

## # A tibble: 6 x 6
##   treatment time  variable     n  mean    sd
##   
## 1 ctr       t1    score       12  88    8.08
## 2 ctr       t2    scor