good

0th

Percentile

Approximate number of contingency tables

Approximate number of contingency tables with specified marginal totals.

Keywords
array
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.

Note

Method “A” is the exact answer, given by no.of.boards(). Do not use this on large matrices! Methods “B”, “C”, and “D” give (according to Good) increasingly accurate approximations to the exact answer.

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

no.of.boards

• good
Examples
a <- matrix(1,15,15)

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


Documentation reproduced from package aylmer, version 1.0-11, License: GPL-2

Community examples

Looks like there are no examples yet.