conditional.rareMETALS.range: Perform conditional analysis for gene-level tests

Description Usage Arguments

View source: R/conditional.rareMETALS.range.R

Description

Perform conditional analysis for gene-level tests

Usage

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)

Arguments

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.


dajiangliu/rareGWAMA documentation built on Sept. 13, 2019, 9:14 a.m.