Exact and approximate REML for heteroscedastic regression

G.K. Smyth, A.F. Huele, A.P. Verbyla
2001 Statistical Modelling  
Exact REML for heteroscedastic linear models is compared with a number of approximate REML methods which have been proposed in the literature, especially with the methods proposed by Lee & Nelder (1998) (LN98) and Smyth & Verbyla (1999) (SV99) for simultaneous mean-dispersion modelling in generalized linear models. It is shown that neither of the LN98 or SV99 methods reduces to REML in the normal linear case. Asymptotic variances and efficiencies are obtained for these and other estimators of
more » ... e same general form. A new algorithm is suggested, similar to one suggested by Huele et al. (2000) , which returns the correct REML estimators and an improved approximation to the standard errors. It is possible to obtain REML estimators by alternating between two generalized linear models but the final fitted generalized linear model objects will not return the correct standard errors for the variance coefficients. The true REML likelihood calculations therefore fit only partially into the double generalized linear model framework.
doi:10.1191/147108201128140 fatcat:xd7hwoubrrgfxffl2syttjtetu