Devlog · M1c/M2/M3a 联合批次（2026-10-01 夜间）
Source:vignettes/articles/devlog-2026-10-01-M1c-M2-M3a-night.qmd
ID83 终审（2×P0 → 已修）
- key4 缺
.min(5)封顶——(AC)⁹ 类 ≥7 拷贝串联 debug panic/release 错桶（测试恰好只扫 6 拷贝）；补 4 条 oracle 对拍 golden。 - release profile 编译失败（
label_ofcfg 门 vs debug_assert 类型检查）——CRAN 模式 R CMD INSTALL 会失败；release 门纳入验收电池。 P2：SparseSignatures 归属勘误（LS+LASSO 非 KL-L1L1——更新式是 Hoyer/Cichocki KL-MU，数学正确归属需改）、线搜索参照 divergence 声明、tol 域对齐（R validator 由 ≥0 改 >0 对齐内核）、ARD 早期 β 平台误剪声明、runif 开区间勘误。
ard.rs w_active P0（ID83 终审连带发现，非驯化教训的镜像）
剪除装配按列拼接写，消费方按行主读——k_est≥2 经 R 的 ms_ard 提取全部转置打乱（内核 1.0 vs 经 R 0.354；k_est=1 golden 恰好掩盖）。修复 + 双侧回归（Rust row-major 恒等/双存活提取；R 路径可分 96×12 ARD 恢复 cosine≥0.95）。
U-M3a-01 likelihood.rs（交付，1514 行）
多项/NB LL（零伪计数；NB size 锚 8/50/100）+ χ²₁ LRT（统计量钉 mSigAct v3.0.3-branch 原文；默认 p = ½·erfc(√(D/2)) Self–Liang 边界渐近；mSigAct 全上尾口径保留为 p_chisq1_upper，两口径逐位 =2×）+ Fisher 内点 SE + 清零判据（0.01/0.05/10 三锚）+ 自实现 ln_gamma（Lanczos）/erfc（级数+连分式）。211 测试双工具链：size α=0.05→0.0525（精确 0.0522）、NB→Poisson 极限单调、erfc 对 libm 锚 1e-11、变异复核（删 ½ 因子 size 测红）。偏离：内核签名 Result 化（契约 5）；Δ 等价口径（多项系数项）；ECM/引导采样属 U-M3a-02/03。
M2 关闭与遗留
批次 A/C 终审：4 golden 修复非驯化 ✓、错误协议 P1 已修、round-22 注释双重勘误（监控量=22；全格=7；中间的 17 为我方误改——三个数字三个出处，审核员实测裁决）。M2 验收 sigminer 对拍 PASS（恢复 0.99999+、互 best-match 0.99999+、auto-KL 独立选 K=3）→ ROADMAP M2 勾选就绪。遗留：cv=FALSE 预算级跳过的专用 FFI 面、Wilcoxon 诊断实现、sigminer run_xval 的 CANON 注释勘误（type-major 描述有误，实际 5’-flank-major）。
环境项
macOS 27.0 SDK tapi malformed file（clang 17 链接 msuiter.so 失败）→ SDKROOT=MacOSX26.5.sdk 绕过；1.71 工具链同样需要。CI 不受影响。rustc 1.71 clippy 的 needless_borrow（build.rs）已修。
P0 修复：presence 池化零分布（Binomial(½)-χ² 混合参照）— 2026-10-03
独立审核 P0 裁决：ms_test_presence 的池化 raw p 在池化 D 上评估单一 χ²₁ 半尾，仅 n=1 恰确；H0 下每样本 D_j ~ ½δ₀+½χ²₁（Self–Liang）独立求和 ⇒ 精确零分布 = Σ_{m=0}^{n} C(n,m)/2ⁿ·χ²_m（Binomial(n,½) 混合）。实测 H0 pass_bh 率（n=1/2/8/16）0.035/0.135/0.545/0.860 —— 真实队列（n≥10）会把一切签名判为”存在”，BH 不救。修法 = PI 裁决方案①（保留池化 estimand，修正参照）。
实现（src/rust/src/fit.rs）：新私有 pooled_presence_p(d, n) = 混合严格上尾（m=0 原子项 = [D<0]，对可容许 D≥0 恒 0）；χ²_m 上尾 chisq_survival = Q(m/2, D/2)（复用 engine ln_gamma；df=1 路由 engine p_chisq1_upper erfc 路径 ⇒ n=1 逐位退化为旧 Self–Liang 半尾；df≥2 级数/连分式 NR 路径，1200 迭代帽，quantile 锚 1e-10 钉测）。权重：n≤53 精确 dyadic C(n,m)/2ⁿ（恰表出，n=1 位级锚的载体）；n>53 ln Γ 递推 + 权重下限 e⁻⁷⁰⁰ 跳过（跳过质量 <1e-12，n=1500 钉测）。R 面 R/fit.R MsFit 方法同病同修：纯 R 孪生 .ms_presence_pooled_p（dbinom 权重 × pchisq 上尾；R 的 pchisq(x, df=0, lower.tail=FALSE) 在 x=0 返回 1，与严格上尾约定不符——原子项显式处理）；.ms_bh_adjust 未动；catalog 面走 Rust 自动修复。D13 口径：χ² 参照修正，非新统计声明——estimand、逐样本 lrt_p 列、BH/pass_bh 逻辑均不变。vignettes/articles/ffi-surface.qmd presence 行同步。
H0 校准矩阵（固定种子 0xCA1B，1000 reps，Poisson-limit NB，k=2 BH，缺位签名 pass_bh 率；Rust 内核，同 draws 前后对照）：
| n | 旧（χ²₁ 半尾） | 新（精确混合） |
|---|---|---|
| 1 | 0.058 | 0.058（逐位同） |
| 2 | 0.117 | 0.062 |
| 4 | 0.261 | 0.051 |
| 8 | 0.557 | 0.047 |
| 16 | 0.920 | 0.050 |
（审核实测 0.135/0.545/0.860 @ n=2/8/16 与本表旧列同现象、不同 draws。）R 面经 .ms_test_presence_rust（种子 9000+n，1000 reps）：0.044/0.057/0.056/0.047/0.043 @ n=1/2/4/8/16——两faces全部落 [0.02, 0.08]（α=0.05，±4SE 窗）。Power 对照：弱效应（means 520/480，h₁=50/1000）n=1: 0.346 → n=8: 0.901（池化收益，测试钉死 0xF0CE）。
测试：Rust 新增 5（survival 锚 9 点 + χ²₂ 闭式、n=1 位级锚 + n=2 mpmath 金标 + 单调性、权重和=1−2⁻ⁿ @ n∈{2..1500}、H0 校准矩阵、power 随 n 增长）；既有 presence_single_sample_matches_fit_lrt_columns（位级）与 n=1 smoke 原样通过。R 新增 3（孪生锚/金标/退化阈值、H0 校准矩阵、MsFit 面手推混合公式 1e-10）；p_raw ≤ 0.5 断言修为 ≤ 1（混合尾不再以 ½ 为上界）。
验收：cargo test –workspace 绿（stable + 1.71.0/SDKROOT=MacOSX26.5.sdk 双跑，234+234）；clippy –all-targets 对 fit.rs 零警告（工作树中无关未提交 WIP engine/src/assign.rs 有 3 条预存警告，不在本修复范围）；cargo deny 四节绿；Cargo.lock md5 零变化；devtools::test() 全绿（251 块 / 120,391 expectations / 0 failed 0 error 0 skipped）。遗留：ARCHITECTURE §5 MSU-Fit 行的 “χ²₁ LRT（Self–Liang 边界渐近）” 措辞待随下批次对齐为”逐样本 Self–Liang + 池化 Binomial(½) 混合”（本修复不动 ARCHITECTURE）。