[R-sig-ME] 2 level logit, 2 REs, large sample - log likelihood returns "NaN"

Daniel Adkins deadkins at vcu.edu
Fri Jun 10 10:08:36 CEST 2011

I am fitting a large (j=50K, i=9K) 2-level logit with random intercept
and age slope and 14 covariates. Model estimates become stable at
nAGQ>=150 (large, I know). Based on simpler models (random intercept
only, random
slope only, ordinary logit, etc) the solution looks sound. However,
all the fit indices return a value of "NaN", which naturally stands
for "not a number". Why is this? This model should yield a scalar log
likelihood, no? Any advice would be appreciated.


Daniel E. Adkins, PhD
Assistant Professor
Center for Biomarker Research and Personalized Medicine
School of Pharmacy
Virginia Commonwealth University
McGuire Hall, Room 216B
1112 East Clay Street
Richmond, VA 23298

More information about the R-sig-mixed-models mailing list