Line 140: |
Line 140: |
| | | |
| ==PhoneHome== | | ==PhoneHome== |
− | ===--noPhoneHome===
| |
− | ===--phoneHomeThinning===
| |
| | | |
− | | + | {{PhoneHomeParameters|hdr===|bullet=1}} |
− | | |
− | * --prefix is optional.
| |
− | * If --prefix is not specified, the output file names will be:
| |
− | traitname.singlevar.score.txt
| |
− | traitname.singlevar.cov.txt
| |
− | * Otherwise, the output file names are:
| |
− | prefix.traitname.singlevar.score.txt
| |
− | prefix.traitname.singlevar.cov.txt
| |
− | * --LDwindow specifies the length of the window that LD Matrix should be generated upon each variant. The default is 1MB.
| |
− | * --zip gives users the option of writing compressed files (bgzip compressed) automatically for convenient sharing.
| |
− | * --thin tells RMW to thin points when generating QQ plot and Manhattan plots, so the file size is smaller.
| |
− | * --labelHits tells RMW to to label the hits using pvalue threshold 0.05/(#of variants tested) with gene name, based on human genome build 19.
| |
− | | |
− | ==== VC Options ====
| |
− | * When --vcShared and --vcX are specified, RMW knows that you want to fit shared environment and/or chromosome X variance component together with genetic component and non-shared environment.
| |
− | * When --makeResiduals is specified, RMW understands covariates should be read from PED/DAT file. Covariates are modeled as fixed effects.
| |
− | | |
− | ==== Trait Options ====
| |
− | * --makeResiduals tells RMW to adjust the covariates and analyze residuals instead of the original phenotypes. If either --kinGeno or --kinPedigree option is used, then a variance component model will be fit based on residuals. If the --inverseNormal option is also used, then the residuals will be quantile normalized before fitting variance component model.
| |
− | * --traitName is created for situations when you have many traits saved in your PED and DAT file, but you are interested in one or a few of them. It can read a file ending with .txt with each trait of interest in a separate line, or trait names separated with "/". An example to handle one trait or multiple traits is in the following:
| |
− | --traitName LDL
| |
− | --traitName LDL/HDL/TG
| |
− | --traitName traitsOfInterest.txt
| |
− | * If --traitName is not used, all traits in PED/DAT file will be analyzed.
| |
− | | |
− | ==== Model Options ====
| |
− | * additive model is used in RMW as default.
| |
− | * --recessive allows additional association results (pvalue, effect size, and standard error) generated using recessive model. If VCF file is used, then non-reference allele is considered the recessive allele. If PED/DAT files are used for genotype, then minor allele is considered the recessive allele.
| |
− | * --dominant allows additional association results (pvalue, effect size, and standard error) generated using dominant model. If VCF file is used, then non-reference allele is considered the dominant allele. If PED/DAT files are used for genotype, then minor allele is considered the dominant allele.
| |
− | * --recessive and --dominant options can be used together.
| |
− | * Recessive and dominant results are stored in separate files.
| |
− | | |
− | ==== Kinship Source ====
| |
− | * --kinPedigree allows RMW to generate kinship matrix from pedigree, when pedigree information is available.
| |
− | * --kinGeno informs RMW to generate kinship matrix from all available variants that pass the criteria, specified in --kinMaf and --kinMiss options. The default will take variants with MAF>0.05 and genotype missing rate <0.05.
| |
− | * --kinGeno option can NOT be used with --kinPedigree or --kinFile option. Only one of three options or none of them can be used in the same run.
| |
− | * --kinFile let RMW read in a kinship matrix from a file. The first row of the kinship file has to be the sample IDs included in the kinship file. If a sample of interest is not included in the kinship file, fatal error will occur and the program will be terminated. A sample of interest is a sample that is phenotyped and has all covariates measured when --makeResiduals is specified.
| |
− | * --kinSave allows you to save the kinship matrix.
| |
− | | |
− | ==== Kinship Options ====
| |
− | * --kinMiss and --kinMaf should be used with --kinGeno together.
| |
− | * --kinMiss specifies the maximum genotype missing rate when calculating kinship from genotypes. The default is 0.05.
| |
− | * --kinMaf specifies the minimum minor allele frequency used when calculating kinship from genotypes. The default is 0.05.
| |
− | | |
− | ==== Chromosome X ====
| |
− | * --xLabel should have a value of a string which specifies how variants on chromosome X are coded. The default is "X".
| |
− | * --xStart and --xEnd specifies the start and end of non-pseudo-autosomal regions on chromosome X. These options should be specified when --vcX is used.
| |
− | * The default for --xStart is 2699520 and default for --xEnd is 154931044, according to NCBI genome build 37.
| |
− | | |
− | Please refer to the following for the analysis of X-linked variants [[RAREMETALWORKER_X|'''ANALYZING CHROMOSOME X''']].
| |
− | | |
− | {{PhoneHomeParameters|hdr=====|bullet=1}} | |