I have found the problem and a solution.
The problem comes from the functions addterm.glm() and dropterm.glm()
from MASS package. They use extractAIC() without transmitting the ...
argument:
aic <- aic + (extractAIC(object, k = k)[2L] - aic[1L])
I replace the call with:
aic <- aic + (extract
I would like test AICc as a criteria for model selection for a glm using
stepAIC() from MASS package.
Based on various information available in WEB, stepAIC() use
extractAIC() to get the criteria used for model selection.
I have created a new extractAIC() function (and extractAIC.glm() and
e
2 matches
Mail list logo