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

Description Usage Arguments Value Examples

View source: R/wgsPowerCalc_helper.R

Description

Tests the power to detect a burden difference using case control data varying sample size.

Usage

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

Arguments

R

Relative risk (numeric).

N

Maximum 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 sample size is varied.

Examples

1
2
3
4
5
6
7
8
9
R <- 2.5
N <- 20000
r <- 1
q <- 126
f <- 0.2
p_thres_cc_burden <- 0.05 / 1000
name <- "Test1"
col <- "#e41a1c"
plotCcBurdenBySampleSize()

stephansanders/wgsPowerTest documentation built on Sept. 10, 2017, 12:25 a.m.