总览 · ← 返回运行 20261003-043412-search-t2-heart-interp-g24-D-s1

节点 n20

在 T2HI-01 插值上把 log-线性 target RMS 乘压实因子 K^(4t(1−t))（K=0.62，端点为1、中点最大），修正心脏窗口尺度过冲；shape_scale 56.2→71.9，A半 62.6→66.5。

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

来自 programs.git 6bdac83dcc:solution/METHOD.md

target_rms = log_interp(rms_a, rms_b, t, SCALE_DAMP)
contract = CONTRACT_K ** (4.0 * t * (1.0 - t))   # 1 at endpoints, K at t=0.5
target_rms = target_rms * contract

调研员的计划

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

对比：父节点版本 88dd0156db。改动的文件：solution/METHOD.md +48 −29、solution/run.py +10 −0

diff --git a/solution/METHOD.md b/solution/METHOD.mdindex 15ccc50..3476dca 100644--- a/solution/METHOD.md+++ b/solution/METHOD.md@@ -1,49 +1,68 @@-按 PLAN 在 T2HI-01 输出坐标端实现确定性 PCA 规范化（减质心→PCA 旋转→det=+1→三阶矩定向）并做 ON/OFF 对照；对照显示形状对齐已旋转不变、规范化仅微伤 occupancy，按 PLAN 判据关闭提交，产物与父节点12逐位一致。+在 T2HI-01 插值上把 log-线性 target RMS 乘压实因子 K^(4t(1−t))（K=0.62，端点为1、中点最大），修正心脏窗口尺度过冲；shape_scale 56.2→71.9，A半 62.6→66.5。  ## 方法族与实现 -PLAN 指定家族 **T2HI-01**（对齐 → logit 组成插值 → 共有类型内 k-NN 簇抽样 + 逐细胞表达凸组合，坐标 hold），family 机制全部保持打开（`T2HI_INTERP=1`, `T2HI_EXPR=1`, `T2HI_COORD=hold`, K=25/C=10，父节点12 已验证的 matched-align 保持 OFF）。+PLAN 家族 **T2HI-05**（压实尺度校正）移植到 **T2HI-01** 基座（父节点14：对齐 → logit 组成插值 → 共有类型内 k-NN 簇抽样 + 逐细胞表达凸组合，坐标 hold）。其余路径逐位不变。 -本节点新增 PLAN 机制 **canonicalise**（`canonicalize_coords`，`solution/run.py`）：在最终写文件之前对输出坐标做一次-1. 减质心；-2. 对称 `eigh` 求协方差特征分解，旋转到 PCA 主轴（特征值降序）；-3. 若旋转矩阵 det=-1，将最后一轴乘 -1 强制 det=+1；-4. 逐轴按第三中心矩符号定向（正偏度指向 +轴）。+改动仅一处（`solution/run.py`）： -只在三个最终输出点各应用一次（interp / control-mix / 单阶段 fallback），不影响对齐、插值、抽样等任何中间步骤。开关 `T2_CANONICAL`，**提交默认 0（关）**；设 1 可复现机制对照。+```python+target_rms = log_interp(rms_a, rms_b, t, SCALE_DAMP)+contract = CONTRACT_K ** (4.0 * t * (1.0 - t))   # 1 at endpoints, K at t=0.5+target_rms = target_rms * contract+```++- `CONTRACT_K` 由环境变量 `T2_CONTRACT_K` 控制，**提交默认 0.62（机制打开）**。+- bump 形因子只依赖视图内数据：t 来自 `interp_bracket`（相对时间差），rms_a/rms_b 现场从两个输入阶段坐标计算；无阶段名/细胞数/绝对时间硬编码。+- 单阶段退路：fallback 分支在 target_rms 计算之前返回，机制自然关闭；即便走到公式，t=0 → 因子=1。+- 该因子对 ca/cb 两个对齐云和最终输出统一生效（hold 坐标模式下的输出云被 `scale_to_rms` 拉到同一 target_rms）。  ## 机制生效证据（实际改变了什么） -机制确实激活并改变输出（非空转），且为纯刚体变换：+ON（K=0.62）vs OFF（K=1.0），proxy 视图 seed 0，官方评分器 A 半： -- **坐标全部改变**：ON vs OFF `spatial_3D` 不相等（整体旋转+反射作用于每个细胞，位移随位置而异）；规范化后质心 ≈ (0,0,0)（~1e-5），协方差非对角元 ≤1.3e-5（对角化），det=+1，三轴第三矩均为正（305k / 624k / 84k）——PLAN mechanism_evidence 的四条几何性质全部实测成立。-- **表达完全不变**：ON vs OFF 的 `.X` `array_equal=True`；RMS 不变（346.33）；k=5 近邻索引集合 `array_equal=True`（距离结构逐位保留）。-- **四组分变化**：只有 shape_scale 动了（56.21 → 55.97）；cell_state 65.73、expression_change 64.77、local_spatial 63.66 三组逐项不变——符合刚体变换预期。-- 不是常数位移（每细胞位移不同、质心被减掉后旋转作用于全云），也不是组成重加权（细胞数、类型配额、表达值完全不变）。+| 量 | OFF（=父节点14） | ON | 变化 |+|---|---:|---:|---|+| 输出坐标 RMS | 346.3 | 218.9 | 缩到 proxy 真值 ~217 量级 |+| scale_log_ratio | 0.4669 | 0.0080 | 尺度偏差几乎消除 |+| shape_scale | 56.21 | 71.87 | **+15.66** |+| cell_state | 65.73 | 65.73 | 不变 |+| expression_change | 64.77 | 64.77 | 不变 |+| local_spatial | 63.66 | 63.66 | 不变（neighborhood_mmd 0.0549 逐项相同） |+| board (A半) | 62.59 | 66.51 | +3.92 | -## 对照结果（A 半，seed 0，官方评分器）+符合 PLAN mechanism_evidence 全部四条：RMS 降到真值量级、scale_log_ratio <0.05、shape_scale +15.7、其余三组逐位不变（均匀缩放不改表达、细胞配额和邻域拓扑）。不是常数位移（全局乘性缩放，各细胞位移∝其位置），也不是组成重加权（细胞数/类型配额/表达值与 OFF 完全一致）。 -| 配置 | board | shape_scale | local_spatial | cell_state | expr_change |-|---|---:|---:|---:|---:|---:|-| canonical ON | 62.53 | 55.97 | 63.66 | 65.73 | 64.77 |-| **canonical OFF（提交，=父节点12）** | **62.59** | **56.21** | 63.66 | 65.73 | 64.77 |+## 对照结果（mechanism_off_control） -**关闭对照**：`T2_CANONICAL=0`（提交默认）与父节点12输出**逐位一致**（`.X` 与 `spatial_3D` 均 `array_equal=True`），seed 0 两次运行逐位相同。+`T2_CONTRACT_K=1.0` 时因子对任意 t 恒为 1，公式退化为纯 log-线性：输出 RMS 346.3、评分与父节点14 提交态的 A 半分（62.59，见父 METHOD.md）逐项一致，四组分数与全部 metrics 相同。ON 显著优于 OFF（+3.92 A半），按 PLAN 判据（shape_scale 提升 >1 分）**提交 ON（K=0.62）**。 -**结论：机制生效但无收益，按 PLAN 判据回退关闭。** 细分指标给出了机理：ON/OFF 的 `d2_shape`（0.02542）与 `scale_log_ratio`（0.4669）**完全相同**——评分器的形状对齐本身已对旋转/手性不变，k026 假设的"旋转方差源"在 shape 分量上不存在；规范化唯一改变的是 `occupancy_dice`（0.8247 → 0.8207，占据栅格依赖绝对朝向，PCA 规范帧不如 Procrustes 对齐到 a 阶段的原始帧），使 shape_scale 净降 0.24（噪声内但方向为负）。PLAN 成功判据「shape_scale 提升 > 1」未满足 → 提交 OFF。+## K 网格（A半，seed 0） -## 已验证 / 未验证+| K | board | shape_scale | scale_log_ratio |+|---:|---:|---:|---:|+| 0.50 | 64.27 | 62.93 | -0.199 |+| 0.55 | 65.18 | 66.58 | -0.107 |+| 0.60 | 66.30 | 71.05 | -0.023 |+| **0.62** | **66.51** | **71.87** | **0.008** |+| 0.64 | 66.04 | 69.99 | 0.039 |+| 0.70 | 64.96 | 65.68 | 0.125 |+| 1.00（OFF） | 62.59 | 56.21 | 0.467 | -**已验证**：(1) proxy 全流程 ~2s 纯 CPU，`vec-check` ok；(2) 提交默认与父节点12逐位一致、对 seed 确定（`np.random.default_rng(seed)`）；(3) ON/OFF 对照如上，机制的四条几何性质逐条实测；(4) 视图无关：`canonicalize_coords` 只用输出坐标本身（数据驱动、无绝对时间、无路径/board/manifest 字段依赖），且提交默认 OFF 路径 = 父节点12（已证视图无关）。+峰值在 0.62，与节点13（mix 分支，独立代码系）验证值一致——两条谱系给出同一最优 K，降低单点过拟合嫌疑。occupancy_dice/d2_shape 全程不变（0.8247/0.02542）：缩放不改变形状与占据结构，shape_scale 的提升全部来自尺度项。 -**未验证 / 风险**：(1) occupancy_dice 对朝向的敏感性提示"换一个非 PCA 的规范朝向"可能有 ≤0.1 量级的波动，但父帧（Procrustes 到 a 阶段）已是实测更优朝向，无证据存在更好选择，未做搜索（避免过拟合 A 半噪声）；(2) 未在伪装视图实跑，仅代码层面论证；(3) shape_scale 的真正瓶颈仍是 target_rms（log-线性 346 vs proxy 真值 217，`scale_log_ratio`≈0.467），任何刚体变换都无法触及。+## 验证过 / 没验证 -## 生物学知识来源+验证过：+- 确定性：同 seed 重跑 `.X` 与 `spatial_3D` 逐位相同。+- 视图无关：伪装视图（全阶段时间 +1 天、manifest 键序反转、随机路径）输出与真实视图逐位相同；因子只依赖相对时间差 t。+- `vec-check` 通过（n=17616，在 [1000, 17616] 内）。+- 完整视图运行 ~2s / <1GB，远低于 limits。 -未引入外部文献/数据库/检索知识，未用保留阶段（E8.5/E10.5/E12.5、禁窗 (8.25,8.75) 与 9.5<E≤13.5）或保留基因型的任何信息；输入为视图给定的 E8.25_late 与 E9.5（均在允许范围）。PCA、特征分解、第三矩定向均为通用统计/几何方法。未读 `external/` 与 `prior/`。+没验证：+- **final 迁移**（真实括号 E8.25+E8.75→E8.5，t=0.5，因子=0.62）：proxy 只在 t=0.4 处校准过。若真实轨迹是平滑压实谷（用已发布 E9.5 尺度对 log-RMS 做二次拟合外推 E8.5≈256），K=0.62 给出的 172 会过压；若谷在 E8.5 附近则合适。保留阶段真值不可知，此为 PLAN 已登记的风险，本节点按 proxy 证据提交。+- B 半分数（本表全为 A 半）；+3.92 远超 T2 ~1 分噪声，方向应可重现。 -## 下一步+## 生物学知识来源 -1. shape_scale 已被证明对**姿态类**改法完全不敏感（matched-align、PCA 规范化双双证伪，d2_shape 恒为 0.02542）——不要再做任何对齐/朝向类尝试；唯一剩余杠杆是 target_rms 估计（如 RMS~log 细胞数 / 阶段间隔的视图内回归），但需警惕 proxy 非单调 dip（354→217→335）导致的过拟合。-2. occupancy_dice 是唯一随朝向变的形状子指标且现帧更优；若未来改动细胞空间分布（如簇参数），需连同 occupancy 一起复测。-3. cell_state/expr/local 前沿仍被父配置占满，突破需换机制（容量约束双向互匹配的表达插值等，见父 ANALYSIS）。+心脏发育 E8.25–E9.5 窗口内心管环化与心肌致密化使空间尺度非单调（先压实后生长）——通用形态发生机制知识（不针对保留阶段的定量测量）；修正量的校准通过官方允许的 vec-score 代理查分完成，程序运行时所有尺度量均从视图输入现场计算。diff --git a/solution/run.py b/solution/run.pyindex 46e31de..f82697e 100644--- a/solution/run.py+++ b/solution/run.py@@ -55,6 +55,13 @@ MATCH_WEIGHT = os.environ.get("T2_MATCH_WEIGHT", "count")  # count | corr | corr # shape_scale 56.21 -> 55.97. Submitted OFF per the PLAN criterion; set # T2_CANONICAL=1 to re-enable the implemented mechanism for comparison. CANONICAL = os.environ.get("T2_CANONICAL", "0") != "0"+# contract: multiplicative compaction correction on the log-linear target RMS,+# K^(4*t*(1-t)) — bump-shaped in the interpolation parameter t (1 at the+# endpoints, maximal at the midpoint). Corrects the systematic over-estimate of+# log-linear scale interpolation inside windows where the organ compacts+# non-monotonically (validated on node 13, mix branch). K=1 disables the+# mechanism (control). Set T2_CONTRACT_K to override.+CONTRACT_K = float(os.environ.get("T2_CONTRACT_K", "0.62")) MIN_TYPE_CELLS = 5 KNN_K = int(os.environ.get("T2_KNN_K", "25")) N_CLUSTERS = int(os.environ.get("T2_N_CLUSTERS", "10"))@@ -323,6 +330,9 @@ def run(view: str, out: str, seed: int) -> None:         aligned_a, aligned_b, info = align_pair(sa.coords, sb.coords, sa.labels, sb.labels, "procrustes")     rms_a, rms_b = rms_radius(sa.coords), rms_radius(sb.coords)     target_rms = log_interp(rms_a, rms_b, t, SCALE_DAMP)+    # compaction correction: bump factor K^(4t(1-t)), 1 at endpoints, K at t=0.5+    contract = CONTRACT_K ** (4.0 * t * (1.0 - t))+    target_rms = target_rms * contract     ca = scale_to_rms(aligned_a, target_rms)     cb = scale_to_rms(aligned_b, target_rms)     n_total = interp_count(sa.n, sb.n, t, lo, hi, COUNT_DAMP)

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

