[Rd] MASS:::dropterm.glm() and MASS:::addterm.glm() should use ... for extractAIC()
marc_grt at yahoo.fr
Wed Aug 23 14:12:53 CEST 2017
I have sent this message to this list the July, 7th. It was about a
problem in MASS package.
Until now there is no change in the devel version.
As the problem occurs in a package and not in the R-core, I don't know
if the message should have been sent here. Anyway, I have added a copy
to Pr Ripley.
I hope it could have been fixed.
Le 09/07/2017 à 16:05, Marc Girondot via R-devel a écrit :
> Here is a change required from MASS:::dropterm.glm() and
> The stepAIC() function from package MASS uses extractAIC() to get the
> AIC from a model.
> Several methods exist:
> extractAIC.glm() for example, some in MASS packages and some in stats
> The parameters for extractAIC() are:
> fit, scale, k = 2, ...
> The ... are not used in most of the extractAIC.xxx() methods, from
> example in stats:::extractAIC.glm() or MASS:::extractAIC.loglm() but
> its presence could be necessary if extractAIC() is changed to use for
> example to use AICc rather than AIC.
> Within stepAIC(), extractAIC() uses always the ... parameter. So all
> is ok for that.
> However, stepAIC() uses dropterm() or addterm().
> Within MASS:::dropterm.glm() and MASS:::addterm.glm(), extractAIC() is
> also used but without the ... parameter.
> It prevents to use new version of extractAIC() that could use this
> The solution is simple:
> In MASS:::dropterm.glm(), line 60 and MASS:::addterm.glm(), line 84:
> aic <- aic + (extractAIC(object, k = k)[2L] - aic[1L])
> must be changed to
> aic <- aic + (extractAIC(object, k = k, ...)[2L] - aic[1L])
> Other dropterm.xxx() and addterm.xxx() do not suffer this problem.
> R-devel at r-project.org mailing list
More information about the R-devel