Блог

Когда я призываю emmeans к очень простым повторным измерениям anova, подобным этому:

 library(tidyverse)
library(emmeans)

df <- tibble(patient = rep(1:5, each=4),
                       drug = rep(1:4, times=5),
                       response = c(30, 28, 16, 34,
                                    14, 18, 10, 22,
                                    24, 20, 18, 30,
                                    38, 34, 20, 44,
                                    26, 28, 14, 30)) %>%  
  mutate(drug = as.factor(drug), patient = as.factor(patient)) 


my_aov <- aov(response ~ drug   Error(patient / drug), data = df) 

emmeans(my_aov, "drug")
 

Я получаю следующую ошибку, и я не знаю, почему:

Note: re-fitting model with sum-to-zero contrasts Error in terms.formula(formula, "Error", data = data) : 'data' argument is of the wrong type