Advanced Search

Journal Navigation

Journal Home

Subscriptions

Archive

Contact Us

Table of Contents

Click here to sign up for SAGE Journal Email Alerts today!

Sign In to gain access to subscriptions and/or personal tools.
Statistical Modelling
This Article
Right arrow Full Text (PDF)
Right arrow References
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to Saved Citations
Right arrow Download to citation manager
Right arrow Add to My Marked Citations
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Right arrow Citing Articles via Scopus
Google Scholar
Right arrow Articles by Smyth, G. K
Right arrow Articles by Verbyla, A. P
Right arrow Search for Related Content
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?

Exact and approximate REML for heteroscedastic regression

Gordon K Smyth

Department of Mathematics, University of Queensland, Australia, smyth{at}wehi.edu.au

A Frederik Huele

Centre for Quantitative Methods BV, Eindhoven, The Netherlands

Arunas P Verbyla

Biometrics SA, University of Adelaide and South Australian Research and Development Institute, Australia

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 and Nelder (LN98) and Smyth and Verbyla (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 the same general form. A new algorithm is suggested, similar to one suggested by Huele et al., 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.

Key Words: maximum likelihood • modified profile likelihood • residual maximum likelihood • restricted maximum likelihood • sandwich estimator

Statistical Modelling, Vol. 1, No. 3, 161-175 (2001)
DOI: 10.1177/1471082X0100100301


Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Reddit Reddit   Add to Technorati Technorati   Add to Twitter Twitter    What's this?


This article has been cited by other articles:


Home page
Statistical ModellingHome page
S. Kuhnt and M. Erdbrugge
A strategy of robust parameter design for multiple responses
Statistical Modeling, December 1, 2004; 4(4): 249 - 264.
[Abstract] [PDF]


Home page
Statistical ModellingHome page
J. Taylor and A. Verbyla
Joint modelling of location and scale parameters of the t distribution
Statistical Modeling, July 1, 2004; 4(2): 91 - 112.
[Abstract] [PDF]