Display loglikihood value in first order MQL

Welcome to the forum for R2MLwiN users. Feel free to post your question about R2MLwiN here. The Centre for Multilevel Modelling take no responsibility for the accuracy of these posts, we are unable to monitor them closely. Do go ahead and post your question and thank you in advance if you find the time to post any answers!

Go to R2MLwiN: Running MLwiN from within R >> http://www.bris.ac.uk/cmm/software/r2mlwin/
Post Reply
Posts: 11
Joined: Wed May 06, 2015 9:22 am

Display loglikihood value in first order MQL

Post by pimverbunt »

Dear all,

I am modelling a single level unordered multinomial logistic model in r2mlwin. Because my model is single level, I have to estimate it in first order MQL. In MQL, however, the likelihood value is not shown. To compare my results I would like to obtain a likelihood value (on the website I found that is possible to display the likelihood value for IGLS estimation).

Is it possible to obtain such a measure (or an approximation) in MQL? And if so, what code do i need to enter?

This is an example of the simple code I am currently using:

model_a = logit(dep, denom) ~ 1 + var1 + var2
model_a<- runMLwiN(model_a, D = "Unordered Multinomial", data =mydata)

Kind regards,

Pim Verbunt
Posts: 1248
Joined: Mon Oct 19, 2009 10:34 am

Re: Display loglikihood value in first order MQL

Post by ChrisCharlton »

Unfortunately as discrete response models estimated in MLwiN using (R)IGlS use quasi-likelihood it is not possible to obtain an accurate log-likelihood. For an explanation in the context of the -runmlwin- command see https://www.cmm.bristol.ac.uk/forum/viewtopic.php?t=279.
Post Reply