Closed:Extracting SNP status using QCTOOL from an impute2 output
1
0
Entering edit mode
6.4 years ago
Mardze • 0

Hi everyone!

I am doing a GWAS with imputed genotypes (phased with Eagle, imputed with Impute2), that resulted in a format similar to this:

--- rs17149628:126006965:C:T 126006965 C T

After imputation, an analysis using SNPTEST went absolutely fine, without producing errors.

However, I am now trying to extract certain SNP status for all individuals using QCTOOL and it seems there's an issue when I try to specify the rsID (second field, i.e. rs17149628:126006965:C:T) as no output is produced (however no errors are listed s). My guess is QCTOOL reads ":" as a field separator and isn't able to treat it as a SNP. Would anyone know how to overcome the problem? The whole log is as follows:

 qctool -g data.impute2 -s sample.txt -incl-rsids snp.txt -og extract.txt
     (C) 2009-2011 University of Oxford
    Opening genotype files                       :
    (1/1,0.0s,33.6/s)
    ===================================================================
    Input SAMPLE file(s):           "sample.txt" Output
    SAMPLE file:             "(n/a)". Sample statistic output file:  
    "(n/a)". Sample exclusion output file:   "(n/a)".
    Input GEN file(s):
   (not computed)  "snp-id-data-filtered:data.impute2"  (total 1 sources, number of snps not computed).
     Number of samples: 653 Output GEN file(s):  "extract.txt" Output SNP position file(s):    (n/a) SNP statistic
    output file(s):  Sample filter:                  (none). SNP filter:
    (none).

    # of samples in input files:    653.
    # of samples after filtering:   653 (0 filtered out).
    ===============================================================
    Processing SNPs                              :  (1/?,21.2s,0.0/s)
    ===============================================================

    Number of SNPs:
                         -- in input file(s):                 (not computed).  -- in output file(s):                1

    Number of samples in input file(s):   653.

    Output GEN files:                     (1      snps) 
    "extract.txt"
                                          (total 1 snps).
    ========================================================================
    Thank you for using qctool.

So, all seems fine, but the extract.txt is empty. I am becoming desperate here, would anyone know what/why is happening and how to overcome it? I'll be thankful for any comments! Nad

SNP GWAS Impute2 SNPTEST QCTOOL • 1.8k views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 1488 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6