我有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