Information methods for model selection in linear mixed effects models with application to HCV data

被引:12
作者
Dimova, Rositsa B. [2 ]
Markatou, Marianthi [1 ]
Talal, Andrew H. [3 ]
机构
[1] IBM TJ Watson Res Ctr, New York, NY 10532 USA
[2] Columbia Univ, Mailman Sch Publ Hlth, Dept Biostat, New York, NY 10032 USA
[3] Weill Cornell Med Coll, Div Gastroenterol & Hepatol, New York, NY 10065 USA
关键词
Model selection; Linear mixed effects models; AIC; REML information criteria; CHRONIC HEPATITIS-C; ALPHA-2A PLUS RIBAVIRIN; PEGYLATED INTERFERON-ALPHA-2B; AKAIKE INFORMATION; INFECTED PATIENTS; VIRUS; LIKELIHOOD; CRITERION; PHARMACODYNAMICS; ASSOCIATION;
D O I
10.1016/j.csda.2010.10.031
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
In this paper, we derive a small sample Akaike information criterion, based on the maximized loglikelihood, and a small sample information criterion based on the maximized restricted loglikelihood in the linear mixed effects model when the covariance matrix of the random effects is known. Small sample corrected information criteria are proposed for a special case of linear mixed effects models, the balanced random-coefficient model, without assuming the random coefficients covariance matrix to be known. A simulation study comparing the derived criteria and several others for model selection in the linear mixed effects models is presented. We illustrate the behavior of the studied information criteria on real data from a study of subjects coinfected with HIV and Hepatitis C virus. Robustness of the criteria, in terms of the error distributed as a mixture of normal distributions, is also studied. Special attention is given to the behavior of the conditional AIC by Vaida and Blanchard (2005). Among the studied criteria, GIC performs best, while cAIC exhibits poor performance. Because of its inferior performance, as demonstrated in this work, we do not recommend its use for model selection in linear mixed effects models. (C) 2011 Elsevier B.V. All rights reserved.
引用
收藏
页码:2677 / 2697
页数:21
相关论文
共 41 条
[11]  
Demidenko E., 2004, MIXED MODELS THEORY
[12]   Association of pretreatment serum interferon γ inducible protein 10 levels with sustained virological response to peginterferon plus ribavirin therapy in genotype 1 infected patients with chronic hepatitis C [J].
Diago, M ;
Castellano, G ;
García-Samaniego, J ;
Pérez, C ;
Fernández, I ;
Romero, M ;
Iacono, OL ;
García-Monzón, C .
GUT, 2006, 55 (03) :374-379
[13]   Pegylated interferon-α2b:: Pharmacokinetics, pharmacodynamics, safety, and preliminary efficacy data [J].
Glue, P ;
Fang, JWS ;
Rouzier-Panis, R ;
Raffanel, C ;
Sabo, R ;
Gupta, SK ;
Salfi, M ;
Jacobs, S .
CLINICAL PHARMACOLOGY & THERAPEUTICS, 2000, 68 (05) :556-567
[14]   Selecting the best linear mixed model under REML [J].
Gurka, MJ .
AMERICAN STATISTICIAN, 2006, 60 (01) :19-26
[15]  
HARVILLE DA, 1977, J AM STAT ASSOC, V72, P320, DOI 10.2307/2286796
[16]   Counting degrees of freedom in hierarchical and other richly-parameterised models [J].
Hodges, JS ;
Sargent, DJ .
BIOMETRIKA, 2001, 88 (02) :367-379
[17]   REGRESSION AND TIME-SERIES MODEL SELECTION IN SMALL SAMPLES [J].
HURVICH, CM ;
TSAI, CL .
BIOMETRIKA, 1989, 76 (02) :297-307
[18]   Bootstrapping log likelihood and EIC, an extension of AIC [J].
Ishiguro, M ;
Sakamoto, Y ;
Kitagawa, G .
ANNALS OF THE INSTITUTE OF STATISTICAL MATHEMATICS, 1997, 49 (03) :411-434
[19]   Robustness of the linear mixed model to misspecified error distribution [J].
Jacqmin-Gadda, Helene ;
Sibillot, Solenne ;
Proust, Cecile ;
Molina, Jean-Michel ;
Thiebaut, Rodolphe .
COMPUTATIONAL STATISTICS & DATA ANALYSIS, 2007, 51 (10) :5142-5154
[20]  
Jiang J., 2003, Sankhya, V65, P23