| good {aylmer} | R Documentation |
Approximate number of contingency tables with specified marginal totals.
good(x, method = "D", ...)
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() |
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.
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.
Robin K. S. Hankin
a <- matrix(1,15,15) good(a,"B") good(a,"C") good(a,"D") #identical to answer given by method "C"