Volume 116, Number 1, January-February 2011 Journal of Research of the National Institute of Standards and Technology
[J. Res. Natl. Inst. Stand. Technol. 116, 539-556 (2011)] Maximum Likelihood and Restricted Likelihood Solutions in Multiple-Method Studies
Volume 116 Number 1 January-February 2011
Andrew L. Rukhin A formulation of the problem of considered to be known an upper bound combining data from several sources is on the between-method variance is National Institute of Standards discussed in terms of random effects obtained. The relationship between and Technology, models. The unknown measurement likelihood equations and moment-type Gaithersburg, MD 20899-8980 precision is assumed not to be the same equations is also discussed. for all methods. We investigate maximum likelihood solutions in this model. By representing the likelihood equations as Key words: DerSimonian-Laird simultaneous polynomial equations, estimator; Groebner basis; [email protected] the exact form of the Groebner basis for heteroscedasticity; interlaboratory studies; their stationary points is derived when iteration scheme; Mandel-Paule algorithm; there are two methods. A parametrization meta-analysis; parametrized solutions; of these solutions which allows their polynomial equations; random effects comparison is suggested. A numerical model. method for solving likelihood equations is outlined, and an alternative to the maximum likelihood method, the restricted Accepted: August 20, 2010 maximum likelihood, is studied. In the situation when methods variances are Available online: http://www.nist.gov/jres
1. Introduction: Meta Analysis and difficult when the unknown measurement precision Interlaboratory Studies varies among methods whose summary results may not seem to conform to the same measured property. This article is concerned with mathematical aspects Reference [6] provides some practical suggestions for of an ubiquitous problem in applied statistics: how to dealing with the problem. combine measurement data nominally on the same In this paper we investigate the ML solutions of the property by methods, instruments, studies, medical random effects model which is formulated in Sec. 2. By centers, clusters or laboratories of different precision. representing the likelihood equations as simultaneous One of the first approaches to this problem was sug- polynomial equations, the so-called Groebner basis gested by Cochran [1], who investigated maximum for them is derived when there are two sources. A likelihood (ML) estimates for one-way, unbalanced, parametrization of such solutions is suggested in heteroscedastic random-effects model. Cochran re- Sec. 2.1. The maxima of the likelihood function are turned to this problem throughout his career [2, 3], and compared for positive and zero between-labs variance. Ref. [4] reviews this work. Reference [5] discusses A numerical method for solving likelihood equations applications in metrology and gives more references. by reducing them to an optimization problem for a The problem of combining data from several sources homogeneous objective function is given in Sec. 2.2. is central to the broad field of meta-analysis. It is most An alternative to the ML method, the restricted maxi-
539 Volume 116, Number 1, January-February 2011 Journal of Research of the National Institute of Standards and Technology
σ^ 2 mum likelihood is considered in Sec. 3. An explicit i = 0. In order to find these estimates one can replace formula for the restricted likelihood estimator is μ in (2) by discovered in Sec. 3.2 in the case of two methods. −1 Section 4 deals with the situation when methods ∑ x ()σσ22+ μ = i ii, (3) variances are considered to be known, and an upper σσ221+ − ∑ ()i bound on the between-method variance is obtained. i The Sec. 5 discusses the relationship between likeli- which reduces the number of parameters from p +2 to hood equations and moment-type equations, and p +1. Sec. 6 gives some conclusions. All auxiliary material Our goal is to represent the set of all stationary points related to an optimization problem and to elementary of the likelihood equations as solutions to simultaneous symmetric polynomials is collected in the Appendix. polynomial equations. To that end, note that
−=μ 22 − 2 2. ML Method and Polynomial Equations ∑∑∑∑wi() x i wx ii ( wx ii )/() w i iiii(4) =−2 To model interlaboratory testing situation, denote by ∑∑wwij()/(). x i x j w i 1≤