chisq.theo.multcomp function

Pairwise comparisons after a chi-squared test for given probabilities

Performs pairwise comparisons after a global chi-squared test for given probabilities.

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

Arguments

• x: numeric vector (counts).
• p: theoretical proportions.
• p.method: method for p-values correction. See help of p.adjust.

Details

Since a chi-squared test is an approximate test, an exact test is preferable when the number of individuals is small (200 is a reasonable minimum). See multinomial.theo.multcomp in that case.

Returns

• 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.

• statistic: statistics of each test.

• p.value2: corrected p-values.

• p.value: data frame of results.

Author(s)

Maxime HERVE maxime.herve@univ-rennes1.fr

See Also

chisq.test, multinomial.test, multinomial.theo.multcomp

Examples

counts <- c(49,30,63,59)
p.theo <- c(0.2,0.1,0.45,0.25)
chisq.test(counts,p=p.theo)
chisq.theo.multcomp(counts,p=p.theo)