Possible bug in confint.profile.gnm()

[dwoll]

I believe there is a subsetting bug in confint.profile.gnm() which is triggered when there is only 1 covariate. Consider

library(gnm)
d <- as.data.frame(cautres)
fit <- gnm(Freq ~ vote, eliminate=class, family=poisson, data=d)
confint(fit)

For me, this returns an error message which appears to be due to lines 30-31 in file confint.profile.gnm.R:

std.err <- attr(object, "summary")$coefficients[, "Std. Error"]
parm <- parm[!is.na(std.err)[parm]]

After replacing these lines with

std.err <- attr(object, "summary")$coefficients[, "Std. Error", drop=FALSE]
parm <- parm[!is.na(std.err)[parm, ]]

The error disappears.

Many thanks for developing gnm and making it freely available! Conditional Poisson regression is very useful in epidemiology.

[hturner]

Thanks for the comprehensive bug report! This should be fixed in the GitHub version ready for the next release.