【发布时间】:2015-01-14 11:47:11
【问题描述】:
我在 R 中工作。我有三个矩阵:d1、d2、d3,并希望从矩阵的每一行生成一个箱线图。像这样:boxplot(d1[1,], d2[1,], d3[1,])。
如何为矩阵的每一行自动生成图。我需要一个循环吗?应用功能可以帮助我吗?此外,在一个绘图窗口中绘制 12 个箱线图会很棒。 感谢您的帮助!
【问题讨论】:
-
为什么不先“整理”数据呢?您为每个变量使用一行。
我在 R 中工作。我有三个矩阵:d1、d2、d3,并希望从矩阵的每一行生成一个箱线图。像这样:boxplot(d1[1,], d2[1,], d3[1,])。
如何为矩阵的每一行自动生成图。我需要一个循环吗?应用功能可以帮助我吗?此外,在一个绘图窗口中绘制 12 个箱线图会很棒。 感谢您的帮助!
【问题讨论】:
这是我认为您想要的一个可重复的示例。假设你有三个 4x40 矩阵:
set.seed(1)
d1 <- matrix(rnorm(40), nrow = 4)
d2 <- matrix(rnorm(40), nrow = 4)
d3 <- matrix(rnorm(40), nrow = 4)
您应该首先组织您的数据,例如堆叠它们并添加一列:
d <- rbind(t(d1), t(d2), t(d3))
d <- cbind(d, rep(1:3, each = 10))
现在您可以绘制前三个列的箱线图:
boxplot(d[, 1] ~ d[, 5])
如果您希望每个矩阵有四个,您可以使用循环:
par(mfrow=c(4,1))
for (i in 1:4){
boxplot(d[, i] ~ d[, 5])
}
par(mfrow=c(1,1))
但是,如果您想要更令人印象深刻的图形,请尝试使用 lattice 或 ggplot2 包。
希望对你有帮助
【讨论】:
你可以在转置矩阵上调用boxplot:
# Create some faux data
x <- matrix(rnorm(50), 5, 10, dimnames = list(LETTERS[1:5], letters[1:10]))
# Transpose and plot
boxplot(t(x), col = rainbow(5))
【讨论】:
boxplot(d1[1,], d2[1,], d3[1,]) 替换为 @ 987654325@) 我最终会得到 100 个绘图,这就是为什么我更愿意在一个绘图窗口中绘制 12 个绘图。
# your matrices (you have 5 variables and 20 observations per matrix, so you'll have to transpose)
m1 <- matrix(rnorm(100), 5, 20)
m2 <- matrix(rnorm(100), 5, 20)
m3 <- matrix(rnorm(100), 5, 20)
# tidy data
library(reshape2)
Df <- melt(list(data.frame(t(m1)), data.frame(t(m2)), data.frame(t(m3))))
# boxplot
library(ggplot2)
ggplot(data = Df, aes(L1, value, color = factor(L1))) + geom_boxplot() + facet_wrap(~variable)
【讨论】: