균일하게 분포 된 다변량 변수 사용
Taeke는 구체적으로 2-norm 및 1-norm 사례를 설명함으로써 아래 텍스트가보다 직관적 인 기사에 대한 링크를 제공합니다.
2- 노름 ∥x∥2≤r
샘플 방향
이 결과를 사용할 수 있습니다 http://mathworld.wolfram.com/HyperspherePointPicking.html
다변량 가우시안 분포 변수 (ID 공분산 행렬 포함)는 거리 또는 제곱합에만 의존합니다.X
f(X1,X2,...,Xn)=∏1≤i≤n12π−−√e12x2i=12π−−√e12∑1≤i≤nx2i
따라서 는 n- 차원 대구면에 균일하게 분포되어있다.X∥X∥2
샘플 거리
완료하려면 구의 균일 분포를 공의 균일 분포로 변경하기 위해 거리 만 샘플링하면됩니다. (디스크 포인트 선택에 대한 링크 된 예와 다소 비슷합니다)
단순히 균일 한 분포로 을 샘플링 하는 경우 중심 근처에서 상대적으로 밀도가 높아집니다 (볼륨은 r n으로 스케일 되므로 점 의 일부 r 은 볼륨 r n 에 도달합니다. 중심이며 균일 한 분포를 의미하지 않습니다)rrnrrn
대신 균일 분포에서 샘플링 된 변수 의 번째 근 을 사용하면 균일 한 분포를 얻게됩니다.n
1- 노름 ∥x∥1≤r
방향
이 경우 가우스 분포 대신 라플라스 분포에서 를 샘플링 하고 1- 노름으로 나눕니다. XX 은 n- 차원 1- 노름 구체에 균일하게 분포된다.X|X|1
공식적인 증거는없고 직감 만 있습니다
(pdf가 위치와 독립적이므로, 1- 노름이 동일한 무한대 영역 / 볼륨은 같은 확률 를 가질 것이며, 이것을 단위 표면으로 접을 때 동일한 f ( x )를 기대합니다 d Af(x)dVf(x)dA )
그러나 시뮬레이션 테스트는 좋아 보인다.
library(rmutil)
x <- abs(rlaplace(20000))
y <- abs(rlaplace(20000))
z <- abs(rlaplace(20000))
rn <- abs(x)+abs(y)+abs(z)
xi <- (x/rn)
yi <- (y/rn)
zi <- (z/rn)
plot(sqrt(0.5)*(xi-yi),
sqrt((0.5-0.5*(xi+yi))^2+zi^2),
pc=21,bg=rgb(0,0,0,0.02), col=rgb(0,0,0,0),cex=1)
거리
거리는 2- 노름 경우와 유사합니다 (볼륨은 여전히 rn ).
p- 노름 ∥x∥p≤r
이 같은 원리를 따라하고자하는 경우이 경우, 당신이 분포에서 샘플링 할 필요가 (가설). 이것들은 일반 정규 분포 이며 아마도 분포 G 와 관련이 있습니다.f(x)∝e|x|pG()