Glm2

Fitting generalized linear models with convergence problems

Ian C. Marschner*

*Corresponding author for this work

    Research output: Contribution to journalArticle

    64 Citations (Scopus)

    Abstract

    The R function glm uses step-halving to deal with certain types of convergence problems when using iteratively reweighted least squares to fit a generalized linear model. This works well in some circumstances but non-convergence remains a possibility, particularly with a nonstandard link function. In some cases this is because step-halving is never invoked, despite a lack of convergence. In other cases step-halving is invoked but is unable to induce convergence. One remedy is to impose a stricter form of stephalving than is currently available in glm, so that the deviance is forced to decrease in every iteration. This has been implemented in the glm2 function available in the glm2 package. Aside from a modified computational algorithm, glm2 operates in exactly the same way as glm and provides improved convergence properties. These improvements are illustrated here with an identity link Poisson model, but are also relevant in other contexts.

    Original languageEnglish
    Pages (from-to)12-15
    Number of pages4
    JournalR Journal
    Volume3
    Issue number2
    Publication statusPublished - 2011

    Fingerprint Dive into the research topics of 'Glm2: Fitting generalized linear models with convergence problems'. Together they form a unique fingerprint.

  • Cite this