Skip to content

模型Model

均匀样本的随机间距

Uniform spacings · Dirichlet spacings

把均匀样本连同两端点之间的间距合起来,得到单纯形上的均匀随机向量,也可由独立指数变量归一化生成。

形式陈述 ​

在单位区间上随机放置 n 个点,留下的 n+1 段长度怎样共同变化?设整数 n≥1,U1,…,Un 是来自 (0,1) 均匀分布的独立同分布样本。将其排序并加入固定端点

U(0)=0<U(1)<⋯<U(n)<U(n+1)=1.

定义间距 Di=U(i)−U(i−1),1≤i≤n+1。它们非负、总和为一,其联合分布是参数全为一的 Dirichlet 分布,即一类支持在单纯形上的联合概率律。精确地说,前 n 个坐标的密度是

f(d1,…,dn)=n!1{di>0, ∑i=1ndi<1},Dn+1=1−∑i=1nDi.

这是相对于前 n 个自由坐标的 n 维 Lebesgue 密度;整个向量位于 Rn+1 的和为一平面上,没有普通的 (n+1) 维密度。

还有一个等价构造:若 E1,…,En+1 独立服从率一指数分布,S=∑iEi,则

(D1,…,Dn+1)=d(E1/S,…,En+1/S).
直觉

随机切点在区间上没有偏爱的位置,但切出来的段长不独立。一段占得多,其他段就必须少一些,因为总长度固定为一。

指数归一化构造先让各段独立生成一个正“份额”,再一起除以总份额。独立性存在于归一化之前;共享同一个分母之后,长度约束与负相关就出现了。

例子与边界

两个切点留下三段 ​

取 n=2。前两个间距 (D1,D2) 在三角形 d1>0,d2>0,d1+d2<1 上密度为二,第三段为 1−d1−d2。

每一段的边缘密度都是 2(1−d),0<d<1,所以

EDi=13,EDi2=16,Var(Di)=118.

对不同两段直接积分,

E(D1D2)=2∫01∫01−d1d1d2dd2dd1=112,

因此 Cov(D1,D2)=1/12−1/9=−1/36,相关系数为 −1/2。两段不可能都很长,联合几何给出可量化的负依赖。

最长一段超过 1/2 的事件也容易算。任意一段超过 1/2 的概率为 (1/2)2=1/4,而两段不能同时严格超过 1/2,所以

P(max(D1,D2,D3)>12)=34.

均匀切点不等于独立均匀段长 ​

若独立抽三个 Vi∼Unif(0,1) 再归一化,虽然结果也正且和为一,却一般不是均匀单纯形分布。归一化后的密度取决于允许总长度的范围,偏向某些比例。上面的指数构造之所以正确,是因为其联合密度在分离出总和后恰好不再依赖比例。

若原数据来自一般连续总体,应用真实分布函数 F 后的 F(Xi) 才均匀;原始数值间距 X(i)−X(i−1) 通常没有本页分布。概率坐标间距与原始单位间距必须分开。

推论与应用

排序与间距变换的 Jacobian ​

独立均匀点的联合密度为一。除去概率零的并列情形,每个有序向量恰对应 n! 种原始标签顺序,所以排序后的联合密度在 0<u1<⋯<un<1 上为 n!。

从 d 回到 u 的变换是 uj=d1+⋯+dj。其Jacobian矩阵为下三角、对角元全为一,行列式为一。于是密度保持 n!,支持从有序区域变成间距单纯形。

指数总和与比例怎样分离 ​

令 Ei=SDi(i≤n)、En+1=S(1−∑i=1nDi)。从 (D1,…,Dn,S) 到 (E1,…,En+1) 的 Jacobian 绝对值为 Sn。指数样本的联合密度 e−∑Ei 因而变成

e−ssn1{s>0, di>0,∑di<1}=(sne−sn!)(n!1{di>0,∑di<1}).

它分解为总和密度与比例密度,既证明间距表示,也证明比例向量与指数总和独立。对 s 积分使用 ∫0∞sne−sds=n!,归一化可直接核验。

一般边缘与最长空隙 ​

由单纯形切片,Di 的密度为 n(1−d)n−1。因此

EDi=1n+1,Var(Di)=n(n+1)2(n+2),Cov(Di,Dj)=−1(n+1)2(n+2) (i≠j).

最后一式也可用交换对称性与 Var(∑iDi)=0 得出。

若指定 j 段都超过 x,先从每段减去 x,剩余单纯形的总长度为 1−jx,体积按其 n 次方缩放。对这些事件应用有限概率形式的容斥原理,得到最长间距 Dmax 的完整分布:

P(Dmax≤x)=∑j=0n+1(−1)j(n+1j)(1−jx)+n,x>0.

它能量化随机布点留下大空隙的概率,而不是只报告平均间距 1/(n+1)。

参考资料
关系图谱19 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
分类位置

上位 / 更一般

下位 / 直接特例

暂未标注直接特例。

类型化关系