Objectives: The statistical analysis of fossil remains is essential to understand the evolution of the genus Homo. Unfortunately, the human fossil record is straight away scarce and plagued with severe loss of information caused by taphonomic processes. The recently developed field of Virtual Anthropology helps to ameliorate this situation by using digital techniques to restore damaged and incomplete fossils.Materials and methods: We present the package Arothron, an R software suite meant to process and analyze digital models of skeletal elements. Arothron includes tools to digitally extract virtual cavities such as cranial endocasts, to statistically align disarticulated or broken bony elements, and to visualize local variations between surface meshes and landmark configurations.Results: We describe the main functionalities of Arothron and illustrate their usage through reproducible case studies. We describe a tool for segmentation of skeletal cavities by showing its application on a malleus bone, a Neanderthal tooth, and a modern human cranium, reproducing their shape and calculating their volume. We illustrate how to digitally align a disarticulated model of a modern human cranium, and how to combine piecemeal shape information on individual specimens into one. In addition, we present useful visualization tools by comparing the morphological differences between the right hemisphere of the Neanderthal and the modern human brain.Conclusions: The Arothron R package is designed to study digital models of fossil specimens. By using Arothron, scientists can handle digital models with ease, investigate the inner morphology of 3D skeletal models, gain a full representation of the original shapes of damaged specimens, and compare shapes across specimens.

Objectives: The advent of micro‐computed tomography (μCT) made cancellous bone more accessible than ever before. Nevertheless, the characterization of cancellous bone is made difficult by its inherent complexity and the difficulties in defining homology across datasets. Here we propose novel virtual methodological approaches to overcome those issues and complement existing methods.Materials and methods: We present a protocol for the isolation of the whole cancellous region within a μCT scanned bone. This method overcomes the subsampling issues and allows studying cancellous bone as a single unit. We test the protocol on a set of primate bones. In addition, we describe a set of morphological indices calculated on the topological skeleton of the cancellous bone: node density, node connectivity, trabecular angle, trabecular tortuosity, and fractal dimension. The usage of the indices is shown on a small comparative sample of primate femoral heads.Results: The isolation protocol proves reliable in isolating cancellous structures from several different bones, regardless of their shape. The indices seem to detect some functional differences, although further testing on comparative samples is needed to clarify their potential for the study of cancellous architecture.Conclusions: The approaches presented overcome some of the difficulties of trabecular bone studies. The methods presented here represent an alternative or supporting method to the existing tools available to address the biomechanics of cancellous bone.