也许这将有所帮助 – 它是从我写的一段时间后的代码改编而成。它现在也绘制箭头。
PCbiplot
# PC being a prcomp object
data
plot
plot
datapc
mult
(max(data[,y]) - min(data[,y])/(max(datapc[,y])-min(datapc[,y]))),
(max(data[,x]) - min(data[,x])/(max(datapc[,x])-min(datapc[,x])))
)
datapc
v1 = .7 * mult * (get(x)),
v2 = .7 * mult * (get(y))
)
plot
plot
plot
}
fit
PCbiplot(fit)
您可能想要改变文字的大小,以及透明度和颜色来品尝;这将很容易使它们的功能的参数。注意:对我来说,这适用于prcomp,但您的示例是与princomp。您可能还需要相应地调整代码。注2:从注释到OP链接的邮件列表中借用geom_segment()的代码。