multinomial.theo.multcomp: Pairwise comparisons after an exact multinomial test for given probabilities

Description

Performs pairwise comparisons after a global exact multinomial test for given probabilities. These comparisons are performed using exact binomial tests.

Usage

multinomial.theo.multcomp(x, p = rep(1/length(x), length(x)), prop = FALSE,
  p.method = "fdr")

Value

method: name of the test.
data.name: a character string giving the name(s) of the data.
observed: observed counts.
expected: expected counts.
p.adjust.method: method for p-values correction.
p.value2: corrected p-values.
p.value: data frame of results.

Arguments

x: numeric vector (counts). Can also be a factor; in that case table(x) is used as counts.
p: theoretical proportions.
prop: logical indicating if results should be printed as counts (FALSE) or proportions (TRUE).
p.method: method for p-values correction. See help of p.adjust.

Author

Maxime HERVE <maxime.herve@univ-rennes1.fr>

Details

Since chi-squared and G tests are approximate tests, exact tests are preferable when the number of individuals is small (200 is a reasonable minimum).

An exact multinomial test with two groups is strictly the same than an exact binomial test.

Examples

Run this code

counts <- c(5,15,23)
p.theo <- c(0.2,0.5,0.3)
multinomial.test(counts,p=p.theo)
multinomial.theo.multcomp(counts,p=p.theo)

Run the code above in your browser using DataLab