# \donttest{
# Load and adjust data.
data(epiDat)
simResults <- prep_data(data = epiDat, type = "counts", Sest.method = "average",
cases = 5, N = 100, sites = 10,
n = 5, m = 5, k = 30,
transformation = "none", method = "bray",
dummy = FALSE, useParallel = FALSE,
model = "single.factor")
# }
simResults
# Computing components of variation
compVar <- scompvar(data = simResults)
compVar
# Cost-benefit optimization
cboResult <- sim_cbo(comp.var = compVar, ct = 20000, ck = 100, cj = 2500)
cboResult
# Determination of statistical power
epiBetaR <- sim_beta(simResults, alpha = 0.05)
epiBetaR
# Visualization of statistical power
plot_power(data = epiBetaR, n = NULL, m = 3, method = "both")
Run the code above in your browser using DataLab