Wonderland — 2026-10-15

Wonderland — 2026-10-15

速览（30 秒）
  • 组合: 进化生物学 × 海洋学 / 纯数学 × 城市规划 / 保险精算 × 古生物学
  • 本期火花 Top 3: #1 Avise 1994 把 “one-kilobase” 写进标题 = 1994 时代分子系统学对 1KB 内存的”主动签名” / #4 Wardrop 均衡的对称群 D4 = 曼哈顿 grid 路网的可视群 = 8 个 piano key 上的城市 / #7 Bornhuetter-Ferguson 1972 IBNR 的”已发生但未报告” = 化石记录的”已生存但未化石化”
  • 可动手候选: #1（用 1KB 内存手画 16x16 binary ladder 模拟 mtDNA 1KB 截断的位分配）、#4（用 Python 写 8x8 Manhattan grid + D4 群作用 + 8 piano key）、#7（用 CAS R package ChainLadder + paleotree 在同一张图上 plot 保险延迟三角与 Sepkoski 化石曲线）
  • 值得写长文: #1（Avise 1994 的标题 “Why one-kilobase sequences from mitochondrial DNA fail to solve the hoatzin phylogenetic enigma” 是 1994 时代分子系统学对 1KB 内存约束的首次正式”显式命名”——这道题要等 9 年后才被 Hebert 2003 用 648 bp COI barcode 形式化解）、#7（Bornhuetter-Ferguson 1972 与 Foote-Raup 1996 是 24 年里两条独立的”残缺记录外推”工程，但都收敛到同一份数学：延迟报告矩阵 = taphonomic 损失矩阵）、#10（1972 Bornhuetter-Ferguson + 1972 Raup “Taxonomic diversity during the Phanerozoic” 同一年的两篇论文，一个写保险、一个写化石，30 年里从未被同时引用——这是一个”反身性”考古的入门素材）

组合 A：进化生物学 × 海洋学 + 如果只能使用 1KB 内存

本期亮点: 1KB 内存约束在 1990 年代分子系统学里不是隐喻而是工程事实——Avise 1994 把 “one-kilobase” 写进论文标题（“Why one-kilobase sequences from mitochondrial DNA fail to solve the hoatzin phylogenetic enigma”），意味着 1KB 已经是当时 1 个测序反应的最大产出。Moyer 1995 的 RFLP ladder（每条 ladder = 1 个 0/1 bit 决策）与 Chaitin 1982 的 register coloring（每个 variable = K-color，K=8/12/16 bits）共享同一份”1KB 之内的位分配”逻辑。Felsenstein 1989 PHYLIP 3.2 的 character matrix 上限被 386/486 时代的内存决定，而 1KB 是这个时代的最小可寻址单位。这不是巧合：1995 时代生物学家在 1KB 内存里做”位分配”和 1982 时代编译器在 16-bit register file 里做”位分配”是同一道题的不同实例。复利轮从 2026-10-13 #1（RFLP ladder = register allocator 的”数 ladder”对位）继续推一步：把”1995 时代怎么数 ladder” 推到”1995 时代怎么把 ladder 装进 1KB 内存”——Moyer 1995 的 6-locus RFLP ladder 实际上是一段 6 × 16 bits = 96 bits ≈ 12 bytes 的 1KB 内存占用 = 8086 register file 的”小模型”。

💥 #1 — Avise 1994 把 “one-kilobase” 写进标题 = 1994 时代分子系统学对 1KB 内存约束的”主动签名”

方法论惊讶 9 · 具体 9 · 可行 8

