비 중심 지수 분포의 예상 로그 값


9

가 위치 k 및 rate \ lambda 와 함께 비 중심적으로 지수 분포되어 있다고 가정 합니다. 그런 다음 E (\ log (X)) 는 무엇입니까 ?XkλE(log(X))

난에 대해 알고 k=0 , 응답이 log(λ)γ 여기서 γ 오일러 - 마스케 로니 상수이다. k> 0 일k>0어떻습니까?


Mathematica에 통합을 시도 했습니까?

4
I는 가정 k>0 (밀도는 다음과 같이 기입된다 λexp{λ(xk)} ), 그렇지 않으면 x<0 대 무서운 결과를 초래할 가능성> 0으로 Elogx .
jbowman

2
I있어 E[log(X)]=ekλΓ(0,kλ)+log(k) . Assumptions매개 변수 공간 지정 명령 을 사용하면 Mathematica가 더 빠릅니다 .

4
불완전한 상단 감마 함수가 닫힌 형태 로 계산 됩니까? (나에게는 그렇지 않습니다.) 이것은 단순히 표기법을 통해 정수를 숨기는 것입니다.
추기경

2
@NeilG 이것은 Mathematica 코드 Integrate[Log[x + k]*\[Lambda]*Exp[-\[Lambda]*x], {x, 0, \[Infinity]}, Assumptions -> k > 0 && \[Lambda] > 0]입니다. 복사하여 .nb 파일에 붙여 넣기 만하면됩니다. Wolfram Alpha가 제한을 포함 할 수 있는지 확실하지 않습니다.

답변:


11

원하는 적분은 무차별 대입 조작으로 제출 될 수 있습니다. 여기서 우리는 대신 약간 더 확률적인 맛을 가진 대체 유도를 제공하려고 시도합니다.

하자 의 위치 파라미터와 noncentral 지수 확률 변수 일 및 레이트 파라미터 . 그런 다음 여기서 입니다.XExp(k,λ)k>0λX=Z+kZExp(λ)

참고 그렇게하는 사용 음수 확률 변수의 기대를 계산하기위한 표준 사실 , 그러나 에서 부터 그래서 마지막 평등은 치환 따른다log(X/k)0

Elog(X/k)=0P(log(X/k)>z)dz=0P(Z>k(ez1))dz.
P(Z>k(ez1))=exp(λk(ez1))z0ZExp(λ)
Elog(X/k)=eλk0exp(λkez)dz=eλkλkt1etdt,
t=λkez, 입니다.dz=dt/t

최근 디스플레이의 오른쪽 크기 적분 그냥 정의에 의해 그렇게 질문에 대한 의견에서 @ Procrastinator의 Mathematica 계산 에 의해 확인 된 .Γ(0,λk)

ElogX=eλkΓ(0,λk)+logk,

NB : 대신 등가 표기법 도 종종 사용됩니다 .E1(x)Γ(0,x)


4
+1 @Michael Chernick 모든 사람이 게으른 것은 아닌 것 같습니다;).

정말 대단합니다. 불완전한 감마 함수의 많은 구현이 첫 번째 매개 변수를 엄격하게 양수로 제한한다는 것을이 구현하는 사람에게 지적하고 싶습니다. ID 는 그 사소한 문제를 해결합니다. Γ(0,z)=Ei(z)
Neil G
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.