0

我正在使用带有rstatix包的 r 来分析我的数据，但是rstatix当我尝试在第三个变量的每个级别确定简单的双向交互时，我收到了来自包的警告消息。

我的实验设计是 3 路混合方差分析，其中 2 和 1 在主题因素内。主体因子之间的 2 是rpm和ratio，而因子之间的 1 是time。是在几个时期val测量的因变量。time

我已将我的数据简化如下：

library(tidyverse)
library(rstatix)

df <- tribble(
  ~id, ~time, ~rpm, ~ratio, ~val,
  1,   0,     50,   1,      6.49,
  2,   0,     50,   1,      6.51,
  3,   0,     80,   1,      6.49,
  4,   0,     80,   1,      6.49,
  5,   0,     50,   3,      6.56,
  6,   0,     50,   3,      6.56,
  7,   0,     80,   3,      6.50,
  8,   0,     80,   3,      6.44,
  9,   0,     50,   6,      6.52,
  10 , 0,     50,   6,      6.59,
  11,  0,     80,   6,      6.65,
  12,  0,     80,   6,      6.60,
  1,   3,     50,   1,      4.61,
  2,   3,     50,   1,      4.80,
  3,   3,     80,   1,      4.33,
  4,   3,     80,   1,      4.36,
  5,   3,     50,   3,      4.81,
  6,   3,     50,   3,      4.70,
  7,   3,     80,   3,      4.80,
  8,   3,     80,   3,      4.82,
  9,   3,     50,   6,      4.60,
  10 , 3,     50,   6,      4.63,
  11,  3,     80,   6,      4.74,
  12,  3,     80,   6,      4.72
)

df <-df %>% mutate(across(id:ratio, as.factor))

twoway_1 <- df %>%
  group_by(time) %>%
  anova_test(dv = val, wid = id, between = c(rpm, ratio))  # OK
twoway_1

twoway_2 <- df %>%
  group_by(rpm) %>%
  anova_test(dv = val, wid = id, between = c(time, ratio)) # Warning messages
twoway_2

我收到的警告信息是：

Warning: Problem with `mutate()` column `data`.
i `data = map(.data$data, .f, ...)`.
i The 'wid' column contains duplicate ids across between-subjects variables. Automatic unique id will be created.

我认为警告是因为该id列有两次数字 1 到 12 的集合。该组重复两次，因为每个受试者id在 = 0 时测量两次，time然后在 3 时测量。

更新：

我认为这是rstatix. 以下等效项不会产生任何问题：

model1 <- aov(val ~ (time * ratio), data = subset(df, rpm=="50"))
model2 <- aov(val ~ (time * ratio), data = subset(df, rpm=="80"))
4

0 回答 0