我有一个包含 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