R PCA图使用Hotelling的T2表示置信区间

Question

对于我的测试数据，我使用prcomp函数计算了前两个主成分。然后，我使用绘图函数绘制了前两个组件。

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’