모든 합리성을지지하는 이산 rv 구성


19

이것이이 질문 의 구성주의 속편이다 .

구간 에서 모든 합리성을 지원하는 불연속 균일 랜덤 변수를 가질 수 없다면 다음으로 가장 좋은 것은 다음과 같습니다. [0,1]

확률 변수 구조체 이 지원 갖는다 , 그것은 다음 것을 일부 분포. 그리고 저의 장인은이 랜덤 변수가 우리가 얻고 자하는 것을 추상적으로 정의하여 생성되는 것이 아니라 기존 분포 로 구성 해야합니다 .QQQ[0,1]

그래서 나는 다음을 생각해 냈습니다.

하자 파라미터와 기하 분포 변종 II는 다음의 이산 랜덤 변수 일 , 즉X0<p<1

X{0,1,2,...},P(X=k)=(1p)kp,FX(X)=1(1p)k+1

또한하자 동일한 매개 변수로 기하 분포 변종 I는 다음의 이산 확률 변수 일 수 , 즉,Yp

Y{1,2,...},P(Y=k)=(1p)k1p,FY(Y)=1(1p)k

XY 는 독립적입니다. 이제 랜덤 변수를 정의하십시오

Q=XY

조건부 분포를 고려하십시오.

P(Qq{XY})

즉 이른바 "조건부 Q 의 비율 X 이상 Y 조건부 X 작거나 같 인 Y ." 이 조건부 분포는 {0,1,1/2,1/3,...,1/k,1/(k+1),...,2/3,2/4,...}=Q[0,1] .

"질문"은 다음과 같습니다. 누군가 조건부 확률 질량 함수를 제공 할 수 있습니까?

코멘트가 "폐쇄해야합니까?" 오늘날 닫힌 양식을 구성하는 것은 그렇게 분명하지 않으므로 다음과 같이하겠습니다. 우리는 [0,1] 에서 유리수를 입력 [0,1]하고 확률을 얻을 수있는 함수형을 찾고 있습니다 (일부 경우 물론 매개 변수 p 의 지정된 값 ), pmf를 나타내는 그래프 로 이어집니다. 그런 다음 p 를 변화시켜 그래프가 어떻게 변하는 지 확인하십시오.

그것이 도움이된다면, 지원의 한 가지 또는 두 가지 경계를 열어 둘 수는 있지만, 이러한 변형은 pmf 의 상한 및 / 또는 하한 값을 확실히 그래프로 나타낼 수있는 능력을 박탈합니다 . 또한 상한을 열면 컨디셔닝 이벤트 고려해야합니다 .{X<Y}

또는 pmf와 함께 사용되는 한이 지원이있는 다른 rv도 환영 합니다.

지오메트리 분포는 지지점에 0을 포함하지 않는 두 가지 변형을 쉽게 사용할 수 있기 때문에 사용했습니다 (따라서 0으로 나누는 것을 피할 수 있습니다). 분명히, 일부 잘림을 사용하여 다른 이산 rv를 사용할 수 있습니다.

나는 확실히이 질문에 현상금을 넣을 것이지만, 시스템은 이것을 즉시 허용하지 않습니다.


1
당신은 의미합니까 ? (무조건 변수를 조건부로 정의하는 것은 말이되지 않습니다. 이런 식으로 만 분포를 정의 할 수 있습니다)Q=XY1{XY}
Stéphane Laurent

1
Q는 셀 수 있습니다. N = {1, 2, ...}와 Q 사이에 1-1 대응 관계가 있다는 것을 알고 있습니다. 그러한 대응 관계를 찾을 수 있으면 해결책은 N에 대한 분포를 선택하여 사용하는 것입니다. Q의 해당 요소를 선택합니다.
Adrian

