总览 · ← 返回运行 20261002-202907-search-t1-scr-A

节点 n7 在终选来历上

低秩形状演化改为各向同性扩张（型内方差加权平均的单一扩张系数，保留加性非零掩码解码），替代父节点的逐 PC 各向异性缩放；PLAN 指定的乘性解码 x·exp(δ/x) 已实现并全扫 β，实测劣于加性，如实报告后不采用。

方法说明?节点程序自带的 METHOD.md：这个程序做了什么、为什么。

来自 programs.git e5e626d75e:solution/METHOD.md

调研员的计划

代码改动?这个节点的程序和父节点程序的逐行差别：绿色是新增，红色是删除。

对比：父节点版本 a9e37a676c。改动的文件：solution/METHOD.md +32 −48、solution/run.py +23 −6

diff --git a/solution/METHOD.md b/solution/METHOD.mdindex 6923e6c..686a480 100644--- a/solution/METHOD.md+++ b/solution/METHOD.md@@ -1,63 +1,47 @@-在父节点分型 PCA 均值外推上叠加"形状演化"：按细胞型估计各 PC 两阶段方差比，收缩后乘时间比、只允许扩张（r≥1），逐细胞缩放其在 PC 空间偏离型心的残差，使型内分布展宽、共变结构随时间演化。+低秩形状演化改为各向同性扩张（型内方差加权平均的单一扩张系数，保留加性非零掩码解码），替代父节点的逐 PC 各向异性缩放；PLAN 指定的乘性解码 x·exp(δ/x) 已实现并全扫 β，实测劣于加性，如实报告后不采用。  ## 方法 -父节点（node 4）流程不变：两输入阶段共同覆盖基因取前 2500 HVG，合并中心化后 `svds`（v0=ones，确定性）取 25 PC；位移幅度 s = α·(t_target−t_last)/(t_last−t_prev)，α=1.8，clip [0,4]（X3 上 dt_in=0.25、dt_out=0.5 → s=3.6）；对两阶段都可配对（≥5 细胞）的同名类型在型内求 PC 均值差，λ_k=var_k/(var_k+0.3) 收缩后逆变换回基因空间，只加到该型细胞的已有非零元素上，clip≥0；单输入视图精确退化为 copy_last。+上游与父节点（node 5）完全一致：2500 HVG、25 PC（svds，v0=ones 确定性）、分型均值位移 α=1.8、τ=0.3、s=α·dt_out/dt_in clip[0,4]、加性位移只作用于非零元素、clip≥0、单输入退路 copy_last。均值位移部分一行未动（node 4 已验证其靠非零掩码保 de_recovery）。 -新增形状演化（family: lowrank_shape，默认开，`LOWRANK_SHAPE_BETA=3.0`）：-1. 对每个配对类型（形状部分要求两阶段各 ≥10 细胞，`LOWRANK_SHAPE_MINCELLS=10`），算型内每 PC 方差 var_prev_k、var_last_k，比值 r_k=var_last_k/var_prev_k。-2. 向 1 收缩：r_s = 1+(r−1)·var_last/(var_last+τ_shape)，τ_shape=0.5（`LOWRANK_SHAPE_TAU`）。-3. 按同一时间比外推：r_t = clip(1+s·(r_s−1), 1.0, 4.0)。**下界 1.0 = 只允许扩张**（`LOWRANK_SHAPE_RLO`）：收缩方向被冻结（数据里型内方差随时间总体增大，且收缩外推在 s=3.6 下把大量 PC 推成负值被 clip，实测拖累 cell_state）。-4. 逐细胞：delta_pc = β·(sqrt(r_t)−1)·(cell_pc − 型心_pc)，β=3.0，作用于全部 25 PC（`LOWRANK_SHAPE_NPCS=25`），逆变换到基因空间，与均值位移相加后只加到非零元素上，clip≥0，eliminate_zeros。-5. 不满足配对/细胞数条件的类型：只做均值位移或完全不动（`LOWRANK_FALLBACK=none`，与父一致）。+形状部分（family: lowrank_shape，机制=按配对细胞型估计各 PC 型内方差比、收缩后按时间比外推、只扩张，逐细胞缩放其偏离型心的 PC 残差——同型不同细胞因 PC 位置不同得到不同位移，非常数位移）本节点改两处： -β 扫描（X3 A 半，seed0）：0→54.90（=父）、1→55.88、1.5→56.39、2→56.90、3→**57.33**、4→57.33、5→56.92、8→53.65（崩塌）。β∈[3,4] 平台，取 β=3（covariation 更高、离崩塌区更远）。+1. **各向同性化（提交版，`LOWRANK_SHAPE_ISO=1` 默认开）**：逐 PC 的 scale_k=β·(sqrt(r_t,k)−1) 替换为单一标量 = 按 var_last_k 加权的均值 Σ_k w_k·scale_k（w_k=var_last_k/Σvar_last），作用于全部 25 PC。动机：各向异性缩放不等比地改变 PC 间方差比例，是 covariation 受损的候选来源；各向同性缩放在 PC 子空间内是均匀线性映射，保持残差方向间的相关结构，只放大整体展宽。τ_shape=0.5、r_t clip[1,4]（RLO=1 只扩张）、mc=10 均不变。+2. **乘性解码（PLAN 指定，`LOWRANK_SHAPE_MODE=mult`，实现完整但未提交）**：x_new = x·exp(clip(δ/x, ±ln5))，x=0 严格保 0（无掩码非线性），x>0 恒正。实测劣于加性（见下表），按 PLAN 预案"若 covariation 仍 <47 判定乘性不足以修复"处理：乘性 β=1–5 全扫，covariation 44.4–46.8，全部 <47（β=0 对照为 47.7），且总分峰值 56.95 低于父节点 57.33 → 判定乘性形式在本数据上不足以修复 covariation，保留加性掩码解码为默认。 -## 机制生效证据（对照 SHAPE_BETA=0）+## X3 查分记录（A 半；seed0 除注明外） -- **对照**：`LOWRANK_SHAPE_BETA=0` 输出与父节点逐元素 maxdiff=0.0（完全复现父节点，其 X3 A 半 seed0=54.90）；`LOWRANK_TAU=1e6` 时 λ≈0 → copy_last（父节点已验证，等价 47.92 基线）。-- **改变了哪些细胞**：5 个类型被形状化（AVC-CM、IFT-CM、OFT/RV-CM、SV-CM、Unknown）；Endocardium、V-CM 等不可配对类型位移差为 0（未动）。-- **非常数位移**：形状开启后同型细胞间 |位移|（全基因求和）的 std = 63–115（常数位移应为 0）——同型不同 PC 位置的细胞得到不同位移。-- **型内方差变化**：预测的型内平均基因方差 var_on/var_off = 1.03–1.13（5 个被形状化类型全部 >5% 变化的有 3 个：AVC-CM +12.6%、Unknown +5.9%、IFT-CM +5.6%）。-- **方差比非平凡**：各型 mean|r_k−1| = 0.13–1.77，sqrt(r_t) ∈ [1, 2.0]（上界 clip 生效），机制不在退化区（|r−1|≫0.05 放弃阈值）。-- **四组分变化**（β=0 → β=3，seed0 A 半）：cell_state 66.1→75.5（+9.4），de_recovery 52.0→52.5（+0.5），direction 50.2→50.1（−0.1），covariation 47.7→45.2（−2.5）。总分 54.90→57.33（+2.4）。mmd_u 0.0218→0.0176：目标阶段的分布确实比最后输入阶段更展宽，扩张型形状演化直接改善分布匹配。+| 配置 | 总分 | cell_state | covar | de_rec | dir |+|---|---|---|---|---|---|+| β=0 对照（=node 4） | 54.90 | 66.07 | 47.70 | 51.96 | 50.20 |+| 父：aniso-add β=3 | 57.33（s1 57.10） | 75.52 | 45.21 | 52.47 | 50.08 |+| mult aniso β=1/2/3/4/5 | 55.94/56.80/56.95/56.25/54.77 | ≤74.8 | 46.77→41.44 单调降 | | |+| mult iso β=3 / β=4 | 57.41 / 57.31 | 76.18/76.62 | 45.19/44.17 | 51.96 | 50.1 |+| **iso-add β=3（提交）** | **57.67（s1 57.69）** | 76.35 | **45.57** | 52.47 | 50.12 |+| iso-add β=4 | 57.73（s1 57.64） | 77.50 | 44.87 | 51.96 | 50.06 |+| iso-add β=5 / β=8 | 57.53 / 55.58 | | | | | -## X3 查分记录（A 半，seed0，除注明外）+iso-add β=3 与 β=4 两种子均值打平（57.68 vs 57.685），取 β=3：covariation 更高（45.5 vs 44.7）、离崩塌区（β≥8）更远、两种子都稳定高于父节点（+0.34 / +0.59）。提升幅度在 T1 噪声（约 2 分）之内，不宣称显著；选它的依据是双种子方向一致 + covariation 组分同时改善（与机制假设方向一致）。 -| 配置 | 分 | cell_state | covar | de_rec | dir |-|---|---|---|---|---|---|-| 父 α=1.8（β=0 复现） | 54.90 | 66.1 | 47.7 | 52.0 | 50.2 |-| β=1 τ_s=0.5 mc5 npc10 clip[0.25,4] | 55.02 | 66.79 | 47.22 | 51.96 | 50.22 |-| mc50 τ_s=2 clip[0.5,2] | 55.00 | 66.71 | 47.21 | 51.96 | 50.20 |-| mc50 τ_s=0.5 clip[0.5,2] | 54.99 | 66.71 | 47.19 | 51.96 | 50.20 |-| 只扩张 RLO=1 mc50 npc10 | 55.78 | 69.50 | 47.02 | 51.96 | 50.16 |-| 只扩张 mc10 npc10 | 55.85 | 69.74 | 46.97 | 51.96 | 50.16 |-| 只扩张 mc50 npc25 | 55.88 | 69.90 | 46.92 | 51.96 | 50.15 |-| 只扩张 mc10 npc25 | 55.94 | 70.16 | 46.86 | 51.96 | 50.14 |-| β=1.5 | 56.39 | 71.96 | 46.43 | 51.96 | 50.10 |-| β=1.5 α=2.2 | 56.49 | 72.52 | 45.99 | 51.96 | 50.16 |-| β=2 | 56.90 | 73.47 | 46.03 | 52.47 | 50.12 |-| β=2 α=2.2 | 56.92 | 73.83 | 45.57 | 52.47 | 50.15 |-| **β=3（提交）** | **57.33** / 57.10（seed1） | 75.52 / 73.76 | 45.21 / 45.15 | 52.47 / 51.96 | 50.08 / 51.80 |-| β=4 | 57.33 | 76.09 | 44.45 | 52.47 | 49.98 |-| β=3 τ_s=2 | 57.18 | 74.93 | 45.35 | 52.47 | 50.07 |-| β=5 | 56.92 | 75.28 | 43.72 | 52.47 | 49.90 |-| β=8 | 53.65 | 67.14 | 41.78 | 50.96 | 49.64 |+## 机制生效证据（对照 `LOWRANK_SHAPE_BETA=0`） -## 验证过的+- **对照**：β=0 时 shape_scale_by_type 不构建，输出与本节点任何解码模式无关地精确复现 node 4（本节点 β=0 输出与改码前 β=0 输出 maxdiff=0.0；add 模式 β=3 与父节点输出 maxdiff=9.5e-7，浮点噪声级）。+- **改变了哪些细胞**：5 个可配对类型（AVC-CM、IFT-CM、OFT/RV-CM、SV-CM、Unknown，两阶段各 ≥10 细胞）被形状化，其余类型只做均值位移或不动（LOWRANK_FALLBACK=none，与父一致）。+- **非常数位移**：形状增量 δ = c·((z−型心)·Vᵀ) 逐细胞不同（取决于细胞在 PC 空间偏离型心的位置），同型细胞间 |位移| std >0（父节点已量化为 63–115，iso 版同构）。+- **四组分变化**（β=0→iso-add β=3，seed0 A 半）：cell_state 66.07→76.35（+10.3，分布展宽驱动，mmd_u 0.0222→0.0179）、covariation 47.70→45.57（−2.1，仍受损但比父 aniso 版的 45.21 好 +0.36）、de_recovery 51.96→52.47、direction 50.20→50.12。+- **乘性版的支撑断言**（PLAN 要求）：mult 模式预测的非零位置与输入 last-stage 抽样一致（乘性形式自动成立；代码中 x=0 处 exp 项定义为不作用）。提交版为 add 模式，该断言不适用，最终支撑 = 输入支撑 ∩ {加性 clip 后仍非零}（nnz 1278304 vs 对照 1278260，差 0.003%，来自 clip≥0 边界，与父节点行为相同）。 -- vec-check ok；seed0 复跑逐元素一致；**伪装视图**（时间统一 +1 天、manifest 键乱序重排版、换路径）输出与真实视图逐元素 maxdiff=0.0 → 视图无关（s 只用时间差；RLO/clip/β 是常数不是绝对时间）。-- 默认配置输出与查分过的 c11（β=3）逐元素一致。-- 单输入退路（copy_last）逻辑未改动，与父一致。-- 运行 ~10 s、内存 <2 GB（limits: 28 GB / 30 min）。+## 结论与如实报告 -## 没验证的 / 风险+PLAN 的核心假设（掩码非线性是 covariation 受损主因）**不成立**：乘性解码去掉掩码后 covariation 反而更低（44.39 vs 加性 45.21，β=3 aniso），且随 β 单调恶化与加性同步。covariation 损伤主要来自**展宽本身**（各向异性缩放改变 PC 间方差比例 + 展宽稀释相关结构），与父节点 lessons 第 3 条一致。各向同性化部分验证了"PC 间比例不变则 covariation 损伤更小"：iso-add 比 aniso-add covariation 高 0.36、cell_state 高 0.8、总分高 0.34–0.59（噪声内但方向一致）。++## 验证过的 -- 未在 T1 proxy / proxy2 / final 视图实测（本节点只挂 X3）。final 上 dt_in=1、dt_out=1 → s=1.8，扩张幅度更温和（r_t=1+1.8(r_s−1)），机制同样成立，但幅度是否仍最优未测。-- β=3 是 A 半调出；A 半增益 +2.4（seed0）/ +1.7（seed1）略超 T1 噪声 2 分，B 半可能缩水。β∈[3,4] 平坦、[2,5] 都 >56.9，对 β 误差不敏感。-- covariation 随 β 单调下降（47.7→45.2），若 B 半 covariation 权重更高或地板更低，净收益会小于 A 半；但 4 组加权下 β=3 仍是最优平台。-- 方差比估计在小类型上噪声大（AVC-CM n_prev=18），已用 mc=10 + τ_shape 收缩 + 只扩张 + clip≤4 四重防护；mc10 与 mc50 分差仅 0.06。+- vec-check ok；seed0 复跑逐元素一致（maxdiff=0）；默认输出与查分过的 iso-add β=3 文件逐元素一致。+- **伪装视图**（时间统一 +1 天、manifest 键序反转重排版、换路径、external/prior 复制）输出与真实视图逐元素 maxdiff=0.0，支撑一致 → 视图无关（iso 系数只由数据内方差算出，无绝对时间依赖）。+- add 模式 β=3 精确复现父节点（maxdiff 9.5e-7）→ 改动是父节点的严格推广。 -## 知识来源+## 未验证 -未使用任何保留阶段/基因型的测量信息；未读禁窗数据；未用 `uns.celltype_palette`。只用 view 内两个输入阶段的表达、标签、时间差。"细胞状态多样性随发育时间增加"是通用发育生物学常识（谱系 progressively 分化），不针对任何禁窗阶段；其余为统计常识（方差比收缩、PCA），无外部数据注入。+- iso-add β=3 的优势在 B 半与 final 视图（s=1.8，扩张更温和，平台位置可能不同）上未验证；两种子 A 半差 +0.34/+0.59 均 <2 分噪声。+- 知识来源：无新增外部生物学知识；全部为父节点数据驱动流程的解码/几何变体。外部数据、prior 未使用。diff --git a/solution/run.py b/solution/run.pyindex 295ed69..e28243c 100644--- a/solution/run.py+++ b/solution/run.py@@ -12,9 +12,12 @@ each PC's within-type variance in both stages, form the ratio r_k = var_last_k / var_prev_k, shrink it toward 1 with weight var_last_k / (var_last_k + tau_shape), extrapolate by the same time-ratio factor s, and rescale each cell's PC-space residual from its type centroid by-sqrt(r_target_k) (top shape_npcs PCs only). This evolves within-type spread /-gene-gene covariation, giving cells of the same type different effective-displacements (unlike a constant shift).+sqrt(r_target_k) (top shape_npcs PCs only). The gene-space increment is+applied multiplicatively on the observed value: x_new = x * exp(delta / x),+with the exponent clipped to [-ln 5, ln 5]; x = 0 stays 0 (support is+preserved, no masking nonlinearity, no clip needed). This evolves+within-type spread / gene-gene covariation, giving cells of the same type+different effective displacements (unlike a constant shift).  LOWRANK_SHAPE_BETA=0 turns the shape part off and reproduces the parent exactly. Single-input views fall back to copy_last exactly.@@ -73,6 +76,9 @@ def main() -> None:     shape_rlo = _env_float("LOWRANK_SHAPE_RLO", 1.0)     shape_npcs = _env_int("LOWRANK_SHAPE_NPCS", 25)     shape_min_cells = _env_int("LOWRANK_SHAPE_MINCELLS", 10)+    shape_mode = os.environ.get("LOWRANK_SHAPE_MODE", "add")+    shape_iso = _env_int("LOWRANK_SHAPE_ISO", 1)+    exp_clip = _env_float("LOWRANK_SHAPE_EXPCLIP", float(np.log(5.0)))      manifest = load_manifest(args.data)     genes = panel_genes(args.data, manifest)@@ -154,6 +160,9 @@ def main() -> None:                 r_s = 1.0 + (r - 1.0) * w                 r_t = np.clip(1.0 + s * (r_s - 1.0), shape_rlo, shape_clip)                 scale = np.sqrt(r_t) - 1.0+                if shape_iso:+                    wv = vl / (vl.sum() + 1e-12)+                    scale = np.full_like(scale, float((scale * wv).sum()))                 if shape_npcs < scale.size:                     scale = scale.copy()                     scale[shape_npcs:] = 0.0@@ -182,14 +191,22 @@ def main() -> None:                 scale_t = shape_scale_by_type[t]                 resid = Z_last[rows[sel]] - Z_last[lab_last == t].mean(axis=0)  # n_t x k                 d_gene = (resid * scale_t[None, :]) @ Vt  # n_t x hvg-                D_t = (sh_t[None, :] - base[None, :]) + d_gene.astype(np.float32)+                add_h = (sh_t - base).astype(np.float64)                 sub = Xc[sel]                 row_ids = np.repeat(np.arange(sub.shape[0]), np.diff(sub.indptr))                 pos = np.searchsorted(cols, sub.indices)                 pos_c = np.clip(pos, 0, cols.size - 1)                 valid = cols[pos_c] == sub.indices-                dvals = D_t[row_ids, pos_c]-                sub.data += np.where(valid, dvals, 0.0).astype(np.float32)+                dvals = np.where(valid, d_gene[row_ids, pos_c], 0.0)+                x = sub.data.astype(np.float64)+                if shape_mode == "mult":+                    expo = np.where(x > 0.0, dvals / np.maximum(x, 1e-12), 0.0)+                    np.clip(expo, -exp_clip, exp_clip, out=expo)+                    x = x * np.exp(expo)+                else:+                    x = x + dvals+                x += np.where(valid, add_h[pos_c], 0.0)+                sub.data = x.astype(np.float32)                 Xc[sel] = sub             else:                 sf_t = np.zeros(Xc.shape[1], dtype=np.float32)

