用R语言计算Hotelling T2 检验T2 和 f 统计量对应的P值

在 R 语言中，我们可以使用 `Hotelling.test()` 函数来进行 Hotelling T2 检验，并计算 T2 和 f 统计量对应的 P 值。下面是一个示例代码：

# 生成两个样本数据
set.seed(123)
x1 <- matrix(rnorm(100, mean = 5, sd = 1), ncol = 5)
x2 <- matrix(rnorm(100, mean = 6, sd = 1), ncol = 5)

# 进行 Hotelling T2 检验
library(multtest)
res <- Hotelling.test(x1, x2)

# 输出 T2 和 f 统计量，以及对应的 P 值
cat("T2 statistic:", res$T2, "\n")
cat("f statistic:", res$f, "\n")
cat("P-value (T2):", res$p.value.T2, "\n")
cat("P-value (f):", res$p.value.f, "\n")

在上面的示例代码中，我们生成了两个样本数据 `x1` 和 `x2`，然后使用 `Hotelling.test()` 函数进行 Hotelling T2 检验。最后输出了 T2 和 f 统计量，以及对应的 P 值。注意，这里我们使用了 `multtest` 包中的 `Hotelling.test()` 函数，该函数还可以进行多组样本的 Hotelling T2 检验。