Skip to content

Report analytical (w/o J) uncertainty separately from full σ #1

Description

@FormosaRes

問題

AgeCalc 每個年代只報一個 σ_age(含 σ_J 的全誤差,目前以 2σ 顯示)。但 σ_J 是整管照射共用的系統誤差:J 偏一點,全管樣品年代一起等比例平移,不改變樣品之間的差。

所以拿全誤差去做同管樣品互比(例如 transect 沿線年代梯度 resolvable 與否)會把門檻灌水,可能把真實存在的年代差判成「分不開」。

範例(NO.65 9.77 Ma vs 假想同管 B 10.05 Ma,差 0.28)

  • 全誤差 ±0.28(含 J):B−A 的合成誤差 ~0.40 → 0.28 < 0.40,結論「沒有梯度」
  • 分析誤差 ±0.12(w/o J):合成 ~0.17 → 0.28 > 0.17,結論「B 確實較老」

同一組數據結論相反。馬遠溪 transect 全同管,梯度判斷應該用 w/o J 誤差。

該做什麼

同時報兩個誤差:age ± σ_analytical (w/o J) / ± σ_full (含 J),例如 9.77 ± 0.12 / 0.28 Ma (2σ)

  • 同管互比 / transect 梯度 → 用 analytical(w/o J)
  • 絕對年代、跨照射管、跟文獻比 → 用 full(含 J)
  • 跟 U-Pb 等其他定年系統比 → 再加 σ_λ、monitor age(之後再說)

實作

σ_J 的貢獻在規劃中的 uncertainty-budget 功能本來就要拆出來,所以這項併進去最省事:
σ_analytical = sqrt(σ_full² − σ_J貢獻²),其中 J 貢獻 = |∂t/∂J|·σ_J = (F/(λ(1+J·F)))·σ_J

動機 / 參考

對照 Pychron pipeline 發現它把 age error 分「w/o J」「含 J」「含 decay-constant」分開報,是 Ar/Ar 標準做法。McDougall & Harrison (1999)、Schaen et al. (2021) 都強調同源樣品互比要用 analytical-only error。

狀態

使用者 2026-06-12 決定先不動,先記錄。實作時可與 uncertainty-budget panel 一起做。

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions