Как создать гистограмму с динамическими переменными для каждой группы и включенными столбцами ошибок? - PullRequest
Новая
       15

Как создать гистограмму с динамическими переменными для каждой группы и включенными столбцами ошибок?

0 голосов
/ 29 апреля 2018

Я пытаюсь создать динамическую гистограмму в том же формате, что и этот график модели здесь, используя R. К сожалению, мой учитель не предоставил нам R-коды для графика модели и, следовательно, эту динамическую Гистограмма отнимает у меня часы и часы, и я до сих пор не могу сделать ее точно такой же, как график модели. Данные, которые я использовал для создания динамической гистограммы, находятся здесь: 

tidyadhd <- structure(list(voltage = c(9.5, 10, 9.2, 6, 9, 11, 9, 9.2, 8, 
9, 11, 9, 8.2, 8, 9, 9, 6, NA, 9.5, 9, 7.1, 9.3, 9, 9.5, 9, 9.2, 
9.3, 9, 8.5, 9, 9.2, 9.3, 8.2, 9, 9.3, NA, 9.5, 9, 9.2, 9.3, 
9, 9.5, 9, 9.2, 9.3, 9, 9.5, 9, 9, 9, 9.3, 8, 7, 8, 5.1, 4.3, 
4.5, 5.3, 4.5, 5.1, 4.3, 4.5, 5.3, 4.5, 5.1, 4.3, 4.3, 4.3, 5.3, 
6.4, 5, 8), state = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L
), .Label = c("Rest", "Cue"), class = "factor"), treatment = c("placebo", 
"placebo", "placebo", "placebo", "placebo", "placebo", "placebo", 
"placebo", "placebo", "placebo", "placebo", "placebo", "placebo", 
"placebo", "placebo", "placebo", "placebo", "placebo", "placebo", 
"placebo", "placebo", "placebo", "placebo", "placebo", "placebo", 
"placebo", "placebo", "placebo", "placebo", "placebo", "placebo", 
"placebo", "placebo", "placebo", "placebo", "placebo", "ritalin", 
"ritalin", "ritalin", "ritalin", "ritalin", "ritalin", "ritalin", 
"ritalin", "ritalin", "ritalin", "ritalin", "ritalin", "ritalin", 
"ritalin", "ritalin", "ritalin", "ritalin", "ritalin", "ritalin", 
"ritalin", "ritalin", "ritalin", "ritalin", "ritalin", "ritalin", 
"ritalin", "ritalin", "ritalin", "ritalin", "ritalin", "ritalin", 
"ritalin", "ritalin", "ritalin", "ritalin", "ritalin"), patient = c(1L, 
2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 
16L, 17L, 18L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 
12L, 13L, 14L, 15L, 16L, 17L, 18L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 
8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 1L, 2L, 
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 
17L, 18L)), .Names = c("voltage", "state", "treatment", "patient"
), row.names = c(NA, -72L), class = "data.frame")

и изображение части таблицы вывода показано на этом рисунке здесь: two-way ANOVA table

Может ли кто-нибудь научить меня, как создать динамический линейный график с точно таким же форматом, как у этой модели, пожалуйста? Спасибо.

1 Ответ

0 голосов
/ 29 апреля 2018

Вы можете начать с этого сюжета. 

library(tidyverse)
tidyadhd %>%
  group_by(treatment, state) %>%
  mutate(n = n()) %>%
  ungroup() %>%
  group_by(treatment, state, n) %>%
  summarise_at(vars(voltage), funs(mean, sd = sd), na.rm = TRUE) %>%
  mutate(se = sd / sqrt(n)) %>%
  ggplot(aes(treatment, mean, group = state, fill = state)) + geom_bar(stat = "identity", position = "dodge", color = "black") + 
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se), width = 0.2, position = position_dodge(0.9)) +
  theme_classic() + scale_fill_manual(values = c("gray", "white")) +
  theme(
    legend.position = "bottom",
    panel.grid.minor = element_blank(),
    panel.grid.major.y = element_blank(),
    legend.title=element_blank())

enter image description here

...