Description Usage Arguments Value
View source: R/rareMETALS.estimatePVE.R
Estimate the variance explained for a list of variants (possibly in LD);
1 2 3 | rareMETALS.estimatePVE(score.stat.file, cov.file, range, range.name, refaltList,
maf.cutoff, out.digits, callrate.cutoff, hwe.cutoff, correctFlip = TRUE,
lambda = 0)
|
score.stat.file |
files of score statistics |
cov.file |
Covariance matrix files |
range |
tabix range for each gene/region; Tabix range needs to be in the format of "1:123-1234". Quotation marks are necessary. |
range.name |
The name of the range,e.g. gene names can be used; |
refaltList |
A list of reference, alternative allele, a vector of alternative allele frequencies; Specifically, the list should consist of |
maf.cutoff |
MAF cutoff used to analyze variants; Default value is 1, i.e. no cutoffs are applied. MAFs are based upon the sample MAFs. |
out.digits |
Number of digits used in the output |
callrate.cutoff |
Cutoffs of call rate, lower than which will NOT be analyzed (labelled as missing) |
hwe.cutoff |
Cutoffs of HWE p-values |
correctFlip |
Correct for flipped sites for score statistics and their covariance matrices |
lambda |
Tuning parameter to the variant the variance correlation matrix to avoid near singular LD matrix; Assume the orignal LD matrix is R2, R2.adj=R2+lambda*Id; |
a list consisting of results;
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.