多変量ガウスのコンターのハイパーボリューム


8

から抽出されたサイズサンプルで、原点までのユークリッド距離が最小の観測値の%の(行列式の対数)の共分散の()値の値を 探しています、二変量標準ガウス。α Nnαn

-楕円のハイパーボリュームは、その共分散行列の行列式に比例するため、タイトルに比例します-

標準変量ガウス--By、Iは平均長さ2の0のベクトルであり、ランク2単位行列であるが.--- 0 2 IN2(02,II2)02II2

、数値が前後の 場合よりも、シミュレーションで簡単に確認でき。α=52/701.28

library(MASS)
n<-10000
p<-2
x<-mvrnorm(n,rep(0,p),diag(2))
h<-ceiling(0.714286*n)
p<-ncol(x)
w<-mahalanobis(x,rep(0,p),diag(p),inverted=TRUE) #These are eucledian distances, because the covariance used is the identity matrix
s<-(1:n)[order(w)][1:h]
log(det(cov(x[s,])))

しかし、これについて正確な式を取得する方法(または失敗した場合、より良い近似)を思い出しません。


1
あなたのテキストでは、二変量分布のパラメーターについては何も言われていません。また、ユークリッドdではなく、マハラノビスdに関するコードのようです。
ttnphns 2014

1
標準ガウスとは、原点を中心とし、恒等共分散を持つものを意味します(これは編集します)。アイデンティティ共分散行列に対するマハラノビス距離==ユークリッド距離。
user603 14

1
コードを使用している場合、またはコードのヘルプを求めている場合は、使用している言語またはプログラムを明記してください。
ウルフィー14

回答:


7

さて、この質問は時々出てくるようですので、一概にお答えしたいと思います。

[1]で、著者は、 with symmetric positive definite、andΣ S αxxiNp(μμ,ΣΣ),i=1,,nΣSα

(0)Sα={i:(xxiμμ)Σ1(xxiμμ)qα}

以下のための及びqα=χp2(α),0<α1

(1)Cα=coviSαxxi

次に、漸近的に、はに収束し。CαlαΣ

(2)lα=Fχp+22(qα)α

この近似は本当に良いです(ここではalpha = 60/70の場合):

library(MASS)
alpha<-60/70
p<-2
n<-1000000

radius<-sqrt(qchisq(alpha,df=p))
x0<-mvrnorm(n,rep(0,p),diag(p),empirical=TRUE)
Id<-which(rowSums(x0*x0)<=radius**2)
cov(x0[Id,])

qalpa<-qchisq(alpha,p)
diag(1/(alpha/(pchisq(qalpa,p+2))),p)

したがって、最後に、質問に答えるために、起点に対する最小ユークリッドノルムを持つ観測の共分散行列の行列式(これは、および)は、log[αn]Σ=IIpμμ=00p

(3)plogFχp+22(qα)plogα
  1. Croux C.、Haesbroeck G.(1999)。最小共分散行列式散布行列推定器の影響関数と効率。多変量解析のジャーナル。71. 161--190。
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.