这个问题问关于
=1−∏i=1n(1−Fi)(1)
随着增长和均匀收缩,使得(a)全部为非负数,并且(b)它们加总为1。(这些来自的构造和概率公理。)nFi Fi
根据定义,该乘积是其对数的指数:
∏i=1n(1−Fi)=exp(∑i=1nlog(1−Fi)).
应用于泰勒定理(具有余数的拉格朗日形式)确定log
log(1−Fi)=−Fi−12ϕ2i≥−Fi−12F2i
对于一些在区间。换句话说,这些对数等于至多是一些术语至多倍。但是,当足够大以确保所有都小于给定的(由的均匀收缩确保的条件)时,则(b)表示因此ϕi[0,Fi]−Fi 1/2F2inFiϵ>0Finϵ>∑Fi=1
∑i=1nF2i≤∑i=1nϵ2<∑i=1n(1n)2=1n.
所以
−1=−∑i=1nFi≥∑i=1nlog(1−Fi)≥−∑i=1nFi−121n=−1−12n
压缩两个收敛到序列之间的对数。由于是连续的,因此乘积收敛到该限制的指数。所以−1exp∏ni=1(1−Fi)exp(−1)
limn→∞(1−∏i=1n(1−Fi))=1−exp(−1)≈0.632,
QED。
仔细查看此分析,该近似值中的误差(始终为下界)不会大于 例如,将标准正态分布划分为和之间的切片会在模式附近产生一个最大,它将近似等于该矩形的面积。前述界限确定式的值将在其极限值的以内。实际误差要小一个数量级,
(exp((n/2)max(F2i))−1)exp(−1).
n=400−44Fi0exp(−1/2)/50≈0.012(1)0.0110.001041。这是其中的计算
R
(我们可以信任,因为相对于都不小):
fi1
f <- diff(pnorm(seq(-4, 4, length.out=401))) # The normal "slices".
f <- f / sum(f) # Make them sum to unity.
exp(-1) - prod(1 - f) # Compute the error.
实际上,1 - prod(1-f)
是而是。0.6331615…1−exp(−1)0.6321206…