In multi-center randomized clinical trial the treatment eect may be changed over centers. It is thus important to investigate the heterogeneity in treatment eect between centers. For this, uncorrelated random-eect models assuming independence between random-eect terms have been often used, which may be a strong assumption. In this paper we propose a correlated frailty modelling approach of investigating such heterogeneity using the hierarchical-likelihood method when the outcome is time-to-event. In particular, we show how to construct a proper prediction interval for frailty, which explores graphically the potential heterogeneity for a treatment-by-center interaction term. The proposed method is illustrated via numerical studies based on data from the design of a multi-center clinical trial.