タグ付けされた質問 「multivariate-normal」

1
結合信頼区間を計算するためのガウス相関不等式の結果
Quanta Magazineのこの非常に興味深い記事によると、「長い間求められていた証拠、発見され、ほとんど失われた」、- 多変量ガウス分布を持つベクトルが与えられたことが証明されました。そして間隔所与I 1、... 、Iはn個の対応する構成要素の手段を中心Xを、次いで、x=(x1,…,xn)x=(x1,…,xn)\mathbf{x}=(x_1,\dots,x_n)I1,…,InI1,…,InI_1,\dots,I_n xx\mathbf{x} p(x1∈I1,…,xn∈In)≥∏i=1np(xi∈Ii)p(x1∈I1,…,xn∈In)≥∏i=1np(xi∈Ii)p(x_1\in I_1, \dots, x_n\in I_n)\geq \prod_{i=1}^n p(x_i\in I_i) (ガウス相関不等式またはGCI。より一般的な定式化については、https: //arxiv.org/pdf/1512.08776.pdfを参照してください)。 これは本当に素晴らしく簡単に思えますが、記事は、それが共同信頼区間に結果をもたらすと述べています。しかし、それに関しては私にはまったく役に立たないようです。我々はパラメータ推定されていると仮定 、我々は推定した^ θ 1、... 、^ θ n個ある(多分漸近的に)共同ノーマル(例えば、MLE推定)。次に、各パラメーターの95%信頼区間を計算すると、GCIはハイパーキューブI 1 × … I nが(θ1,…,θnθ1,…,θn\theta_1,\dots,\theta_nθ1^,…,θn^θ1^,…,θn^\hat{\theta_1},\dots,\hat{\theta_n}I1×…InI1×…InI_1\times\dots I_n ...これは、適度な nでもかなり低いカバレッジです。(0.95)n(0.95)n(0.95)^n nnn したがって、共信頼領域を見つける賢い方法ではないようです。多変量ガウス、つまり超楕円体の通常の信頼領域は、共分散行列が既知で、よりシャープであるかどうかを見つけるのは難しくありません。共分散行列が不明な場合に信頼領域を見つけることが役立つかもしれませんか?GCIと共同信頼領域の計算との関連性の例を教えてください。

4
多変量正規分布の分位数(アイソライン?)を決定する方法
多変量分布の分位数を計算する方法に興味があります。図では、特定の単変量正規分布の5%および95%の分位点を描画しました(左)。適切な多変量正規分布の場合、アナログは密度関数の基底を囲む等値線になると想像しています。以下は、パッケージを使用してこれを計算する試みの例ですが、mvtnorm成功しません。多変量密度関数の結果の等高線を計算することでこれを行うことができると思いますが、別の選択肢(たとえばの類似体qnorm)があるかどうか疑問に思っていました。ご協力いただきありがとうございます。 例: mu <- 5 sigma <- 2 vals <- seq(-2,12,,100) ds <- dnorm(vals, mean=mu, sd=sigma) plot(vals, ds, t="l") qs <- qnorm(c(0.05, 0.95), mean=mu, sd=sigma) abline(v=qs, col=2, lty=2) #install.packages("mvtnorm") require(mvtnorm) n <- 2 mmu <- rep(mu, n) msigma <- rep(sigma, n) mcov <- diag(msigma^2) mvals <- expand.grid(seq(-2,12,,100), seq(-2,12,,100)) mvds <- …

2
最尤推定量-多変量ガウス
環境 多変量ガウス分布は機械学習で頻繁に使用され、次の結果は多くのMLブックおよび派生物なしのコースで使用されます。 次元行列の 形式のデータが与えられ、データが 平均()および共分散行列(変量ガウス分布に従うと仮定した場合)最尤推定量は次によって与えられます:XX\mathbf{X} m×pm×p m \times ppppμμ\mup×1p×1p \times 1 ΣΣ\Sigmap×pp×pp \times p μ^=1m∑mi=1x(i)=x¯μ^=1m∑i=1mx(i)=x¯\hat \mu = \frac{1}{m} \sum_{i=1}^m \mathbf{ x^{(i)} } = \mathbf{\bar{x}} Σ^=1m∑mi=1(x(i)−μ^)(x(i)−μ^)TΣ^=1m∑i=1m(x(i)−μ^)(x(i)−μ^)T\hat \Sigma = \frac{1}{m} \sum_{i=1}^m \mathbf{(x^{(i)} - \hat \mu) (x^{(i)} -\hat \mu)}^T 多変量ガウスの知識は多くのMLコースの前提条件であることを理解していますが、多くの自己学習者が統計を跳ね回っていると感じているので、自己完結型の回答に完全に由来することが役立つと思います。 stackexchangeおよびmath.stackexchange Webサイトで回答を探しています。 質問 多変量ガウスの最尤推定量の完全な導出は何ですか 例: これらの線形判別分析の講義ノート(11ページ)、またはこれらのものは結果を利用すると、以前の知識を前提としています。 また、部分的に回答またはクローズされている投稿もいくつかあります。 多変量正規分布の最尤推定器 多変量正規分布の最尤推定を理解するのに助けが必要ですか?

2
共同分布が多変量正規分布である場合、ピアソンのρは関連性の網羅的な尺度にすぎないのはなぜですか?
この主張はこの質問への一番の回答で提起されました。「なぜ」という質問は、新しいスレッドを保証するほど十分に異なると思います。グーグルの「関連性の徹底的な尺度」はヒットを生み出さず、そのフレーズが何を意味するのか分かりません。


1
多変量ガウス分布から値を生成する
現在、平均ベクトルと共分散行列多変量正規分布を持つ次元のランダム変数値をシミュレートしようとしています。X μ = (μ 1、。。。、μ N )T SNNNXXXμ=(μ1,...,μN)Tμ=(μ1,...,μN)T\mu = (\mu_1,...,\mu_N)^TSSS 逆CDF法に似た手順を使用したいと考えています。つまり、最初にNNN次元の一様なランダム変数UUUを生成し、次にこの分布の逆CDFにプラグインして、値Xを生成しますXXX。 手順が十分に文書化されておらず、MATLABのmvnrnd関数とウィキペディアで見つけた説明にわずかな違いがあるため、問題が発生しています。 私の場合、分布のパラメーターもランダムに選択しています。特に、平均分布μiμi\mu_i一様分布U(20,40)U(20,40)U(20,40)ます。次に、次の手順を使用して共分散行列SSSを作成します。 下三角行列作成L L(I、I)= 1のための私は1..N =及びLを(I、J)= U(-1,1)のために 、I &lt;JLLLL(i,i)=1L(i,i)=1L(i,i) = 1i=1..Ni=1..Ni=1..NL(i,j)=U(−1,1)L(i,j)=U(−1,1)L(i,j) = U(-1,1)i&lt;ji&lt;ji < j してみましょうS = LL ^ T L ^ Tはの転置表すLを。S=LLTS=LLTS = LL^TLTLTL^TLLL この手順により、SSSが対称かつ正定であることを確認できます。また、S = LL ^ Tになるように下三角行列Lを提供します。これは、分布から値を生成するために必要だと思います。LLLS=LLTS=LLTS = LL^T ウィキペディアのガイドラインを使用すると、次のようにN次元のユニフォームを使用してXの値を生成できるはずです。XXXNNN X=μ+L∗Φ−1(U)X=μ+L∗Φ−1(U)X = \mu + L * …

2
共同正規性は、正常なランダム変数の合計が正常であるための必要条件ですか?
関連する質問に対する私のこの回答に続くコメントで、ユーザーssdecontrolとGlen_b は、合計正規性を主張するためにと共同正規性が必要かどうかを尋ねました。ジョイントの正規性が十分であることは、もちろんよく知られています。この補足的な質問はそこでは取り上げられておらず、おそらくそれ自体で検討する価値があります。XXXYYYX+YX+YX+Y 共同正規性は限界正規性を意味するので、私は尋ねます が通常のランダム変数であるが、とが 一緒に通常のランダム変数ではないような 通常のランダム変数とが存在しますか?XXXYYYX+YX+YX+YXXXYYY 場合はと正規分布を持つ必要はありません、正常な確率変数を簡単に見つけることができます。1つの例は、以前の回答にあります(リンクは上記のとおりです)。上記のハイライトされた質問に対する答えは「はい」であると信じており、この質問に対する答えとして例を(私が思うに)掲載しています。XXXYYY



2
通常、標準シンプレックスからのサンプリング
与えられた平均と共分散行列を使用して、合計が1になるように範囲に切り捨てられた次元多変量ガウス分布から値を生成できるようにしたいと考えています。nnn[0,1][0,1][0, 1] これはガウス分布による標準 -simplex からのサンプリングと同じだと思いますが、どうすればこれを実行できますか?んnn
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.