chisq {broman}

Chi-square test by simulation for a two-way table

Description

Calculate a p-value for a chi-square test by Monte Carlo simulation.

Usage

chisq(tab, n.sim = 1000)

Arguments

tab

A matrix of counts.

n.sim

Number of samples of permuted tables to consider.

Details

This is like the function stats::chisq.test(), but calculates an approximate P-value rather than refering to asymptotics. This will be better for large, sparse tables.

Value

A single number: the P-value testing independence of rows and columns in the table.

See Also

stats::chisq.test(), stats::fisher.test(), fisher()

Examples

TeaTasting <- matrix(c(3,1,1,3),nrow=2)
chisq(TeaTasting,1000)