Motivated by a study of muscular dystrophy in MD STARnet,a joint model of bivariate survival times and longitudinal data is developed.We propose to analyze correlated bivariate survival responses associated with a longitudinal biomarker in the Frequentist paradigm.A Gamma frailty variable is used to account for the correlation between the two correlated survival outcomes in addition to the random variables that account for the correlation between the survival times and longitudinal maker.The EM algorithm is adopted to compute the maximum profile likelihood estimate.The bootstrap method is applied to estimate the standard error of estimated model parameters.The simulation study is conducted to demonstrate the validity of the proposed methodology.Finally the method is applied to the MD STARnet for illustration.