plotDnLocusByRelativeRisk: Tests the power to detect a gene locus using de novo data...

In stephansanders/wgsPowerTest: Whole Genome Sequencing Power Calculations

Description

Tests the power to detect a gene locus using de novo data varying relative risk.

Usage

plotDnLocusByRelativeRisk(R, N, r = 1, q, f, p_thres_denovo_locus = 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_denovo_locus	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 <- 25
N <- 5000
r <- 1
q <- 0.0249
f <- 0.2
p_thres_denovo_locus <- 0.05 / 20000
name <- "Test1"
col <- "#e41a1c"
plotDnLocusByRelativeRisk()

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