ガウス過程とウィシャート分布の共分散行列


11

一般化されたウィシャートプロセス(GWP)に関するこのペーパーを読んでいます。この論文では、2乗指数共分散関数、つまりを使用して、さまざまな確率変数(ガウスプロセスに従って)間の共分散を計算します。次に、この共分散行列はGWPに従います。K(x,x)=exp(|(xx)|22l2)

私は、線形共分散関数()K(x,x)=xTxから計算された共分散行列は、適切なパラメーターを使用してウィシャート分布に従うと考えていました。

私の質問は、二乗指数共分散関数を使用してウィシャート分布に従う共分散をまだどのように仮定できるかです。また、一般的に、Wishart分散共分散行列を生成するための共分散関数に必要な条件は何ですか?

回答:


8

混同されているのは、ガウス過程が定義されている周囲空間に関する共分散仕様と、有限次元のガウス確率変数を変換してウィシャート分布を生成する演算です。

もしあるとの次元のガウス確率変数(列ベクトル)は、平均0及び共分散行列、の分布は、ウィシャート分布です。は行列であることに注意してください。これは、2次形式の がガウス分布をウィシャート分布に変換する方法に関する一般的な結果 です。これは、任意の正定共分散行列に当てはまります。iid観測がある場合は、XN(0,Σ)pΣW=XXTWp(Σ,1)Wp×p

xxxT
ΣX1,,Xn次に、て、 はウィシャート分布。割る、我々は経験的共分散行列を取得の推定値。Wi=XiXiT
W1++Wn
Wp(Σ,n)nΣ

ガウスプロセスの場合、周囲空間が存在します。例として、であるとしましょう。考慮される確率変数は、周囲空間の要素によってインデックスが付けられます。つまり、プロセスを検討します。有限次元の周辺分布がガウスの場合、つまり for all。OPで言及されているように、共分散関数の選択により、共分散行列、つまり R(X(x))xR

X(x1,,xp):=(X(x1),,X(xp))TN(0,Σ(x1,,xp))
x1,,xpR
cov(X(xi),X(xj))=Σ(x1,,xp)i,j=K(xi,xj).
の選択を無視すると、 はWishart -distribution。K
X(x1,,xp)X(x1,,xp)T
Wp(Σ(x1,,xp),1)

これに答えてくれてありがとう。いくつか質問があります。あなたの答え-Gaussian distをWishart distに変換する変換が+ veの明確なcov行列の任意の選択に対して成立すると言うとき、このcov行列にはどのようなdiff選択肢がありますか?また、明確にするために-cov関数で定義されたcovマトリックスについて、iとjはガウスプロセスの周囲空間の要素を示しています(たとえば、時間プロセスの場合、時刻t_1およびt_2)。
ステディフィッシュ2013

@steadyfish、はい、インデックスとは、環境空間内の点とを参照し、時間的なプロセスでは2つの時点を参照します。共分散行列は常に正の(半)確定です。製剤は、どのような方法で結果を制限するのではなく、それは、任意の選択のために保持していることを強調することを意図していなかった限り共分散行列です。私がいる可能性取り残さ単数正規分布などの無関係な問題に答えを乱雑に避けるためにちょうど半正定値ことができますijxixjΣ ΣΣ
NRH

@NRHに感謝します。アンビエントスペースについてポイントを取得します。共分散行列について、私の質問は、以外の共分散行列を定義する方法が他にあるかどうかでした(正定性または正の半定特性についてではありません)。(今回は問題が明確であることを願っています!)xTx
ステディフィッシュ

@steadyfish、ああ、なるほど。実際、転置や、ベクトルが行ベクトルであるか列ベクトルであるかについては、ずさんでした。これを正確にして、経験的共分散行列と理論的共分散行列の関係について少し追加しました。理論は観測値の観点から定義されていません。
NRH 2013
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.