Next Article in Journal
Autopsy by Imaging: The Last 10 Years
Next Article in Special Issue
Subadult Age Estimation Using the Mixed Cumulative Probit and a Contemporary United States Population
Previous Article in Journal
Enhanced Efficiency of PrepFiler™ BTA DNA Extraction Kit with PrepFiler™ Lysis Buffer Top-Up over Sample Remnants
Previous Article in Special Issue
Dental Age Estimation Standards for Hispanic Children and Adolescents in California
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Bayesian Approach to Estimating Age from the Auricular Surface of the Ilium in Modern American Skeletal Samples

1
Department of Anthropology, California State University Sacramento, Sacramento, CA 95819, USA
2
Sociology and Anthropology Department, University of La Verne, La Verne, CA 91750, USA
*
Author to whom correspondence should be addressed.
Submission received: 15 September 2022 / Revised: 4 October 2022 / Accepted: 7 October 2022 / Published: 19 October 2022
(This article belongs to the Special Issue Estimating Age in Forensic Anthropology)

Abstract

:
Age estimation from human skeletal remains is a critical component of the biological profile for unidentified decedents. Using a Bayesian approach, we examine two popular methods (Lovejoy–LJ, and Buckberry zand Chamberlain–BC) for estimating age from the auricular surface of the ilium. Ages of transition are generated from a modern Portuguese skeletal sample (n = 466) and are coupled with an informative prior from historic Spitalfields, London (n = 179) to estimate age in a sample of modern Americans from the Bass Donated collection (n = 639). The Bass collection was challenging to statistically model, potentially due to higher morbidity and mortality characteristics of the central southern United States. The highest posterior density ranges provide a realized accuracy between 84–89% for males and 85–91% for females using the LJ method, and a realized accuracy between 79–82% for males and 65–71% for females using the BC method. Both methods worked well for older individuals. Cumulative binomials showed that both methods significantly underperformed; however, results were better for the LJ method, which also showed lower bias. Reference tables for aging modern American samples are provided, and the data meet Daubert guidelines, i.e., legal criteria for acceptable scientific evidence in a court of law in the United States.

1. Introduction

An important component in building the biological profile for an unidentified decedent is the estimation of age-at-death. Adult age estimation is particularly challenging because it is based on the degeneration of the skeletal tissue, which is further affected by an individual’s life history. These biological and environmental processes result in considerable variation in age indicators, especially for middle to older adults.
Traditionally, skeletal age is estimated by scoring features on a reference sample where morphological changes in the skeleton are linked with known age and then applying these estimated ages to an unknown target sample or individual. In a severe critique, Bocquet-Appel and Masset [1] outlined the bias and inaccuracy inherent in this approach, where the age distribution of the target sample will reflect the age distribution of the reference sample, a condition which later came to be known as age mimicry. Boldsen and colleagues [2] proposed transition analysis, a method that reduces or eliminates the issue of age mimicry. When transition analysis parameters are combined with Bayesian statistics using a prior age distribution, the resulting age estimates more accurately reflect the senescence changes in the target sample while minimizing bias [2]. Bayesian statistics with transition analysis have been successfully applied to pelvic age indicators, including the pubic symphysis [3,4,5,6] and the auricular surface [7,8,9].
Lovejoy and colleagues [10] developed the traditional, standardized method for estimating age-at-death based on morphological changes to the auricular surface of the ilium. Their method (referred to here as the Lovejoy method, LJ), based on an eight-phase system, describes changes to several criteria (e.g., transverse organization, texture, porosity, etc.) and is considered advantageous because the changes manifest well into advanced years. While widely used across various geographic regions, most authors reported broad agreement when applying the Lovejoy method [11,12,13,14,15,16,17,18,19]. However, as one might expect, this traditional approach suffered from typical issues of age mimicry and some researchers struggled to assign a myriad of morphological changes to one stage. This latter issue was specifically addressed by Lovejoy and colleagues, who proposed that scholars should focus on the critical age indicators that best represent the aging process and use auxiliary features to adjust age estimates.
Buckberry and Chamberlain [20] used a different approach to revise the auricular surface aging method (referred to here as Buckberry and Chamberlain method, BC). They argued that the auricular surface features aged independently from one another and should be scored separately and not grouped together into a single phase. Thus, instead of combining various indicators into a single phase, they created a composite scoring system where five features (transverse organization, texture, microporosity, macroporosity, and apical changes) were independently scored and then combined to estimate a final composite age for the individual. They reported increased replicability, due in part to lower intra- and inter-observer error. Importantly, while results were presented using a traditional statistical approach, the authors also used Bayesian statistics to provide posterior probabilities for ten-year age ranges. Some authors [21,22] compared the LJ and BC methods and found BC markedly easier to apply and more accurate than the LJ method. In 2016, Hens and Godde [7] applied Bayesian methodology to both auricular surface methods, statistically modeling age in Portuguese males. They found that the application of transition analysis coupled with Bayesian statistics significantly improved age estimates for both methods, work that was later independently confirmed by Kim and Algee-Hewitt [23].
The main objective of this research is to use a Bayesian analysis to apply an informative prior probability distribution with transition analysis to estimate age-at-death in a modern American skeletal sample of both males and females. We compare the Lovejoy and Buckberry and Chamberlain methods. The resulting age estimates are tested for performance against a hold-out sample of Americans, an issue considered critical by Braga et al. [24]. This Bayesian approach allows for the estimation of accuracy and precision (i.e., bias) for each method. Our second objective is to provide the highest posterior density ages and highest posterior density age ranges at 75%, 90%, and 95% accuracy (i.e., coverages) that may be used by forensic anthropologists practicing in a legal setting who need to estimate the age of unknown individuals. These results meet the criteria outlined in Daubert v. Merrell Dow Pharmaceuticals, Inc. [25] which required rigorous testing and evaluation of forensic methods. Referred to as the Daubert guidelines, this case law outlines the criteria needed for the presentation of scientific evidence in a court of law. Specifically, the United States Supreme Court determined that in order for evidence to be considered reliable, it must be grounded in scientific methods and procedures, i.e., the technique has been subjected to peer review, error rates are reported, and the technique is generally accepted within the relevant scientific community. The 75% coverage would be acceptable for skeletal remains in archaeological contexts, where personal identification of decedents is highly unlikely or impossible and age estimation is provided as a foundation for further demographic or pathological interpretations. However, the 90% and 95% coverages meet the standards outlined for reliability in the Daubert guidelines where accuracy and precision in age estimates are essential.