어쨌든 모든 돌이킬 수없는 분수 대해 를 계산해야 하며 이것은 . p / q Pr ( X = p , X = 2 p , ) × Pr ( Y = q , Y =Pr(X/Y=p/q)p/qPr(X=p,X=2p,)×Pr(Y=q,Y=2q,)
Stéphane Laurent

1
pmf를 제공해야한다는 요구가 폐쇄 양식이 필요함을 의미합니까? 또는 @ StéphaneLaurent의 무한 합이 조건을 충족시키기에 충분합니까?
Juho Kokkala

1
하자 하고 Y 게시물의 RV. P r [ Q = q ] = P r [ Y = f 1 ( q ) ]f:NQ[0,1]Pr[Q=q]=Pr[Y=f1(q)]
Adrian

답변:


19

확률 질량 을 갖는 세트 를 지원 하는 이산 분포 를 고려하십시오.{ (F{(p,q)|qp1}N2

F(p,q)=321+p+q.

이것은 쉽게 합산되어 (관련된 모든 계열은 기하학적 임) 실제로 분포임을 보여줍니다 (총 확률은 1입니다).

제로가 아닌 유리수를 들어 하게 최저 관점에서 표현 될 :이고, 과 .a / b = x b > 0 gcd ( a , b ) = 1xa/b=xb>0gcd(a,b)=1

G [ 0 , 1 ] QF 는 규칙을 통해 에서 불연속 분포 를 유도합니다.G[0,1]Q

G(x)=G(ab)=n=1F(an,bn)=321+a+b2.

(및 ). 모든 유리수 는 0이 아닌 확률을 갖습니다 ( 양의 확률로 값 중 포함해야하는 경우 과 같은 다른 숫자에서 일부 확률을 빼서 할당 ).(0)=00 1 0(0,1]010

이 구성을 이해하려면 다음 묘사를보십시오 .에프

[F의 그림]

p , q F p / q p q 0 1 G G G ( 1 ) 1 F ( 1 , 1 ) + F ( 2 , 2 ) + F ( 3 , 3 ) +에프 는 양의 적분 좌표 로 모든 점 에서 확률 질량을 제공 합니다. 값은 원형 심볼의 컬러 영역으로 표시됩니다. 선은 플롯에 나타나는 가능한 모든 좌표 및 조합에 대해 경사 를 갖습니다 . 그들은 원형 기호와 같은 방식으로 색이 지정됩니다 : 경사에 따라. 따라서, (명확 내지 기울기 통해 받는)와 색에 대응 인수 의 및 값 각 라인에 누워 모든 원의 면적을 합산하여 얻어진다. 예를 들어,p,에프p/01G(1) = 주어진 기울기 의 주 대각선을 따라 모든 (빨간색) 원의 면적을 합하여 구합니다. 입니다.13 / 8 + 3 / 32 + 3 / 128 + = 1 / 2F(1,1)+F(2,2)+F(3,3)+3/8+3/32+3/128+=1/2

그림

이 그림은 을 제한 하여 달성 한 대한 근사값을 보여줍니다 . 에서 까지의 합리적인 숫자로 값을 표시합니다 . 최대 확률 질량은 입니다.q 100 3044Gq1003044(1) (1)1/100112,314,110,362,362,142,

다음은 의 전체 CDF입니다 (이미지의 해상도에 정확함). 방금 나열된 6 개의 숫자는 눈에 띄는 점프의 크기를 나타내지 만 CDF의 모든 부분은 예외없이 점프로 구성됩니다.G

그림 2


1
감사! 나는 건축을 이해하는 과정에 있습니다. a) 는 이변 량이지만 연결하는 표현 에서는 단 변량으로 나타납니다. 뭔가 빠졌습니까? b) 는 단 변량 이므로 인상적으로 보이는 첫 번째 그래프의 모든 점이 가로 축의 다른 값을 나타내는 것 같습니다 (물론 이러한 척도로 정확하게 표현할 수는 없지만) 맞습니까? FGGG
Alecos Papadopoulos 2016

나는 당신의 의견, Alecos를 다룰 수있는 그림을 완성하고 답변에 추가했습니다. 내가 함께 시작했습니다 수 있다는 것을 참고 있는 이산 분포 및 건설 같은 방식으로; 이 특정 분포는 계산을 쉽게하기 위해 선택되었습니다. GFG
whuber

이전 주석의 첫 번째 질문과 같이 대신 합니까? 즉, 이고 입니까? F(ab,n)p=a/bq=nF(abn)p=a/bq=n
Alecos Papadopoulos 2016

이것은 내 것보다 더 나은 대답입니다! 나는 두 가지 작은 것을 발견했다. 나는 당신의 F (p, q)가 4로 쓰여진 것으로 생각한다. 또한 아래의 식에서 "F는 불연속 분포 G를 유도합니다"F (na, nb) no?
Adrian

@Adrian, Alecos 오타를 잡아 주셔서 감사합니다. 은 이어야 하고 의 표기는 분명히 틀립니다. 내가 바로 고쳐 줄게 - (1) F11F
whuber

8

의견을 정리하고 명확하게하기 위해 답변으로 게시하겠습니다. 그러나 내가하는 일은 문제를 다른 문제로 줄이는 것이므로 매우 만족스럽지 않을 것으로 기대합니다.

내 표기법 :

QQ 그의 지원하는 RV입니다 - 내 있다 되지 는 AS 같은 자신의 영업 구조 . 우리는 와 사용 하여이 를 정의 할 것 입니다. Q Q XQ[0,1]QQ QYfXYQYf

Y 누구의 지원이 어떤 RV입니다 더 - 영업 이익에 의해 주어진 일하는 것이, 예를 들어. YN{1,2,}Y

f : Nf 는 일대일 대응입니다. 이고 은 그 역수입니다. 우리는 이것이 존재한다는 것을 알고 있습니다.f:NQ[0,1]f1

이제 와 찾는 것으로 문제를 줄일 수 있다고 주장합니다 .ff1

그냥하자 하면 모든 작업이 완료된다. 의 PMF 는 입니다. Q Pr [ Q = q ] = Pr [ Y = f 1 (Q=f(Y)QPr[Q=q]=Pr[Y=f1(q)]

편집하다:

다음은 일대일 대응이 아니더라도 (중복 때문에) 의 역할을하는 함수 g입니다 .f

g <- function(y) {
    y <- as.integer(y)
    stopifnot(y >= 1)
    b <- 0
    a <- 0
    for (unused_index in seq(1, y)) {
        if (a >= b) {
            b <- b+1
            a <- 0
        } else {
            a <- a+1
        }
    }
    return(sprintf("q = %s / %s", a, b))
    ## return(a / b)
}

(+1) 아니요, 귀하의 접근 방식은 매우 적용 가능한 결과와 알고리즘 에 도달하기 위해 추상적 접근 방법을 생각하고 사용하는 방법의 훌륭한 예라고 생각 합니다. 내가 지금 이해하고있는 주요 요점은, 를 지원 하는 이산 형 분포 의 pmf를 기능적인 형태로 사용하여 원하는 구성을 얻을 수 있다는 것입니다 . 물론 와 을 찾아야 합니다. 당신이 나 f f 1N{1,2,}ff1
보다이

jcu.edu/math/vignettes/infinity.htm 참조 : 유사한 "대각선 패턴"을 사용할 수 있습니다. 어려운 부분은 대한 표현을 얻는 것 입니다. 어떻게 해야할지 모르겠지만 math.stackexchange.com에 요청할 수 있습니다 (또는 먼저 인터넷 검색을 더하십시오). f1
Adrian

당신이 제공 한 링크에서 그것은 어느 시점에서 다음과 같이 말합니다 . "여기서 실제로 공식을 나타 내기보다는 어떤 일이 발생했거나 존재한다는 것을 보여 주어야한다." 글쎄, 내 질문의 요점은 실제로 공식을 제시하는 것입니다 . 나는이 질문을 "구성 주의자"라고 불렀습니다.
Alecos Papadopoulos 2016

1
작동하는 알고리즘을 제공 할 수 있다고 생각합니다. 조금 더 생각해 보겠습니다.
Adrian

나는 무언가를 올렸다 .Q를 시뮬레이션 할 수는 있지만 PMF 문제는 해결하지 못한다.
Adrian
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.