In geanes/kidstats: Subadult age and sex estimations based on long bone measurements

library(pander)
library(caret)
library(dplyr)

Generated by KidStats on r Sys.Date()

---

Input

• Reference sample: r switch(input$refsamp, 'za' = 'South Africa')
• Transformation on Age: r switch(input$transform, 'none' = 'None', 'sqrt' = 'Square root', 'cbrt' = 'Cube root')

Input measurements

| | Diaphyseal Length | Proximal Breadth | Midshaft Breadth | Distal Breadth | | :-------- | :---------------: | :--------------: | :--------------: | :------------: | | Femur | r input$fdl | | r input$fmsb | r input$fdb | | Tibia | r input$tdl | r input$tpb | r input$tmsb | r input$tdb | | Humerus | r input$hdl | r input$hpb | r input$hmsb | r input$hdb | | Radius | r input$rdl | r input$rpb | r input$rmsb | r input$rdb | | Ulna | r input$udl | | r input$umsb | |

---

Age

Estimation

| Lower Prediction Interval | Estimated Age | Upper Prediction Interval | | :-----------------------: | :--------------------------: | :------------------------: | | r earth_mod()[[2]][2] | r earth_mod()[[2]][1] | r earth_mod()[[2]][3] |

Model

• Sample size used to model: r nrow(age_samp())
• Age model Generalized R^2: r round(earth_mod()[[1]]$grsq, digits = 4)
• Excluded measurements: r if (length(input$ex_age) == 0) pander("none") else pander(input$ex_age)

hinges <- add_rownames(data.frame(earth_mod()[[1]]$cuts), "Coefficients")
hinges <- hinges[earth_mod()[[1]]$selected.terms,]
pander(hinges)

---

SEX

Sex Estimation

| Posterior Female | Posterior Male | | :------------------: | :------------------: | | r fda_mod()[[2]]$F | r fda_mod()[[2]]$M |

Model

• Sample size used to model: r nrow(sex_samp())
• Classification Accuracy: r if (input$bstrap_ca) paste(round(fda_mod()[[3]]$t0, digits = 4)," (bootstrapped)") else round(fda_mod()[[3]], digits = 4)
• Excluded measurements: r if (length(input$ex_sex) == 0) pander("none") else pander(input$ex_sex)

fda_coef <- data.frame(coef(fda_mod()[[1]]))
pander(fda_coef)

geanes/kidstats documentation built on May 17, 2019, 12:15 a.m.