A multivariate random-effects model with restricted parameters: Application to assessing radiation therapy for brain tumours

Abstract

In clinical studies, multiple endpoints are often measured for each patient longitudinally. The multivariate random-effects or random coefficient model has been a useful method for analysis. However, medical research problems may impose restrictions on the model parameters of interests. For example, in a paediatric brain tumour study on radiation therapy, there is a natural ordering in the white matter relaxation time of brain tissues among different regions surrounding the primary tumour, i.e. the closer a specific region of brain tissues is to the centre of primary tumour, the shorter is the relaxation time. Such parameter constraints should be accounted for in the analysis. This article proposes a class of multivariate random coefficient models with restricted parameters and derives its maximum likelihood estimates (MLE). We propose a modified EM algorithm for the quadratic optimalization with linear inequality constraints necessary in deriving the MLE. The method is applied to analysing the paediatric brain tumour study. Copyright © 2005 John Wiley & Sons, Ltd.