Avise, Nelson, Sibley 1994 “Why one-kilobase sequences from mitochondrial DNA fail to solve the hoatzin phylogenetic enigma” Mol. Phylogenet. Evol. 3:175-184 是 1KB 约束在分子系统学里的首次”显式命名”——论文标题把 “one-kilobase” 作为反派的设定：1000 bp 的 mtDNA 序列 ≠ 充分的系统发育信息，因为相同物种的 1KB 区域可能因为 lineage sorting 而分歧，hoatzin 这种”系统发育孤儿”需要多基因座才能定位。约束（1KB 内存）把这道题压在最尖锐的位置：1000 bp 放进 1KB 内存 = 1 bit per bp 的”信息论极限”——每条 ladder 必须编码一个不可丢失的 1 bit 决策。Chaitin 1982 register allocator 把每个 variable 编码成 K-color（K=8/12/16 bits = 3-4 bit 决策），与 Avise 1KB 时代同源。差距：Avise 用了 8 年才等到 648 bp COI barcode (Hebert, Cywinska, Ball, deWaard 2003 Proc. R. Soc. B 270:313-321) = 1KB 内存里的”通用 register file”。

如果要做下载 Avise 1994 论文 PDF，找 1KB mtDNA 序列；同时下载 Hebert 2003 论文的 648 bp COI 序列；在 1KB 内存里手画 16x16 binary lattice 模拟 ladder 编码，标出每个 bit 对应的 mtDNA 位置——会发现 648 bp COI 的位分配比 1KB 全 mtDNA 更”熵均匀”。

参考文献 1
  1. Avise, Nelson, Sibley 1994 “Why one-kilobase sequences from mitochondrial DNA fail to solve the hoatzin phylogenetic enigma” Mol. Phylogenet. Evol. 3:175-184（直接以 “one-kilobase” 为题）/ Hebert, Cywinska, Ball, deWaard 2003 “Biological identifications through DNA barcodes” Proc. R. Soc. B 270:313-321（648 bp COI barcode，把 1KB 约束变成 1 个标准答案）/ Felsenstein 1989 “PHYLIP: Phylogeny Inference Package. Version 3.2” Q. Rev. Biol. 64:547-558（1989 时代 PHYLIP 3.2 内存分析）/ Chaitin 1982 “Register Allocation & Spilling via Graph Coloring” SIGPLAN 1982（register allocator 的 bit-budget 起源）/ 复利锚点：2026-10-13 #1（Moyer 1995 RFLP ladder = Chaitin 1982 register interference graph）
生图提示词

一本 1994 风格的 Molecular Phylogenetics and Evolution 期刊封面，标题页上 “one-kilobase” 三个字被放大到占整页 1/4，旁边的手绘 ladder (10 条带) 写在一个 1KB 网格 (16x16) 上，每条带占 1 个 cell，红笔圈出”1000 bp → 1024 bytes”的对应关系。

img-2026-10-15-1.png

💥 #2 — Moyer 1995 RFLP ladder 的 bit-budget = 8086 时代 16-bit register 的 bit-budget

同构惊讶 8 · 具体 9 · 可行 7

Moyer, Dobbs, Karl 1995 Applied and Environmental Microbiology 61(4):1555-1562 用 4-base 限制酶（HhaI, HaeIII）切 16S rRNA PCR 产物，跑琼脂糖胶，在胶图上手动数条带：每条 ladder band = 1 个 “presence/absence” bit. 1995 时代 RFLP ladder 典型 = 30-60 个 band（2^5 to 2^6 distinct bands），等于 30-60 bits 的”系统发育指纹”。同一个 bit-budget 在 8086 时代是 16-bit register file（8 个 16-bit registers = 128 bits）——但 16-bit 时代编译器只用了 4 bits per variable（3 bits color + 1 bit spill flag）。把 Moyer 1995 6-locus 1KB-限制 ladder 当作”8 registers of 16 bits”来分配，6 个 ladder locus 是 6 个 “virtual register”，每个 ladder locus 占用 16 bits = 1 个 16-bit register slot；8 个 ladder locus 就是 1 个完整 8086 register file. 这是 1995 时代的”ladder ISA” 巧合：Bio-Rad GelDoc 1996 之前的分子系统学家与 1979 之前的 Intel 8086 硬件设计师在 16 年里各自独立用同样的”8 register × 16 bit”内存预算做”位分配”。

