NSM3 (version 1.1)

cSkilMack: Computes a critical value for the Skillings-Mack SM distribution.

Description

This function computes the critical value for the Skillings-Mack SM distribution at (or typically in the "Exact" and "Monte Carlo" cases, close to) the given alpha level.

Usage

cSkilMack(alpha, obs.mat, method = NA, n.mc = 10000)

Value

Returns a list with "NSM3Ch7c" class containing the following components:knumber of treatmentsnnumber of blocksssnumber of treatments per blockcutoff.Uupper tail cutoff at or below user-specified alphatrue.alpha.Utrue alpha level corresponding to cutoff.U (if method="Exact" or "Monte Carlo")

Details

The incidence matrix, obs.mat, will be an n x k matrix of ones and zeroes, which indicate where the data are observed and unobserved, respectively.

Examples

##Hollander, Wolfe, Chicken Example 7.8 Effect of Rhythmicity of a Metronome on Speech Fluency
obs.mat<-matrix(c(rep(1,10),0,rep(1,13)),ncol=3,byrow=TRUE)
#cSkilMack(.01,obs.mat)
cSkilMack(.01,obs.mat,n.mc=5000)