Scale parameter estimation in
gam depends on the type of
extended families then the RE/ML estimate is used. For conventional exponential families, estimated by the default
outer iteration, the scale estimator can be controlled using argument
gam.control. The options are
The Pearson estimator is the (weighted) sum of squares of the pearson residuals, divided by the
effective residual degrees of freedom. The Fletcher (2012) estimator is an improved version of the Pearson estimator.
The deviance estimator simply substitutes deviance residuals for Pearson residuals.
Usually the Pearson estimator is recommended for GLMs, since it is asymptotically unbiased. However, it can also be unstable at finite sample sizes, if a few Pearson residuals are very large. For example, a very low Poisson mean with a non zero count can give a huge Pearson residual, even though the deviance residual is much more modest. The Fletcher (2012) estimator is designed to reduce these problems.
For performance iteration the Pearson estimator is always used.
gamm uses the estimate of the scale parameter from the underlying call to
bam uses the REML estimator if the method is
"fREML". Otherwise the estimator is a Pearson estimator.
Simon N. Wood email@example.com with help from Mark Bravington and David Peel
Fletcher, David J. (2012) Estimating overdispersion when fitting a generalized linear model to sparse data. Biometrika 99(1), 230-237.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.