mv.2way.test: Randomized Complete Block Design. In MNM: Multivariate Nonparametric Methods. An Approach Based on Spatial Signs and Ranks

Description

Multivariate tests for testing the null hypothesis that there is no treatment effect in a randomized complete block design using different scores.

Usage

 ```1 2 3 4 5``` ```mv.2way.test(x, block, treatment, score = c("identity", "sign", "rank"), stand = c("outer", "inner"), method = c("approximation", "permutation"), n.simu = 1000, eps=1.0e-10, n.iter=10000, na.action = na.fail) ```

Arguments

 `x` a numeric data frame or matrix of response variables. `block` a factor with at least two levels. `treatment` a factor with at least two levels. `score` the score to be used. Possible choices are `identity`, `sign` and `rank`. `stand` the standardization method used. Possible choices are `outer` and `inner`. `method` method for the computation of the p-value for the spatial sign and spatial rank tests. Possible choices are `approximation` and `permutation`. `n.simu` number of simulated permutations if `method="permutation"`. `eps` convergence criterion. `n.iter` maximum number of iterations. `na.action` a function which indicates what should happen when the data contain 'NA's. Default is to fail.

Details

This implements the tests described in chapter 12 of the MNM book.

Value

A list with class 'htest' containing the following components:

 `statistic ` the value of the test statistic. `parameter` the degrees of freedom for the test statistic or the number of replications in the simulation. `p.value` the p-value for the test. `null.value` the specified null hypothesis value of the location. `alternative` a character string with the value 'two.sided'. `method` a character string indicating what type of test was performed. `data.name` a character string giving the name of the data set and of the grouping vector.

Author(s)

Jyrki M<f6>tt<f6>nen [email protected]

References

Oja, H. (2010), Multivariate Nonparametric Methods with R, Springer.

`mv.1sample.test`, `mv.Csample.test` , `mv.2way.est`
 ```1 2 3 4``` ```blocks <- gl(10, 5) treatments <- factor(rep(1:5, 10)) X <- rmvnorm(n = 50, mean = c(1,2,3), sigma = diag(3)) mv.2way.test(X, blocks, treatments, score="r", stand="i", method="a") ```