mvnormalTest (version 1.0.0)

power.mswV: Power Calculation using the SWT-based Villasenor-Alva and Gonzalez-Estrada (VAGE) Test Statistic

Description

Empirical power calculation using VAGE test statistic.

Usage

power.mswV(a, n, p, B = 1000, FUN, ...)

Arguments

a

significance level ($\alpha$).

n

number of rows (observations).

p

number of columns (variables), $n>p$.

B

number of Monte Carlo simulations, default is 1000 (can increase B to increase the precision).

FUN

self-defined function for generate multivariate distribution. See example.

...

optional arguments passed to FUN.

Value

Returns a numeric value of the estimated empirical power (value between 0 and 1).

References

Villasenor Alva, J. A., & Estrada, E. G. (2009). A generalization of Shapiro<U+2013>Wilk's test for multivariate normality. Communications in Statistics<U+2014>Theory and Methods, 38(11), 1870-1883.

Examples

# NOT RUN {
set.seed(12345)

## Power calculation against bivariate (p=2) independent Beta(1, 1) distribution ##
## at sample size n=50 at one-sided alpha = 0.05 ##

power.mswV(a = 0.05, n = 50, p = 2,  B = 100, FUN=IMMV, D1=runif)

# }