Computing likelihood ratio pseudo-R2

Question

Andreas 2014-1-17

编辑： Arseny 2017-9-4

I wish to compute the likelihood ratio pseudo-R2 from a logistic regression using glmfit.

According to the Wikipedia definition it is defined as

R2 = (DevianceNull-DevianceModel)/(DevianceNull).

Does DevianceModel correspond to dev returned from glmfit?

[~, dev] = glmfit(X, Y, 'binomial');

In that case, is the pseudo-R2 computed as?

[~, devModel] = glmfit(X, Y, 'binomial'); 
[~, devNull]  = glmfit(zeros(length(Y), 1), Y, 'binomial');
R2 = (DevNull-DevModel)/(DevNull);

Many thanks

Andreas 2014-1-20

Nobody?

Answer 1

Arseny 2017-9-4

编辑：Arseny 2017-9-4

It would be your code, but with

[~, devNull] = glmfit(X, Y(randperm(length(Y))), 'binomial');

Or better yet, this row repeated several times, and devNull estimation averaged over these several attempts.