Error in if (abs(dev - devold)/(0.1 + abs(dev)) < control$epsilon)

[monicawang-i]

Hi,

I'm currently using gnm for constructing conditional poisson models. I had the model specified as
gnm(event_count ~ factor(exposure) + cov1 + ... + cov n, data = data, family = 'poisson', offset = logtime, eliminate = factor(analytical_unit))

I encountered the error message as Error in if (abs(dev - devold)/(0.1 + abs(dev)) < control$epsilon) { : missing value where TRUE/FALSE needed

However, when I tried the formula specified differently, using a subset of covariates, it would run. I also tried to run the same formula on a dataset with same structure, the function would run as well. Additionally, I checked the data that I feed into the function had no missing values.

I could not find documentation about this error on stackoverflow or CRAN. Do you know where the problem is? Thanks you so much!

Best,
Monica

[Alessa19]

Exactly same issue here. Did you find out what was causing the error?