I have run glm (for a logistic regression) and as result I got that the algorithm diverges.
I want to try the function logistf to check whether Firth's approach can improve my data analysis.
Let dat denote my dataframe and y ~ f1 + f2 be my model. Trying
fit <- logistf(y ~ f1 + f2, dat)
or
fit <- logistf(y ~ f1 + f2, data=dat)
I always get the error
Error in logistf.fit(x = x, y = y, weight = weight, offset = offset, firth, :
no loop for break/next, jumping to top level
As glm runs smoothly, I do not see how my data structure could be the source of error (this is why I have not described dat explicitly). Do you see any potential conflict? I thank you very much,
Avitus