# \donttest{
# Load and adjust data.
data(epiDat)
simResults <- prep_data(data = epiDat, type = "counts", Sest.method = "average",
cases = 5, N = 100, M = 10,
n = 5, m = 6, k = 20,
transformation = "none", method = "bray",
dummy = TRUE, useParallel = FALSE,
model = "single.factor")
# }
simResults
# Computing components of variation
compVar <- scompvar(data = simResults)
compVar
# Determination of statistical power
epiBetaR <- sim_beta(simResults, alpha = 0.05)
epiBetaR
# Cost-benefit optimization
cboResult <- sim_cbo(epiBetaR, cn = 75)
cboResult
# Visualization of statistical power
plot_power(data = epiBetaR, method = "power")
Run the code above in your browser using DataLab