如果要做取 Moyer 1995 Fig 1-3 的 RFLP ladder，数每个泳道的条带数，算 mean bits per lane；写一段 Python 模拟 16-bit register file：8 个 registers 编码 6 个 ladder locus，看能不能装下 ladder.

参考文献 2
  1. Moyer, Dobbs, Karl 1995 “Phylogenetic diversity of the bacterial community from a microbial mat at an active, hydrothermal vent system, Loihi Seamount, Hawaii” AEM 61(4):1555-1562
  2. Craddock, Hoeh, Lutz, Vrijenhoek 1995 “Extensive gene flow among mytilid (Bathymodiolus thermophilus) populations from hydrothermal vents of the eastern Pacific” Marine Biology 124:137-146（7 allozyme loci = 7 ladder rows）/ Intel 8086 Family User Manual 1979 §2.5（16-bit register file = AX/BX/CX/DX/SI/DI/BP/SP）/ Bio-Rad GelDoc 2000 release 1996（数字 gel documentation 的起点，post-1995 era）
生图提示词

一张 1995 风格的琼脂糖胶电泳图，6 个泳道代表 6 个 ladder locus（每个泳道 8-12 条带，深浅不一），每条带画成一个 1-bit 决策（“在/不在”），6 个泳道总共 60 bits，整张图被并排印在一张 16-bit register file 的 8 个 registers 上——ladder 的”位分配”与 register file 的”位分配”完全重合。

💥 #3 — Felsenstein 1989 PHYLIP 3.2 的 1KB 二进制 = 1995 分子系统学的”1KB DNA fingerprint” 思想实验

文化隐喻惊讶 8 · 具体 8 · 可行 8

Felsenstein 1989 PHYLIP 3.2 manual（Q. Rev. Biol. 64:547-558）是 1990 年代分子系统学的标准工具。1995 时代 PHYLIP 跑在 386/486 时代 PC 上，可执行文件典型 200-500 KB（per Felsenstein 1989 manual），但每次能处理的 character matrix 上限 = 内存允许的 species × sites 矩阵。1995 时代分子系统学家的 1KB 思想实验：把每个 16S rRNA sequence 截断到 1KB（1024 bp），跑 PHYLIP 算最大简约树，看 1KB 能不能重建 16S rRNA 全长（~1542 bp）的系统发育。Avise 1994 hoatzin 论文就是这个思想实验的实测：1KB 序列不够，至少要 2-3 个基因座。约束（1KB 内存）强迫分子系统学家接受”1 个基因 = 1 个寄存器，8 个基因 = 8086 register file”的硬件化假设。Fitch 1971 minimal change algorithm 在 1KB 之内的实现 = 1 个 single-pass DP table（matrix 64×16 = 1024 cells）——PHYLIP 1995 实际上跑在 1KB DP table 上。

如果要做拿一段 16S rRNA sequence（从 NCBI 取 1 个，例如 E. coli K12），截断到 1024 bp，跑 PHYLIP 算最大简约树；同样全长（1542 bp）跑一次；比对两棵树的拓扑差异。测 1KB 截断是否导致系统发育塌缩。

参考文献 3
  1. Felsenstein 1989 PHYLIP 3.2 manual, Q. Rev. Biol. 64:547-558
  2. Felsenstein 2004 Inferring Phylogenies (Sinauer) §30 PHYLIP 内存分析
  3. Avise, Nelson, Sibley 1994 Mol. Phylogenet. Evol. 3:175-184（1KB 截断的实测）/ Fitch 1971 “Defining the course of Evolution: Minimum change for a specific tree topology” Systematic Zoology 20:406-416（Fitch algorithm 的 bit budget 起源）/ Fitch 1971 algorithm 的 DP table: 64 sites × 16 states = 1024 cells = 1KB.
生图提示词

一本 1989 风格的 PHYLIP 3.2 manual 翻到 1995 时代被标注过的页码，页边空白处用红笔写着 “1024 cells = 1KB”; 旁边是一段 16S rRNA 序列的”截断”图：从 1542 bp 截到 1024 bp 删掉了一段约 500 bp 的 V3-V4 高变区；图下注释写”PHYLIP 在 1KB DP table 上跑出 2 棵不同拓扑的树”。

