upper and lower bound of CI are equal, but no failure flag

[louisahsmith]

Sometimes the confint() function returns upper and lower CI bounds that are identical but the fit$return_flag is “SUCCESS”. Other times, on the same data, it doesn’t – in the latter case I don’t know whether they or the coefficient itself is wrong. Running the function twice in a row on the same object can lead to different results (unpredictable, but sometimes happens with the example below).

library(Cyclops)
test_dat <- data.frame(exposure = c(0, 1), 
                       person_time = c(13775, 10115), 
                       outcomes = c(0, 3))
cyclops_dat <- createCyclopsData(outcomes ~ exposure, data = test_dat,
                                 time = log(test_dat$person_time), modelType = "pr")
cyclops_fit <- fitCyclopsModel(cyclops_dat)
coef(cyclops_fit)
confint(cyclops_fit, "exposure")