chisq(x, ...)## S3 method for class 'htest':
chisq(x, ...)
## S3 method for class 'table':
chisq(x, correct = FALSE, ...)
## S3 method for class 'default':
chisq(x, correct = FALSE, ...)
"htest"
a coming from a Chi-squared test,
an object of class "table"
, or
the inputs to tally
.tally
or chisq.test
.stat
if(require(mosaicData)) {
Mites.table <- tally( ~ outcome + treatment, data=Mites )
Mites.table
chisq.test(Mites.table)
chisq(Mites.table)
chisq(chisq.test(Mites.table))
## Randomization test. Increase replications to decrease Monte Carlo error.
do(3) * chisq( tally( ~ outcome + shuffle(treatment), data=Mites ) )
Mites.rand <- do(1000) * chisq( tally( ~ outcome + shuffle(treatment), data=Mites ) )
tally( ~(X.squared >= chisq(Mites.table)), data=Mites.rand, format="proportion")
}
Run the code above in your browser using DataLab