Несколько коробочных участков для 4 групп на одном участке - PullRequest
0 голосов
/ 03 июля 2018

enter image description here Я хотел бы построить коробчатые графики для набора данных, который включает четыре категориальные, такие как: Good, Bad, VeryGood и VeryBad, и четыре нормальных распределения. Мой вопрос, как сделать четыре категориальных графика с четырьмя различными нормальными распределениями на одном графике отдельно друг от друга, я пытался (см. Ниже), но там выглядит беспорядок. Я использовал пример, который нашел здесь и внес некоторые изменения. Я добавил еще один отредактированный график, который мне бы хотелось, чтобы каждый график в виде квадрата был похож на этот, который более четкий, и каждые четыре категориальных (синий, желтый, красный и зеленый) являются четкими. enter image description here

par(mfrow=c(2,2))
df <- data.frame(id = c(rep("Good",200), rep("Bad", 200),
                        rep("VeryGood",200), rep("VeryBad",200)),
           F1 = c(rnorm(200,10,2), rnorm(200,8,1), rnorm(200,5,2),rnorm(200,7,3)),
           F2 = c(rnorm(200,7,1),  rnorm(200,6,1), rnorm(200,8,1),rnorm(200,12,4)),
           F3 = c(rnorm(200,6,2),  rnorm(200,9,3),rnorm(200,12,3),rnorm(200,15,2)),
           F4 = c(rnorm(200,12,3), rnorm(200,8,2),rnorm(200,8,5),rnorm(200,5,1)))

boxplot(df[,-1], xlim = c(0.5, ncol(df[,-1])+0.9), 
        boxfill=rgb(1, 1, 1, alpha=1), border=rgb(1, 1, 1, alpha=1)) #invisible boxes
boxplot(df[which(df$id=="Good"), -1], xaxt = "n", add = TRUE, boxfill="red", boxwex=0.25, 
        at = 1:ncol(df[,-1]) - 0.15) #shift these left by -0.15
boxplot(df[which(df$id=="Bad"), -1], xaxt = "n", add = TRUE, boxfill="blue", boxwex=0.25,
        at = 1:ncol(df[,-1]) + 0.15) #shift these right by +0.15
boxplot(df[which(df$id=="VeryBad"), -1], xaxt = "n", add = TRUE, boxfill="green", boxwex=0.25,
        at = 1:ncol(df[,-1]) + 0.25) #shift these right by +0.15
boxplot(df[which(df$id=="VeryGood"), -1], xaxt = "n", add = TRUE, boxfill="yellow", boxwex=0.25,
        at = 1:ncol(df[,-1]) + 0.45) #shift these right by +0.15

1 Ответ

0 голосов
/ 04 июля 2018

Если вы не настроены на использование графики Base R и смотрите на новый сюжет, который вы добавили в вопрос, я думаю, это то, что вы ищете:

library(dplyr)
library(tidyr)
library(ggplot2)

df <- data.frame(id = c(rep("Good",200), rep("Bad", 200),
                        rep("VeryGood",200), rep("VeryBad",200)),
                 F1 = c(rnorm(200,10,2), rnorm(200,8,1), rnorm(200,5,2),rnorm(200,7,3)),
                 F2 = c(rnorm(200,7,1),  rnorm(200,6,1), rnorm(200,8,1),rnorm(200,12,4)),
                 F3 = c(rnorm(200,6,2),  rnorm(200,9,3),rnorm(200,12,3),rnorm(200,15,2)),
                 F4 = c(rnorm(200,12,3), rnorm(200,8,2),rnorm(200,8,5),rnorm(200,5,1)))

df2 <- tidyr::gather(df, key = "FVar", value = "value", F1:F4)

df2 %>% 
  ggplot(aes(id, value, fill = id)) + 
  geom_boxplot() + 
  facet_grid(. ~ FVar) + 
  theme(axis.text.x = element_text(angle = 90, hjust = 0.5, vjust = 0.5))

enter image description here

...