到boxplot

Question

我有代码来生成好的barplot，并且我尝试用相同的数据创建一个盒图。条形图显示所有人员(id)的“响应”计数。我想为每种类型的“响应”创建一个方格图来替换这3条。应从每个参与者的具体“响应”计数中计算出“方框图”。到目前为止，没有运气，因为我被困在如何计算每个参与者的反应。

当前代码：

df %>%
  ggplot(position = dodge) + 
  labs(title= "question") +
  geom_bar(aes(x = response), fill="red") +
  labs(y = "count", x = "responses") +
  scale_y_continuous(breaks=seq(0,100,20), limits = c(0,100)) 

产出：

​

数据样本：

structure(list(id = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 
3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 
4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5), 
    response = c(0, 1, 1, 0, 0, 0, 1, -1, 1, -1, 0, 1, -1, 1, 
    0, 0, 0, 0, 1, 1, 1, -1, 0, 1, 0, 1, 1, -1, 0, 1, 1, 1, 0, 
    1, 0, 0, 1, -1, 0, 1, 1, 1, -1, 1, 1, 1, 0, 0, -1, 1, 1, 
    -1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 
    1, 1, 0, 0, 0), iscorrect = c(0, 1, 1, 0, 0, 0, 1, 0, 1, 
    0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 1, 1, 0, 
    0, 1, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 
    0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 
    0, 1, 0, 0, 1, 1, 0, 0, 0), min = c(100, 150, 150, 
    50, 50, 50, 150, 100, 100, 100, 50, 100, 50, 150, 150, 150, 
    50, 100, 100, 100, 150, 150, 50, 50, 50, 150, 150, 100, 50, 
    100, 100, 150, 150, 50, 50, 50, 150, 100, 100, 100, 50, 100, 
    50, 150, 150, 150, 50, 100, 100, 100, 150, 150, 50, 50, 50, 
    150, 150, 100, 50, 100, 100, 150, 100, 50, 50, 50, 150, 100, 
    100, 50, 150, 100, 50, 150, 150), max = c(125.4, 180.8, 
    180.8, 62.4, 62.4, 62.4, 180.8, 125.4, 125.4, 125.4, 62.4, 
    125.4, 62.4, 180.8, 180.8, 180.8, 62.4, 125.4, 125.4, 125.4, 
    180.8, 180.8, 62.4, 62.4, 62.4, 180.8, 180.8, 125.4, 62.4, 
    125.4, 125.4, 180.8, 180.8, 62.4, 62.4, 62.4, 180.8, 125.4, 
    125.4, 125.4, 62.4, 125.4, 62.4, 180.8, 180.8, 180.8, 62.4, 
    125.4, 125.4, 125.4, 180.8, 180.8, 62.4, 62.4, 62.4, 180.8, 
    180.8, 125.4, 62.4, 125.4, 125.4, 180.8, 125.4, 62.4, 62.4, 
    62.4, 180.8, 125.4, 125.4, 62.4, 180.8, 125.4, 62.4, 180.8, 
    180.8), time = c(5, 7, 9, 5, 1, 7, 1, 1, 7, 3, 9, 9, 
    3, 5, 3, 1, 9, 5, 1, 7, 9, 3, 5, 7, 1, 5, 7, 3, 3, 9, 5, 
    7, 9, 5, 1, 7, 1, 1, 7, 3, 9, 9, 3, 5, 3, 1, 9, 5, 1, 7, 
    9, 3, 5, 7, 1, 5, 7, 3, 3, 9, 9, 7, 5, 7, 5, 9, 5, 3, 1, 
    1, 9, 7, 3, 3, 1)), row.names = c(NA, -75L), class = c("tbl_df", 
"tbl", "data.frame"))

Answer 1

您可以使用以下代码：

data %>% 
  group_by(id) %>%
  count(response) %>%
  mutate(response = as.factor(response)) %>%
  ggplot(aes(x = response, y = n)) +
  geom_boxplot(fill = "red") +
  labs(y = "count", x = "responses")

输出：

​

​

Answer 2

你可以试试：

library(dplyr)
library(ggplot2)

df %>% 
  group_by(id, response) %>% 
  count() %>%
  mutate(id = factor(id), response = factor(response)) %>%
  ggplot(aes(response, n)) + 
  geom_boxplot(fill = "red") +
  scale_y_continuous(name = "Number of responses per participant")

​

​

请注意，对于小计数之类的离散数据，方格图不能很好地工作(除非实际数据有更多的参与者，每个响应的计数要高得多)