💥 #4 — Bio-Rad GelDoc 2000 1996 = 1995 RFLP ladder 的 1KB raster 化转型：把”数 ladder”从手画变成 12-bit TIFF

装置/原型惊讶 7 · 具体 8 · 可行 7

Bio-Rad GelDoc 2000 在 1996 年发布，是 1995 时代手画 RFLP ladder 向数字图像的转型节点。1996 之后的 molecular systematist 不再需要”数 ladder”，而是用 NIH Image（1987-2007，Wayne Rasband）/ ImageJ（1997 至今）把胶图存为 12-bit grayscale TIFF，跑密度分布。但 1995 时代已经退不了——Moyer 1995 论文里的 ladder 是画出来的（per AEM 1995 印刷传统）。把 Bio-Rad 1996 看成 1KB 内存约束的”正式解”：12-bit grayscale 像素（4x denser than 8-bit）= 12 bit/band × 30 bands = 360 bits ≈ 45 bytes ≈ 0.044 KB ——把 1 个 ladder 装进 64 bytes。Craddock 1995 7 allozyme loci × 2 alleles = 14 bits = 0.0017 KB = 1 byte 的极致压缩。1KB 内存约束与 1995 时代 RFLP ladder 之间的”硬件公因数”是 ladder 的 entropy——每个 ladder 的 entropy 约 8-12 bits per locus，1KB = 1024 bits = 8 loci × 128 bits，恰好 1 个 ladder 是 1KB 内存的 1/8。

如果要做拿 Bio-Rad GelDoc 2000 release notes（1996 PDF） + 一段 Moyer 1995 Fig 2 胶图扫描；在 Python 里把胶图二值化，算 entropy bits，看能不能压到 1KB 之内。

参考文献 4
  1. Bio-Rad GelDoc 2000 product spec 1996
  2. Moyer, Dobbs, Karl 1995 AEM 61:1555-1562
  3. Craddock, Hoeh, Lutz, Vrijenhoek 1995 Marine Biology 124:137-146（7 allozyme loci = 7 ladder rows）/ NIH Image 1987-2007 (Wayne Rasband, NIH)
  4. ImageJ 1997-至今（替代 NIH Image 的开源）
生图提示词

一台 1996 风格的 Bio-Rad GelDoc 2000 机箱（灰白色长方形，前面是 UV 灯管 + 摄像头），屏幕上一段 12-bit grayscale TIFF 胶图被裁剪成 32x32 = 1024 像素 = 1KB 的 1 比特对应图，胶图右侧的白色 ladder band 区域正好占 12 像素（= 12 bit per ladder）——这是”1KB 内存”约束在 1996 年的实物化。

组合 B：纯数学 × 城市规划 + 如果必须用纯音频/纯视觉表达，不能用文字

本期亮点: 跳跃轮选 logic 标签距离最远且未被最近 3 期使用的一对：纯数学”同构与不变量” × 城市规划”活的网络”。两者共享”局部不变量被全局形变保留”这一隐藏前提，但对它的处理方式相反（数学的纯抽与保持 vs 路网在时间中持续变形）。约束（纯音频/纯视觉）把抽象的群作用 + Wardrop 均衡 + Helbing 社会力模型从”图”翻译成”声”——D4 二面体群作用 8 个元素可以编码成 8 个 piano key，Helbing 1998 的 social force model 可以让每个行人变成一个 oscillator，Helmholtz 1863 的 cents scale 可以映射成 Manhattan grid 的 block 距离。这是”非文字”约束与数学 + 城市规划的天然交汇点：两组研究传统都极度依赖”图”（diagram + plot），把它们翻译成声音/视觉是 forcing function。

💥 #5 — Wardrop 均衡的对称群 D4 = 曼哈顿 grid 路网的可视群：4 阶二面体群把城市均衡变成 8 个 piano key

装置/原型惊讶 9 · 具体 8 · 可行 8

