Rでロジスティック回帰を行ったときのdeviance残差の計算方法

Question

Rでロジスティック回帰を行った際のresidual devianceを取得したいと考えています。
glm.fit関数を用いたところ、以下のような完全分離の警告が出ました。

glm.fit: fitted probabilities numerically 0 or 1 occurred 

そこで、こちらを参考にlogistf関数を利用してみました。

cont <- logistf.control(maxit=1000,maxstep = 100)
plcont <- logistpl.control(maxit=1000)
res <- logistf(y ~ x,firth = TRUE, control = cont,plcontrol = plcont)

summary(res)の中身は以下のとおりです。

> summary(res)
logistf(formula = dat[, j] ~ Dat_res[, i], control = cont, plcontrol = plcont, firth = TRUE)

Model fitted by Penalized ML
Coefficients:
                   coef   se(coef)  lower 0.95 upper 0.95    Chisq        p
(Intercept)   0.1375199 0.11492692 -0.08696857  0.3642894 1.439322 0.230249
Dat_res[, i] -0.7204668 0.05044381 -0.82565105 -0.6269034      Inf 0.000000
             method
(Intercept)       2
Dat_res[, i]      2

Method: 1-Wald, 2-Profile penalized log-likelihood, 3-None

Likelihood ratio test=895.6586 on 1 df, p=0, n=1000
Wald test = 203.9916 on 1 df, p = 0

この結果から、こちらのサイトでも説明されているdeviance残差(residual deviance)を計算しようと考えていますが、導出する方法がよく分かっていません。
理解不足で申し訳ありませんが、ご教授よろしくお願いします。