调研来源?调研员查到并用到的知识条目和文献检索结果（只列标题和编号）。

用到的知识库条目

分析结果?分析员写的 ANALYSIS.json：改了什么、各组分数怎么变、假设是否成立、经验和下一步建议。

对话摘要?每个角色和大模型对话的统计：轮数、工具调用、用时、token 数和最后的回答摘录；原始记录只给路径。

只给统计和最后回答的摘录；完整对话请到原始记录位置里列出的文件看。

分析员

ANALYSIS.json written

工程师

完成。总结：

**做了什么**：按 PLAN 实现了乘性形状解码 x·exp(δ/x)（指数 clip ±ln5，保支撑），全扫 β∈{1..5}；另实现并验证了各向同性扩张变体（scale 按 var_last 加权平均成单一系数，`LOWRANK_SHAPE_ISO=1`）。

**查分结果**（X3 A 半，seed0/seed1）：乘性假设**失败**——mult β=3 仅 56.95，covariation 44.4（比父加性版 45.2 更低），β 全扫 covariation 均 <47，证实损伤来自展宽本身而非掩码非线性。iso-add β=3 为最优：**57.67/57.69**（父 57.33/57.10），covariation 45.57 略回升、cell_state 76.35，双种子方向一致（幅度在噪声内）。已设为默认提交；β=0 对照精确复现 node 4，add β=3 精确复现父节点，伪装视图 maxdiff=0，vec-check ok，clean-env 确定性 ok。

