주어진


16

n 번째 누적이 1 로 주어진 분포에 대한 정보가 있습니까?1n ? 누적 생성 함수는

κ(t)=01etx1x dx.
나는 임의의 변수의 제한 분포로 그것을 가로 질러 보았지만 그것에 대한 정보를 찾을 수 없었습니다.

당신이 부여 한이 함수 에 소유권이 주장 된 속성 이 있음을 알 수 없습니다 ! yoiur 작업을 수정해야합니다. 1 + t x 로 0에 가까운 정수 n에 지수를 근사하면 , 0에 가까운 정수는 t / x 가되어 분기됩니다. 따라서 적분은 누적 생성 함수를 나타낼 수 없습니다. κ(t)1+txt/x
kjetil b halvorsen

@kjetilbhalvorsen 내가 확실하지 않습니다. 1 + t x로 를 근사 하면 t x가됩니다.etx1+tx정수의 경우 x =t입니다. 도에있어서,I 준 쌍곡선 함수는 사인 및 코사인 적분 측면에서 알려진 적분있다. 그 표시하려면κ(t을)청구 된 속성이 단지 주변 전체 테일러 시리즈을0에 대한전자tX와의 테일러 시리즈 얻을 합을 통해 통합을 밀어κ(t)의 주위에공을. txx=tκ(t)0etxκ(t)0
guy

sympy는 적분이 (자신의 편심적인 방식으로) 분 기적이라고 말합니다. 그러나 sympy는 틀렸어 야합니다. 지금은 수치 적분을 실험하여 잘 작동합니다. 다시 시도합니다.
kjetil b halvorsen

Wolphram 알파 결과를 살펴보면 t도 0에 가까울 때 0이 아닌 한계를 가지지 만 정확하지 않습니다. κ(0)=0
kjetil b halvorsen

2
나는 그것이 에 절대적으로 연속적이라고 생각합니다 . 그것은 포아송 임의 변수의 한계로 실현된다. 같은 N 레이트 화합물과 포아송 1 1 / N (1)(0,)n및 점프 분포 밀도fn(x)11/n11x dx는이 분포에 약하게 수렴합니다. fn(x)1xI(1/n<x<1)
guy

답변:


8

cumulants의 값을 알면이 확률 분포의 그래프가 어떻게 보이는지 알 수 있습니다. 분포의 평균과 분산은

E[Y]=κ1=1,Var[Y]=κ2=12

왜도 및 초과 첨도 계수는

γ1=κ3(κ2)3/2=(1/3)(1/2)3/2=223

γ2=κ4(κ2)2=(1/4)(1/2)2=1

{0,1,...,m}{p0,p1,...,pm},k=0mpk=1, and then use the cumulants to calculate the raw moments, with the purpose of forming a system of linear equations with the probabilities being the unknowns. Cumulants are related to raw moments by

κn=μni=1n1(n1i1)κiμni
Solved for the first five raw moments this gives (the numerical value at the end is specific to the cumulants in our case)
μ1=κ1=1μ2=κ2+κ12=3/2μ3=κ3+3κ2κ1+κ13=17/6μ4=κ4+4κ3κ1+3κ22+6κ2κ12+κ14=19/3μ5=κ5+5κ4κ1+10κ3κ2+10κ3κ12+15κ22κ1+10κ2κ13+κ15=243/15
If we (momentarily) set m=5 we have the system of equations

k=05pk=1,k=05pkk=1k=05pkk2=3/2,k=05pkk3=17/6k=05pkk4=19/3,k=05pkk5=243/15s.t.pk0k

Of course we do not want m to be equal to 5. But increasing gradually m (and obtaining the value of the subsequent moments), we should eventually reach a point where the solution for the probabilities stabilizes. Such an approach cannot be done by hand -but I have neither the software access, nor the programming skills necessary to perform such a task.


This is cool. Maybe I could do some kind of Edgeworth expansion as well? Actually, I have an idea of what the density looks like already (assuming it exists) since I can simulate directly from it. It is very strange - it looks uniform over some range (0,a) and then on (a,) it decays with something like an exponential tail (it's been a long time since I did the simulation).
guy

Thanks. Of course you can always perform an Edgworth expansion based on the cumulants, but I wonder how well it will perform, given the strange shape you describe. It would be interesting to contrast the two.Can you tell me the value for a?
Alecos Papadopoulos

Dug up my old code and found a1. If Yκ(t) then [YY<1] is approximatey U(0,1) and [Y1Y>1] is approximately gamma distributed with shape 1.4 and mean 0.64.
guy

What do you mean by Yκ(t)?
Alecos Papadopoulos

1
So what does the pdf look like then? As for fitting by moments, is the fit 'robust' and 'stable' as one increases the number of moments used (4, 5, 6, 7 or 8 etc), or is it all over the place?
wolfies
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.