R 获取多元数据分位数的代码

我有一个由365个数据值组成的向量，每个值都遵循其自身的正态分布。我有一个平均向量和一个标准偏差的对角矩阵。我想从这些数据中得到分位数（0.1,0.2，…0.9）。我使用了mvtnorm包中的qmvnorm，但没有成功。我想知道是否有其他方法，也许使用循环

编辑说明：粘贴来自OP中随后删除的其他材料

对不起。为了澄清，我有一个平均值向量（总共365个值）和365x365个标准偏差矩阵。我想计算每个平均值的分位数。我使用

qmvnorm（0.1，mean=sun.base，sigma=sun.std）

得到第10个分位数。我希望结果是365个值的向量。对于我的数据，它返回如下，这不是我需要的输出

$quantile [1] 9.634588

$f.quantile [1] -6.997342e-06 attr(,"error") [1] 0 attr(,"msg") [1] "Normal Completion"

$iter [1] 10

$estim.prec [1] 6.103516e-05

---

什么不起作用，以什么方式失败？请看，这样我们可以重新创建问题。您想得到什么的分位数？数据？平均数是多少？SD？我支持一个可复制示例的请求。也许这会有所帮助：