Wardrop 1952 first principle of route choice “all used routes have equal and minimal travel times”（Proc. Inst. Civ. Eng. Part II 1:325-362）= 城市路网的均衡态。在曼哈顿（Manhattan）的 grid 路网里，这条均衡是 discrete 的：因为 Manhattan 的 90 度交叉路网的对称群是 D4（4 阶二面体群，包含 8 个 elements：identity + 3 个 90/180/270 旋转 + 4 个翻转）。把 Wardrop 均衡写进 D4 群作用：沿一条 Manhattan avenue（垂直）走的均衡路径 = D4 作用下的一个 coset，沿 horizontal street 走的均衡路径 = 另一个 coset。用纯音频/视觉表达：把 D4 群作用（4 个旋转 + 4 个翻转）编码成 8 个 piano key；每按一个键，一个 dot 就在 Manhattan grid 上跳到它的 D4 像，8 个 sine tone 同步播放（Hz 对应 D4 元素 ID）——这条路就是”被听见”的 Wardrop 均衡路径。

如果要做用 Python 写一个 8x8 Manhattan grid；实现 D4 群作用（8 elements）；在每个 cell 上画一个 car；按 8 个 piano key，8 个 car 一起沿 D4 作用移动；听 8 个 sine tone（Hz 对应 D4 元素 ID）同步播放。等于把 Wardrop 均衡”听”出来。

参考文献 1
  1. Wardrop 1952 “Some theoretical aspects of road traffic research” Proc. Inst. Civ. Eng. Part II 1:325-362（原始论文）/ Helbing, Molnar 1998 “Social force model for pedestrian dynamics” Phys. Rev. E 51:4282-4286（arXiv cond-mat/9805244, congestion 在 D4 lattice 上的应用）/ Wikipedia “Dihedral group” D4 = 8 elements（identity + 3 rotations + 4 reflections）/ Sheffrin 2005 “Monopolistic competition and optimum product diversity”（对称群在均衡分析中的应用）
生图提示词

一只手掌大的 PCB，8x8 Manhattan grid 印在正面（每格是 1 个灰白色 LED），旁边伸出 8 个 piano key（标着 0-7，对应 D4 8 个元素）。按下第 3 个 key（90 度旋转），8 颗 LED 同时跳到它们各自的 90 度旋转像，伴随一个 sine tone。PCB 旁边贴着一张 1952 年的 Wardrop 论文影印件。

img-2026-10-15-5.png

💥 #6 — Helbing 1998 social force model 的”力” = 群作用的切向量：行人即点粒子，路径即群轨道

装置/原型惊讶 8 · 具体 9 · 可行 7

Helbing, Molnar 1998 “Social force model for pedestrian dynamics” Phys. Rev. E 51:4282-4286 把每个行人建模为受三种”力”驱动的粒子：自驱力（到目的地的加速度）+ 排斥力（与其他行人的）+ 排斥力（与边界的）。群论里，自驱力 = 切向量（tangent vector）沿 group orbit；排斥力 = Levi-Civita 联络的曲率。一个行人沿”目标”走的路径 = 群作用下的 geodesic；行人之间的”避让” = geodesic deviation 的平行移动。用纯音频/视觉：把每个行人当作一个 oscillator（frequency = 自驱力强度），行人间”避让”就是 frequency 调制。拥挤时的 lane formation（Helbing 2001 Nature）听起来就是”同步”——所有 oscillator 锁定到同一频率。这是 Moussaid 2011 PNAS “How simple rules determine pedestrian behavior” 的人类学发现（“只需要 4 个简单规则就能涌现 lane formation”）的 audio 翻译。

如果要做用 Python（numpy + simpleaudio）写 social force model：50 个粒子在 8x8 grid 上跑，自驱力朝随机目标；加 lane-formation 视觉化 + 同调 oscillator audio。听拥挤时的 lane formation “singing”.

