wanova_pwr.test {WAnova}R Documentation

wanova_pwr.test

Description

This function computes the approximate sample size required to achieve a desired power level for Welch's F-test in a one-way heteroscedastic ANOVA. The function takes the initial sample size, group means, and standard deviations, and iteratively determines the necessary sample size to meet the specified power and significance level.

Usage

wanova_pwr.test(n, means, sd, power = 0.9, alpha = 0.05)

Arguments

n

A numeric vector of initial sample sizes for each group.

means

A numeric vector of the means for each group.

sd

A numeric vector of the standard deviations for each group.

power

Desired power of the test. Default is 0.90.

alpha

Significance level for the test. Default is 0.05.

Details

The function adjusts the sample sizes iteratively until the desired power level is achieved. It uses the F-distribution to determine critical values and compute the power of the test. The output includes the sample size required for each group to achieve the specified power.

Reference: Levy, K. J. (1978a). Some empirical power results associated with Welch's robust analysis of variance technique. Journal of Statistical Computation and Simulation, 8, 43-48.

Value

A list containing:

n

The adjusted sample sizes for each group. If all adjusted sizes are the same, it returns a single value.

alpha

The significance level used in the computation.

power

The calculated power of the test, rounded to two decimal places.

See Also

print.wsize

Examples


# Example usage of the wsize function
n <- c(10, 10, 10, 10)
means <- c(1, 0, 0, -1)
sd <- c(1, 1, 1, 1)
result <- wanova_pwr.test(n, means, sd, power = 0.90, alpha = 0.05)
print(result)



[Package WAnova version 0.4.0 Index]