2. Materials and Methods

2.1. Samples and Scoring

This research uses two contemporary, modern skeletal samples: (1) American males and females of European (i.e., White) ancestry derived from the William M. Bass Donated Skeletal Collection at the University of Tennessee Knoxville (Bass Collection) and (2) Portuguese males and females from the Luis Lopes Collection of the National Museum of Natural History and Science in Lisbon, Portugal (Lopes Collection). Both collections are associated with documented demographic characteristics. Background on collection history and composition are presented by Shirley et al. [26] and Cardoso [27], and the samples have been thoroughly described in the literature [3,4,7,8,28]. A third sample from Christ Church Spitalfields (Natural History Museum, London) dates to post-medieval London residents and is dated to 1646–1859 AD. This sample also has documented age, sex, and occupation for a portion of the collection [29] and served as the sample for Buckberry and Chamberlain’s [20] work. The Spitalfields sample has a similar age profile to our two modern samples. Table 1 depicts a breakdown of the samples by sex and age for each collection.
The American and Portuguese samples were scored by one experienced researcher (SMH) for the two auricular surface methods, LJ and BC, using protocols to maintain a blind study. Pathological specimens and individuals of unknown sex were not scored. All Portuguese pelves were scored across a period of two weeks, while the Americans were scored during one week. To account for intraobserver error, ten pelves from each collection were scored three times each on different days and there was perfect agreement. Additionally, due to the challenging nature of the LJ method, wherein one must combine numerous indicators into one stage, scoring followed guidelines outlined by Lovejoy et al. [10] where individuals exhibiting traits indicative of an advanced stage were placed into a higher phase/stage even when youthful characteristics were retained. We lament cases where authors have not reported their scoring techniques thoroughly in the literature.

2.2. Statistical Methodology

