aylmer (version 1.0-10)

good: Approximate number of contingency tables

Description

Approximate number of contingency tables with specified marginal totals.

Usage

good(x, method = "D", ...)

Arguments

x
An integer-valued matrix (with no NA entries)
method
The method to use; one of A, B, C, D. See details section
...
Further arguments (notably n), passed to no.of.boards()

Details

All formulae and terminology taken from Good 1976. The letters A-D are from the approximations given on pages 1167 and 1168. Note This function will accept matrices with any NA entries, but a warning is given. The number of permissable boards will be less than the number of permissible contingency tables considered by Good.

The approximations are intended to provide rough-and-ready upper bounds for boards that have a few NAs.

References

  • I. J. Good 1976.On the Application of Symmetric Dirichlet Distributions and Their Mixtures to Contingency Tables. The Annals of Statistics 4(6):1159--1189

See Also

no.of.boards

Examples

Run this code
a <- matrix(1,15,15)

good(a,"B")
good(a,"C")
good(a,"D")  #identical to answer given by method "C"

Run the code above in your browser using DataCamp Workspace