| wp.anova.binary {WebPower} | R Documentation |
Statistical Power Analysis for One-way ANOVA with Binary Data
Description
The power analysis procedure for one-way ANOVA with binary data is introduced by Mai and Zhang (2017). One-way ANOVA with binary data is used for comparing means of three or more groups of binary data. Its outcome variable is supposed to follow Bernoulli distribution. And its overall test uses a likelihood ratio test statistics.
Usage
wp.anova.binary(k = NULL, n = NULL, V = NULL, alpha = 0.05,
power = NULL)
Arguments
k |
Number of groups. |
n |
Sample size. |
V |
Effect size. See the research by Mai and Zhang (2017) for details. |
alpha |
Significance level chosed for the test. It equals 0.05 by default. |
power |
Statistical power. |
Value
An object of the power analysis.
References
Mai, Y., & Zhang, Z. (2017). Statistical Power Analysis for Comparing Means with Binary or Count Data Based on Analogous ANOVA. In L. A. van der Ark, M. Wiberg, S. A. Culpepper, J. A. Douglas, & W.-C. Wang (Eds.), Quantitative Psychology - The 81st Annual Meeting of the Psychometric Society, Asheville, North Carolina, 2016: Springer.
Zhang, Z., & Yuan, K.-H. (2018). Practical Statistical Power Analysis Using Webpower and R (Eds). Granger, IN: ISDSA Press.
Examples
#To calculate the statistical power for one-way ANOVA (overall test) with binary data:
wp.anova.binary(k=4,n=100,V=0.15,alpha=0.05)
# One-way Analogous ANOVA with Binary Data
#
# k n V alpha power
# 4 100 0.15 0.05 0.5723443
#
# NOTE: n is the total sample size
# URL: http://psychstat.org/anovabinary
#To generate a power curve given a sequence of sample sizes:
res <- wp.anova.binary(k=4,n=seq(100,200,10),V=0.15,alpha=0.05,power=NULL)
res
# One-way Analogous ANOVA with Binary Data
#
# k n V alpha power
# 4 100 0.15 0.05 0.5723443
# 4 110 0.15 0.05 0.6179014
# 4 120 0.15 0.05 0.6601594
# 4 130 0.15 0.05 0.6990429
# 4 140 0.15 0.05 0.7345606
# 4 150 0.15 0.05 0.7667880
# 4 160 0.15 0.05 0.7958511
# 4 170 0.15 0.05 0.8219126
# 4 180 0.15 0.05 0.8451603
# 4 190 0.15 0.05 0.8657970
# 4 200 0.15 0.05 0.8840327
#
# NOTE: n is the total sample size
# URL: http://psychstat.org/anovabinary
#To plot the power curve:
plot(res)
#To calculate the required sample size for one-way ANOVA (overall test) with binary data:
wp.anova.binary(k=4,n=NULL,V=0.15,power=0.8, alpha=0.05)
# One-way Analogous ANOVA with Binary Data
#
# k n V alpha power
# 4 161.5195 0.15 0.05 0.8
#
# NOTE: n is the total sample size
# URL: http://psychstat.org/anovabinary
#To calculate the minimum detectable effect size for one-way ANOVA (overall test) with binary data:
wp.anova.binary(k=4,n=100,V=NULL,power=0.8, alpha=0.05)
# One-way Analogous ANOVA with Binary Data
#
# k n V alpha power
# 4 100 0.1906373 0.05 0.8
#
# NOTE: n is the total sample size
# URL: http://psychstat.org/anovabinary
#To generate a power curve given a sequence of effect sizes:
wp.anova.binary(k=4,n=100,V=seq(0.1,0.5,0.05),alpha=0.05,power=NULL)
# One-way Analogous ANOVA with Binary Data
#
# k n V alpha power
# 4 100 0.10 0.05 0.2746396
# 4 100 0.15 0.05 0.5723443
# 4 100 0.20 0.05 0.8402271
# 4 100 0.25 0.05 0.9659434
# 4 100 0.30 0.05 0.9961203
# 4 100 0.35 0.05 0.9997729
# 4 100 0.40 0.05 0.9999933
# 4 100 0.45 0.05 0.9999999
# 4 100 0.50 0.05 1.0000000
#
# NOTE: n is the total sample size
# URL: http://psychstat.org/anovabinary