Rで解析：ggplot2でPCA結果をBiplotで表示「ggbiplot」パッケージ

ggplot2パッケージを利用して主成分分析の結果をBiplotで表示するパッケージの紹介です。

パッケージバージョンは0.55。実行コマンドはR version 4.2.2で確認しています。

＜おすすめのRに関する書籍です＞

言葉と数式で理解する多変量解析入門 | 小杉 考司

Amazonで小杉 考司の言葉と数式で理解する多変量解析入門。アマゾンならポイント還元本が多数。小杉 考司作品ほか、お急ぎ便対象商品は当日お届けも可能。

www.amazon.co.jp

パッケージのインストール

下記コマンドを実行してください。

#パッケージのインストール
install.packages("devtools")
devtools::install_github("vqv/ggbiplot")

実行コマンド

詳細はコマンド、パッケージのヘルプを確認してください。

#パッケージの読込
library("ggbiplot")

###データ例の作成#####
set.seed(1234)
TestData <- matrix(rnorm(1000), 200)
TestData <- as.data.frame(TestData)
#データの調整
TestData[1:67,] <- TestData[1:67,] + 1
TestData[68:135,] <- TestData[68:135,] + 2
TestData[136:200,] <- TestData[136:200,] + 3
#グループ情報を付与
TestData <- cbind(c(rep("Group1", 67),
                    rep("Group2", 67),
                    rep("Group3", 66)),
                  TestData)
colnames(TestData) <- c("Group", paste0("ColName", seq(5)))
########

#statパッケージ:prcompコマンドで主成分分析
TestPrc <- prcomp(TestData[, 2:6], scale. = FALSE)

#biplotの作成:ggbiplotコマンド
#prcomp()またはprincomp()の結果を指定:pcobjオプション
#プロットする主成分を指定:choicesオプション
#グループ情報を指定:groupsオプション
#グループ毎の確率楕円を描写:ellipseオプション
#相関円を描写:circleオプション
#ggplot2のコマンドが利用可能です
ggbiplot(pcobj = TestPrc, choices = 1:2, obs.scale = 1, var.scale = 1,
         groups = TestData[, 1], ellipse = TRUE, circle = TRUE) +
  scale_colour_manual(values = c("#FF0000", "black", "#00FF00")) +
  theme(legend.direction = "horizontal", legend.position = "top")

出力例

・ggbiplotコマンド

少しでも、あなたの解析が楽になりますように！！

＜おすすめのRに関する書籍です＞

RとRStudioによる教育テストデータの分析 | 堀 一輝, 福原 弘岳, 山田 剛史

Amazonで堀 一輝, 福原 弘岳, 山田 剛史のRとRStudioによる教育テストデータの分析。アマゾンならポイント還元本が多数。

www.amazon.co.jp