American J Phys Anthropol - 2012 - Gilmore - Estimating Age From Adult Occlusal Wear A Modification of The Miles Method
American J Phys Anthropol - 2012 - Gilmore - Estimating Age From Adult Occlusal Wear A Modification of The Miles Method
ABSTRACT The Miles method of age estimation with age categories estimated from traditional methods;
relies on molar wear to estimate age and is widely used in this suggests that the modifications have not hampered
bioarcheological contexts. However, because the method the ability of the Miles method to estimate age even in
requires physical seriation and a sample of subadults to heterogeneous samples. As expected, the error increases
estimate wear rates it cannot be applied to many samples. with age and in populations with smaller sample sizes.
Here, we modify the Miles method by scoring occlusal These modifications allow the Miles method to be applied
wear and estimating molar wear rates from adult wear to skeletal samples of adult crania that were previously
gradients in 311 hunter-gatherers and provide formulae only amenable to cranial suture age estimation, and
to estimate the error associated with each age estimate. A importantly, provide a measure of uncertainty for each
check of the modified method in a subsample (n 5 22) age estimate. Am J Phys Anthropol 149:181–192, 2012.
shows that interval estimates overlap in all but one case C 2012
V Wiley Periodicals, Inc.
The ability to estimate adult age-at-death from skele- 1971; Griffin et al., 1979; Ruff, 1981; Kieser et al., 1983;
tal remains is critical to bioarcheological studies (e.g., Levers and Darling, 1983; Lovejoy, 1985; Whittaker
paleodemography and paleopathology) and in forensic et al., 1985, 1987; Harper, 1994; Visser, 1994; Borrman
contexts. As such, a staggering volume of literature is et al., 1996; Sullivan, 2004, 2005; Papageorgopoulou
devoted to methods of age estimation; commonly used et al., 2009; Sagona et al., 2010; Cole and Waldron, 2011)
methods include scoring the degeneration of the pubic and tests of the method routinely suggest that it per-
symphysis (e.g., Acsádi and Nemeskéri, 1970; Meindl et forms as well or better than other methods of skeletal
al., 1985; Brooks and Suchey, 1990) and auricular sur- age estimation (Nowell, 1978; Kieser et al., 1983; Love-
face (e.g., Lovejoy et al., 1985b; Buckberry and Chamber- joy, 1985; Lovejoy et al., 1985a; but see Santini et al.,
lain, 2002; Igarashi et al., 2005), cranial suture fusion 1990). Practically, the Miles method requires that the
(e.g., Meindl and Lovejoy, 1985), and dental wear (e.g., entire collection be seriated by tooth wear (Helm and
Miles, 1962; Brothwell, 1981; Walker et al., 1991). Given Prydso, 1979) and that it include roughly 20 juveniles to
the superior preservation of dental tissues and the his- make reliable wear rate estimates (Nowell, 1978).
toric preference of early anthropologists to collect skulls Unfortunately, many skeletal collections are curated
over postcranial elements (Walker et al.1991), the ability by separate institutions and therefore cannot be seriated
to estimate age from the dentition alone can be particu- together according to wear. Furthermore, juvenile
larly useful for researchers studying past populations. remains are more fragile than those of adults and are of-
Age estimation from dental wear is further distinguished ten less well preserved in archaeological assemblages.
from other methods in that it is based on an indicator This study was necessitated by just such a situation:
(enamel) that does not have the potential to remodel while collecting data to compare antemortem tooth loss
(Boldsen, 2005), and that it can be calibrated to the pop- frequencies in hunter-gatherers, one of us (CCG) was
ulation in question, rather than requiring the applica- obliged to use collections composed primarily (sometimes
tion of aging standards from a reference population exclusively) of adult crania and curated by diverse repo-
(Lovejoy et al., 1985a). This combination of qualities
likely explains why dental wear was found to be the best
single age indicator when compared with pubic symphy-
sis, auricular surface, femoral radiographs, and cranial Grant sponsors: The L.S.B. Leakey Foundation, the University of
suture closure in a known-age sample (Lovejoy et al., California, Davis, Department of Anthropology Summer Fellow-
1985a). ships, and the UC Davis Institute of Governmental Affairs.
Though researchers have long noted the value of den-
tal wear to classify individuals into broad age categories, *Correspondence to: Cassandra Gilmore, Department of Anthro-
Miles (1962, 1963, 1969, 2001) was the first to use molar pology, University of California, Davis, One Shields Avenue, Davis,
wear to estimate age systematically by progressively CA 95616, USA.
extrapolating wear rates from younger individuals—for E-mail: [email protected]
whom age can be estimated with relative certainty—to
older individuals where the determination of age is more Received 1 December 2011; accepted 22 May 2012
problematic. The Miles method, or the aging chart
accompanying the original publication, is commonly used DOI 10.1002/ajpa.22106
Published online 5 July 2012 in Wiley Online Library
to estimate age in bioarcheological contexts (Wolpoff, (wileyonlinelibrary.com).
C 2012
V WILEY PERIODICALS, INC.
10968644, 2012, 2, Downloaded from https://onlinelibrary.wiley.com/doi/10.1002/ajpa.22106 by Egyptian National Sti. Network (Enstinet), Wiley Online Library on [22/09/2025]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
182 C.C. GILMORE AND M.N. GROTE
TABLE 1. Sample composition molar had been in functional occlusion. For example, if
Population N the second molars are just erupting then the individual is
roughly 12 years old, and the first molars have been in
Afalou 23 occlusion for about 6 years. The wear on the first molars is
Alaska 41 thus equivalent to 6 years of wear. Miles seriated the sam-
Australian/Tasmanian 70 ple of adults and projected the wear rates observed in the
Buriat 27
Chukchi 16 subadult sample to estimate the functional ages of young
Sami 13 adult molars. The functional age based on wear stage plus
Fuegian 9 the age of eruption results in an estimate of the age of an
Khoi-San 38 individual. Using the estimated functional ages of young
Native Californian 74 adult molars, Miles was then able to identify stages of
Total 311 wear at older functional ages. This process repeats, build-
ing upon successively older individuals to estimate func-
tional ages of progressively more worn molars. The accu-
racy of the method is expected to diminish with each
sitories across North America and Europe. The ages of extrapolation. During the course of his investigation,
individuals in these collections are unknown, but a com- Miles observed a consistent relationship between wear
parison of antemortem tooth loss between populations rates of the first, second, and third molars; he estimated
without considering sample age structures would ignore that it would take the second molar 6.5 years and the
a potentially large source of variation. third molar 7 years to reach a wear stage that took the
The purpose of the present study is to suggest a prac- first molar only 6 years to achieve.
tical modification of the Miles method using a widely Miles (1962) states two primary assumptions of this
accepted method of estimating wear rates from wear gra- method: 1) individuals belonging to a single biological
dients in adults (e.g., Smith, 1972; Lunt, 1978; Benfer population have relatively consistent wear rates; 2)
and Edwards, 1991; Constandse-Westermann, 1997; Lev- ‘‘. . .the rates of wear of the molars remain[ed] constant
Tov and Smith, 2006; Bulbeck and O’Connor, 2011). We throughout the life of the dentition’’ (884). In practice,
will call it the modified Miles method. The modification Miles did not require assumption 2, because he updated
relies on data routinely collected during skeletal analysis molar-specific wear rates in a step-wise manner as he
and can be applied to samples where juveniles and post- extrapolated to successively older individuals. However,
cranial elements are under-represented. Unlike many Miles does require a slightly weaker assumption: that
other methods of dental age estimation, it does not the ratio of wear rates among molar regions M1, M2,
require destruction or disarticulation of the remains and and M3 remains constant through the lifespan. We will
is applicable to archaeological collections. Here, we call this assumption 3. Please see Appendix A for a
describe the application of the modified Miles method to mathematical treatment of this assumption and its rela-
a sample of hunter-gatherers and further develop a tionship to the assumptions of the modified method pro-
method of quantifying the error associated with each age posed here.
estimate. The modified Miles method makes two changes to the
Miles method. First, we score occlusal wear following
MATERIALS Scott (1979) in lieu of physically seriating the sample.
Second, we use adult wear gradients and eruption sched-
The sample is composed of individuals from hunter- ules from nonindustrialized populations to estimate wear
gatherer and pastoralist populations from a variety of rate. We require all three assumptions, 1, 2, and 3, dis-
latitudes with diverse population histories and diets cussed in the context of Miles’ method above, along with
(Table 1). The Khoi-San, Australian aborigine, Fuegian, one additional assumption: 4) that molars in the same
and North Eurasian (Buriat, Chukchi, Sami) collections individual wear at the same rate. This is essentially a
all derive from modern exploratory expeditions, while more restrictive version of assumption 3; in addition to
the Afalou, Alaskan, and Native Californian samples are requiring molars to wear at a constant ratio, the modi-
archaeological. Afalou is a Late Pleistocene cave site in fied Miles method requires that this ratio be 1:1:1. We
Algeria, dating to roughly 13,000 years ago (Hachi, further describe each change and discuss each assump-
1996), while the Alaskan and Californian remains are tion below.
significantly younger. The Alaskan sample comes from
Paleo-Eskimo remains from the Point Hope site in
Alaska, dating from about 2,360 to 1,619 BP (Rainey and Modification 1: Scoring occlusal wear
Ralph, 1959), and the Californian remains are from a
We chose to seriate individuals by quantifying molar
collection of sites in Central California ranging from
wear with Scott’s (1979) method of scoring occlusal wear.
about 5,000 to 500 BP (Milliken et al., 2007). Each skele-
Ordinal scoring methods like Scott’s represent the least
tal population conforms to a general geographic area and
dramatic departure from Miles’ assessments of occlusal
broadly similar diet and lifestyle, but none represent the
wear and are routinely used by researchers to estimate
ideal small biological populations existing at a single
wear rates from molar gradients (see section on Modifi-
time and place for which Miles created his method of
cation 2, below). They are quickly and easily applied and
age estimation.
have the advantage of being exceptionally nondestruc-
tive in that they do not require sectioning, disarticulat-
METHODS ing, or even coming into physical contact with a speci-
The Miles method and its assumptions men. Scott’s method is recommended by Standards of
Data Collection from Human Skeletal Remains (Buikstra
Miles (1962, 1963, 1969, 1978) relied on molar eruption and Ubelaker, 1994) and is commonly recorded by
schedules to calculate the number of years that a child’s researchers documenting human remains, suggesting
Fig. 1. Average molar wear score for each molar class plotted against the others for every individual in every population in the
sample. A) The relationship between the M1 and M2; B) the relationship between the M1 and M3; C) the relationship between
the M2 and M3. The M1 and the M2 are more closely correlated with each other than either is with the M3, demonstrating the
increased variance in the wear scores of the latter molar class.
late wear rates for a skeletal population. When calcu- TABLE 2. Eruption ages and standard deviations from studies
lated separately, our data suggested a slightly higher of nonindustrialized populations that reported sample sizea
wear rate in the maxillary than the mandibular molars; Tooth Eruption age (yr) Standard deviation (yr) Covariance
results from some previous work agree (Murphy, 1959;
Lunt, 1978; Kieser et al., 1983; Mo and Peng, 1983; M1 6.02 0.85 –
McKee and Molnar, 1988), though others have found the M2 11.31 1.13 –
M3 17.87 1.16 –
opposite pattern (Molnar et al., 1983; Lovejoy, 1985;
M1, M2 – – 0.717
Mays, 2002), or no difference at all (Dreier, 1994; Li and M1, M3 – – 1.018
Ji, 1995; Mays et al., 1995). Reports on the effect of sex M2, M3 – – 1.362
on wear rate are similar; sexual dimorphism seems to be
a
significant in some populations (Tomenchuck and May- Data are averaged from Steggerda and Hill (1942), Fried-
hall, 1979; e.g., Benfer and Edwards, 1991), and unim- lander and Bailit (1969), Brown (1978), and Mayhall et al.
portant in others (Lunt, 1978; Richards and Brown, (1978).
1981; Lovejoy, 1985; Walker et al., 1991; Mays, 2002).
Ultimately, we elected to follow Miles’ procedure and
combine maxillary and mandibular arcades as well as approach zero and will not be representative of the true
the sexes. By combining the maxilla and mandible of wear rate. Therefore, if any quadrant of individual i has
each individual to estimate a wear rate, we increase the a wear score of 10 (the maximum), the average wear
molar sample size used for each estimate of individual score for that molar class in individual i is excluded
wear rate and decrease the influence of ‘‘noise’’ from from calculating the population average wear rate.
intraindividual variation. For example, antemortem After creating average molar class wear scores for
tooth loss and dental disease have been predicted to ei- each individual, we calculated individual wear score dif-
ther reduce the wear rate of the occluding molar (Mays, ferences: Di 5 wearM1i 2 wearM2i. Across the sample,
2002), to increase wear rate in non-occlusal partners most such differences are positive. The difference in
(Miles, 1962), or both, though Mays (2002) found no evi- wear of the M1-M3 and M2-M3 could also be used to cal-
dence of a change in wear rate of the occlusal partner in culate wear rate. Though Mays (2002) finds that the
a low wear rate sample. One could also form distinct third molar is equally well correlated with age as the
prediction equations for males and females, but the loss first and second molars, Figure 1 demonstrates that M3
of precision (resulting from splitting the sample by sex) wear is less closely correlated with M1 or M2 wear than
is compounded with the uncertainty of assigning sex cor- M1 and M2 are with each other in these data, probably
rectly. Our procedures arguably introduce bias, because because third molars tend to be more variable in erup-
distinct categories of teeth and sex are combined, but tion time (see Table 2) and morphology than first or sec-
also decrease variance by taking full advantage of sam- ond molars. For this reason, we use only the M1-M2 dif-
ple degrees of freedom. ference to calculate wear rate.
To calculate wear rates we average (equally weighted) Figure 2 plots M1-M2 differences, D, against M1 wear
the wear scores of the maxillary and mandibular molars scores, wearM1i, and shows that variation in D appears
of individual i to create a first molar wear score in Scott to increase with wear. This pattern likely results from
units (wearM1i), a second molar wear score (wearM2i), the accumulation of interindividual variation over a life-
and a third molar wear score (wearM3i). We refer to time of wear and supports previous observations that
these as average molar class wear scores. At some point, error in age estimation based on wear increases with age
each tooth will reach a maximum wear stage and cease (Miles, 1962; Kieser et al., 1983; Mays, 2002). It is possi-
to change with time. The M1 is likely to be the first ble that the lack of variation at very low wearM1i is in
molar class to reach this maximum stage, while the M2 part an artifact of the method—maximum differences in
and M3 continue to wear (Benfer and Edwards, 1991). wear are constrained by the relation Di wearM1i 2 1,
As the M1 ceases to wear and the other molar classes and this constraint is especially severe when wearM1i is
continue to wear, the difference in wear between the small. The finite property of the scoring system suggests
first and second molars or first and third molars will that very lightly worn teeth could skew estimates of the
Brooks and Suchey, 1990). The Miles method and this We therefore suppose that molar wear scores are real-
modification are also able to estimate age without refer- valued, differentiable functions of t 5 age, denoted as
ence to standards from a separate population and there- M1(t), M2(t), M3(t). Our development is concerned with
fore avoid mimicking the age distribution of the refer- M1(t) and M2(t)—the wear scores on which the modified
ence population (Bocquet-Appel and Masset, 1982; Aiello Miles method focuses—but the other pairs (M1(t) and
and Molleson, 1993) or making incorrect assumptions M3(t), M2(t) and M3(t)) could be treated analogously.
about the variation in the target population. Though the Taking ‘‘rate’’ as a synonym for derivative with respect
modification does not appear to improve upon the tend- to t, assumption #3 translates as
ency of commonly applied methods of skeletal age esti-
mation to underestimate age in older adults (Molleson dM1ðtÞ=dt
et al., 1993; Aykroyd et al., 1999), a large, known-age ¼ constant: ð1Þ
dM2ðtÞ=dt
sample is needed to confidently characterize the biases.
Methods of age estimation from molar wear are expected
to perform more poorly in recent, low wear rate popula- A necessary and sufficient condition for Eq. (1) is that
tions than in prehistoric populations, as interindividual there exist real numbers l and b (b = 0) such that
variation in wear may overwhelm signals of age when
wear rate is low (Molnar et al., 1983; Walker et al., M2ðtÞ ¼ l þ bM1ðtÞ: ð2Þ
1991; Mays, 2002). Nevertheless, we encourage further
testing to explore the utility of the modified Miles
method, particularly in cases where age is known and The precise form of M1(t) in Eq. (2) does not need to be
wear rates are high (rarely available, as ages are typi- specified in order for Eq. (1) to hold. Equation (2) implies
cally only known for relatively recent samples) or per- that the difference between first and second molar wear
haps in a high wear population with ages estimated by scores at age t is
tooth cementum annuli methods.
Although it still suffers from some of the same limita- M1ðtÞ M2ðtÞ ¼ ð1 bÞM1ðtÞ l: ð3Þ
tions as other methods of age estimation, the modifica-
tion of the Miles method we present here is a viable We understand Eqs. (1)–(3) to be consequences of Miles’
option for estimating age in skeletal samples that might assumption #3, for numerical wear scores that vary con-
otherwise be aged only by cranial suture closure. This tinuously with t.
modification is an improvement over cranial suture age In the modified Miles method, a population average
estimation in that aging standards are determined from wear score difference D, and subsequently a wear rate R,
the population in question and that point-age estimates are obtained by assuming that
with error, effectively ranking individuals, are produced.
M1ðtÞ M2ðtÞ ¼ constant; ð4Þ
ACKNOWLEDGMENTS which is equivalent to requiring b 5 1 in Eq. (3). Thus,
The first author offers gratitude to curators and col- where Miles (1962) finds that 6 years of wear on the M1
leagues at the following museums for their permission to is roughly equivalent to 6.5 years of wear on the M2, the
study their collections: American Museum of Natural modified Miles method assumes the ratio is 1:1 (assump-
History, Institut de Paléontolgie Humaine, Musée de tion #4). The wear rate R ¼ D=ðTM2 TM1Þ, where TM1
l’Homme, National Museum of Natural History (Smith- and TM2 are respectively the mean ages of M1 and M2
sonian), Peabody Museum at Harvard University, eruption, is understood to approximate the average wear
Phoebe A. Hearst Museum of Anthropology at the Uni- rate over the portion of the lifespan for which the modi-
versity of California, Berkeley, and the Anthropology fied method is intended (broadly speaking, early-adult
Department at the Universität Wien. The authors thank through mid-adult, with decreasing accuracy for older
Dr. Eric Bartelink for his willingness to share his age adults).
estimates of the Native Californians at the Phoebe An estimated age (based on M1 wear) is obtained
Hearst Museum and Drs. Eric Bartelink, Henry under the modified method using the linear equation
McHenry, Teresa Steele, and Tim Weaver for their com- wearM1
ments on this manuscript. Finally, the authors thank Age ¼ TM1 þ ; ð5Þ
R
the Editor-in-Chief, Associate Editor, and two anony-
mous reviewers for valuable suggestions on an earlier
version that led to substantial revisions. where wearM1 is the M1 wear score of an individual of
unknown age. Thus assumption #2, that molar wear
rates remain constant throughout the life of the denti-
APPENDIX A: CONNECTING THE ASSUMPTIONS tion, is implied by the use of Eq. (5). Although Miles
OF THE MILES AND MODIFIED MILES METHODS (1962) states assumption #2, he does not seem to use it
in practice; in contrast, the modified method requires
The relationship between the Miles and modified Miles that assumption #2, and consequently Eq. (5), hold to a
methods rests on assumption #3: that the ratio of wear reasonable approximation across the middle lifespan.
rates among molar regions M1, M2, and M3 remains Figure A1 presents a scheme that is compatible with
constant through the lifespan. Although Miles (1962) the Miles method and shows the linear approximations
scored molar wear qualitatively and cautioned against implicit in the modified method. This scheme has wear
overt mathematization of his method, a mathematical scores increasing with age, though other schemes
translation seems necessary here for concreteness, and (including those for which ‘‘scores’’ decrease monotoni-
for comparison with the numerical approach of the modi- cally with age) are compatible with the two methods and
fied method. could be depicted analogously.
Step 5
Calculate wear rate per year for population p, where
Dp is defined by the equation in Step 4, and TM1, TM2,
and TM3 are the mean ages of M1, M2, and M3 erup-
tion, respectively, in years (Table 2):
Dp
Rp ¼
TM1 TM2
Step 6
Estimate age of each individual i in population p:
Fig. A1. A molar wear scoring scheme compatible with the wearM1i
Miles and modified Miles methods. M1(t) and M2(t) are hypotheti- AgeM1i ¼ TM1 þ
cal wear score functions compatible with Miles’ assumption #3. The Rp
tangent lines indicated by dashes have slopes conforming to the wearM2i
6:6.5 relationship noted by Miles, and tangent lines at all common AgeM2i ¼ TM2 þ
time-points of M1(t) and M2(t) have this property. The modified
Rp
Miles method assumes that the ratio of slopes is 1:1 at all common wearM3i
points, and further approximates wear score functions by straight AgeM3i ¼ TM3 þ
Rp
lines (linear M1 and M2 in the figure). Given a particular molar
wear score, age is estimated by reflecting the score across a corre-
sponding wear function to a value on the ‘‘age’’ axis. The modified APPENDIX C: FORMULAE FOR VARIANCE
method is expected to perform best where M1(t) and linear M1 are ESTIMATION
close in value, shown here in the mid-adult years. The modified
method may notably underestimate the ages of older adults. TM1 The function for estimation based on M1 wear can be
and TM2 are eruption times of M1 and M2, respectively. expressed as
APPENDIX B: FORMULAE FOR AGE ESTIMATION TM2 TM1
AgeM1i ¼ TM1 þ wearM1i ¼ f ðTM1; TM2; Dp Þ:
Each step in estimating age with this method is detailed Dp
below, including the steps involved in estimating popula-
tion wear rates discussed in the methods section above. Estimates based on M2 or M3 have similar forms. The ap-
proximate variance of AgeM1i by the Delta-Method is then:
Step 1
2
For each individual i, subtract the average M2 class @f
varðAgeM1i Þ varðTM1Þ
wear score from the average M1 class wear score: @TM1
2 !2
Di ¼ wearM1i wearM2i @f @f
þ varðTM2Þ þ varðDp Þ
@TM2 @Dp
Step 2
@f @f
Average the differences for all individuals in popula- þ2 covðTM1; TM2Þ;
@TM1 @TM2
tion p and calculate the standard deviation, where Np is
uncorrected sample size and Di is as in Step 1:
in which sample statistics, variances, and covariances from
PNp sPffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi Tables 2 and 3 are substituted at appropriate places. In deriv-
Np 2
i¼1 Di i¼1 ðDi Dp Þ ing the approximate variance, we have assumed that wearM1i
Dp ¼ sdðDp Þ ¼
Np Np 1 is measured without error, and that covðTM1; Dp Þ ¼ cov
ðTM2; Dp Þ ¼ 0. Delta-Method approximations are under-
Step 3 stood to be more accurate in large samples. Similar calcu-
lations produce variances for age estimates based on M2
Remove lightly worn individuals from the dataset by or M3, though in the latter case, the function f depends on
keeping only specimens that fit the requirements of the four variables: TM1, TM2, TM3, and Dp . After substitu-
following equation: tion, the Delta-Method variances are: