gwas_power: Calculate power for GWAS of a binary trait

View source: R/functions.R

gwas_powerR Documentation

Calculate power for GWAS of a binary trait

Description

Calculate power for GWAS of a binary trait

Usage

gwas_power(
  eaf = 0.5,
  beta_alt = 0.5,
  n_case = 1000,
  n_control = 1000,
  alpha = 5e-08,
  n_perm = 1000
)

Arguments

eaf

effect allele frequency

beta_alt

beta under the alternate hypothesis (i.e. per-allele log odds ratio in your discovery GWAS)

n_case

number of cases

n_control

number of controls

alpha

desired alpha

n_perm

number of permutations for power calculations

Value

the estimated power for each SNP given the allele frequency, n, alpha, effect size, and number of permutations. Note that this is for a one-tailed test.


benjacobs123456/popPoweR documentation built on March 28, 2022, 4:32 a.m.