In computer-assisted orthopaedic surgery, recovering three-dimensional patient-specific anatomy from incomplete information has been focus of interest due to several factors such as less invasive surgical procedures, reduced radiation doses, and rapid intra-operative updates of the anatomy. The aim of this paper is to report results obtained combining statistical shape modeling and multivariate regression techniques for predicting bone shape from clinically and surgically relevant predictors, including sparse observations of the bone surface but also morphometric and anthropometric information. Different state of the art methods such as partial least square regression, principal component regression, canonical correlation analysis, and non-parametric kernel-based regression are compared. Clinically relevant surrogate variables and combinations are investigated on a database of 142 femur and 154 tibia shapes obtained from CT images. The results are evaluated using cross validation to quantify the prediction error. The proposed approach enables to characterize the added value of different predictors in a quantitative and localized fashion. Results indicate that complementary sources of information can be efficiently exploited to improve the accuracy of shape prediction.