私のテスト データでは、関数を使用して最初の 2 つの主成分を計算しましたprcomp
。次に、プロット関数を使用して最初の 2 つのコンポーネントをプロットしました。
df <- data.frame(x1 = runif(10), x2 = runif(10), x3 = runif(10), x4 = runif(10))
rownames(df) <- letters[1:10]
res <- prcomp(df, center = TRUE, scale = TRUE, na.action = na.omit)
plot(res$x, xlim = c(-5, 5), ylim = c(-5, 5))
ここで、ホテリングの T2 統計を使用して 95% の信頼区間を表すデータの周りに楕円を配置したいと思います。
SensoMineR
パッケージ、特にpanellipse
機能は、私が探しているもののようです。ただし、私の R の知識はまだかなり限られているため、そのような信頼区間を構築するために正しい方向に向けていただけないでしょうか。
どうもありがとう!ジャッキー