用到的知识库条目

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

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

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

分析员

ANALYSIS.json written

工程师

**做了什么**：按 PLAN 把节点13 的压实尺度校正 K^(4t(1−t)) 移植到 T2HI-01 插值分支（`solution/run.py`，`T2_CONTRACT_K` 默认 0.62，端点=1、中点最大，t 与 RMS 均从视图现场计算），其余路径逐位不变。

**查分结果**（A半，proxy）：K 网格 {0.50,0.55,0.60,0.62,0.64,0.70} 峰值在 0.62 → **66.51**（seed 1 复测 66.59），shape_scale 56.21→71.87（+15.66），scale_log_ratio 0.467→0.008，其余三组逐位不变；OFF 对照（K=1.0）= 62.59，与父节点14 一致，符合 PLAN 全部四条机制证据。已验证确定性、伪装视图不变（时间+1/键序打乱输出逐位相同）、vec-check 通过、运行 ~2s。

**下一步最值得试**：final 迁移风险——真实括号 E8.25+E8.75→E8.5 的 t=0.5 处因子=0.62 未经校准（二次拟合外推暗示可能过压至 172 vs 估计真值 ~256）；后续节点可用已发布 E9.5 尺度做视图内谷形拟合来替代全局常数 K。另一方向是父 ANALYSIS 建议的 expression_change/cell_state 机制（容量约束双向互匹配表达插值）。

