logo
Free, unlimited AI code reviews that run on commit
git-lrc git-lrc GitHub Install Now We'd appreciate a star git-lrc - Free, unlimited AI code reviews that run on commit | Product Hunt git-lrc - Free, unlimited AI code reviews that run on commit | Product Hunt

QTLtools gwas - Genome-wide association tests

Authors

       Olivier Delaneau (olivier.delaneau@gmail.com), Halit Ongen (halitongen@gmail.com)

QTLtools-v1.3                                      06 May 2020                                  QTLtools-gwas(1)

Bugs

       Versions up to and including 1.2, suffer from a bug in reading missing genotypes in VCF/BCF files.   This
       bug  affects variants with a DS field in their genotype's FORMAT and have a missing genotype (DS field is
       .) in one of the samples, in which case genotypes for all the samples are  set  to  missing,  effectively
       removing this variant from the analyses.

       Please submit bugs to <https://github.com/qtltools/qtltools>

Citation

       Delaneau  O., Ongen H., Brown A. A., et al. A complete tool set for molecular QTL discovery and analysis.
       NatCommun8, 15452 (2017).  <https://doi.org/10.1038/ncomms15452>

Description

       This mode runs a linear regression with every genotype vs. every phenotype.

Example

       o Run the association between all the variants and phenotypes, correcting the phenotypes  with  technical
         covariates and normal transforming them:

         QTLtools  gwas  --bed  quantifications.bed.gz  --vcf  genotypes.bcf --cov covariates.txt --normal --out
         gwas_results.txt

Name

       QTLtools gwas - Genome-wide association tests

Options

--vcf[in.vcf|in.bcf|in.vcf.gz|in.bed.gz]...
              Genotypes in VCF/BCF/BED format.  The DS field in the genotype FORMAT of a variant (dosage of  the
              genotype calculated from genotype probabilities, e.g. after imputation) is REQUIRED.  REQUIRED.

       --bedquantifications.bed.gz ...
              Molecular phenotype quantifications in BED format.  REQUIRED.

       --outfilename.txt
              The output file name.  REQUIRED.

       --covcovariates.txt
              Covariates to correct the phenotype data with.

       --normal
              Rank  normal  transform  the  phenotype  data  so  that  each  phenotype  is normally distributed.
              RECOMMENDED.

Output File

outputfile
        Space separated output file with the following columns.
        1     The phenotype ID
        2     The variant's chromosome
        3     The variants's start position
        3.1   The variant's end position.  Only printed if a BED file is supplied to --vcf
        4     The variant's ID
        5     The p-value
        6     The correlation coefficient

See Also

QTLtools(1)

       QTLtools website: <https://qtltools.github.io/qtltools>

Synopsis

QTLtoolsgwas--vcf[in.vcf|in.vcf.gz|in.bcf|in.bed.gz]--bedquantifications.bed.gz--outoutput.txt[OPTIONS]

See Also