Firth’s model in r ” missing values in resampled performance measures” Here is an example with sample data