Continued fraction representation for expected cell counts of a 2 x 2 table: a rapid and exact method for conditional maximum likelihood estimation.
The expected cell count for a 2 x 2 contingency table, governed by the noncentral (extended) hypergeometric distribution, is expressed as a terminating continued fraction. The coefficients in the continued fraction are better behaved than the multinomial coefficients required for the usual moment calculation. The expected cell count must be calculated repeatedly in a conditional maximum likelihood analysis of K2 x 2 contingency tables. Since the continued fraction can be easily evaluated, a rapid and numerically stable computational algorithm results. Once this first moment is known, higher moments can be obtained as shown by Harkness (1965, Annals of Mathematical Statistics 36, 938-945). A BASIC program to implement the continued fraction algorithm is given in an appendix.
Satten
G A
GA
Department of Biostatistics, School of Public Health, University of North Carolina, Chapel Hill 27599-7400.
Kupper
L L
LL
D000465:Algorithms; D002452:Cell Count; D016013:Likelihood Functions; D015233:Models, Statistical; D016017:Odds Ratio; D012984:Software