[论文解读] WARN: an R package for quantitative reconstruction of weaning ages in archaeological populations using bone collagen nitrogen isotope ratios

Takumi Tsutaya, Minoru Yoneda|arXiv (Cornell University)|Apr 9, 2013

Forensic Anthropology and Bioarchaeology Studies参考文献 56被引用 3

一句话总结

本文介紹了 WARN，一個開源 R 套件，可利用骨膠原氮同位素比值（δ¹⁵N）對考古學中的青少年個體進行定量、機率性的斷奶年齡重構。透過整合新估算的青少年骨膠原更新速率與近似貝葉斯計算（ABC）框架，WARN 提供關鍵參數（如斷奶開始/結束年齡、母體至嬰兒 δ¹⁵N 增強值、斷奶食物 δ¹⁵N 值）的後驗機率與可信區間，克服了先前視覺評估方法的主觀性。

ABSTRACT

Nitrogen isotope analysis of bone collagen has been used to reconstruct the breastfeeding practices of archaeological human populations. However, weaning ages have been estimated subjectively because of a lack of both information on subadult bone collagen turnover rates and appropriate analytical models. Here, we present a model for analyzing cross-sectional delta-15N data of subadult bone collagen, which incorporates newly estimated bone collagen turnover rates and a framework of approximate Bayesian computation. Temporal changes in human subadult bone collagen turnover rates were estimated anew from data on tissue-level bone metabolism reported in previous studies. A model for reconstructing precise weaning ages was then developed and incorporating the estimated turnover rates. The model is presented as a new open source R package, WARN (Weaning Age Reconstruction with Nitrogen isotope analysis), which computes the age at the start and end of weaning, 15N-enrichment through maternal to infant tissue, and delta-15N value of collagen synthesized entirely from weaning foods with their posterior probabilities. A precise reconstruction of past breastfeeding and weaning practices over a wide range of time periods and geographic regions could make it possible to understand this unique feature of human life history and cultural diversity in infant feeding practices.

研究动机与目标

• 解決缺乏客觀、定量方法來估算考古學青少年個體骨膠原 δ¹⁵N 資料斷奶年齡的問題。
• 克服先前模型忽略青少年骨膠原更新速率，並將關鍵參數視為固定值的限制。
• 發展一個統計框架，透過估算斷奶參數的後驗機率與可信區間來納入不確定性。
• 提供一個可重現、開源的工具（WARN），供研究人員分析考古學人群的橫斷面 δ¹⁵N 資料。

提出的方法

• 利用組織層級骨代謝研究資料，重新估算青少年骨膠原更新速率的時間變化。
• 建立一個考慮年齡依賴更新速率的青少年骨膠原 δ¹⁵N 變化動態模型。
• 納入可變參數：母體至嬰兒的 δ¹⁵N 增強值與斷奶食物的 δ¹⁵N 值作為未知數以進行估算。
• 應用近似貝葉斯計算（ABC）以計算斷奶開始/結束年齡、增強值與斷奶食物 δ¹⁵N 的後驗分配與可信區間。
• 將模型實作為使用者友善的 R 套件（WARN），並配備結果解讀的視覺化工具。
• 使用斯皮塔爾菲爾茲考古學人群資料集驗證模型，顯示關鍵參數具有穩定的後驗機率分配。

实验结果

研究问题

• RQ1 斯皮塔爾菲爾茲人群的斷奶開始與結束年齡範圍為何？這些估計的機率為何？
• RQ2 母體至嬰兒組織間的 δ¹⁵N 增強量為多少？此增強因子的可信區間為何？
• RQ3 在所研究人群中，完全由斷奶食物合成的膠原蛋白 δ¹⁵N 值估計為何？
• RQ4 年齡依賴的骨膠原更新速率如何影響斷奶年齡估算的準確性？
• RQ5 此模型在重構斷奶習俗方面，相較於先前的主觀或點估計方法，改善程度為何？

主要发现

• 模型成功估算出，在斯皮塔爾菲爾茲人群中，斷奶開始於 0.0 至 1.2 歲之間，結束於 1.2 至 2.0 歲之間的聯合後驗機率為 95.6%。
• 母體至嬰兒組織間 δ¹⁵N 增強值落在 1.6–2.4‰ 範圍的邊際後驗機率為 0.967。
• 由斷奶食物合成的膠原蛋白 δ¹⁵N 值落在 12.4–13.0‰ 範圍的邊際後驗機率為 0.975。
• 模型納入年齡特定的骨膠原更新速率，顯著提升了斷奶年齡估算的準確性，相較於假設更新速率恆定的模型。
• 該套件透過後驗機率分配提供完整的不確定性量化，使研究人員能評估重構斷奶參數的可靠性。
• WARN R 套件使跨多樣考古學與古人類學情境的青少年 δ¹⁵N 資料分析具備可重現性、透明性與統計穩健性。