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