View source: R/conditional.rareMETALS.range.R
Perform conditional analysis for gene-level tests
1 2 3 4 | conditional.rareMETALS.range(range.name = NULL, score.stat.file, cov.file,
candidate.variant.vec, known.variant.vec, test = "GRANVIL", maf.cutoff,
alternative = c("two.sided", "greater", "less"), ix.gold = 1,
out.digits = 4, callrate.cutoff = 0, hwe.cutoff = 0, max.VT = NULL)
|
range.name |
name of the range to be analyzed (for example, it can be a gene name e.g. APOE) |
score.stat.file |
files of score statistics |
cov.file |
covariance matrix files |
candidate.variant.vec |
Vectors of candidate variants: e.g. c("1:123","1:1234"). The quotation is necessary!! |
known.variant.vec |
Vectors of known variants: e.g. c("1:123","1:1234"). The quotation is necessary!! |
test |
test of rare variant tests |
maf.cutoff |
Cutoffs of MAF used for determining rare variants |
alternative |
Alternative hypothesis to be tested |
ix.gold |
Index of the gold standard population: Used for flipping alleles |
out.digits |
The number of digits used in the output |
callrate.cutoff |
Cutoff of call rates. Sites with callrates lower than the cutoff will be labeled as missing |
hwe.cutoff |
Cutoff of HWE p-values. Sites with HWE pvalues lower than the cutoff will be labeled as missing |
max.VT |
The maximum number of thresholds used in VT; Setting max.VT to 10 can improve the speed for calculation without affecting the power too much. The default parameter is NULL, which does not set upper limit on the number of variable frequency threhsold. |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.