# E0 校正的逐特征恢复估计量

> 正式口径：`E0_adjusted_Pstar_Qstar_v1`。旧的 raw-E8 P 结果仅保留用于历史审计，不能再作为正式恢复调用。

## 估计量

- `B = ZMF_E0 − GMF_E0`；
- `P* = (ZMF_E8 − GMF_E8) − B`：E0 校正的暴露终点偏离；
- `Q*_r = (ZMF_Rr − GMF_Rr) − B`：E0 校正的恢复期残余偏离；
- `S_r = Q*_r − P*`：E8 到恢复期的来源谱系差变化；
- 正式 IUT 为 `max(p_P*, p_S)` 后按预声明特征族 BH；方向性水平回归要求 `P*×S<0` 且 `|Q*|<|P*|`，节律要求 `dot(P*,S)<0` 且 `||Q*||<||P*||`。
- `R+8−R+4` 虽然代数上仍是原始组间差的阶段变化，但收缩方向必须用 `Q*_R4` 和 `Q*_R8` 判断。

所有阶段一次性进入同一个八群体格模型，含 centered sex 和 cell-specific 24 h sin/cos；标本为随机横断面粪便，不作同鼠恢复解释。

## 调用数

| 数据层 | 恢复阶段 | 特征数 | P*水平 | S水平 | 水平IUT | 方向+效应水平 | P*节律 | S节律 | 节律IUT | 方向节律 |
|---|---|---:|---:|---:|---:|---:|---:|---:|---:|---:|
| Alpha | R+4 | 3 | 1 | 2 | 0 | 0 | 3 | 2 | 2 | 1 |
| Alpha | R+8 | 3 | 1 | 1 | 1 | 1 | 3 | 2 | 2 | 1 |
| Function | R+4 | 2,838 | 1,119 | 1,204 | 558 | 30 | 2,209 | 1,359 | 1,116 | 1,068 |
| Function | R+8 | 2,838 | 1,119 | 1,795 | 711 | 32 | 2,209 | 1,652 | 1,220 | 802 |
| Species | R+4 | 6,673 | 2,926 | 3,315 | 1,495 | 259 | 4,960 | 4,082 | 3,213 | 2,838 |
| Species | R+8 | 6,673 | 2,926 | 4,434 | 1,870 | 322 | 4,960 | 3,924 | 2,816 | 1,892 |
| Metabolite | R+4 | 1,171 | 250 | 323 | 49 | 45 | 237 | 326 | 50 | 41 |
| Metabolite | R+8 | 1,171 | 250 | 249 | 33 | 31 | 237 | 240 | 38 | 31 |

## 正式输出

- 菌群/功能/alpha：`microbiome/outputs/*_recovery_e0_adjusted_*.csv.gz`；
- 代谢组：`metabolomics/outputs/recovery_E0_adjusted_*.csv.gz`；
- 本目录 `outputs/` 保存同一冻结结果与 validation。

旧 `*_recovery_pqs.csv.gz` 和 `recovery_P_S_Q_*` 表使用 raw E8 gap，只能做历史数值审计，不得用于新图、候选排序或正式结论。