세트의 측정 지수의 편견 추정기?


12

우리가 (측정 가능하고 적절하게 동작하는) 세트 을 가지고 있다고 가정하자 . 여기서 는 컴팩트하다. 또한, 우리가 걸쳐 균일 한 분포에서 샘플을 그릴 수 있다고 생각 르 베그 측정 WRT 우리가 측정 알고 . 예를 들어, 아마도 는 포함 하는 상자 입니다 .SBRnBBλ()λ(B)B[c,c]nS

고정 경우 에서 점을 균일하게 샘플링하고 점이 내부 또는 외부에 있는지 확인 하여 를 추정하는 간단한 편견없는 방법이 있습니까?αReαλ(S)BS

제대로 작동하지 않는 것의 예로, 포인트 을 샘플링한다고 가정 합니다. 그런 다음 Monte Carlo 추정값 그러나 는 의 편견 추정치 이지만 가 의 편견 추정치 인 경우는 아닙니다. . 이 알고리즘을 수정하는 방법이 있습니까?kp1,,pkUniform(B)

λ(S)λ^:=#{piS}kλ(B).
λ λ(S)E-α λ E-αλ(S)λ^λ(S)eαλ^eαλ(S)

답변:


11

다음과 같은 자원을 사용할 수 있다고 가정하십시오.

  1. 견적 자 액세스 할 수 있습니다 .λ^
  2. λ^ λ(S) 는 대해 편향되지 않습니다.λ(S)
  3. λ^ C 는 거의 확실하게 묶여 있습니다.C
  4. 당신은 일정을 알고 하고,C
  5. 원하는만큼 를 독립적으로 구현할 수 있습니다 .λ^

이제 인 경우 다음을 유지합니다 (Taylor 확장 ).u>0expx

eαλ(S)=eαCeα(Cλ(S))=eαCk0(α[Cλ(S)])kk!=eαCeuk0eu(α[Cλ(S)])kk!=euαCk0ukeuk!(α[Cλ(S)]u)k

이제 다음을 수행하십시오.

  1. 샘플 .KPoisson(u)
  2. 폼 의 IID 바이어스 추정기로서 .λ^1,,λ^Kλ(S)
  3. 견적 자 반환

Λ^=euαC(αu)Ki=1K{Cλ^i}.

Λ^ 는 의 음이 아닌, 편향되지 않은 추정량입니다 . 이 때문입니다λ(S)

E[Λ^|K]=euαC(αu)KE[i=1K{Cλ^i}|K]=euαC(αu)Ki=1KE[Cλ^i]=euαC(αu)Ki=1K[Cλ(S)]=euαC(αu)K[Cλ(S)]K

따라서

E[Λ^]=EK[E[Λ^|K]]=EK[euαC(αu)K[Cλ(S)]K]=euαCk0P(K=k)(αu)K[Cλ(S)]K=euαCk0ukeuk!(α[Cλ(S)]u)k=eαλ(S)

이전 계산으로.


흥미 롭습니다! 위의 의해 묶여 있기 때문에 질문에 설명 된 대한 추정기는 여기에서 작동 하지 않습니다 . 또한 이것이 어떻게 @whuber의 답변과 모순되지 않습니까? 이것이 편견이없는 이유는 쉬운가? 많은 질문에 대해 유감 λ^λ(B)<
Justin Solomon

1
알고 있기 때문에 설명하는 추정기가 작동합니다 . 나는 이것이 가정 때문에 다른 답변과 모순되지 않는다고 생각한다 . 편견이없는 추정자에 대한 유한 액세스 권한이 주어지면이 구성이 작동하지 않을 것이라고 생각합니다. 불편 함은 에 대한 기대치를 위의 거듭 제곱 과 비교함으로써 비롯됩니다 . 나는 대답을 더 명확하게 할 것입니다. λ(B)5Λ^
πr8

편견의 증거의 두 번째 줄에서 제품과 기대를 교환 할 수 있습니까?
jbowman

2
iid로 계산 되었기 때문에 괜찮은 것 같습니다.
저스틴 솔로몬

2
+1 이것이 흥미롭고 교훈적인 예라고 생각합니다. 그것은 내 대답에 암묵적인 가정을하지 않음으로써 성공합니다 : 샘플 크기가 지정되었거나 최소한 경계가되어 있습니다.
whuber

10

대답은 부정적입니다.

균일 한 표본에 대한 충분한 통계량 은 에있는 것으로 관측 된 점 의 개수 입니다 이 개수에는 이항 분포 분포가 있습니다. 및 씁니다XS.(n,λ(S)/λ(B))p=λ(S)/λ(B)α=αλ(B).

샘플 크기를 들어 하자 임의 (unrandomized) 추정기 수 기대는n,tnexp(αλ(S))=exp((αλ(B))p)=exp(αp).

E[tn(X)]=x=0n(nx)px(1p)nxtn(x),

이는 대부분의 학위 다항식 동일 에서 그러나 지수 는 에서 다항식으로 표현할 수 없습니다 (한가지 증명 : 도함수를 취 하십시오. 기대 결과는 0이지만 지수 의 지수 인 도함수는 0 수 없습니다.)np.αp0,exp(αp)p.n+1p,

무작위 추정 추정자에 대한 시연은 거의 동일하다. 기대에 부딪 치면 다시p.

결과적으로 편향 추정기가 존재하지 않습니다.


1
아, 다우 너! 좋은 증거 주셔서 감사합니다. 그러나 대한 Taylor 시리즈는 상당히 빠르게 수렴됩니다. "아마도"편의되지 않은 "추정기가 있습니까? 확실하지 무엇을 (나는별로 통계 학자 :-)의이야) 수단 있음exp(t)
저스틴 솔로몬

정확히 얼마나 빨리요? 대답은 의 값에 달려 있으며 그 값이 무엇인지 모르기 때문에 문제가 있습니다. 과 사이에 있다는 것만 알고 있습니다 원하는 경우이를 사용하여 바이어스에 대한 경계를 설정할 수 있습니다. αp0α.
whuber

내 응용 프로그램에서 는 의 많은 부분을 차지할 것으로 기대 합니다. 의사-마지막 메트로폴리스-해 스팅 수용 비율에서이 값을 사용하고 싶습니다.이 방법이 제어 가능한 수준의 편향도 처리 할 수 ​​있는지 확실하지 않습니다 ...SB
Justin Solomon

4
BTW이 질문에 대한 다른 답변에 대한 귀하의 의견에 진심으로 감사드립니다!
저스틴 솔로몬
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.