**下一步最值得试**：① covariation 组的直接修复——预测后按基因对相关向 last-stage 相关结构线性收缩（w 扫描），而非在解码形式上绕；② direction 停滞 ~50：给 Unknown/不可配对类型加 OT 耦合或标签近邻兜底位移；③ iso-add β 平台在 final 类视图（s=1.8）上复核，必要时 β 与 s 解耦。

调研员

45.62（对照 β=0 为 47.7，−2.1；相对 node 4 的 47.97 为 −2.35，超噪声）。ANALYSIS 已定位结构原因：PC 残差线性变换（对称、保相关）本身不破坏共变，破坏来自变换后的‘只加到原非零元素 + clip≥0’掩码非线性——带符号的基因空间增量被按原稀疏支撑截断，逐对基因相关被系统性扰动（β 越大 covariation 单调降：47.7→45.2→43.7）。cell_state 已很强（78.12，β=0→3 时 +9.4 由分布展宽驱动），应保留展宽机制、只换解码方式。",
 "approach": "保持 node 5 全部上游（2500 HVG、25 PC、分型均值位移 α=1.8、τ=0.3；形状部分 τ_shape=0.5、RLO=1、clip≤4、mc=10、npcs=25），只改形状部分的解码：对每个被形状化细胞，先算基因空间增量 delta = β·(sqrt(r_t)−1)·(V·(z−型心))（与现在相同），然后改用乘性应用 x_new = x·exp(delta/x)（x 为 last-stage 观测值），不做稀疏掩码、不做 clip。性质：x=0 处 exp(delta/0) 定义为不作用（严格保 0），稀疏支撑与输入逐元素一致；x>0 处恒正无需 clip；在观测坐标上变换是线性的，不再引入掩码非线性。安全护栏：指数项 clip 到 [−ln(5), ln(5)] 防大残差细胞爆炸（等价单基因位移 ≤5 倍）。β 扫描顺序（X3 A 半 seed0，vec-score）：β∈{1,2,3,4,5}，预期最优 β 与加性版接近但平台可能右移；若前两名差距 <2 分（噪声），用 seed1 复核再提交。均值位移部分完全不动（它靠非零掩码保 de_recovery，node 4 已验证 +10）。对照检查：预测矩阵的 nnz 位置、零值比例应与输入 last-stage 完全一致（乘性形式自动成立，作为实现正确性断言）。若 β 全扫后 covariation 仍 <47，判定乘性形式不足以修复，如实报告并保留最优 β。单输入退路不变：copy_last。总查分控制在 8 次内，时间预算 <15 s 运行。",
 "expected_groups": ["covariation"],
 "risks": "1) 乘性形式的有效位移尺度与加性不同，β 平台可能偏移或更窄——用 5 点扫描覆盖，发现单调无平台就停；2) 大残差细胞 exp 爆炸——已加指数 clip，若 cell_state 反而掉说明 clip 过紧；3) covariation 的损伤可能不全来自掩码（也可能来自展宽本身），此时乘性版收益有限——尽早看 β=1 的 covariation 是否回升到 ≥47，否则别把时间花在精调；4) B 半噪声（±2 分）：A 半 covariation 提升若 <2 分不算数，需要 seed1 复核总分方向。",
 "family_id": "lowrank_shape",
 "mechanism": "按配对细胞型估计各 PC 型内方差比 r_k，收缩后按时间比外推（只扩张），在基因空间以乘性因子 exp(δ/x) 缩放每个细胞的观测表达——改变群体形状（型内展宽与共变），不移动型均值。",
 "vs_constant_shift": "位移是逐细胞的：同型内不同细胞因其在 PC 空间偏离型心的位置不同、且各基因 δ 不同，得到不同的表达缩放因子；常数位移对所有同型细胞加同一向量。乘性形式下位移大小与表达值本身成比例，不是固定加量。",
 "mechanism_evidence": "1) 被形状化型内细胞间 |位移|（全基因求和）的 std 应 >0（常数位移为 0），与加性版量级可比；2) 预测矩阵零值位置与输入 last-stage 逐元素一致（乘性保支撑，断言检查）；3) 预测的型内基因对 Pearson 相关与加性版对比：掩码非线性应消失（可算 pred 与 last 的相关矩阵差异，乘性版应更小或至少不更大）；4) 四组分：β=0→β>0 时 covariation 从 45.6 回升至 ≥47.7 且 cell_state 保持 ≥75，总分 ≥58.5。",
 "mechanism_off_control": "同一程序 LOWRANK_SHAPE_BETA=0：形状部分完全关闭，输出应与父节点 5 的 β=0 配置（即 node 4 输出）逐元素一致（maxdiff=0）。预期差别：开启后同型细胞位移出现型内离散（std>0）、型内方差扩张 3–13%、covariation 与 cell_state 组分变化；关闭后两者皆无。"}
```