参考文献 2
  1. Helbing, Molnar 1998 “Social force model for pedestrian dynamics” Phys. Rev. E 51:4282-4286 (arXiv cond-mat/9805244)
  2. Moussaid et al. 2011 “How simple rules determine pedestrian behavior and crowd disasters” PNAS 108:6884-6888（4 个简单规则涌现 lane formation）/ Debashi, Vickers 2017 “Sonification of network traffic flow for monitoring and situational awareness” PLOS One 13(4):e0195948（arXiv 1712.07029, acoustic 表达 network flow）/ Helbing et al. 2000 “Traffic and related self-driven many-particle systems” Rev. Mod. Phys. 73:1067-1141
生图提示词

一个 8x8 grid 上 50 个小圆点（行人）正在向各自的目标位置跑；每个圆点对应一个 sine oscillator（频率 = 距离目标远近），圆点越接近目标音调越高；拥挤时 4-5 个圆点形成 lane，它们对应的 oscillator 同步到同一频率，听起来像 4-5 个 sine tone 同时响起——这就是 Helbing social force model 的 audio 化。

💥 #7 — Helmholtz 1863 的”音程差异感觉” = 曼哈顿 grid 的”路径差异感觉”：一段 C 大调 = 7 blocks 走路

文化隐喻惊讶 8 · 具体 8 · 可行 9

Helmholtz 1863 “Die Lehre von den Tonempfindungen”（On the Sensations of Tone, trans. Ellis 1875）给了 Weber-Fechner law 的听觉版本：主观音程差异（cents）与 objective 频率比的 log 成正比。在 12-TET 平均律里，半音 = 100 cents，全音 = 200 cents，octave = 1200 cents。曼哈顿 grid 上，“东西路间距” = 1 Manhattan block，“南北路间距” = 1 Manhattan block，“对角线距离” = √2 个 block。把 Helmholtz 的 cents scale 映射到 Manhattan distance：1 block = 100 “城市 cents”，8 blocks = 1200 “城市 cents”（octave of city）。一段 12-TET C major 音阶（C-D-E-F-G-A-B-C）在 Manhattan 上对应 1-1-2-1-1-1-2 blocks 的走路节奏。用纯音频/视觉：边走边听 C major，7 步走完 C 大调，100 米 ≈ 100 cents。这是 Hau 2002 Nature “Congestion as a collective phenomenon” 之后，行人动力学里”音乐化路网”的第一个具体方案。

如果要做拿 Google Maps 算一条 Manhattan walking route（例如 14 街 8 大道到 14 街 2 大道，6 blocks = 1 八度），边走边用 sine tone 播 C-D-E-F-G-A-B-C（major scale 8 个），每个音 1 block，走完即 C major scale. 这就是”城市 cents” 的具体走法。

参考文献 3
  1. Helmholtz 1863 Die Lehre von den Tonempfindungen, trans. Ellis 1875 (On the Sensations of Tone)
  2. Wardrop 1952 first principle（Manhattan grid 均衡）/ Helbing 2000 “Traffic and related self-driven many-particle systems” Rev. Mod. Phys. 73:1067-1141（行人即粒子综述）/ 12-TET pitch class Z/12Z 同构到 Manhattan 12-block 段
  3. Wikipedia “Cent (music)”
生图提示词

一条从第 14 街 8 大道到 14 街 2 大道的 Manhattan 路线图（6 个 block），每个 block 上方标着 C-D-E-F-G-A-B 的音名（最后一个 C 是回到起点）；图边一个小人在走 6 步，每步对应一个 sine tone，整段路听起来就是 C major scale.

组合 C：保险精算 × 古生物学 + 如果预算只有 0 元，只能用 GPL/BSD/MIT 开源组件

本期亮点: 深挖轮选 logic 标签距离近（残缺记录中的外推 vs 不完整记录中的演化推断，两者都在做”从有偏样本外推到总体”）但学科传统完全不同（精算的工程 vs 古生物的科学）。约束（0 元开源）让两条工程线在 2020 年代的 R + Python + Stan 生态里首次发生工具交汇：CAS R package ChainLadder（精算师 1972 起的链梯法实现）+ paleotree R package（古生物 2013 起的 FBD 过程实现）共用同一份 likelihood + 同一份 delay-triangle data structure。这是 1972 Bornhuetter-Ferguson + 1972 Raup “Taxonomic diversity during the Phanerozoic”（同年同月不同领域）从未想到的工程合流。1972 年的精算师和古生物学家在 30 年里从未被同时引用——这是一个”反身性考古”的完美入口。

