Интересно, видимо, устроен алгоритм оценивания параметров в Firth's penalized maximum likelihood (с profile penalized log likelihood выводом).
В коде тестовая статистика, соответствующая p-значению, вычисляется из p-значения, а не наоборот. Наверное, какие-то занятные вычислительные аспекты есть.