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

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

SensoMineR 很好。

还有一个包叫Package ‘Hotelling’这很容易在数据周围放置一个椭圆。

Package ‘Hotelling’