调研员

Looking at the experiment table and parent analysis:

- Parent 14 (=12=10) has shape_scale 56.09 as weakest group
- Node 13 proved contraction factor K^(4t(1-t)) gives +15.65 shape_scale on the mix branch
- T2HI-01 branch still uses plain log-linear RMS (346 vs truth ~217, scale_log_ratio 0.467)
- All pose/orientation approaches exhausted (matched-align, PCA canonical both rejected)

The clear move: port node 13's proven contraction correction to this branch.

```json
{"name": "T2HI-01 移植 K^(4t(1-t)) 压实尺度校正",
 "motivation": "父节点14 shape_scale=56.09 是最弱组，瓶颈为 target_rms 过估（log-线性 346 vs proxy 真值 217，scale_log_ratio≈0.467）。节点13 在 mix 分支用 K^(4t(1-t))（K=0.62）将 shape_scale 从 55.93 提到 71.58（+15.65），但该修正尚未移植到 T2HI-01 插值分支。父 ANALYSIS 明确建议此方向。",
 "approach": "1) 在 run.py 的 target RMS 计算处（当前为 log_interp 或 log-线性），乘以压实因子 K^(4*t*(1-t))，K 由环境变量 T2_CONTRACT_K 控制，默认 0.62（节点13 已验证值）。t 为当前预测在两括号阶段间的插值参数（单阶段时 t=0，因子=1，退路自然成立）。2) 搜索范围：K∈{0.50,0.55,0.62,0.70}，以 vec-score 查 shape_scale 选最优；预期 0.62 附近最优。3) 若 proxy 查分 shape_scale 提升 >1 分则提交，否则回退 K=1（=父节点）。4) 全流程 <5s，无需小样本。",
 "expected_groups": ["shape_scale"],
 "risks": "1) T2HI-01 坐标为 hold 模式（真实对齐位置），压实因子通过缩放输入阶段间接影响输出坐标分布，效果可能略小于 mix 分支；若 shape_scale 提升 <1 分则回退。2) proxy 心脏尺度非单调（354→217→335），K=0.62 是针对中间时间点的，若 final 的 t 不同需确认公式仍适用。3) 尽早在第一次 vec-score 查分确认 shape_scale 是否 >57，否则立即停止。",
 "family_id": "T2HI-05",
 "mechanism": "心脏发育中间阶段存在非单调压实，log-线性 RMS 插值系统性过估尺度；乘以 bump 形因子 K^(4t(1-t)) 将中间时间点的目标 RMS 收缩至真实量级。",
 "vs_constant_shift": "不是每细胞加常数位移，也不是全局组成重加权；它是对目标 RMS 的乘性尺度修正，改变所有坐标的绝对散布范围（缩放），且修正量随插值参数 t 非线性变化（端点为1、中点最大），反映心脏发育的时变压实生物学。",
 "mechanism_evidence": "Engineer 对比 ON（K=0.62）vs OFF（K=1）：1) 输出坐标的 RMS 应从 ~346 降至 ~217 量级；2) scale_log_ratio 从 0.467 降至 <0.05；3) shape_scale 分组分应提升 >10 分；4) cell_state/expression_change/local_spatial 三组应不变（均匀缩放不改表达和邻域拓扑）。",
 "mechanism_off_control": "设 T2_CONTRACT_K=1.0（即 K^x=1 对任意 x），此时公式退化为纯 log-线性 RMS，输出应与父节点14逐位一致（.X 和 spatial_3D array_equal）。预期差别：ON 时 shape_scale 显著提升，OFF 时回到 56.09。"}
```

审查员

REVIEW.json written