# Superpower package
library(Superpower)
## Warning: package 'Superpower' was built under R version 4.0.3
## Registered S3 methods overwritten by 'lme4':
## method from
## cooks.distance.influence.merMod car
## influence.merMod car
## dfbeta.influence.merMod car
## dfbetas.influence.merMod car
n <- 240
sd <- .8
r <- 0.84
alpha = 0.05
f <- 0.125
mu <- c(-.125,.125,-.125,.125,-.125,(.125+.125))
design = "3b*2w"
labelnames <- c("condition","control","equivalent","intervention","time", "pre", "post")
design_result <- ANOVA_design(design = design,
n = n, mu = mu, sd = sd, r = r,
labelnames = labelnames, plot = T)
design_result
## The design is specified as: 3b*2w
## Summary of means (mu) and standard deviations (SD)
## mu SD condition time
## 1 -0.125 0.8 control pre
## 2 0.125 0.8 control post
## 3 -0.125 0.8 equivalent pre
## 4 0.125 0.8 equivalent post
## 5 -0.125 0.8 intervention pre
## 6 0.250 0.8 intervention post
##
## Correlation Matrix
## control_pre control_post equivalent_pre equivalent_post
## control_pre 1.00 0.84 0.00 0.00
## control_post 0.84 1.00 0.00 0.00
## equivalent_pre 0.00 0.00 1.00 0.84
## equivalent_post 0.00 0.00 0.84 1.00
## intervention_pre 0.00 0.00 0.00 0.00
## intervention_post 0.00 0.00 0.00 0.00
## intervention_pre intervention_post
## control_pre 0.00 0.00
## control_post 0.00 0.00
## equivalent_pre 0.00 0.00
## equivalent_post 0.00 0.00
## intervention_pre 1.00 0.84
## intervention_post 0.84 1.00
alpha_level <- 0.05
nsims <- 500
simulation_result <- ANOVA_power(design_result,
alpha_level = alpha_level,
nsims = nsims, verbose = F)
simulation_result
## Power and Effect sizes for ANOVA tests
## power effect_size
## anova_condition 14.4 0.004371
## anova_time 100.0 0.296223
## anova_condition:time 88.6 0.019133
##
## Power and Effect sizes for pairwise comparisons (t-tests)
## power
## p_condition_control_time_pre_condition_control_time_post 100.0
## p_condition_control_time_pre_condition_equivalent_time_pre 4.8
## p_condition_control_time_pre_condition_equivalent_time_post 93.4
## p_condition_control_time_pre_condition_intervention_time_pre 6.0
## p_condition_control_time_pre_condition_intervention_time_post 100.0
## p_condition_control_time_post_condition_equivalent_time_pre 93.2
## p_condition_control_time_post_condition_equivalent_time_post 5.8
## p_condition_control_time_post_condition_intervention_time_pre 92.4
## p_condition_control_time_post_condition_intervention_time_post 41.4
## p_condition_equivalent_time_pre_condition_equivalent_time_post 100.0
## p_condition_equivalent_time_pre_condition_intervention_time_pre 4.4
## p_condition_equivalent_time_pre_condition_intervention_time_post 100.0
## p_condition_equivalent_time_post_condition_intervention_time_pre 94.0
## p_condition_equivalent_time_post_condition_intervention_time_post 39.6
## p_condition_intervention_time_pre_condition_intervention_time_post 100.0
## effect_size
## p_condition_control_time_pre_condition_control_time_post 0.5579678
## p_condition_control_time_pre_condition_equivalent_time_pre 0.0030674
## p_condition_control_time_pre_condition_equivalent_time_post 0.3174970
## p_condition_control_time_pre_condition_intervention_time_pre 0.0024219
## p_condition_control_time_pre_condition_intervention_time_post 0.4731168
## p_condition_control_time_post_condition_equivalent_time_pre -0.3108922
## p_condition_control_time_post_condition_equivalent_time_post 0.0029781
## p_condition_control_time_post_condition_intervention_time_pre -0.3122625
## p_condition_control_time_post_condition_intervention_time_post 0.1579784
## p_condition_equivalent_time_pre_condition_equivalent_time_post 0.5584144
## p_condition_equivalent_time_pre_condition_intervention_time_pre -0.0005692
## p_condition_equivalent_time_pre_condition_intervention_time_post 0.4696139
## p_condition_equivalent_time_post_condition_intervention_time_pre -0.3155467
## p_condition_equivalent_time_post_condition_intervention_time_post 0.1550900
## p_condition_intervention_time_pre_condition_intervention_time_post 0.8306050