【发布时间】:2011-09-17 05:35:00
【问题描述】:
我使用以下代码制作了多个图:
set.seed(12345)
a <- data.frame(Glabel=LETTERS[1:7], A=rnorm(7, mean = 0, sd = 1), B=rnorm(7, mean = 0, sd = 1), C=rnorm(7, mean = 0, sd = 1))
T <- data.frame(Tlabel=LETTERS[11:20], A=rnorm(10, mean = 0, sd = 1), B=rnorm(10, mean = 0, sd = 1), C=rnorm(10, mean = 0, sd = 1))
for(i in 2:(ncol(a)-1))
{
for(j in (i+1):ncol(a))
{
r <- 0.08
win.graph(width=10, height=10, pointsize=12)
plot(a[, i],a[, j], pch=19, cex=1, panel.first=grid(col="gray", lty="dotted"))
points(T[, i],T[, j],pch=19, col="white", cex=1)
text(a[, i],a[, j], rownames(a), cex=1,col="black", pos=1)
text(T[, i],T[, j], rownames(T), cex=1,col="blue", pos=1)
arrows(x0=0, y0=0, x1=T[, i], y1=T[, j], length=0.1,col=2,lty=1)
}
}
使用win.graph 我可以维护基地中的所有地块。但我想在 ggplot2 中有我的情节。以下是我的代码 ggplot2。
library(ggplot2)
for(i in 2:(ncol(a)-1))
{
for(j in (i+1):ncol(a))
{
r <- 0.08
p <- ggplot(data=a, mapping=aes(x=a[, i], y=a[, j])) + geom_point() + theme_bw()
p <- p + geom_text(data=a, mapping=aes(x=a[, i], y=a[, j], label=Glabel),
size=3, vjust=1.35, colour="black")
p <- p + geom_segment(data = T, aes(xend = T[ ,i], yend=T[ ,j]),
x=0, y=0, colour="black",
arrow=arrow(angle=25, length=unit(0.25, "cm")))
p <- p + geom_text(data=T, aes(x=T[ ,i], y=T[ ,j], label=Tlabel), size=3, vjust=0, colour="red")
print(p)
}
}
上面的代码只给了我最后一张图。我想知道如何在 ggplot2 中获得所有可能的图表。谢谢
【问题讨论】:
-
尝试在
print(p)之前添加dev.new() -
另一方面,不建议您使用
aes(x=a[, i], ...)映射变量,您可能应该使用melt()ed(= 长格式)data.frames。 -
@baptiste:感谢您的评论。有效。再次感谢。
-
@baptiste:你是对的。我现在对
aes(x=a[, i], ...)有疑问,但可以弄清楚如何制作这样的图表。如果您使用melt()方法回答这个问题,我将不胜感激。谢谢