Methods
run(u, v, w, mafs) → {Array.<Number>}
This code corresponds roughly to: https://github.com/statgen/raremetal/blob/2c82cfc5710dbd9fd56ef67a7ca5f74772d4e70d/raremetal/src/Meta.cpp#L3456
Parameters:
| Name | Type | Description |
|---|---|---|
u |
Array.<Number> | Vector of score statistics (length m, number of variants) |
v |
Array.<Number> | Covariance matrix of score statistics |
w |
This parameter is ignored for VT. Weights are calculated automatically from mafs. |
|
mafs |
Minor allele frequency of each variant |
Returns:
Array of: [largest test statistic over all burden tests across allele frequency thresholds, VT analytical p-value, effect size of largest test statistic, standard error of effect size] Note that this method returns the largest test statistic and corresponding effect/se over all burden tests. Remember that these do not correspond to the VT p-value, because VT uses the MVN and covariance of all test statistics to compute the p-value.
- Type
- Array.<Number>