| Title: | Estimation of Developmental Age |
|---|---|
| Description: | Bayesian methods for estimating developmental age from ordinal dental data. For an explanation of the model used, see Konigsberg (2015) <doi:10.3109/03014460.2015.1045430>. For details on the conditional correlation correction, see Sgheiza (2022) <doi:10.1016/j.forsciint.2021.111135>. Dental scoring is based on Moorrees, Fanning, and Hunt (1963) <doi:10.1177/00220345630420062701>. |
| Authors: | Valerie Sgheiza [aut, cre] |
| Maintainer: | Valerie Sgheiza <[email protected]> |
| License: | GPL |
| Version: | 0.1.0 |
| Built: | 2024-12-01 08:03:31 UTC |
| Source: | CRAN |
A function for finding and correcting fuzzy posteriors produced by mvcp_est.f
find_fuzzies.f(cases, prior)find_fuzzies.f(cases, prior)
cases |
as input data.table of rows of collapsed dental development scores, or vector of scores from single individual |
prior |
as string argument for desired prior |
for one case, list object with interpretation and data.table of true age, lower and upper bounds of HDR, mode, and posterior. If multiple cases are entered, only the data.table is returned.
find_fuzzies.f(c(NA,NA,9,10,11,14,15,10,15,11),prior="jeff")find_fuzzies.f(c(NA,NA,9,10,11,14,15,10,15,11),prior="jeff")
WARNING: this function is based on scoring with an additional crypt stage. Pre-collapse staging: 0-no crypt, 1-crypt, 2-Ci, 3-Cco, 4-Coc, 5-Cr1/2, 6-Cr3/4, 7-Crc, 8-Ri, 9-Cli, 10-R1/4, 11-R1/2, 12-R3/4, 13-Rc, 14-A1/2, 15-Ac
mfh_collapse(cases)mfh_collapse(cases)
cases |
as input data.table of Moorrees et al. dental development scores, or vector of scores from single individual |
returns data.table of collapsed scores
mfh_collapse(c(NA,NA,15,15,14,12,11,15,11,8))mfh_collapse(c(NA,NA,15,15,14,12,11,15,11,8))
Age estimation from a single case, called by mvcp_est.f
mv.probit(case, prior, drop)mv.probit(case, prior, drop)
case |
as input vector from object dat |
prior |
as string argument for desired prior |
drop |
as column to be dropped based on outcome of find_fuzzies.f |
returns probs as a vector of posterior probabilities
A function for estimating age from 8 teeth via multivariate cumulative probit and a user-determined prior. Call directly iff you do NOT want to correct for fuzzy posteriors. Call find_fuzzies.f otherwise.
mvcp_est.f(prior, dat)mvcp_est.f(prior, dat)
prior |
as string argument for desired prior |
dat |
as input data.table of rows of dental development scores, may have only 1 row |
returns postm as a data.table of posterior probability distributions
Collapsed London Atlas scores of dental development from NMDID images by Stull. Do not apply mfh_collapse before using. Stages are already collapsed!
nmdid.testnmdid.test
## 'nmdid.test' A data.table with 188 rows and 10 columns:
Decedent record number from NMDID
age in decimal years
Collapsed London atlas score of left permanent mandibular teeth I1-M3