3

我正在创建一个带有多变量数据的 PCA 双图。

有没有办法指定线段的颜色/透明度/位置ggbiplot？此命令的所有参数均不提供此选项。

我知道ggbiplot是基于ggplot- 它可能接受aes论点吗？或者可以在创建的图上使用一层颜色/透明度/位置来覆盖默认值？

具体来说，关于位置，如果可能的话，我想抖动线段（尽管使它们更透明可能已经解决了问题）。

工作示例，使用iris数据：

#load required packages
library(ggplot2)
library(devtools)
library(ggbiplot)

#load dataset
data(iris)

#perform principal component analysis
pca = prcomp(iris[ , 1:4], scale=T)

#define classes, generate & view PCA biplot
class = iris$Species
ggbiplot(pca, obs.scale = 1, var.scale = 1, groups = class, circle = FALSE, 
         varname.size = 1, varname.adjust = 6)

非常感谢 - 任何帮助表示赞赏！

亲切的问候。

4

2 回答 2

7

看来您需要ggbiplot稍微更改功能。在控制台中输入ggbiplot，将代码复制到编辑器。在arglist中function，为箭头添加颜色、线型和透明度（“alpha”）的“name = expression”术语。

ggbiplot2 <- function (pcobj, choices = 1:2, scale = 1, pc.biplot = TRUE, 
          obs.scale = 1 - scale, var.scale = scale, groups = NULL, 
          ellipse = FALSE, ellipse.prob = 0.68, labels = NULL, labels.size = 3, 
          alpha = 1, var.axes = TRUE, circle = FALSE, circle.prob = 0.69, 
          varname.size = 3, varname.adjust = 1.5, varname.abbrev = FALSE, 
          color = muted("red"), # <- add new arguments to the function
          linetype = "solid",
          alpha_arrow = 1)

然后搜索该部分，并为和geom_segment添加参数：colorlinetypealpha

g <- g + geom_segment(data = df.v, aes(x = 0, y = 0, xend = xvar, yend = yvar),
                      arrow = arrow(length = unit(1/2, "picas")),
                      color = color, linetype = linetype, alpha = alpha_arrow)

将已编辑的功能分配给一个新名称，例如ggbiplot2。试试看，在其中设置箭头的默认值以外的值：

ggbiplot2(pca, obs.scale = 1, var.scale = 1, groups = class, circle = F, varname.size = 1, varname.adjust = 6,
color = "blue", linetype = "dashed", alpha_arrow = 0.5) # <- use the new arguments

在此处输入图像描述

于 2014-09-23T14:27:53.450 回答
2

谢谢你的建议。

我写了这个函数，以便它可以很容易地使用。您还可以指定要选择的列（=变量）以及分配给每个选定变量的颜色。到目前为止，它仅基于 princomp()，因为我还没有使用 PCA 的其他函数。随意贡献：

https://github.com/EhrmannS/r-snippets/blob/master/graphics/ggplot/ggbiplot_more_graphics_options.R

只需指定：

g <- ggbiplot2(pcobj, 
               coi <- list(c("variables", "with", "first", "colour"), 
                           c("variables", "with", "second", "colour")), 
               arrow.alpha = c(0.2, 1, 1),
               arrow.color = c(muted("red"), "black", "red"))

连同你的其他论点。

于 2015-06-11T14:49:51.687 回答