R: Хорошее соответствие с logitfmx

Question

Я оценил некоторые бинарные логистические модели с предельными эффектами, используя mfx пакет в R logitmfx(), Теперь мне нужна добрая подгонка, например, в Нагелкерке R^2. С использованием NagelkerkeR2() из fmsb Пакет возвращает ошибку:

"Ошибка в -rr$null: неверный аргумент для унарного оператора"

Как я могу получить хорошее соответствие для моих моделей?

Спасибо!