Subject-specific musculoskeletal models can predict accurate joint and muscle biomechanics thereby helping clinicians and surgeons. Current modeling strategies do not incorporate accurate subject-specific muscle parameters. This study reports a statistical shape model (SSM) based method to predict subject-specific muscle attachment regions on shoulder bones and illustrates the concurrent validity of the predictions. Augmented SSMs of scapula and humerus bones were built using bone meshes and five muscle attachment (origin/insertion) regions which play important role in the shoulder motion and function. Muscle attachments included Subscapularis, Supraspinatus, Infraspinatus, Teres Major and Teres Minor on both the bones. The regions were represented by subset of vertices on the bone meshes and were tracked using vertex identifiers. Subject-specific muscle attachment regions were predicted using external set of bones not used in building the SSMs. Validity of predictions was determined by visual inspection and also by using four similarity measures between predicted and manually segmented regions. Excellent concurrent validity was found indicating the higher accuracy of predictions. This method can be effectively employed in modeling pipelines or in automatic segmentation of medical images. Further validations are warranted on all the muscles of the shoulder complex.