В R & dabestr, как правильно получить сгруппированные различия?

Question

Использование dabestr пакета Я пытаюсь получить различия между двумя наборами контрольных и тестовых данных.Немного изменив пример из файла справки, который я попробовал:

library(dabestr)

N <- 70
c1 <- rnorm(N, mean = 50, sd = 20)
t1 <- rnorm(N, mean = 200, sd = 20)
ID <- seq(1:N)
long.data <- tibble::tibble(ID = ID, Control1 = c1, Test1 = t1)

meandiff1 <- long.data %>%
  tidyr::gather(key = Group, value = Measurement, Control1:Test1)

ID <- seq(1:N) + N
c2 <- rnorm(N, mean = 100, sd = 70)
t2 <- rnorm(N, mean = 100, sd = 70)
long.data <- tibble::tibble(ID = ID, Control2 = c2, Test2 = t2)

meandiff2 <- long.data %>%
  tidyr::gather(key = Group, value = Measurement, Control2:Test2)

meandiff <- dplyr::bind_rows(meandiff1, meandiff2)

paired_mean_diff <- 
  dabest(meandiff, x = Group, y = Measurement,
         idx = c("Control1", "Test1", "Control2", "Test2"),
         paired = TRUE,
         id.col = ID)

plot(paired_mean_diff)

Я получаю следующие результаты:

Так что сравнивается не только вседо Control1, но также опция paired = TRUE, похоже, не имеет никакого эффекта.Я надеялся получить что-то похожее на примеры со страницы пакета:

Есть ли какие-нибудь указатели о том, как этого добиться?

Answer 1

Для парного графика вы хотите вложить опцию ключевого слова idx следующим образом:

paired_mean_diff <- 
  dabest(meandiff, x = Group, y = Measurement,
    idx = list(c("Control1", "Test1"), 
               c("Control2", "Test2")),
    paired = TRUE,
    id.col = ID)