This work developed, validated, and compared statistical shape, statistical intensity, and statistical shape and intensity models (SSMs, SIMs, SSIMs) of scapulae from a clinical population. SSMs efficiently describe bone shape variation while SIMs describe bone material property variation, and SSIM's combine description of both variables. This work establishes these models' efficacy and whether they can be used in surgical planning. Models were developed using shoulder arthroplasty data of patients with bone erosion, which is challenging to treat and would benefit from improved surgical planning. Models were created using previously validated nonrigid registration and material property assignment processes that were optimized for scapula characteristics. The models were assessed using standard metrics, anatomical measurements, and correlation analyses. The SSM and SIM specificity and generalization error metrics were 3.4 mm and <1 mm and 184 HU and 156 HU, respectively. The SSIM did not achieve the same level of performance as the SSM and SIM in this study (e.g., shape generalization: SSIM—2.2 mm versus SSM—<1 mm). Anatomical correlation analysis showed that the SSM more effectively and efficiently described shape variation compared to the SSIM. The SSM and SIM modes of variation were not strongly correlated (e.g., rmax = 0.56 for modes explaining ≤2.1% of variance). The SSIM is outperformed by the SSM and SIM and the latter two are not strongly correlated; therefore, using the SSM and SIM in conjunction will generate synthetic bone models with realistic characteristics and thus can be used for biomechanical surgical planning applications.