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 rep - Replicate QTL associations in an independent dataset

Authors

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

QTLtools-v1.3                                      06 May 2020                                   QTLtools-rep(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 reads phenotype-genotype pairs that were identified in another dataset and  checks  association
       between the same genotype and phenotype in this dataset.

Example

       o Replicate the  results  of  dataset1  in  dataset2,  correcting  dataset2  phenotype's  with  technical
         covariates and normal transforming them:

         QTLtools   rep   --bed   dataset2.bed.gz   --vcf   dataset2.bcf   --cov  dataset2.covariates.txt  --qtl
         significant_qtls_dataset1.txt --normal --out rep_results.txt

Name

       QTLtools rep - Replicate QTL associations in an independent dataset

Options

--vcf[in.vcf|in.bcf|in.vcf.gz|in.bed.gz]...
              Genotypes in VCF/BCF format, or another molecular phenotype in BED format.  If there is a DS field
              in   the  genotype  FORMAT  of  a  variant  (dosage  of  the  genotype  calculated  from  genotype
              probabilities, e.g. after imputation), then this is used as the genotype.  If there is only the GT
              field in the genotype FORMAT then this is used and it is converted to a dosage.  REQUIRED.

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

       --qtlexternal_qtls.txt ...
              The file listing the QTLs identified  in  the  external  dataset.   One  phenotype  genotype  pair
              (separated by space, phenotype first) per line.  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 phenotype's chromosome
         3   The phenotype's start position
         4   The phenotype's end position
         5   The phenotype's strand
         6   The genotype ID
         7   The genotype's chromosome
         8   The genotype's start position
         9   The genotype's end position
        10   The p-value of the association
        11   The slope of the association

See Also

QTLtools(1)

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

Synopsis

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

See Also