원의 확률 분포
확률 이론 과 방향 통계 에서 폰 미제스 분포 (순환 정규 분포 또는 티코노프 분포 라고도 한다)는 원 의 연속 확률 분포 다.이것은 정상 분포의 원형 아날로그인 포장된 정규 분포에 가까운 근사치 입니다.원 위의 자유롭게 확산되는 각도 {{\displaystyle \theta }는 정규 분포를 따르는 랜덤 변수로서, 시간에 따라 선형으로 성장하는 포장 되지 않은 분산이다. 한편, 폰 미제스 분포는 조화 전위, 즉 선호하는 방향을 갖는 원상의 표류 및 확산 과정을 고정적으로 분포하는 것이다.[1] 폰 미제스 분포는 첫 번째 원형 모멘트 의 실제 부분과 가상 부분을 지정할 때 원형 데이터의 최대 엔트로피 분포 다. 폰 미제스 분포는 N차원 구체에 폰 미제스-피셔 분포 의 특수한 경우다. null
정의 각도 x 에 대한 폰 미제스 확률 밀도 함수는 다음과 같이 주어진다.[2]
f ( x ∣ μ , κ ) = e κ cas ( x − μ ) 2 π I 0 ( κ ) {\displaystyle f(x\mid \mu ,\kappa )={\frac {e^{\kappa \cos(x-\mu )}}{2\pi I_{0}(\kappa )}}}}}}} 여기서 I 0 ( κ {\displaystyle \kappa } )는 순서 0의 수정된 베셀 함수 다.
매개변수 μ 및 1/ μ {\displaystyle \kappa } 은 정규 분포에서 μ 및 μ 2 (평균 및 분산)와 유사하다 .
μ 는 위치의 척도(분포가 μ 주위에 군집화되어 있음)이며, κ {\displaystyle \kappa } 은 (는) 농도의 척도(분산 의 역수 척도이므로 1/ κ {\displaystyle \kappa } 은 σ 과2 유사하다.null κ {\displaystyle \kappa } 이 (가) 0이면 분포가 균일하고, 작은 κ {\displaystyle \kappa } 의 경우 균일성에 가깝다. κ {\displaystyle \kappa } 이 (가) 크면 분포가 각도 μ 에 매우 집중되며, with {\displaystyle \kappa } 은 농도의 측정값이다 .실제로 κ {\displaystyle \kappa } 이(가) 증가함에 따라 분포는 평균 μs 및 분산 1/2 κ {\displaystyle \kappa } 의 정규 분포에 접근한다. 확률밀도는 일련의 Besel 함수로[3] 표현될 수 있다.
f ( x ∣ μ , κ ) = 1 2 π ( 1 + 2 I 0 ( κ ) ∑ j = 1 ∞ I j ( κ ) cas [ j ( x − μ ) ] ) {\displaystyle f(x\mid \mu ,\kappa )={\frac {1}{2\pi }}\왼쪽(1+{\frac {2}){ I_{0}(\kappa )}}\sum _{j=1}^{\\infit }I_{j}(\kappa )\cos[j(x-\mu )]\right)}} 여기서 I j (x )는 순서 j 의 수정된 베셀 함수 다.
누적분포함수는 분석적이지 않으며 위의 시리즈를 통합하여 가장 잘 찾을 수 있다. 확률밀도의 무한정 적분은 다음과 같다.
Φ ( x ∣ μ , κ ) = ∫ f ( t ∣ μ , κ ) d t = 1 2 π ( x + 2 I 0 ( κ ) ∑ j = 1 ∞ I j ( κ ) 죄를 짓다 [ j ( x − μ ) ] j ) . {\displaystyle \Phi(x\mid \mu ,\kappa )=\int f(t\mid \mid \mu ,\kappa )\\,dt={\frac {1}{1}{2\pi }}}\flt(x+{\frac {2}{2}{2}{}{}{}}}{}{}}}}}}}{}{}{}{}{} I_{0}(\kappa )}}}\sum _{j=1}^{\\inflat }{{j}(\kappa ){\sin[j(x-\mu )]}{j}\right). } 누적 분포 함수는 통합 x:의0 하한값의 함수가 될 것이다.
F ( x ∣ μ , κ ) = Φ ( x ∣ μ , κ ) − Φ ( x 0 ∣ μ , κ ) . {\displaystyle F(x\mid \mu ,\kappa )=\Phi(x\mid \mood \mu ,\kappa )-\Phi(x_{0}\mid \mid \mu ,\kappa \,},},} 순간 폰 미제스 분포의 모멘트는 대개 각도 x 자체보다는 복합 지수 z = e 의ix 모멘트로 계산된다. 이러한 순간들을 순환적 인 순간이라고 한다. 이러한 순간에서 계산된 분산을 순환 분산 이라고 한다. 이에 대한 한 가지 예외는 "평균"이 보통 복합 평균의 주장 을 가리킨다는 점이다. null
z 의 n번째 원시 모멘트는 다음과 같다.
m n = ⟨ z n ⟩ = ∫ Γ z n f ( x μ , κ ) d x {\displaystyle m_{n}=\langle z^{n}\langle z^\\int_{\\\\Gamma }z^{n}\,f(x \mu ,\kappa )\,dx} = I n ( κ ) I 0 ( κ ) e i n μ {\displaystyle ={\frac {I_{n }(\kappa )}{ I_{0}(\kappa )}e^{in\mu }}}} 여기서 적분은 길이 2㎝의 γ{\displaystyle \Gamma } 간격에 걸쳐 있다. 위 의n 적분을 계산할 때 [4] z = cos(nx ) + i sin(nx ) 및 Besel 함수 ID:
I n ( κ ) = 1 π ∫ 0 π e κ cas ( x ) cas ( n x ) d x . {\displaystyle I_{n}(\kappa )={\frac {1}{\pi }\int_{0}^{\pi }e^{\cappa \cos(x)\cos(nx)\,dx.} 복합 지수 z 의 평균은 다음과 같다.
m 1 = I 1 ( κ ) I 0 ( κ ) e i μ {\displaystyle m_{1}={\frac {I_{1}(\kappa )}{{}}{{}}{ I_{0}(\kappa )}e^{i\mu }}} 그리고 각도 x 의 원형 평균 값은 인수 μ 로 간주된다. 이것은 각 랜덤 변수의 예상 방향 또는 선호 방향이다. z 의 분산 또는 x 의 원형 분산은 다음과 같다.
시합을 하다 ( x ) = 1 − E [ cas ( x − μ ) ] = 1 − I 1 ( κ ) I 0 ( κ ) . {\displaystyle {\textrm {var}(x)=1-E[\cos(x-\mu )]=1-{\frac {I_{1}(\kappa )}{{{}}{{}}}{{}}}}} I_{0}(\kappa )}}. }
제한행동 κ {\displaystyle \kappa } 이(가) 크면 분포가 정규 분포 와 유사하다.좀 더 구체적으로, 큰 양의 실수 ▼ {\displaystyle \kappa },
f ( x ∣ μ , κ ) ≈ 1 σ 2 π 생략하다 [ − ( x − μ ) 2 2 σ 2 ] {\displaystyle f(x\mid \mu ,\kappa )\약 {\frac {1}{\frac {1}{\sqrt {2\pi }}}}\exp \efrac[{\dfrac {-}-(x-\mu )^{2}}\}\오른쪽]}} 여기서 σ2 = 1/ κ {\displaystyle \kappa }, and {\displaystyle \kappa } 이(가) 무한대로 가면서 근사치의 왼쪽과 오른쪽의 차이가 균일 하게 0으로 수렴된다. 또한 κ {\displaystyle \kappa } 이(가) 작을 경우 확률밀도함수는 동일한 분포 와 유사하다.
임이 있는 κ → 0 f ( x ∣ μ , κ ) = U ( x ) {\displaystyle \lim _{\kappa \rightarrow 0}f(x\mid \mu ,\kappa )=\mathrm {U}(x)} where the interval for the uniform distribution U ( x ) {\displaystyle \mathrm {U} (x)} is the chosen interval of length 2 π {\displaystyle 2\pi } (i.e. U ( x ) = 1 / ( 2 π ) {\displaystyle \mathrm {U} (x)=1/(2\pi )} when x {\displaystyle x} is in the interval and U ( x ) = 0 {\displaystyle \mathrm {U} (x )= 0}( x {\displaystyle x}이(가) 간격에 없는 경우).null
모수 추정 폰 미제스 분포에서 도출된 일련 의 N 측정 z n = e i in {\ displaystyle z_{n}=e^{i\theta_{n}}} 을(를) 사용하여 분포의 특정 모수를 추정할 수 있다.[5] 영상 시리즈 z'{\ displaystyle {\overline{z}} 의 평균은 다음과 같이 정의된다 .
z ¯ = 1 N ∑ n = 1 N z n {\displaystyle {\overline{z}={\frac {1}{{N}\sum _{n=1}^{N}z_{n}}}} 그리고 그 기대 가치는 단지 첫 순간일 것이다.
⟨ z ¯ ⟩ = I 1 ( κ ) I 0 ( κ ) e i μ . {\displaystyle \langle {\overline{z}\rangele ={\frac {I_{1}(\kappa )}{{}}{}}{ I_{0}(\kappa )}e^{i\mu }} 즉, z 의 {\ displaystyle {\overline{z}} 은(는) 첫 순간 의 편견 없는 추정기 라는 것이다 . 평균 μs {\displaystyle \mu } 이(가) 간격 [ - π , π ] {\displaystyle [-\pi ,\pi ]}} 에 있다고 가정하면 Arg(z ){\displaystyle({\overline {z})} 은 평균 μs 의 (편향상)가 될 것이다.
z n {\displaystyle z_ { n }}을( 를) 복잡한 평면에서 벡터 집합으로 보면 R 의 2 {\displaystyle {\bar{R}^{2 }} 통계량은 평균 벡터 길이의 제곱이다.
R ¯ 2 = z ¯ z ∗ ¯ = ( 1 N ∑ n = 1 N cas θ n ) 2 + ( 1 N ∑ n = 1 N 죄를 짓다 θ n ) 2 {\displaystyle {\bar {R}}^{2}={\overline {z}}\,{\overline {z^{*}}}=\left({\frac {1}{N}}\sum _{n=1}^{N}\cos \theta _{n}\right)^{2}+\left({\frac {1}{N}}\sum _{n=1}^{N}\sin \theta _{n}\right)^{2}} 그리고 그 기대치는
⟨ R ¯ 2 ⟩ = 1 N + N − 1 N I 1 ( κ ) 2 I 0 ( κ ) 2 . {\displaystyle \langle{\bar{R}^{2}\angle ={\frac {1}{{N}+{\frac {{N:1}}\,{\frac {I_{1}{1}(\kappa )^{2}}:{\}{\frac {I_{1}}}}}}{2}}:{{{{{{}}}}}}}}}}}}}}}} I_{0}(\kappa )^{2}}. } 즉, 통계는
R e 2 = N N − 1 ( R ¯ 2 − 1 N ) {\displaystyle R_{e}^{2}={\frac {N}{{N:1}}\왼쪽({\bar{R}}^{{\frac {1}{N}\오른쪽)}} I 1 ( κ ) 2 I 0 ( κ ) 2 {\ displaystyle {\frac {I_{1}(\kappa )^{2}}:{{}} 의 편향되지 않은 추정기가 될 것이다.I_{0}(\kappa )^{2}}}\,} 등식을 풀고 R e = I 1 ( κ ) I 0 ( κ ) {\ displaystyle R_{e}={\frac {I_{1}(\kappa )}}{{}}}}}}}}{{}}}}}}}}}}}}}{{ I_{0}(\kappa )}}\,} for κ {\displaystyle \kappa \,} will yield a (biased) estimator of κ {\displaystyle \kappa \,} . In analogy to the linear case, the solution to the equation R ¯ = I 1 ( κ ) I 0 ( κ ) {\displaystyle {\bar {R}}={\frac {I_{1}(\kappa )}{ I_{0}(\kappa )}}\,} 은(는) 최대우도 추정치 κ{\ displaystyle \kappa \,} 을(를) 산출하며 , 큰 N의 한계 에서는 둘 다 동일할 것이다.κ {\ displaystyle \kappa \,} 에 대한 대략적인 용액은 von Mises-Fisher 배포 를 참조하십시오 .null
평균 분포 von Mises 분포에 대한 표본 평균 z' = R' e' i 'i' i' i'i'i'i'i'i'i'i'i'i' i'i'i'{\ thea'}}} 는 다음과 같다.[7]
P ( R ¯ , θ ¯ ) d R ¯ d θ ¯ = 1 ( 2 π I 0 ( κ ) ) N ∫ Γ ∏ n = 1 N ( e κ cas ( θ n − μ ) d θ n ) = e κ N R ¯ cas ( θ ¯ − μ ) I 0 ( κ ) N ( 1 ( 2 π ) N ∫ Γ ∏ n = 1 N d θ n ) {\displaystyle P({\bar{R},{\bar {\theta}}})\,d{\bar {R}\,d{\bar {\theta}}}={\frac {1}{{{{{{{pi I_{0}(\cappa )^{{{{{{{{{}}}}}}}}}}^{{{{{}}}}}}}}}}}}}})^{{{{{{}}} N}}}\int _{\Gamma }\prod _{n=1}^{N}\left(e^{\kappa \cos(\theta _{n}-\mu )}d\theta _{n}\right)={\frac {e^{\kappa N{\bar {R}}\cos({\bar {\theta }}-\mu )}}{ I_{0}(\kappa )^{{N}}}\왼쪽({\frac {1}{{{{{{{n}}}}\int_{\Gamma }\prod _{n=1}^{N}d_{n}d_{n}}\right)}} where N is the number of measurements and Γ {\displaystyle \Gamma \,} consists of intervals of 2 π {\displaystyle 2\pi } in the variables, subject to the constraint that R ¯ {\displaystyle {\bar {R}}} and θ ¯ {\displaystyle {\bar {\theta }}} are constant, where R ¯ {\displaystyle {\bar {R}}} 평균 결과:
R ¯ 2 = z ¯ 2 = ( 1 N ∑ n = 1 N cas ( θ n ) ) 2 + ( 1 N ∑ n = 1 N 죄를 짓다 ( θ n ) ) 2 {\displaystyle {\bar {R}}^{2}= {\bar {z}} ^{2}=\left({\frac {1}{N}}\sum _{n=1}^{N}\cos(\theta _{n})\right)^{2}+\left({\frac {1}{N}}\sum _{n=1}^{N}\sin(\theta _{n})\right)^{2}} 그리고 θ{\ displaystyle {\overline {\theta}} 은 (는) 평균 각도:
θ ¯ = A r g ( z ¯ ) . {\displaystyle {\theta}}=\mathrm {Arg}({\overline {z}) \,} 괄호 안의 제품 항은 원형 균일 분포 에 대한 평균의 분포일 뿐이라는 점에 유의하십시오.[7] null
This means that the distribution of the mean direction μ {\displaystyle \mu } of a von Mises distribution V M ( μ , κ ) {\displaystyle VM(\mu ,\kappa )} is a von Mises distribution V M ( μ , R ¯ N κ ) {\displaystyle VM(\mu ,{\bar {R}}N\kappa )} , or, equivalently, V M ( μ , R κ ) {\displaystyle VM(\mu ,R\kappa )} .
엔트로피 정의상 폰 미제스 분포의 정보 엔트로피는 [2]
H = − ∫ Γ f ( θ ; μ , κ ) ln ( f ( θ ; μ , κ ) ) d θ (\displaystyle H=-\int _{\Gamma }f(\theta ;\mu ,\kappa )\,\ln(\theta ;\mu ,\cappa )\,d\theta \,} 여기서 γ {\displaystyle \Gamma} 은 (는) 길이 2 π [\displaystyle 2\pi } 의 임의 간격이다. Von Mises 분포의 밀도에 대한 로그는 간단하다.
ln ( f ( θ ; μ , κ ) ) = − ln ( 2 π I 0 ( κ ) ) + κ cas ( θ ) #\displaystyle \ln(f(\theta ;\mu ,\kappa )=-\ln(2\pi I_{0}(\kappa )+\cappa \cos(\theta )\,},} Von Mises 분포의 특성 함수 표현은 다음과 같다.
f ( θ ; μ , κ ) = 1 2 π ( 1 + 2 ∑ n = 1 ∞ ϕ n cas ( n θ ) ) {\displaystyle f(\theta ;\mu ,\kappa )={\frac {1}{2\pi }}\좌측(1+2\sum _{n=1}^{n}\inflty _{n}\cos(n\tea )\right)}} 여기서 ϕn = I n ( κ ) / I 0 ( κ ) {\displaystyle \phi _{n}= I_{n }(\kappa )/I_{0}(\kappa )}. 이러한 표현을 엔트로피 적분으로 대체하고, 통합과 합산의 순서를 교환하며, 코사인의 직교성을 사용하여 엔트로피를 다음과 같이 쓸 수 있다.
H = ln ( 2 π I 0 ( κ ) ) − κ ϕ 1 = ln ( 2 π I 0 ( κ ) ) − κ I 1 ( κ ) I 0 ( κ ) {\displaystyle H=\ln(2\pi I_{0}(\kappa )-\cappa \pi \pi \{1}=\ln(2\pi I_{0})(\kappa )-\papa {\frac {I_{1}(\kappa )}{}}}}{}}}}}}}}}}}}}}}}}{{{}}}}}}}}}}}}}}}}}} I_{0}(\kappa )}} κ = 0 {\displaystyle \kappa =0} 의 경우, 폰 미제스 분포는 순환 균일 분포 가 되고 엔트로피는 최대값 ln ( 2 π ){\displaystyle \ln (2\pi )} 에 도달한다.
폰 미제스 분포는 첫 번째 원형 모멘트 의 실제와 가상의 부분이 지정되거나[8] 동등하게 원형 평균 과 원형 분산 이 지정될 때 엔트로피를 최대화 한다는 점에 유의한다. null
참고 항목 참조 ^ Risken, H. (1989). The Fokker–Planck Equation . Springer. ISBN 978-3-540-61530-9 . ^ a b Mardia, Kantilal ; Jupp, Peter E. (1999). Directional Statistics . Wiley. ISBN 978-0-471-95333-3 . ^ 아브라모위츠 및 스테건 제9.6.34조 참조 ^ 아브라모위츠 및 스테건 §9.6.19 참조 ^ Borradaile, G. J. (2003). Statistics of earth science data : their distribution in time, space, and orientation . Springer. ISBN 978-3-662-05223-5 . ^ Kutil, Rade (August 2012). "Biased and unbiased estimation of the circular mean resultant length and its variance" . Statistics: A Journal of Theoretical and Applied Statistics . 46 (4): 549–561. doi :10.1080/02331888.2010.543463} . ^ a b Jammalamadaka, S. Rao; Sengupta, A. (2001). Topics in Circular Statistics . World Scientific Publishing Company. ISBN 978-981-02-3778-3 . ^ Jammalamadaka, S. Rao; SenGupta, A. (2001). Topics in circular statistics . New Jersey: World Scientific. ISBN 981-02-3778-2 . Retrieved 2011-05-15 .
추가 읽기 아브라모위츠, M., 스테건, I. A. (edd.), 수학적 기능 핸드북 , 국가표준국, 1964; 도버 출판물 , 1965. ISBN 0-486-61272-4 "알고리즘 AS 86: 폰 미제스 분포 함수", 마디아, 적용 통계, 24, 1975(pp. 268–272) "알고리즘 518, 불완전한 베셀 함수 I0: 폰 미제스 분포", 힐, ACM 수학 소프트웨어 거래, 제3권, 1977년 9월, 페이지 279–284. 베스트, D.와 피셔, N.(1979년). 폰 미제스 분포의 효율적인 시뮬레이션. 적용통계, 28, 152–157. 에반스, M, 헤이스팅스, N, 피콕, B, "본 미세스 배급" 통계분포의 41장, 3차 개정. 뉴욕 와일리 2000. 피셔, 니콜라스 1세, 순환 데이터의 통계 분석. 뉴욕 케임브리지 1993. "통계 분포", 2번째. 에디션, 에반스, 헤이스팅스, 피콕, 존 와일리 앤 선즈 , 1993년 (제39장) ISBN 0-471-55951-2 Borradaile, Graham (2003). Statistics of Earth Science Data . Springer . ISBN 978-3-540-43603-4 . Retrieved 31 Dec 2009 .
이산형 일변도의
연속 일변도의
의 지지를 받고 있는. 경계 간격 의 지지를 받고 있는. 반무한 간격을 두고 지지의 대체로 실선 지지하여 누구의 타입이 다른가.
혼합 일변도의
다변량 (공동) 방향 퇴보하다 그리고 단수 가족들