skills/statistical-power/scripts/simulate_power.py

• PowerEstimate · class · L27-L36 — class PowerEstimate
• str · method · L34-L36 — def str(self)
• _wilson_ci · function · L39-L47 — def _wilson_ci(k, n, z=1.96)
• simulate_power · function · L50-L69 — def simulate_power(gen_and_test, n, n_sims=2000, alpha=0.05, seed=0)
• find_sample_size · function · L72-L97 — def find_sample_size(gen_and_test, target_power=0.80, n_sims=2000, alpha=0.05, lo=10, hi=2000, seed=0, verbose=True)
• example_two_group_difference · function · L105-L115 — def example_two_group_difference(effect=0.5, sd=1.0, alpha=0.05)
• gen_and_test · function · L109-L113 — def gen_and_test(n, rng): # n per group
• example_logistic_regression · function · L118-L140 — def example_logistic_regression(beta=0.8, base_rate=0.2, x_sd=1.0, alpha=0.05)
• gen_and_test · function · L128-L138 — def gen_and_test(n, rng)
• example_cluster_randomized · function · L143-L174 — def example_cluster_randomized(effect=0.3, icc=0.05, cluster_size=20, resid_sd=1.0, alpha=0.05)
• gen_and_test · function · L157-L172 — def gen_and_test(n, rng): # n clusters per arm
• example_linear_mixed_repeated · function · L177-L201 — def example_linear_mixed_repeated(effect=0.4, n_timepoints=3, subj_sd=0.7, resid_sd=1.0, alpha=0.05)
• gen_and_test · function · L187-L199 — def gen_and_test(n, rng): # n subjects