R中的boxplot，美学必须是长度1或与数据相同的长度

Question

我正在对UCI网站每加仑汽车数据进行一些分析：

https://archive.ics.uci.edu/ml/datasets/Auto+MPG

我将第一列考虑到高或低里程：

mpg01 = I(auto1$mpg >= median(auto1$mpg))
Auto = data.frame(mpg01, auto1[,-1])
head(Auto)

    mpg01 cylinders displacement horsepower weight acceleration year origin
 1 FALSE         8          307        130   3504         12.0   70      1
 2 FALSE         8          350        165   3693         11.5   70      1
 3 FALSE         8          318        150   3436         11.0   70      1
 4 FALSE         8          304        150   3433         12.0   70      1
 5 FALSE         8          302        140   3449         10.5   70      1
 6 FALSE         8          429        198   4341         10.0   70      1

现在我想为数据框中的每个列制作一个boxplot，由第一列计算。

vars <- c("cylinders", "displacement", "horsepower", "weight", "acceleration", "year", "origin")
ggplot(Auto) + geom_bar(aes(y=vars, fill=factor(mpg01)))

我得到错误“Aes必须是长度1或与数据相同”

“自动”数据帧的维度为392x8

我可以为每列使用boxplot，但想知道是否有办法将它们合并为一个。 谢谢！

Answer 1

更新以解释生成的错误：生成错误是因为需要定义aes(x, y...)来描述数据框变量应如何映射到geoms中。 在您的情况下，没有为geom_boxplot定义x变量。 为了将x变量定义为df的每个列，需要将df重新整形为长格式（例如，使用reshape2::melt或tidyr::gather ）

以下是应该工作的解决方案，它基于mtcars而不是您的数据。 如果没有，我们可以为您dput(Auto)一次后dput(Auto)进行故障排除。 您生成的图表应该与我附加的图表类似。 首先，重塑您的数据。

library(reshape2)
library(ggplot2)
mtcars_melt <- melt(mtcars)

我现在可以在aes定义x 。 注意：注意与facet_wrap使用时，下面两种情况之间的facet_wrap 。

# First with no facet_wrap
ggplot(mtcars_melt, aes(x=variable, y=value, fill=variable)) + geom_boxplot()
# Case 1 with facet_wrap
ggplot(mtcars_melt, aes(x=variable, y=value, fill=variable)) + geom_boxplot() + facet_wrap(~variable)
# Case 2 with facet_wrap
ggplot(mtcars_melt, aes(x="", y=value, fill=variable)) + geom_boxplot() + facet_wrap(~variable)

在情况1中，我在aes定义x=variable ，但是使用facet_wrap会强制每个facet都存在所有x个变量，但是如果我设置x="" ，它允许每个facet只保存1个x变量。

现在允许y轴有独立的比例，我可以设置scales="free_y"

ggplot(mtcars_melt, aes(x="", y=value, fill=variable)) + geom_boxplot() + facet_wrap(~variable, scales="free_y")

或者，我可以设置scales="free"以应用于x和y轴，并将其与x=variable一起使用以获得类似的解决方案。

ggplot(mtcars_melt, aes(x=variable, y=value, fill=variable)) + geom_boxplot() + facet_wrap(~variable, scales="free")

编辑：以下代码适用于您的特定数据集：

library(reshape2)
library(ggplot2)
vars <- c("cylinders", "displacement", "horsepower", "weight", "acceleration", "year", "origin")
Auto_melt <- melt(Auto[, vars])
ggplot(Auto_melt, aes(x="", y=value, fill=variable)) + geom_boxplot() + facet_wrap(~variable, scales="free_y")

编辑代码按要求分隔mpg：通过包含“mpg01”重新定义变量，并通过mpg id融化数据。 使用mpg01作为aes x值。

Auto <- structure(list(mpg01 = structure(c(2L, 1L, 1L, 1L, 1L), .Label = c("FALSE", "TRUE"), class = "factor"), cylinders = c(8L, 8L, 8L, 8L, 8L), displacement = c(307, 350, 318, 304, 302), horsepower = c(130L, 165L, 150L, 150L, 140L), weight = c(3504L, 3693L, 3436L, 3433L, 3449L), acceleration = c(12, 11.5, 11, 12, 10.5), year = c(70L, 70L, 70L, 70L, 70L), origin = c(1L, 1L, 1L, 1L, 1L)), .Names = c("mpg01", "cylinders", "displacement", "horsepower", "weight", "acceleration", "year", "origin"), row.names = c(NA, 5L), class = "data.frame") 

vars <- c("mpg01", "cylinders", "displacement", "horsepower", "weight", "acceleration", "year", "origin")
Auto_melt <- melt(Auto[, vars], id.vars="mpg01")
ggplot(Auto_melt, aes(x=mpg01, y=value, fill=variable)) + geom_boxplot() +    facet_wrap(~variable, scales="free_y")

Answer 2

我想也许你应该整理你的数据，然后绘制boxplot。 我从网站上下载数据：

> head(df)
  mpg01 cylinders displacement horsepower weight acceleration year origin
1    18         8          307        130   3504         12.0   70      1
2    15         8          350        165   3693         11.5   70      1
3    18         8          318        150   3436         11.0   70      1
4    16         8          304        150   3433         12.0   70      1
5    17         8          302        140   3449         10.5   70      1
6    15         8          429        198   4341         10.0   70      1

使用gather {tidyr}来整理数据。

library("tidyr")
library("dplyr")
library("ggplot2")
tidy_df <- df %>% gather("vars","values",-mpg01)

并且tidy_df是：

> head(tidy_df)
  mpg01      vars values
1    18 cylinders      8
2    15 cylinders      8
3    18 cylinders      8
4    16 cylinders      8
5    17 cylinders      8
6    15 cylinders      8

然后你可以绘制boxplot

ggplot(data=tidy_df,aes(vars,values)) + geom_boxplot(aes(fill=vars))

它看起来像这样：