Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

st: xtmelogit: comparing models

From   Luca Campanelli <>
To   "" <>
Subject   st: xtmelogit: comparing models
Date   Fri, 5 Oct 2012 16:50:32 +0100 (BST)

Dear Stata users, 
I’d like to fit and compare mixed effects logistic regression models with crossed random effects using the function xtmelogit (Stata 12IC for Windows). 

For example (“group” has 2 levels[0,1] and “condition” has 3 levels[1,2,3]): 
(1) xtmelogit resp i.condition , || _all: R.item, covariance(id)  || sbj: , covariance(id)
(2) xtmelogit resp i.condition , || _all: R.item, covariance(id)  || sbj: , covariance(id)

In comparing two models, I found a big discrepancy between lrtest on one side, and AIC-BIC on the other side. lrtest was highly significant, indicating that (2) was better than (1), while AIC and BIC values were clearly smaller for model (1). 
Which should I trust? 

Does this apply to my case ; ? 
If yes, how can I do the Wald test? 
Would it be: 

Is it correct? I saw others using testparm or lincom. 
I would appreciate any help to understand what the appropriate thing to do is. 

thank you, 

*   For searches and help try:

© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index