plotCcBurdenByRelativeRisk: Tests the power to detect a burden difference using case...

In stephansanders/wgsPowerTest: Whole Genome Sequencing Power Calculations

Description

Tests the power to detect a burden difference using case control data varying relative risk.

Usage

plotCcBurdenByRelativeRisk(R, N, r = 1, q, f, p_thres_cc_burden = 0.05,
  name = "Test", col = "#e41a1c")

Arguments

R	Maximum relative risk (numeric).
N	Sample size (integer).
r	Case:Control ratio (numeric).
q	Number of functional rare variants per person in seleceted regions (integer).
f	Proportion of functional rare variants that mediate risk for the disorder (numeric, 0-1).
p_thres_cc_burden	p-value threshold after correction for multiple comparisons (numeric, 0-1).
name	Prefix for plot filenames (text).
col	Color of the line (hex).

Value

A PDF showing the estimated statistical power as the relative risk is varied.

Examples

R <- 2.5
N <- 20000
r <- 1
q <- 126
f <- 0.2
p_thres_cc_burden <- 0.05 / 1000
name <- "Test1"
col <- "#e41a1c"
plotCcBurdenByRelativeRisk()

stephansanders/wgsPowerTest documentation built on May 4, 2019, 1:23 p.m.