Multiple outcomes measured repeatedly for the same subject are common in longitudinal observation.If we use the approach by analyzing each outcome separately,it may lead to wrong conclusions due to the failure of accounting for joint evolution of different outcomes.To adequately capture the interdependence among multiple outcomes,we proposed a joint modeling for multivariate longitudinal data by constructing a linear mixed-effects model for each outcome and accommodating the relationship among multiple outcomes through correlation in random effects.Maximum likelihood method was adopted to estimate parameters in this model.The application of this method was demonstrated through the analysis of stroke data.