Support for the emmeans package is available. See below for an example of using emmeans::emmeans() with a logistf object.

Examples


data(sex2)
fit<-logistf(case ~ age+oc+vic+vicl+vis+dia, data=sex2)

emmeans::emmeans(fit, ~age+dia)
#>  age dia emmean    SE  df lower.CL upper.CL
#>    0   0  -0.23 0.184 232  -0.5924    0.133
#>    1   0  -1.34 0.408 232  -2.1387   -0.532
#>    0   1   2.87 1.500 232  -0.0852    5.818
#>    1   1   1.76 1.490 232  -1.1841    4.705
#> 
#> Results are averaged over the levels of: oc, vic, vicl, vis 
#> Results are given on the logit (not the response) scale. 
#> Confidence level used: 0.95