1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
| library(tidyverse)
data("iris")
pca = prcomp(iris[,1:4])
score = pca[["x"]] %>% as.data.frame()
loading = pca[["rotation"]] %>% as.data.frame()
ggplot(score, aes(PC1, PC2)) +
geom_hline(yintercept = 0, linetype = 'dashed') +
geom_vline(xintercept = 0, linetype = 'dashed') +
geom_point(aes(color = iris$Species)) +
geom_segment(data = loading,
aes(x = PC1, y = PC2, xend = 0, yend = 0),
arrow = arrow(length=unit(0.20,"cm"),
ends="first", type = "closed"),
size = 0.5) +
geom_text(data = loading,
aes(x = PC1, y = PC2, label = rownames(loading)),
size = 3) +
theme_bw()
|
