Сезонные температурные боксы в R

Question

В настоящее время я выполняю некоторые исследования переменных и сгенерировал коробочные графики для 3 различных климатических параметров (Tmin, Tmean, Tmax). Я хотел бы знать, как сгруппировать эти переменные в одном одном блок-графике со структурой, подобной этой: Я видел пару учебных пособий онлайн, но все они требуют, чтобы параметр группировки был назначен в строках фрейма данных, а не в заголовке столбца. Я попытался добавить + tmax в качестве аргумента в tmin, но это привело к ошибкам. Код, который я использовал для генерации моего кода, следующий:

tmin<-  ggplot(prism, aes(x = factor(season, levels=c("spring","summer","fall","winter")), 
                                  y = tmin_c)) +
                geom_boxplot(fill = fill, colour = line, alpha = 0.7) +
                theme_bw() +
                scale_y_continuous(name = "Temperature C") +
                scale_x_discrete(name = "Season") +
                ggtitle("MRL Temperature 1980-2013") +
                theme(plot.title = element_text(hjust = 0.5))
tmin

Решено, вот окончательный рабочий результат для будущей ссылки:

#Temperature
dat <- prism
dat <- dat %>%
       select(1,2,4,5,6) #1year,2season,4tmin,5tmean,6tmax

dat <- reshape2::melt(dat, measure.vars=3:5)


ggplot(dat, aes(y = value, 
                x = factor(season, levels=c("spring","summer","fall","winter")), 
                fill=factor(variable))) +
            geom_boxplot() +
            theme_bw() +
            scale_y_continuous(name = "Temperature C") +
            scale_x_discrete(name = "Season") +
            ggtitle("MRL Temperature 1980-2013") +
            theme(plot.title = element_text(hjust = 0.5))

Результат:

Answer 1

Я ссылаюсь на исходную версию вашего кода с использованием airquality данных.

Представьте себе Month - это ваша season переменная. Сначала преобразуйте данные в длинный формат, используя reshape2::melt; ваш measure.vars будет t_min, t_max, ....

dat <- reshape2::melt(airquality, measure.vars=1:3)

summary(dat)
# Temp           Month            Day          variable       value       
# Min.   :56.00   Min.   :5.000   Min.   : 1.0   Ozone  :153   Min.   :  1.00  
# 1st Qu.:72.00   1st Qu.:6.000   1st Qu.: 8.0   Solar.R:153   1st Qu.: 10.30  
# Median :79.00   Median :7.000   Median :16.0   Wind   :153   Median : 24.00  
# Mean   :77.88   Mean   :6.993   Mean   :15.8                 Mean   : 80.86  
# 3rd Qu.:85.00   3rd Qu.:8.000   3rd Qu.:23.0                 3rd Qu.:136.00  
# Max.   :97.00   Max.   :9.000   Max.   :31.0                 Max.   :334.00  
#                                                              NA's   :44 

Во-вторых, используйте boxplot с : для факторов Month (т. Е. season).

at.key <- c(1:3, 5:7, 9:11, 13:15, 17:19)
# at.key <- (1:(5*4))[(1:(5*4)) %% ((5*4)/4-1) != 0]  ## alternatively

b <- boxplot(value ~ variable:Month, border=2:4, col="white",
             data=dat, at=at.key, xaxt="n",
             main="MRL TMin 1980-2013",
             xlab="Season",
             ylab="Tmin (C)")

mtext(b$names, 1, .5, at=at.key, cex=.8, las=2)