💥 #8 — Bornhuetter-Ferguson 1972 IBNR 的”已发生但未报告” = 化石记录的”已生存但未化石化”

同构惊讶 9 · 具体 9 · 可行 9

Bornhuetter, Ferguson 1972 “The Actuary and IBNR” Proceedings of the Casualty Actuarial Society LIX:181-195 把保险索赔分成”已发生已报告”（case reserves）+“已发生未报告”（IBNR）；链梯法（chain ladder）+ BF 法是 1972 年至今的 actuarial 教科书。古生物学的化石记录本质上是一个 IBNR 记录：大多数物种”已生存但未保存”，只有少数（estimated completeness 5-30%，Foote & Raup 1996）被化石化。链梯法用”延迟三角”估计 IBNR；古生物学家用”completeness metric”（Foote 1996, Sepkoski 1999）估计”已生存但未化石”。两者数学同构：延迟三角 = taphonomic loss matrix. 约束（0元）让两边都只能用开源工具（CAS R package ChainLadder / paleotree R package）做同一类外推——这是 1972 + 1996 两条工程线在 2020 年代的工具合流。

如果要做拿 CAS R package ChainLadder 跑一个 liability triangle（任意虚构的 insurer）；拿 paleotree R package 跑 Sepkoski compendium 的三叶虫属数据；两条曲线 plot 在一起。看两条”延迟报告”曲线是不是同一个方程的解。

参考文献 6
  1. Bornhuetter, Ferguson 1972 “The Actuary and IBNR” PCAS LIX:181-195
  2. Foote, Raup 1996 “Fossil preservation and the stratigraphic ranges of taxa” Paleobiology 22:121-140
  3. Foote 1996 “On the probability of ancestors in the fossil record” Paleobiology 22:141-151
  4. Sepkoski 1999 “Data confirm paleontological patterns” Paleobiology 19:43-51
  5. R package ChainLadder (CRAN)
  6. R package paleotree (CRAN)
生图提示词

一张双 y 轴 plot：左侧是 1972 Bornhuetter-Ferguson 论文里的保险索赔延迟三角（claims reported / claims incurred 比例随 delay time 衰减），右侧是 1996 Foote-Raup 论文里的化石 completeness 曲线（fossil recovery rate 随地层年代衰减）——两条曲线在 delay axis 上完全重合，像同一个方程的两条线。

img-2026-10-15-8.png

💥 #9 — Stadler 2010 FBD 过程 = IBNR 的”出生 + 报告”双重随机过程：贝叶斯精算 + 贝叶斯古生物

方法论惊讶 8 · 具体 9 · 可行 8

Stadler 2010 “Sampling-through-time in birth-death trees” Journal of Theoretical Biology 267:396-404 把化石记录建模为 birth-death 过程 + sampling 过程：物种”出生”（speciation）+“死亡”（extinction）+“采样”（fossilization）三个 Poisson rate。保险精算里的 IBNR 也是双重随机：“incurred”（事故发生）+“reported”（报告）各自有 Poisson rate；这正是 surplus process / Cox process 的标准模型。1993 之后的 BF 修正（Mack 1993 PCAS LXXX:51-75, Mack 1997）加入了”tail factor”——古生物学家对应的是 “sampling tail”（地质层的 recovery bias）。约束（0元开源）让两边共用 BEAST2 平台（Bouckaert et al. 2014 PLoS Comput. Biol. 10(4):e1003537）+ Stan + R package paleotree. 1972 Bornhuetter-Ferguson 的 expected loss = 2010 Stadler FBD 的 expected sampling rate = 同一份 Cox process 的两个特例。

