Counting tables using the double-saddlepoint approximation

Vadim Zipunnikov, James G. Booth, Ruriko Yoshida

Research output: Contribution to journalArticle

Abstract

A double-saddlepoint approximation is proposed for the number of contingency tables with counts satisfying certain linear constraints. Computation of the approximation involves fitting a generalized linear model for geometric responses which can be accomplished almost instantaneously using the iterated weighted least squares algorithm. The approximation is far superior to other analytical approximations that have been proposed, and is shown to be highly accurate in a range of examples, including some for which analytical approximations were previously unavailable. A similar approximation is proposed for tables consisting of only zeros and ones based on a logistic regression model. A higher order adjustment to the basic double saddlepoint further improves the accuracy of the approximation in almost all cases. Computer code for implementing methods described in the article is provided as supplemental material.

Original languageEnglish (US)
Pages (from-to)915-929
Number of pages15
JournalJournal of Computational and Graphical Statistics
Volume18
Issue number4
DOIs
StatePublished - Dec 28 2009
Externally publishedYes

Keywords

  • Contingency tables
  • Darwin's finch data
  • Generalized linear model
  • Quasi-independence
  • Uniform association

ASJC Scopus subject areas

  • Statistics and Probability
  • Discrete Mathematics and Combinatorics
  • Statistics, Probability and Uncertainty

Fingerprint Dive into the research topics of 'Counting tables using the double-saddlepoint approximation'. Together they form a unique fingerprint.

  • Cite this