The topic of studying the growth of human brain development has become of increasing interest in the neuroimaging community. Cross-sectional studies may allow comparisons between means of different age groups, but they do not provide a growth model that integrates the continuum of time, nor do they present any information about how individuals/population change over time. Longitudinal data analysis method arises as a strong tool to address these questions. In this paper, we use longitudinal analysis methods to study tissue development in early brain growth. A novel approach of multivariate longitudinal analysis is applied to study the associations between the growth of different brain tissues. In this paper, we present the methodologies to statistically study scalar (univariate) and vector (multivariate) longitudinal data, and demonstrate exploratory results in a neuroimaging study of early brain tissue development. We obtained growth curves as a quadratic function of time for all three tissues. The quadratic terms were tested to be statistically significant, showing that there was indeed a quadratic growth of tissues in early brain development. Moreover, our result shows that there is a positive correlation between repeated measurements of any single tissue, and among those of different tissues. Our approach is generic in natural and thus can be applied to any longitudinal data with multiple outcomes, even brain structures. Also, our joint mixed model is flexible enough to allow incomplete and unbalanced data, i.e. subjects do not need to have the same number of measurements, or be measured at the exact time points.