Description Usage Arguments Value Examples
Tests the power to detect a gene locus (i.e. locus defined by a gene rather than variant) using case control data varying relative risk.
1 2 | plotCcLocusMultiByRelativeRisk(R, N, r = 1, s, f_gene, f, K, AF_bar, N_rep,
p_thres_cc_locus_multi = 0.05, name = "Test", col = "#e41a1c")
|
R |
Maximum relative risk (numeric). |
N |
Sample size (integer). |
r |
Case:Control ratio (numeric). |
s |
Number of possible rare variants per gene (numeric). |
f_gene |
Number of possible functional rare variants per gene (numeric). |
f |
Proportion of functional rare variants that mediate risk for the disorder (numeric, 0-1). |
K |
Prevalence of the disorder (numeric, 0-1). |
AF_bar |
Average minor allele frequency (numeric, 0-1). |
N_rep |
Number of replicate simulations (integer). |
p_thres_cc_locus_multi |
p-value threshold after correction for multiple comparisons (numeric, 0-1). |
name |
Prefix for plot filenames (text). |
col |
Color of the line (hex). |
A PDF showing the estimated statistical power as the relative risk is varied.
1 2 3 4 5 6 7 8 9 10 11 12 13 | R <- 2.5
N <- 20000
r <- 1
s <- 615
f_gene <- 123
f <- 0.2
K <- 0.01
AF_bar <- 0.001
N_rep <- 50
p_thres_cc_locus_multi <- 0.05 / 20000
name <- "Test1"
col <- "#e41a1c"
plotCcLocusMultiByRelativeRisk()
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.