所有可能的变量组合

所有可能的变量组合

我想为所有可能的变量组合绘制图形。我的代码如下:

 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))

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")
dev.new()
  print(p)
}
 }

此代码可以正常工作。但是不建议在此使用此方法(See @baptiste comment),该方法在功能上不起作用。我想知道什么是完成此任务的最佳建议方法。在此先感谢您的帮助。

最佳答案

好吧,这是垃圾,但是我能做的最好。它的效率极低,因为它会通过lapply在每个循环中重新创建部分数据。也许其他人有更好的东西:

MAT <- outer(names(df)[-1], names(df)[-1], paste)
combs <- sapply(MAT[lower.tri(MAT)], function(x) strsplit(x, " "))
ind <- lapply(combs, function(x) match(x, names(df)))

plotter <- function(cn) { #start junky function
    NAMES <- colnames(df)[cn]
    df2 <- df[cn]
    names(df2)<- c('x1', 'x2')
    p <- ggplot(data=df2, aes(x1, x2)) + geom_point() + theme_bw() +
        scale_x_continuous(name=NAMES[1]) +
        scale_y_continuous(name=NAMES[2])
        dev.new()
        print(p)
} #end of junky function

lapply(ind,  function(x) plotter(cn=x))

编辑:这更好一些:
x <- match(names(df)[-1], names(df))
MAT <- outer(x, x, paste)
combs <- t(sapply(MAT[lower.tri(MAT)], function(x) as.numeric(unlist(strsplit(x, " ")))))

plotter <- function(cn) {
    NAMES <- colnames(df)[cn]
    df2 <- df[cn]
    names(df2)<- c('x1', 'x2')
    p <- ggplot(data=df2, aes(x1, x2)) + geom_point() + theme_bw() +
        scale_x_continuous(name=NAMES[1]) +
        scale_y_continuous(name=NAMES[2])
        dev.new()
        print(p)
}

apply(combs,  1, function(x) plotter(cn=x))

关于r - ggplot2:散点图用于所有可能的变量组合,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/10239497/

10-12 19:58