r - 使用 Hotelling 的 T2 作为置信区间的 R PCA 图

Question

prcomp对于我的测试数据，我使用该函数计算了前两个主成分。然后我使用 plot 函数绘制了前 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))

我现在想使用 Hotelling 的 T2 统计量在表示 95% 置信区间的数据周围放置一个椭圆。

包SensoMineR特别是panellipse功能似乎是我正在寻找的。但是，我的 R 知识仍然非常有限，因此我想知道您是否可以为我指出正确的方向来构建这样的置信区间。

非常感谢！杰基

Answer 1