From Genome Analysis Wiki
Jump to navigationJump to search
13 bytes added
, 21:59, 14 June 2015
Line 24: |
Line 24: |
| | | |
| rareMETALS2.single <- function(score.stat.file,range,alternative=c('two.sided','greater','less'),ix.gold=1,callrate.cutoff=0,hwe.cutoff=0,hwe.ctrl.cutoff=0) | | rareMETALS2.single <- function(score.stat.file,range,alternative=c('two.sided','greater','less'),ix.gold=1,callrate.cutoff=0,hwe.cutoff=0,hwe.ctrl.cutoff=0) |
− | * ''score.stat.file'' files of score statistics
| + | '''Relevant Parameters''': |
− | * ''range tabix'' range of variants to be analyzed
| + | * ''score.stat.file'' files of score statistics |
− | * ''alternative'' alternative hypothesis to be specified
| + | * ''range tabix'' range of variants to be analyzed |
− | * ''ix.gold'' Gold standard population to align reference allele to.
| + | * ''alternative'' alternative hypothesis to be specified |
− | * ''callrate.cutoff'' Cutoffs of call rate, lower than which will NOT be analyzed (labelled as missing)
| + | * ''ix.gold'' Gold standard population to align reference allele to. |
− | * ''hwe.cutoff'' Cutoffs of HWE p-values; Variants with HWE p-value smaller than the cutoffs are removed from subsequent analysis and labelled as missing;
| + | * ''callrate.cutoff'' Cutoffs of call rate, lower than which will NOT be analyzed (labelled as missing) |
− | * ''hwe.ctrl.cutoff'' Cutoffs of HWE p-values using controls; Variants with HWE p-value smaller than the cutoffs are removed from subsequent analysis and labelled as missing; In case control studies, it is recommended to use hwe.ctrl.cutoff, since large effect variants may violate HWE.
| + | * ''hwe.cutoff'' Cutoffs of HWE p-values; Variants with HWE p-value smaller than the cutoffs are removed from subsequent analysis and labelled as missing; |
| + | * ''hwe.ctrl.cutoff'' Cutoffs of HWE p-values using controls; Variants with HWE p-value smaller than the cutoffs are removed from subsequent analysis and labelled as missing; In case control studies, it is recommended to use hwe.ctrl.cutoff, since large effect variants may violate HWE. |