假设是 iid 随机变量。预计序列何时会首次减少?X1,X2, … ,XnX1,X2,…,Xn

机器算法验证 可能性 自习 独立同居
2022-03-06 07:07:48

正如标题中所建议的那样。假设的连续 iid 随机变量考虑 ,的事件,因此是序列第一次减少的时候。那么的值是多少?X1,X2,,XnfX1X2XN1>XNN2NE[N]

我首先尝试评估我有 同样,我得到了随着变大,计算变得更加复杂,我找不到模式。谁能建议我应该如何进行?P[N=i]

P[N=2]=f(x)F(x)dx=F(x)22|=12P[N=3]=f(x)xf(y)F(y)dydx=f(x)1F(x)22dx=F(x)F(x)3/32|=13
P[N=4]=18i

4个回答

如果随机变量的可交换 if and on如果因此, 通过对称性。因此,{Xi}i1

N=min{n:Xn1>Xn},
NnX1X2Xn1
Pr(Nn)=Pr(X1X2Xn1)=1(n1)!,()
E[N]=n=1Pr(Nn)=e2.71828

PS 人们询问的证明。由于序列是可交换的,因此对于任何排列,我们有 因为我们有可能的排列,结果如下。()π:{1,,n1}{1,,n1}

Pr(X1X2Xn1)=Pr(Xπ(1)Xπ(2)Xπ(n1)).
(n1)!

正如 Silverfish 所建议的,我在下面发布了解决方案。

P[N=i]=P[X1X2Xi1>Xi]=P[X1X2Xi1]P[X1X2Xi1Xi]=1(i1)!1i!
P[Ni]=1P[N<i]=1(112!+12!13!++1(i2)!1(i1)!)=1(i1)!

因此E[N]=i=1P[Ni]=i=11(i1)!=e

一种论点:的可能排列我们对增加直到倒数第二个位置然后减少的排序感兴趣:这要求最大值位于位置,而之一位于最终位置。由于有种方法可以在我们的有序序列中挑选出前项中的一项并将其移动到最终位置,因此概率为:Xin!X1,,Xnn1n1Xin1n1

Pr(N=n)=n1n!

注意 ,所以这与积分得到的结果是一致的。Pr(N=2)=212!=12Pr(N=3)=313!=13Pr(N=4)=414!=18

的期望值,我们可以使用:N

E(N)=n=2nPr(N=n)=n=2n(n1)n!=n=21(n2)!=k=01k!=e

(为了使总和更明显,我使用了;对于不熟悉这个总和的读者,请使用泰勒级数并替换。)k=n2 ex=k=0xkk!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让我满意。

编辑:我的回答不正确。我把它作为一个例子,说明像这样一个看似简单的问题是多么容易被误解。

我认为您的数学对于的情况不正确。我们可以通过一个简单的模拟来检查这一点: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

因此,如果我们相信我们的模拟结果,看起来模式是但这也是有道理的,因为您真正要问的是,您所有观察的子集中的任何给定观察是最小观察的概率是多少(如果我们假设 iid,那么我们假设可交换性,因此顺序是任意的)。其中一个必须是最小值,所以真正的问题是随机选择的任何观察值是最小值的概率是多少。这只是一个简单的二项式过程。P[N=X]=1x