Intraskeletal Variation in Cortical Bone Quantity in a Medieval Italian Sample: A Multivariate Exploratory Approach

Bioarcheologists interpret skeletal health by examining variability within and between individuals. Studies of bone loss have generated contradictory and conflicting results regarding the onset and severity of age-related bone loss on a global and temporal scale, perhaps due to mismatched methodologies. Intraskeletal comparisons of bone tissue prove challenging precisely because of heterogeneous baselines in quantity and remodeling of cortical bone throughout the skeleton, as well as evolutionary histories and environmental impacts on growth and development. Here we analyze cortical bone indicators from the rib, metacarpal, and femoral cortical bone in a subset of individuals (n = 72) regions from the medieval Italian archaeological site of Pieve di Pava. To facilitate intraskeletal comparisons across elements with different biological baselines, we standardize cortical bone parameters using z-scores. Variation in relative intraskeletal cortical bone was assessed using accessible multivariate methods (principal component analysis and hierarchical cluster analysis). Results suggest an association between femoral and metacarpal cortical bone values, with stochastic trends in metacarpal and femoral relative bone quantity in relation to the rib bone quantity at the sample level. Our study demonstrates that while intraskeletal analyses are challenging, they are made more robust by synthesizing multivariate methods alongside exploratory data analysis (EDA) methods to tack between sample-level and individual-level scales and variability. Ultimately, we advocate for leveraging multivariate techniques not as a final step, but rather as a means of generating new hypotheses and challenging tendencies to a priori establish typological groups in the research process.