Как построить эффект многоуровневой модели лечения в R - PullRequest
0 голосов
/ 11 октября 2019

Я хочу построить Group-Time-Interaction моей многоуровневой модели. Мне удалось построить его без многоуровневой структуры. Как я могу включить вложение? Каждого пациента измеряли 3 раза (до лечения, в середине и после) на одном и том же приборе. Пациенты были рандомизированы на 2 группы (список ожидания и лечение). Это моя модель:

m8 <- lmer(score ~ 1 + +time*Group +(1| subject), data = df2) 

Это то, что я сделал:

score <- data %>%
  filter(!is.na(score)) %>%  
  group_by(ti, Gruppe, ar, ID) %>%
  summarise(Max_Height = max(score)) %>%  # Calculating max score
  ungroup() %>%  # Need to ungroup so that the pipe doesn't get confused
  group_by(ti, Gruppe, ar) %>%
  summarise(score = mean(Max_Height)) 
(basic_mm_scatter <- ggplot(score, aes(ti, score, colour = Gruppe)) +
    geom_point() +
    theme_bw() +stat_smooth(method = "lm",))

Я также хотел бы включить строки ошибок. Не знаю, как, хотя. Помощь очень ценится !!! Большое спасибо!

...