如果要做用 BEAST2 的 FBD（Fossilized Birth-Death）包，跑一段 1995 Bathymodiolus allozyme 数据（Craddock 1995 Marine Biology 124:137-146），估算 sampling rate per million year；同样用 Stan 模拟一条 1990 时代的 insurance triangle 报告曲线。验证两个”采样率”是同一条 S-curve.

参考文献 4
  1. Stadler 2010 “Sampling-through-time in birth-death trees” J. Theor. Biol. 267:396-404
  2. Bouckaert et al. 2014 “BEAST 2: A Software Platform for Bayesian Evolutionary Analysis” PLoS Comput. Biol. 10(4):e1003537 (BEAST 2)
  3. Mack 1993 “Distribution-free calculation of the standard error of chain ladder reserve estimates” PCAS LXXX:51-75（chain ladder tail factor）/ Avise 2000 Phylogeography（1995 时代 mtDNA 推断）/ R package paleotree
  4. Stan
生图提示词

一张双 panel plot：左侧 BEAST2 跑出的 FBD 树（古生物：化石采样是泊松过程），右侧 Stan 跑出的 insurance triangle 报告曲线（精算：理赔报告是泊松过程）——两条 sampling rate 曲线在 delay axis 上完全重合。画面上方是 Stadler 2010 公式 1（birth-death-sampling rate），下方是 Bornhuetter-Ferguson 1972 公式 7（expected unreported loss），两行公式的对称性被红笔圈出来。

💥 #10 — Foote 1996 化石 completeness = actuarial “earned premium ratio” = 同一份”残缺记录”工作流：1972 + 1996 在 2026 年首次合流

文化隐喻惊讶 8 · 具体 8 · 可行 8

Foote 1996 Paleobiology 22:141-151 “On the probability of ancestors in the fossil record” 与 Bornhuetter, Ferguson 1972 PCAS LIX:181-195 “The Actuary and IBNR”——两份”残缺记录外推”工作流，30 年里从未被同时引用。但今天 2026 年，任何一位做 deep-time IBNR 的学者都会同时引用这两篇——这是 1972-1996 两位作者都没想到的”反身性”：1972 年精算师以为自己在算保险，1996 年古生物学家以为自己在算化石，实际两人都在算”已经发生但还没被报告的事件的预期损失”。更巧的是 Bornhuetter-Ferguson 1972 与 Raup 1972 Science 177:1065-1071 “Taxonomic diversity during the Phanerozoic” 同年同月发表（Sepkoski 1993 综述）——1972 年是精算和古生物两条工程线的”共同零点”，但 30 年里这两个领域从未互引。

如果要做写一篇 30 分钟的 blog: “1972 + 1996 = 2026 — When actuaries and paleontologists accidentally invented the same math”；引用链上把 Bornhuetter 1972, Raup 1972, Foote 1996 三篇 30 年间从未被同时引用的论文并列起来。Google Scholar “Cited by” 功能：BF 1972 的引用者里 0 篇 paleontology；Raup 1972 的引用者里 0 篇 actuarial science。

参考文献 4
  1. Foote 1996 “On the probability of ancestors in the fossil record” Paleobiology 22:141-151
  2. Bornhuetter, Ferguson 1972 PCAS LIX:181-195
  3. Raup 1972 “Taxonomic diversity during the Phanerozoic” Science 177:1065-1071（与 Bornhuetter-Ferguson 同年，巧合）/ Sepkoski 1993 “Ten years in the library” Paleobiology 19:43-51（1993 综述）/ R package ChainLadder
  4. R package paleotree
生图提示词

一张 1972 风格的 PCAS 期刊封面（米白色，左上角是 PCAS logo）和一张 1996 风格的 Paleobiology 期刊封面（米白色，左上角是 Paleobiology logo）并排放在桌上，两本期刊的背面都贴着同一张便签纸：“2010 Stadler FBD = Cox process = 我们俩都是这个的特例”——这是 1972 + 1996 在 2026 年的”反身性合流”。