Multilevel Modeling of Cognitive Function in Schizophrenic Patients and Their First Degree Relatives

We describe multilevel modeling of cognitive function in subjects with schizophrenia, their healthy first degree relatives and controls. The purpose of the study was to compare mean cognitive performance between the three groups after adjusting for various covariates, as well as to investigate differences in the variances. Multilevel models were required because subjects were nested within families and some of the measures were repeated several times on the same subject. The following four methodological issues that arose during the analysis of the data are discussed. First, when the random effects distribution was not normal, non-parametric maximum likelihood (NPML) was employed, leading to a different conclusion than the conventional multilevel model regarding one of the main study hypotheses. Second, the between-subject (within-family) variance was allowed to differ between the three groups. This corresponded to the variance at level 1 or level 2 depending on whether repeated measures were analyzed. Third, a positively skewed response was analyzed using a number of different generalized linear mixed models. Finally, penalized quasilikelihood (PQL) estimates for a binomial response were compared with estimates obtained using Gaussian quadrature. A small simulation study was carried out to assess the accuracy of the latter.

[1]  N. Breslow,et al.  Approximate inference in generalized linear mixed models , 1993 .

[2]  S. Faraone,et al.  Neuropsychological functioning among the nonpsychotic relatives of schizophrenic patients: a diagnostic efficiency analysis. , 1995, Journal of abnormal psychology.

[3]  Cheng Hsiao,et al.  Analysis of Panel Data , 1987 .

[4]  K. Kendler,et al.  The Roscommon Family Study. IV. Affective illness, anxiety disorders, and alcoholism in relatives. , 1993, Archives of general psychiatry.

[5]  N. Laird.  Nonparametric Maximum Likelihood Estimation of a Mixing Distribution , 1978 .

[6]  M. Aitkin.  A General Maximum Likelihood Analysis of Variance Components in Generalized Linear Models , 1999, Biometrics.

[7]  S. Rabe-Hesketh,et al.  Generalized linear latent and mixed models , 2000 .

[8]  R. Gur,et al.  Neuropsychological function in schizophrenia. Selective impairment in memory and learning. , 1991, Archives of general psychiatry.

[9]  Diane Lambert,et al.  Generalizing Logistic Regression by Nonparametric Mixing , 1989 .

[10]  Mark Von Tress,et al.  Generalized, Linear, and Mixed Models , 2003, Technometrics.

[11]  A. Farmer,et al.  The Strength of the Genetic Effect , 1994, British Journal of Psychiatry.

[12]  P. Sham,et al.  Statistical methods in psychiatric genetics , 1998, Statistical methods in medical research.

[13]  N. Breslow,et al.  Bias correction in generalised linear mixed models with a single component of dispersion , 1995 .

[14]  N. Breslow,et al.  Bias Correction in Generalized Linear Mixed Models with Multiple Components of Dispersion , 1996 .

[15]  R. Gur,et al.  Neuropsychological deficits in neuroleptic naive patients with first-episode schizophrenia. , 1994, Archives of general psychiatry.

[16]  Pickles A SkrondalA. Rabe-HeskethS.  GLLAMM: A general class of multilevel models and a STATA programme , 2001 .

[17]  S. Faraone,et al.  Neuropsychological functioning among the elderly nonpsychotic relatives of schizophrenic patients , 1996, Schizophrenia Research.

[18]  J. Heckman,et al.  A Method for Minimizing the Impact of Distributional Assumptions in Econometric Models for Duration Data , 1984 .

[19]  Richard B. Davies,et al.  A joint trip timing store-type choice model for grocery shopping, including inventory effects and nonparametric control for omitted variables , 1987 .

[20]  P. Diggle.  Analysis of Longitudinal Data , 1995 .

[21]  S. Faraone,et al.  Neuropsychological risk indicators for schizophrenia: a review of family studies. , 1994, Schizophrenia bulletin.

[22]  P. Diggle.  An approach to the analysis of repeated measurements. , 1988, Biometrics.

[23]  L. Simar.  Maximum Likelihood Estimation of a Compound Poisson Process , 1976 .

[24]  Eric R. Ziegel,et al.  Generalized Linear Models , 2002, Technometrics.

[25]  Sophia Rabe-Hesketh,et al.  GLLAMM - Generalized Linear Latent and Mixed Models , 2000 .

[26]  R. Darrell Bock,et al.  Fitting a response model forn dichotomously scored items , 1970 .

[27]  T. Sigmundsson,et al.  Minor physical anomalies in familial and sporadic schizophrenia: the Maudsley family study , 1998, Journal of neurology, neurosurgery, and psychiatry.

[28]  Harvey Goldstein,et al.  Improved Approximations for Multilevel Models with Binary Responses , 1996 .

[29]  B. Lindsay,et al.  Semiparametric Estimation in the Rasch Model and Related Exponential Response Models, Including a Simple Latent Class Model for Item Analysis , 1991 .

[30]  C. Morris.  Parametric Empirical Bayes Inference: Theory and Applications , 1983 .

[31]  Stephen D. Kachman,et al.  AN INTRODUCTION TO GENERALIZED LINEAR MIXED MODELS , 2001 .

[32]  R. Morris,et al.  Neural correlates of planning ability: Frontal lobe activation during the tower of London test , 1993, Neuropsychologia.

[33]  Noreen Goldman,et al.  An assessment of estimation procedures for multilevel models with binary responses , 1995 .