使用Hotelling的T2进行置信区间的R PCA图

Question

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

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'的软件包 ，可以很容易地在数据周围放置一个椭圆。

套餐'Hotelling'