来自logistf的森林情节

Question

我在R中运行了几个模型用于惩罚后勤模型 logistf包。不过，我希望为数据绘制一些森林图。

sjPlot包：http://www.strengejacke.de/sjPlot/custplot/

为glm输出提供了出色的功能，但没有logistf函数的功能。

有任何帮助吗？

Answer 1

与logistf对象相比，glm个对象的结构不同，但不是太多。我已经添加了对logistf拟合模型的支持，但是，1）无法打印模型摘要，b）logistf模型目前不支持预测概率图。

我今晚会更新the code on GitHub，因此您可以尝试更新的sjp.glm功能......

library(sjPlot)
library(logistf)
data(sex2)
fit<-logistf(case ~ age+oc+vic+vicl+vis+dia, data=sex2)
# for this example, axisLimits need to be specified manually
sjp.glm(fit, axisLimits = c(0.05, 25), transformTicks = T)