For the goal of this paper, the American data (Bass collection) serve as the target sample to derive and test age-at-death estimates from the two auricular surface aging methods. The Portuguese data (Lopes collection) are used to estimate the transition analysis (TA) parameters. The Spitalfields sample serves as the informative prior, from which we calculated hazard parameters using a Gompertz model. These Gompertz parameters are subsequently combined with the Portugal TA parameters in Bayes Theorem to produce point estimates and ranges of ages to be tested on the American sample. Plots of Kaplan–Meier survivorship calculated on American females and males with Gompertz curves from the Portuguese and Spitalfields data confirm the fit of these samples (Figure 1 and Figure 2). Selection of similar mortuary profiles generally ensures a good fit of the informative prior to the age-at-death structure in the TA sample [6]. However, Godde and Hens [3] reported some deviation from fit may be tolerated; a conclusion that was later confirmed by the independent work of other researchers [9,30].
A number of samples were available from which to choose an informative prior, either through published Gompertz parameters (Terry collection, LA County Coroner from Konigsberg et al. [6]), Knox County, Knox Cemetery, U.S. from census data in Godde [8], or the generation of Gompertz parameters from our own data with published Gompertz parameters (Torino, Portugal, Sardinia) [3,4,7,31], unpublished data of which we were granted access (Forensic Databank, St. Brides [32]), publicly available data (post-medieval England and Wales as described in Godde and Hens [3,4,7,33]), and the data published in Buckberry and Chamberlain [20]. Past attempts at modeling the American sample showed the available samples (except those from post-medieval Europe) were not a good fit and did not produce viable age ranges. Gompertz survivorship plots of the three European samples (Figure 3, Figure 4 and Figure 5) showed St. Brides and Spitalfields had promise as their fit to the Bass Collection was good, as was their relationship to the sample from Portugal. The decision to go with Spitalfields over St. Brides was driven by a slightly better fit of TA parameters.
The Gompertz model is expressed as:
h t = α 3 e β 3 t
S t = exp α 3 β 3 1 exp β 3 t
where h is the hazard, t is age-at-death, and S is survivorship; Equation (1) is the hazard and Equation (2) is the survival function. The Gompertz model parameters are α3 and β3. We did not shift our curve to a specific age. The transition analysis parameters were obtained from a cumulative probit model [6], which represent the age at which an individual transitions from, in this case, one auricular surface phase to the next. The probits were run for: (1) Lovejoy method phases, (2) each of the five indicators for Buckberry and Chamberlain, and (3) the composite scores as phases from Buckberry and Chamberlain. The transition analysis and Gompertz model parameters are combined in Bayes’ theorem [2]:
Pr a | c j = Pr   c j | a f a 0 Pr c j | x f x d x    
where f(a) is a probability density function (PDF) used to derive the highest posterior density (HPD), analogous to a frequentist mode, and highest posterior density region (HPDR), which is a range of values calculated at 75%, 90%, and 95% coverages. Coverage represents the “percentage of individuals expected to fall within the specified HPDR” [6]. Due to this property, as the coverages increase, the regions become larger. In this paper, we use the HPDRs as age ranges.
The accuracy of the HPDRs was tested by a cumulative binomial test using the documented age of the American samples as a measure of the accuracy of the HPDRs at the various coverages. A measure of accuracy and whether it is significantly different than the stated coverage is provided through this process. Finally, bias is estimated and reported with the following equation:
Σ   e s t i m a t e d   a g e a c t u a l   a g e N
All analyses were run in R [34] using originally composed scripts and programming written by Dr. Lyle Konigsberg (http://faculty.las.illinois.edu/lylek/, accessed on 25 June 2022).

3. Results

TA parameters for Lovejoy method phases in females and males are found in Table 2. Female ages-at-transition are spread almost across each decade of life until phases 7–8 at which point it jumps to 2 decades. In males, each age-at-transition until phases 6–7 represents a single decade of life. At phases 6–7, this changes to every 2 decades, similar to the females. The Buckberry and Chamberlain method TA parameters in Table 3, Table 4, Table 5, Table 6, Table 7 and Table 8 have relatively similar patterns across the sexes and show several components that distinguish between young and older adults (i.e., apex and macroporosity). Microporosity does not distinguish between major life stages (e.g., young, middle, and older adult), concentrating on discriminating with an approximately 10–12 year age range in both sexes. Texture ages-at-transition are spread from young–middle adult among the sexes, while transverse organization is spread from young–middle until stages 4–5 at which time it increases to the 80s. Finally, for both sexes the phases derived from composite scores extend relatively evenly until older adult.
The HPDR table for the LJ method phases shows similar ranges for females and males (Table 9). The regions increase with phases, which is expected due to increased heterogeneity as one ages. In comparison with the HPDRs generated from the BC method (Table 10), the Lovejoy method’s regions were narrower by phase. Both aging techniques detected older ages.
The cumulative binomials (Table 11) show that neither method’s accuracy met the corresponding coverage level; they all significantly underperformed. However, the realized accuracies for the Lovejoy method coverages were still good–excellent, especially for the 90% coverage, yielding approximately 91% accuracy in females and 89% accuracy in males. Biases indicate that the LJ technique overaged for both sexes, while the BC method underaged females and overaged males. In general, the bias was greater for the BC method. Taken together, the results from the LJ method technique were superior to BC—at least for the American target sample.

4. Discussion

Age estimation from the auricular surface of the ilium remains a popular method utilized by forensic anthropologists, with the LJ method garnering more than twice as many practitioners as the BC method in a recent survey [35]. The BC method is reportedly easier to apply and outperforms the LJ method when no statistical manipulations are applied [7,21,22]. The BC method should be the technique of choice when practitioners are unable to use Bayesian statistics. However, previous work on Portuguese males [7] strongly indicated that while both methods perform relatively equally when Bayesian statistics are applied, the LJ method showed significant improvement under this rigorous methodology compared with a traditional (i.e., non-Bayesian) approach. There is little reason for forensic anthropologists to apply the traditional LJ method without the use of a Bayesian approach. Additionally, the auricular surface of the ilium outperforms the commonly used Suchey–Brooks pubic symphysis method, providing narrower age ranges with better coverages (see ranges in Godde and Hens [3,4] for comparison).
This paper used Bayesian modeling to estimate age in modern Americans from both sexes and supports our previous work modeling age in Portuguese males [7]. When contrasting these two aging methods on Americans of both sexes, the LJ method showed narrower age ranges and lower bias, around 8%, compared with around 11% using BC. Additionally, the LJ method showed higher realized accuracy at 75%, 90%, and 95% coverages, and markedly so for females, where accuracy was between 15–20% higher using LJ. The HPDRs reported here control for age mimicry, allowing age estimates to be applied across different populations. Forensic anthropologists should feel comfortable using the HPDR tables presented here, along with our previously published values for Portuguese males [7] as look-up tables for their own data and avoiding non-Bayesian auricular surface aging, especially for the LJ method. Indeed, Konigsberg et al. [6] argue that population-specific age indicators are not as important as comparable age structures between populations.
Admittedly, we attempted to model age in the American sample previously but were unable to find a comparable sample with a good fit to use as an informative prior. The Bass Donated collection represents a sample of modern Americans drawn primarily from Tennessee, but also other regions in the central south of the United States; one would think comparable populations would abound. In a 2017 paper, Godde [8] noted the Bass Collection Gompertz survivorship was significantly different than a cemetery from the same city and the mortality data from the U.S. Census for the same county, state, and country (also visible in Figure 3, Figure 4 and Figure 5). Thus, we were challenged in finding a sample with a similar age distribution until we accessed the data from Spitalfields and St. Brides. When choosing a prior, it is customary to search for populations that are comparable in terms of time and geography. Why would a sample of post-medieval Londoners have provided such a good fit to modern Americans from the south? In particular, for females, the lower survivorship of around age 83 and over almost mirrors post-medieval Europe. Additionally, the roughly contemporary Portuguese sample, representing urban as well as rural-to-urban migrants, also matched with the modern Americans. An exploration of the literature on health and mortality provides some insights as to why it was so tough to fit a model to the Bass collection.
In general, mortality in the United States is higher than in countries with similar economies [36]. While health and mortality may vary due to numerous factors including socioeconomic status, sex, and ethnicity, geography remains a significant factor. Geographic inequality in adult mortality in the United States is greater than in western Europe [37] and the central-southern section of the United States is one of the most disadvantaged [38], especially since the mid-20th century [39].
The rural southern United States, especially Tennessee, Kentucky, Alabama, and Mississippi, has higher rates of morbidity and mortality compared with other rural and urban areas in the Northeast, Midwest, and Western United States [40,41]. In particular, Appalachia and the Mississippi Delta regions have the lowest life expectancy in the country and the highest mortality rates due to various health issues [42,43,44,45,46]. People from Tennessee and Appalachia are the main contributors to the Bass Donated collection [26].
Mortality rates in the last two decades in many southern states (including Tennessee) were 30–40% higher than other regions of the United States—translating to 3–4 fewer years of life expected at age 50 [36]. These values represent the effects of numerous environmental and health hazards, including: smoking [36]; lower bone mineral density with higher hip fracture rates [47]; nutritional disease [48]; growth stunting [49]; stroke [50,51]; and high helminth load, with up to 55% of people affected by endemic whipworm, despite campaigns to eradicate such infections [52]. Overall, these studies underscore the poorer health and increased mortality risk from structural inequalities in populations from the central south United States, which is the underlying population for the Bass Donated collection. The health and mortality of these U.S. southerners more closely matches that from post-medieval London and urbanizing Lisbon in the 1800s who would have been affected with similar environments [28].
Despite the challenges inherent in modeling the Americans, we are successful in providing highest posterior density ages and highest posterior age ranges for modern American samples. Tests on a hold-out sample show realized accuracy in the 80–90% range. The Bayesian approach is far superior to traditional (i.e., non-Bayesian) methods because it controls for age mimicry, allowing age estimates to be applied across different populations. Additionally, Bayesian methods provide prediction intervals and/or credible intervals. Some forensic anthropologists may find the wide age ranges disappointing and hesitate to apply them to their own casework. We would remind them that these estimates are far superior to the traditional auricular surface age estimates alone and these statistical approaches are necessary to meet Daubert criteria in American forensic science.

Author Contributions

Conceptualization, S.M.H. and K.G.; methodology, K.G.; software, K.G.; validation, S.M.H. and K.G.; formal analysis, K.G.; investigation, S.M.H. and K.G.; resources, S.M.H. and K.G.; data curation, K.G.; writing—original draft preparation, S.M.H.; writing—review and editing, S.M.H. and K.G.; visualization, K.G.; supervision, S.M.H. and K.G.; project administration, S.M.H. and K.G.; and funding acquisition, S.M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding; but was funded in part by an internal campus grant to SMH from the Office of Research and Sponsored Projects, California State University Sacramento.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data from this study was obtained from several sources: (1) published data are available in the works cited, (2) data collected by the authors are available upon reasonable request, (3) data provided by the Museum of London are publicly available at the WORD website (https://www.museumoflondon.org.uk/collections/other-collection-databases-and-libraries/centre-human-bioarchaeology/osteological-database, accessed on 25 July 2019) and by reasonable request to the Museum of London, and 4) other data sources’ availabilities are under the purview of the owners of the data.

Acknowledgments

Many thanks to the curators who allowed access to the William Bass Donated collection and the Luis Lopes skeletal collection. Further, thank you to the faculty and staff who provided access to FDB data and the Museum of London online database (for the St. Brides data). We could not have conducted any of this work without the significant contributions to the literature and valuable insight offered by Lyle Konigsberg, who served as graduate advisor extraordinaire to both authors many years ago.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bocquet-Appel, J.-P.; Masset, C. Farewell to Paleodemography. J. Hum. Evol. 1982, 11, 321–333. [Google Scholar] [CrossRef]
  2. Boldsen, J.; Milner, G.; Konigsberg, L.; Wood, J. Transition Analysis: A New Method for Estimating Age from Skeletons. In Paleodemography: Age Distributions from Skeletal Samples; Hoppa, R.D., Vaupel, J.W., Eds.; Cambridge University Press: Cambridge, UK, 2002; pp. 73–108. [Google Scholar]
  3. Godde, K.; Hens, S.M. Age-at-Death Estimation in an Italian Historical Sample: A Test of the Suchey-Brooks and Transition Analysis Methods. Am. J. Phys. Anthropol. 2012, 149, 259–265. [Google Scholar] [CrossRef] [PubMed]
  4. Godde, K.; Hens, S.M. Modeling Senescence Changes of the Pubic Symphysis in Historic Italian Populations: A Comparison of the Rostock and Forensic Approaches to Aging Using Transition Analysis. Am. J. Phys. Anthropol. 2015, 156, 466–473. [Google Scholar] [CrossRef] [PubMed]
  5. Kimmerle, E.H.; Konigsberg, L.W.; Jantz, R.L.; Baraybar, J.P. Analysis of Age-at-Death Estimation through the Use of Pubic Symphyseal Data. J. Forensic Sci. 2008, 53, 558–568. [Google Scholar] [CrossRef] [PubMed]
  6. Konigsberg, L.W.; Herrmann, N.P.; Wescott, D.J.; Kimmerle, E.H. Estimation and Evidence in Forensic Anthropology: Age-at-Death. J. Forensic Sci. 2008, 53, 541–557. [Google Scholar] [CrossRef] [PubMed]
  7. Hens, S.M.; Godde, K. Auricular Surface Aging: Comparing Two Methods That Assess Morphological Change in the Ilium with Bayesian Analyses. J. Forensic Sci. 2016, 61 (Suppl. S1), S30–S38. [Google Scholar] [CrossRef]
  8. Godde, K. The Use of Informative Priors in Bayesian Modeling Age-at-Death; a Quick Look at Chronological and Biological Age Changes in the Sacroiliac Joint in American Males. AIMS Public Health 2017, 4, 278. [Google Scholar] [CrossRef]
  9. Nikita, E.; Xanthopoulou, P.; Kranioti, E. An Evaluation of Bayesian Age Estimation Using the Auricular Surface in Modern Greek Material. Forensic Sci. Int. 2018, 291, 1–11. [Google Scholar] [CrossRef] [Green Version]
  10. Lovejoy, C.; Meindl, R.; Pryzbeck, T. Chronological Metamorphosis of the Auricular Surface of the Ilium: A New Method for the Determination of Adult Skeletal Age at Death. Am. J. Phys. Anthropol. 1985, 68, 15–28. [Google Scholar] [CrossRef]
  11. Murray, K.; Murray, T. A Test of the Auricular Surface Aging Techniques. J. Forensic Sci. 1991, 36, 1163–1169. [Google Scholar] [CrossRef]
  12. Martrille, L.; Ubelaker, D.H.; Cattaneo, C.; Seguret, F.; Tremblay, M.; Baccino, E. Comparison of Four Skeletal Methods for the Estimation of Age at Death on White and Black Adults. J. Forensic Sci. 2007, 52, 302–307. [Google Scholar] [CrossRef]
  13. Saunders, S.R.; Fitzgerald, C.; Rogers, T.; Dudar, C.; McKillop, H. A Test of Several Methods of Skeletal Age Estimation Using a Documented Archaeological Sample. Can. Soc. Forensic Sci. J. 1992, 25, 97–118. [Google Scholar] [CrossRef]
  14. Bedford, M.E.; Russell, K.F.; Lovejoy, C.O.; Meindl, R.S.; Simpson, S.W.; Stuart-Macadam, P.L. Test of the multifactorial aging method using skeletons with known ages-at-death from the grant collection. Am. J. Phys. Anthropol. 1993, 91, 287–297. [Google Scholar] [CrossRef]
  15. Falys, C.G.; Schutkowski, H.; Weston, D.A. Auricular Surface Aging: Worse than Expected? A Test of the Revised Method on a Documented Historic Skeletal Assemblage. Am. J. Phys. Anthropol. 2006, 130, 508–513. [Google Scholar] [CrossRef]
  16. Storey, R. An elusive paleodemography? A comparison of two methods for estimating the adult age distribution of deaths at late Classic Copan, Honduras. Am. J. Phys. Anthropol. 2007, 132, 40–47. [Google Scholar] [CrossRef]
  17. Rougé-Maillart, C.; Vielle, B.; Jousset, N.; Chappard, D.; Telmon, N.; Cunha, E. Development of a Method to Estimate Skeletal Age at Death in Adults Using the Acetabulum and the Auricular Surface on a Portuguese Population. Forensic Sci. Int. 2009, 188, 91–95. [Google Scholar] [CrossRef] [Green Version]
  18. Schmitt, A. Age-at-Death Assessment Using the Os Pubis and the Auricular Surface of the Ilium: A Test on an Identified Asian Sample. Int. J. Osteoarchaeol. 2004, 14, 1–6. [Google Scholar] [CrossRef]
  19. Hens, S.M.; Rastelli, E.; Belcastro, G. Age Estimation from the Human Os Coxa: A Test on a Documented Italian Collection. J. Forensic Sci. 2008, 53, 1040–1043. [Google Scholar] [CrossRef]
  20. Buckberry, J.L.; Chamberlain, A.T. Age Estimation from the Auricular Surface of the Ilium: A Revised Method. Am. J. Phys. Anthropol. 2002, 119, 231–239. [Google Scholar] [CrossRef]
  21. Hens, S.M.; Belcastro, M.G. Auricular Surface Aging: A Blind Test of the Revised Method on Historic Italians from Sardinia. Forensic Sci. Int. 2012, 214, 209.e1–209.e5. [Google Scholar] [CrossRef]
  22. Mulhern, D.M.; Jones, E.B. Test of revised method of age estimation from the auricular surface of the ilium. Am. J. Phys. Anthropol. 2005, 126, 61–65. [Google Scholar] [CrossRef] [PubMed]
  23. Kim, J.; Algee-Hewitt, B.F.B. Age-at-Death Patterns and Transition Analysis Trends for Three Asian Populations: Implications for [Paleo]Demography. Am. J. Biol. Anthropol. 2022, 177, 207–222. [Google Scholar] [CrossRef]
  24. Braga, J.; Heuze, Y.; Chabadel, O.; Sonan, N.K.; Gueramy, A. Non-Adult Dental Age Assessment: Correspondence Analysis and Linear Regression versus Bayesian Predictions. Int. J. Legal. Med. 2005, 119, 260–274. [Google Scholar] [CrossRef] [PubMed]
  25. Daubert v. Merrell Dow Pharmaceuticals, Inc., 509 U.S. 579; United States Court of Appeals for the Ninth Circuit: San Francisco, CA, USA, 1993.
  26. Shirley, N.R.; Wilson, R.J.; Meadows Jantz, L. Cadaver Use at the University of Tennessee’s Anthropological Research Facility. Clin. Anat. 2011, 24, 372–380. [Google Scholar] [CrossRef]
  27. Cardoso, H.F.V. Brief Communication: The Collection of Identified Human Skeletons Housed at the Bocage Museum (National Museum of Natural History), Lisbon, Portugal. Am. J. Phys. Anthropol. 2006, 129, 173–176. [Google Scholar] [CrossRef]
  28. Hens, S.M.; Godde, K.; Macak, K.M. Iron Deficiency Anemia, Population Health and Frailty in a Modern Portuguese Skeletal Sample. PLoS ONE 2019, 14, e0213369. [Google Scholar] [CrossRef] [Green Version]
  29. Natural History Museum, London. Human Comparative Collection. Available online: https://www.nhm.ac.uk/our-science/collections/palaeontology-collections/human-comparative-collection.html (accessed on 18 July 2022).
  30. Kim, J.; Algee-Hewitt, B.F.B.; Konigsberg, L.W. Inferring Age at Death for Japanese and Thai Skeletal Samples under a Bayesian Framework of Analysis: A Test of Priors and Their Effects on Estimation. Forensic Anthropol. 2019, 2, 273–292. [Google Scholar] [CrossRef]
  31. Hens, S.M.; Godde, K. New Approaches to Age Estimation Using Palatal Suture Fusion. J. Forensic Sci. 2020, 65, 1406–1415. [Google Scholar] [CrossRef]
  32. WORD. Wellcome Osteological Database. Available online: https://www.museumoflondon.org.uk/collections/other-collection-databases-and-libraries/centre-human-bioarchaeology/osteological-database (accessed on 27 January 2020).
  33. Godde, K.; Hens, S.M. An Epidemiological Approach to the Analysis of Cribra Orbitalia as an Indicator of Health Status and Mortality in Medieval and Post-Medieval London under a Model of Parasitic Infection. Am. J. Phys. Anthropol. 2021, 174, 631–645. [Google Scholar] [CrossRef]
  34. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  35. Garvin, H.M.; Passalacqua, N.V. Current Practices by Forensic Anthropologists in Adult Skeletal Age Estimation. J. Forensic Sci. 2012, 57, 427–433. [Google Scholar] [CrossRef]
  36. Fenelon, A. Geographic Divergence in Mortality in the United States. Popul. Dev. Rev. 2013, 39, 611–634. [Google Scholar] [CrossRef] [Green Version]
  37. Wilmoth, J.; Boe, C.; Barbieri, M. Geographic Differences in Life Expectancy at Age 50 in the United States Compared with Other High-Income Countries. In International Differences in Mortality at Older Ages: Dimensions and Sources; Crimmins, E., Preston, S.H., Cohen, B., Eds.; National Academies Press: Washington, DC, USA, 2010; pp. 333–366. [Google Scholar]
  38. Fenelon, A.; Preston, S.H. Estimating Smoking-Attributable Mortality in the United States. Demography 2012, 49, 797–818. [Google Scholar] [CrossRef] [Green Version]
  39. Ezzati, M.; Friedman, A.B.; Kulkarni, S.C.; Murray, C.J.L. The Reversal of Fortunes: Trends in County Mortality and Cross-County Mortality Disparities in the United States. PLoS Med. 2008, 5, e66. [Google Scholar] [CrossRef]
  40. James, W.; Cossman, J.; Wolf, J.K. Persistence of Death in the United States: The Remarkably Different Mortality Patterns between America’s Heartland and Dixieland. Demogr. Res. 2018, 39, 897–910. [Google Scholar] [CrossRef] [Green Version]
  41. Miller, C.E.; Vasan, R.S. The Southern Rural Health and Mortality Penalty: A Review of Regional Health Inequities in the United States. Soc. Sci. Med. 2021, 268, 113443. [Google Scholar] [CrossRef]
  42. Singh, G.K.; Siahpush, M. Widening Rural–Urban Disparities in Life Expectancy, U.S. 1969–2009. Am. J. Prev. Med. 2014, 46, e19–e29. [Google Scholar] [CrossRef]
  43. Wami, W.; Walsh, D.; Hennig, B.D.; McCartney, G.; Dorling, D.; Galea, S.; Sampson, L.; Dundas, R. Spatial and Temporal Inequalities in Mortality in the USA 1968–2016. Health Place 2021, 70, 102586. [Google Scholar] [CrossRef]
  44. Wang, H.; Schumacher, A.E.; Levitz, C.E.; Mokdad, A.H.; Murray, C.J. Left behind: Widening Disparities for Males and Females in US County Life Expectancy 1985–2010. Popul. Health Metr. 2013, 11, 8. [Google Scholar] [CrossRef] [Green Version]
  45. Kulkarni, S.C.; Levin-Rector, A.; Ezzati, M.; Murray, C.J. Falling behind: Life Expectancy in US Counties from 2000 to 2007 in an International Context. Popul. Health Metr. 2011, 9, 16. [Google Scholar] [CrossRef] [Green Version]
  46. Cullen, M.R.; Cummins, C.; Fuchs, V.R. Geographic and Racial Variation in Premature Mortality in the U.S.: Analyzing the Disparities. PLoS ONE 2012, 7, e32930. [Google Scholar] [CrossRef]
  47. Paulozzi, L.J. Does Inadequate Diet during Childhood Explain the Higher High Fracture Rates in the Southern United States? Osteoporos Int. 2010, 21, 417–423. [Google Scholar] [CrossRef] [PubMed]
  48. Park, Y.K.; Sempos, C.T.; Barton, C.N.; Vanderveen, J.E.; Yetley, E.A. Effectiveness of Food Fortification in the United States: The Case of Pellagra. Am. J. Public Health 2000, 90, 727–738. [Google Scholar] [PubMed] [Green Version]
  49. Bethony, J.; Brooker, S.; Albonico, M.; Geiger, S.M.; Loukas, A.; Diemert, D.; Hotez, P.J. Soil-Transmitted Helminth Infections: Ascariasis, Trichuriasis, and Hookworm. Lancet 2006, 367, 1521–1532. [Google Scholar] [CrossRef]
  50. Howard, G. Why Do We Have a Stroke Belt in the Southeastern United States? A Review of Unlikely and Uninvestigated Potential Causes. Am. J. Med. Sci. 1999, 317, 160–167. [Google Scholar] [CrossRef]
  51. Lanska, D.J.; Kuller, L.H. The Geography of Stroke Mortality in the United States and the Concept of a Stroke Belt. Stroke 1995, 26, 1145–1149. [Google Scholar] [CrossRef]
  52. Starr, M.C.; Montgomery, S.P. Soil-transmitted helminthiasis in the United States: A systematic review—1940–2010. Am. J. Trop. Med. Hyg. 2011, 85, 680. [Google Scholar] [CrossRef]
Figure 1. Kaplan–Meier survivorship (stepped lines) showing the probability of survival across adult ages and Gompertz survivorship curves (solid line) representing age-dependent death rates for males. The black line is Spitalfields informative prior showing an excellent fit to the model. The grey line is the Portuguese Lopes collection (used for transition analysis), which shows higher rates of survivorship for middle to older adults in the sample.
Figure 1. Kaplan–Meier survivorship (stepped lines) showing the probability of survival across adult ages and Gompertz survivorship curves (solid line) representing age-dependent death rates for males. The black line is Spitalfields informative prior showing an excellent fit to the model. The grey line is the Portuguese Lopes collection (used for transition analysis), which shows higher rates of survivorship for middle to older adults in the sample.
Forensicsci 02 00051 g001
Figure 2. Kaplan–Meier survivorship (stepped lines) showing the probability of survival across adult ages and Gompertz survivorship curves (solid line) representing age-dependent death rates for females. The black line is Spitalfields informative prior and the grey line is the Portuguese Lopes collection (used for transition analysis). Both Gompertz curves show an excellent fit nearly overlapping with the survivorship curves.
Figure 2. Kaplan–Meier survivorship (stepped lines) showing the probability of survival across adult ages and Gompertz survivorship curves (solid line) representing age-dependent death rates for females. The black line is Spitalfields informative prior and the grey line is the Portuguese Lopes collection (used for transition analysis). Both Gompertz curves show an excellent fit nearly overlapping with the survivorship curves.
Forensicsci 02 00051 g002
Figure 3. A comparison of Gompertz survivorship curves for males from eight separate skeletal collections. Lines closest to the Bass collection (target sample) and Portuguese collection (transition analysis sample) and showing similar curves were considered best fit lines.
Figure 3. A comparison of Gompertz survivorship curves for males from eight separate skeletal collections. Lines closest to the Bass collection (target sample) and Portuguese collection (transition analysis sample) and showing similar curves were considered best fit lines.
Forensicsci 02 00051 g003
Figure 4. A comparison of the Gompertz survivorship curves for males from six modern database (non-skeletal) collections. Lines closest to the Bass collection (target sample) and showing similar curves were considered best fit lines.
Figure 4. A comparison of the Gompertz survivorship curves for males from six modern database (non-skeletal) collections. Lines closest to the Bass collection (target sample) and showing similar curves were considered best fit lines.
Forensicsci 02 00051 g004
Figure 5. A comparison of Gompertz survivorship curves for females from eight separate skeletal collections. Lines closest to the Bass collection (target sample) and Portuguese collection (transition analysis sample) and showing similar curves were considered best fit lines.
Figure 5. A comparison of Gompertz survivorship curves for females from eight separate skeletal collections. Lines closest to the Bass collection (target sample) and Portuguese collection (transition analysis sample) and showing similar curves were considered best fit lines.
Forensicsci 02 00051 g005
Table 1. Sample age distributions.
Table 1. Sample age distributions.
Age IntervalAmericanPortugueseSpitalfields
Males/FemalesMales/FemalesMales/Females
16–191/05/52/2
20–297/318/225/5
30–3944/1220/1014/9
40–4972/3034/128/11
50–5983/5749/3915/22
60–6980/6735/4325/19
70–7950/6033/6613/16
80–8935/3827/432/9
90–990/00/52/0
Total372/267221/24586/93
Table 2. Transition analysis parameters for Lovejoy et al. method on Portuguese samples.
Table 2. Transition analysis parameters for Lovejoy et al. method on Portuguese samples.
ParameterEstimateAge at
Transition
EstimateAge at
Transition
MALESFEMALES
I-II2.9719.453.0521.17
II-III3.2826.473.3729.16
III-IV3.5534.893.6237.38
IV-V3.8145.183.8647.29
V-VI4.0859.284.0356.45
VI-VII4.2972.914.1663.87
VII-VIII4.5392.754.4081.80
ST DEV0.25 0.23
b4.04 4.28
Log-likelihood−303.65 −324.73
Table 3. Transition analysis parameters for transverse organization using Buckberry and Chamberlain method on Portuguese samples. Estimates and slope are on a log scale.
Table 3. Transition analysis parameters for transverse organization using Buckberry and Chamberlain method on Portuguese samples. Estimates and slope are on a log scale.
ParameterEstimateAge at TransitionEstimateAge at Transition
MALESFEMALES
I-II3.0420.963.1423.26
II-III3.5534.863.6136.95
III-IV3.8647.313.9149.82
IV-V4.4081.634.4383.92
ST DEV0.40 0.43
b2.52 2.30
Log-likelihood−258.58 −293.02
Table 4. Transition analysis parameters for texture using Buckberry and Chamberlain method on Portuguese samples. Estimates and slope are on a log scale.
Table 4. Transition analysis parameters for texture using Buckberry and Chamberlain method on Portuguese samples. Estimates and slope are on a log scale.
ParameterEstimateAge at
Transition
EstimateAge at
Transition
MALESFEMALES
I-II2.9318.693.0420.84
II-III3.2926.943.5133.42
III-IV3.7542.553.8345.98
IV-V4.0054.364.0155.08
ST DEV0.32 0.29
b3.12 3.50
Log-likelihood−205.14 −181.16
Table 5. Transition analysis parameters for apical changes using Buckberry and Chamberlain method on Portuguese samples. Estimates and slope are on a log scale.
Table 5. Transition analysis parameters for apical changes using Buckberry and Chamberlain method on Portuguese samples. Estimates and slope are on a log scale.
ParameterEstimateAge at
Transition
EstimateAge at
Transition
MALESFEMALES
I-II3.4732.143.5534.81
II-III4.3275.474.2167.04
ST DEV0.52 0.47
b1.92 2.11
Log-likelihood−192.38 −204.25
Table 6. Transition analysis parameters for microporosity using Buckberry and Chamberlain method on Portuguese samples. Estimates and slope are on a log scale.
Table 6. Transition analysis parameters for microporosity using Buckberry and Chamberlain method on Portuguese samples. Estimates and slope are on a log scale.
ParameterEstimateAge at
Transition
EstimateAge at
Transition
MALESFEMALES
I-II3.4631.733.3327.89
II-III3.7843.713.6337.82
ST DEV0.53 0.47
b1.89 2.14
Log-likelihood−167.89 −119.71
Table 7. Transition analysis parameters for macroporosity using Buckberry and Chamberlain method on Portuguese samples. Estimates and slope are on a log scale.
Table 7. Transition analysis parameters for macroporosity using Buckberry and Chamberlain method on Portuguese samples. Estimates and slope are on a log scale.
ParameterEstimateAge at TransitionEstimateAge at Transition
MALESFEMALES
I-II3.6337.763.5434.53
II-III4.5897.244.2268.01
ST DEV0.74 0.49
b1.35 2.04
Log-likelihood−214.93 −207.04
Table 8. Transition analysis parameters for Buckberry and Chamberlain Phases on Portuguese samples. Estimates and slope are on a log scale.
Table 8. Transition analysis parameters for Buckberry and Chamberlain Phases on Portuguese samples. Estimates and slope are on a log scale.
ParameterEstimateAge at
Transition
EstimateAge at
Transition
MALESFEMALES
I-II2.9619.253.0420.96
II-III3.2525.843.3829.48
III-IV3.5032.983.5534.92
IV-V3.7743.543.7341.74
V-VI4.0054.733.9451.48
VI-VII4.2671.034.2066.90
ST DEV0.30 0.30
b3.30 3.29
Log-likelihood−293.97 −278.58
Table 9. Highest posterior density (HPD) for each Lovejoy method phase with various highest posterior density regions (HPDR) in the American samples using Spitalfields samples as the informative prior.
Table 9. Highest posterior density (HPD) for each Lovejoy method phase with various highest posterior density regions (HPDR) in the American samples using Spitalfields samples as the informative prior.
PhaseHPD75% HPDR90% HPDR95% HPDR
MALES
63.44% Realized Accuracy83.87% Realized Accuracy88.98% Realized Accuracy
I18.0018.00–25.2418.00–29.1518.00–31.91
II25.1818.71–33.9518.00–39.2518.00–43.31
III34.8525.2–48.222.01–55.1220.26–59.78
IV47.0434.45–63.5030.20–71.3927.79–76.45
V61.7946.74–78.5141.19–85.4637.96–89.66
VI73.2857.96–88.0651.75–93.8148.01–97.23
VII81.3666.74–94.4460.41–99.3956.47–102.34
VIII89.1475.72–100.6469.60–104.9565.68–107.50
FEMALES
65.17% Realized Accuracy85.39% Realized Accuracy91.39% Realized Accuracy
I18.0018.00–26.0718.00–30.1518.00–32.99
II27.7720.53–37.6118.26–42.3318.00–46.27
III37.9427.96–51.4924.58–58.4922.68–63.27
IV49.7837.14–66.0332.80–73.7330.32–78.67
V61.3847.06–77.7841.83–84.7738.78–89.04
VI69.6554.76–85.0248.98–91.1845.53–94.88
VII78.5563.52–92.4157.21–97.6853.35–100.81
VIII89.1075.15–100.9168.80–105.2664.76–107.83
Table 10. Highest posterior density (HPD) for each Buckberry and Chamberlain phase with various highest posterior density regions (HPDR) in the American samples using Spitalfields samples as the informative prior.
Table 10. Highest posterior density (HPD) for each Buckberry and Chamberlain phase with various highest posterior density regions (HPDR) in the American samples using Spitalfields samples as the informative prior.
PhaseHPD75% HPDR90% HPDR95% HPDR
MALES
62.10% Realized Accuracy79.03% Realized Accuracy82.26% Realized Accuracy
I18.0018.0–27.918.00–33.5218.00–37.62
II26.0118.04–37.6918.00–45.9918.00–51.78
III35.6324.04–52.8520.56–61.5818.75–67.09
IV48.8933.73–68.7428.79–77.5326.05–82.83
V62.6745.38–80.8438.99–87.9635.30–92.15
VI73.3655.94–89.1248.76–94.9644.46–98.39
VII85.0068.86–98.1961.43–102.9356.73–105.69
FEMALES
50.56% Realized Accuracy64.79% Realized Accuracy71.16% Realized Accuracy
I18.0018.00–29.8918.00–36.418.00–41.12
II30.1720.19–45.3618.00–53.5718.00–60.24
III40.4827.42–59.5723.43–69.0021.29–74.91
IV50.1134.67–70.3829.65–79.2926.85–84.62
V61.8844.29–80.8137.92–88.2234.27–92.54
VI73.3855.31–89.6147.93–95.5543.54–99.01
VII86.4369.86–99.5962.13–104.2357.21–106.92
Table 11. Cumulative binomial tests of the performance of the Lovejoy (LJ) and Buckberry and Chamberlain (BC) methods on American holdout samples using Portuguese transition analysis parameters and Spitalfields samples as the informative prior.
Table 11. Cumulative binomial tests of the performance of the Lovejoy (LJ) and Buckberry and Chamberlain (BC) methods on American holdout samples using Portuguese transition analysis parameters and Spitalfields samples as the informative prior.
MethodCINumber of SuccessesNumber of Failuresp-ValueProbability of SuccessBias
MALES
LJ75%236136<0.00010.637.95
LJ90%312600.00030.848.33
LJ95%33141<0.00010.898.62
BC75%231141<0.00010.6210.85
BC90%29478<0.00010.7910.77
BC95%30666<0.00010.8211.49
FEMALES
LJ75%174930.00040.657.99
LJ90%228390.01830.858.04
LJ95%244230.01090.918.09
BC75%135132<0.00010.51−10.91
BC90%17394<0.00010.65−11.02
BC95%19077<0.00010.71−11.08
Bold indicates significantly lower performance than probability.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hens, S.M.; Godde, K. A Bayesian Approach to Estimating Age from the Auricular Surface of the Ilium in Modern American Skeletal Samples. Forensic Sci. 2022, 2, 682-695. https://0-doi-org.brum.beds.ac.uk/10.3390/forensicsci2040051

AMA Style

Hens SM, Godde K. A Bayesian Approach to Estimating Age from the Auricular Surface of the Ilium in Modern American Skeletal Samples. Forensic Sciences. 2022; 2(4):682-695. https://0-doi-org.brum.beds.ac.uk/10.3390/forensicsci2040051

Chicago/Turabian Style

Hens, Samantha M., and Kanya Godde. 2022. "A Bayesian Approach to Estimating Age from the Auricular Surface of the Ilium in Modern American Skeletal Samples" Forensic Sciences 2, no. 4: 682-695. https://0-doi-org.brum.beds.ac.uk/10.3390/forensicsci2040051

Article Metrics

Back to TopTop