我会在这里尝试（mvnormtest，mvoutlier，mvtnorm）： http ://cran.r-project.org/web/views/Multivariate.html

你为什么不试试我的 cepp 包。它可以为您计算 PC 的分位数！

http://cran.r-project.org/web/packages/cepp/index.html

你需要这个功能$evaluator$在这个包中。它有两个参数，对应于数据中的行数和列数。它返回一个可以优化以获得空间分位数的函数。

有关最小示例，请参阅包文档中的此片段。

x <- rnorm(500)
dim(x) <- c(250,2)
ev <- evaluator(250,2)
##The Spatial Median
trust(ev, parinit=c(median(x[1,]), median(x[2,])), u=c(0,0),
       rinit=0.5, rmax=2e5, samp = x)
##Quantile for vector (0.2,0.3)
trust(ev, parinit=c(median(x[1,]), median(x[2,])), u=c(0.2,0.3),
       rinit=0.5, rmax=2e5, samp = x)

您需要计算联合分布 f(x1,x2,...,xn) 的 CDF。CDF 将为您提供您寻求的分位数。