仮定


10

タイトルのとおり。仮定のPDFに連続しiid確率変数であり、F。そのイベントを考えるX 1 ≤ X 2 ... ≤ X N - 1 > X N、N ≥ 2従って、Nはシーケンスが最初に低下した場合です。次に、E [ N ]の値は何ですか?X1,X2,…,XnfX1≤X2…≤XN−1>XNN≥2NE[N]

最初にを評価しようとしました。私は P [ N = 2 ]P[N=i] 同様に、P[N=4]=1を得た

P[N=2]=∫−∞∞f(x)F(x)dx=F(x)22|−∞∞=12P[N=3]=∫−∞∞f(x)∫x∞f(y)F(y)dydx=∫−∞∞f(x)1−F(x)22dx=F(x)−F(x)3/32|−∞∞=13
。私は大きくなり、計算が複雑になりますし、私はパターンを見つけることができません。誰かが私がどのように進むべきかを提案できますか?P[N=4]=18i

これはコースや教科書の質問ですか?その場合は、[self-study]タグを追加してWikiをお読みください。
— Silverfish、2015

1
n!1,2,…,nXin≥2n−1n−11/21/31/8

(そして、シリーズの結果を推測できない場合は、平均を見つけるために合計します。おそらく、そのシミュレーションを実行する必要があります。小数点以下の最初の2、3桁がわかります。)
— Silverfish

今日受験した問題です。ヒントをありがとう、私はそれを解決する方法を見つけました。
— Hao The Cabbage

2
Xi

回答:


9

{Xi}i≥1

N=min{n:Xn−1>Xn},
N≥nX1≤X2≤⋯≤Xn−1
Pr(N≥n)=Pr(X1≤X2≤⋯≤Xn−1)=1(n−1)!,(∗)
E[N]=∑n=1∞Pr(N≥n)=e≈2.71828…

(∗)π:{1,…,n−1}→{1,…,n−1}

Pr(X1≤X2≤⋯≤Xn−1)=Pr(Xπ(1)≤Xπ(2)≤⋯≤Xπ(n−1)).
(n−1)!
# Monte Carlo
N <- 10^6
dec <- numeric(N)
for (i in 1:N) {
    j <- 1
    x <- rnorm(1, 0, 1) # plug your favorite distribution here!
    repeat {
        j <- j + 1
        y <- rnorm(1, 0, 1)
        if (y < x) {
            dec[i] <- j
            break
        }
        x <- y
    }
}
cat(mean(dec), "\n")
# A good example of how to program C in R!

2
Pr(Y=y)Pr(Y≥y)

Xi

1
Xi

1
実際には、独立は必要ありません。交換性は十分です。結果はより強力です。これを私の答えに追加します。
— Zen

3
連続変数に対して普遍的であることは直感的です。これを明確にする1つの方法は、確率積分変換を適用してもイベントが変わらないことを認識することです。これにより、変数が共通の一様分布を持つ場合に減少します。
— whuber

8

P[N=i]=P[X1≤X2…≤Xi−1>Xi]=P[X1≤X2…≤Xi−1]−P[X1≤X2…≤Xi−1≤Xi]=1(i−1)!−1i!
P[N≥i]=1−P[N<i]=1−(1−12!+12!−13!+⋯+1(i−2)!−1(i−1)!)=1(i−1)!

E[N]=∑i=1∞P[N≥i]=∑i=1∞1(i−1)!=e


7

Xin!X1,…,Xnn−1n−1Xin−1n−1 順序付けられたシーケンスの項を最後の位置に移動すると、確率は次のようになります。

Pr(N=n)=n−1n!

Pr(N=2)=2−12!=12Pr(N=3)=3−13!=13Pr(N=4)=4−14!=18

N

E(N)=∑n=2∞nPr(N=n)=∑n=2∞n(n−1)n!=∑n=2∞1(n−2)!=∑k=0∞1k!=e

k=n−2 ex=∑k=0∞xkk!x=1

シミュレーションで結果を確認できます。Rのコードを次に示します。

firstDecrease <- function(x) {
    counter <- 2
    a <- runif(1)
    b <- runif(1)
    while(a < b){
        counter <- counter + 1
        a <- b
        b <- runif(1)
    }
    return(counter)
}

mean(mapply(firstDecrease, 1:1e7))

これは2.718347、2.71828私を満足させるのに十分近いほど戻ってきました。


-1

編集:私の答えは間違っています。このような一見単​​純な質問が誤解されやすい例として残しておきます。

P[N=4]

n=50000
flag <- rep(NA, n)
order <- 3
for (i in 1:n) {
  x<-rnorm(100)
  flag[i] <- all(x[order] < x[1:(order-1)])==T
}
sum(flag)/n

私たちに与える:

> sum(flag)/n
[1] 0.33326

order用語を4に変更すると、次のようになります。

> sum(flag)/n
[1] 0.25208

そして5:

> sum(flag)/n
[1] 0.2023

P[N=X]=1x


1
Xnn−1Xin−1

それは、ちょっとした誤解ではありません。あなたは正しい、私は間違っている。
— ダルトンハンス2015
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.