假设我有以下数据集1
我想为同一图表上的每个水记录器值制作箱线图。在我检查的每个地方,人们都有一个要使用的因子变量。但是,我不想要因素,我将水记录器编号作为列名。我可以用通用的 boxplot 命令来做到这一点:boxplot(data$colname1, data$colname2, data$colname3, and so on)
但是我怎样才能用更好的图形来做到这一点,比如在ggplot2
.
假设我有以下数据集1
我想为同一图表上的每个水记录器值制作箱线图。在我检查的每个地方,人们都有一个要使用的因子变量。但是,我不想要因素,我将水记录器编号作为列名。我可以用通用的 boxplot 命令来做到这一点:boxplot(data$colname1, data$colname2, data$colname3, and so on)
但是我怎样才能用更好的图形来做到这一点,比如在ggplot2
.
如果没有实际数据,很难向您展示您需要使用的确切代码,但是在瞥见之后png
,我建议您尝试以下几行:
library(reshape2)
library(ggplot2)
df <- melt(your_data)
ggplot(df, aes(x=variable, y=value)) + geom_boxplot()
这段代码可能需要一些调整。如果它不起作用并且调整不明显,请以方便我们使用的方式发布一些示例数据。屏幕截图中的数据意味着我们必须手动复制粘贴每个数字,而很少有人愿意这样做。
为了澄清一般过程:将melt
所有列“堆叠”在一起,并添加一个名为 的变量variable
,它引用旧的列名。你可以把它交给ggplot
并说不同的值variable
应该在x
轴上,这就是你想要的。例如,看看women
:
head(women)
height weight
1 58 115
2 59 117
3 60 120
4 61 123
5 62 126
6 63 129
str(women)
'data.frame': 15 obs. of 2 variables:
$ height: num 58 59 60 61 62 63 64 65 66 67 ...
$ weight: num 115 117 120 123 126 129 132 135 139 142 ...
您会看到这women
是一个包含 15 个观测值和两列的数据框,height
并且weight
.
现在,让melt
他们:
df <- melt(women)
head(df)
variable value
1 height 58
2 height 59
3 height 60
4 height 61
5 height 62
6 height 63
str(df)
'data.frame': 30 obs. of 2 variables:
$ variable: Factor w/ 2 levels "height","weight": 1 1 1 1 1 1 1 1 1 1 ...
$ value : num 58 59 60 61 62 63 64 65 66 67 ...
现在您看到它有 30 个观察值和两列:variable
和value
. variable
标识旧列。
让我们把它交给ggplot
:
ggplot(df, aes(x=variable, y=value)) + geom_boxplot()
产量:
在这里,您有原始women
数据集中两列的箱线图。
这是基于与 coffeinjunky 相同原理的另一个答案,但更具体到您的数据集。由于您没有提供数据集,我创建了一个具有相似列名的虚拟数据集:
d <- data.frame(x=rep(0,8))
d$`Logger 1_Water_Level` <- c(1,2,3,4,5,3,4,5)
d$`Logger 2_Water_Level` <- c(7,9,2,6,8,9,2,3)
您需要重塑数据集,以便获得标识记录器的因子变量。假设您有两个记录器并且来自记录器的数据存储在第 2 列和第 3 列中,您可以使用以下代码从存储数据的宽格式(即每个记录器的单独列)到长绘图所需的格式ggplot2
(即用于水位测量的单列,每个记录器由名为 的列中的数字标识Logger
)
d_long <- reshape(d, varying=2:3, direction="long", timevar="Logger",v.names="Water_Level", times=1:2)
d_long$Logger <- as.factor(d_long$Logger)
现在您可以使用以下方法绘制测量值ggplot2
:
p <- ggplot(d_long, aes(x=Logger, y=Water_Level))
p <- p + geom_boxplot()
p