[R] Changes to stats::glm function between R versions 3.4.0 and 3.5.1

[Purver, Mark]

Dear Martin,

Thanks so much for this response. I have attempted a reproducible example below. The 'real' issue involves a much larger data set, but the very small example below produces very small differences in the standard errors between different R versions, which may have the same root cause. I hope to understand whether such differences indicate any change to the 'glm' function between version 3.4.0 of the 'stats' package and later versions. These differences are only apparent when the model does not converge.

Thanks,
Mark


dat <- data.frame(
  a = c(1, 3, 8, 9, 2, 3, 4, 5),
  b = c(1, 3, 8, 8, 2, 3, 4, 5),
  x = c(1, 0, 0, 1, 0, 0, 0, 0)
)

m <- glm(x~a+b, family=binomial(), data=dat)

print(summary(m)))

# Standard errors of variables 'a' and 'b' are 580356.74 and 638232.37 in R 3.4.0, but 580356.25 and 638231.93 in R 3.4.2
________________________________
 This e-mail and any attachments is intended only for the attention of the addressee(s). Its unauthorised use, disclosure, storage or copying is not permitted. If you are not the intended recipient, please destroy all copies and inform the sender by return e-mail. Internet e-mail is not a secure medium. Any reply to this message could be intercepted and read by someone else. Please bear that in mind when deciding whether to send material in response to this message by e-mail. This e-mail (whether you are the sender or the recipient) may be monitored, recorded and retained by the Ministry of Justice. Monitoring / blocking software may be used, and e-mail content may be read at any time. You have a responsibility to ensure laws are not broken when composing or forwarding e-mails and their contents.