A better understanding of how the shape and density of the shoulder vary among members of a population can help design more effective population-based orthopedic implants. The main objective of this study was to develop statistical shape models (SSMs) and statistical density models (SDMs) of the shoulder to describe the main modes of variability in the shape and density distributions of shoulder bones within a population in terms of principal components (PCs). These PC scores were analyzed, and significant correlations were observed between the shape and density distributions of the shoulder and demographics of the population, such as sex and age. Our results demonstrated that when the overall body sizes of male and female donors were matched, males still had, on average, larger scapulae and thicker humeral cortical bones. Moreover, we concluded that age has a weak but significant inverse effect on the density within the entire shoulder. Weak and moderate, but significant, correlations were also found between many modes of shape and density variations in the shoulder. Our results suggested that donors with bigger humeri have bigger scapulae and higher bone density of humeri corresponds with higher bone density in the scapulae. Finally, asymmetry, to some extent, was noted in the shape and density distributions of the contralateral bones of the shoulder. These results can be used to help guide the designs of population-based prosthesis components and pre-operative surgical planning.