library(tidyverse) dat <- iris %>% select_at(1:4) %>% prcomp(scale = T) biplot(dat, xpd = T) これ自体はいいんですが。 時々、この図の、例えばSepal.Widthの赤い矢印を指して、「平面中のプロットにおいて、このパラメータは、この矢印の方向に行くに従って大きくなります(なると期待されます)。」と解説する方がいらっしゃいます。 これは、アウトです。 確かめてみましょう。 Q. 「この平面内で、Sepal.Widthが大きくなって行く(と期待される)方向(線)」を求めたい。 A. Sepal.Widthに対して、主成分PC1とPC2を重回帰します。 # パラメータの整理 y <- iris$Sepal.Width %>% scale(center = mean(.), scale = sd(.))