-
Notifications
You must be signed in to change notification settings - Fork 101
Description
Dear hakyimlab,
I'm really appreciate the software of TWAS you created. However, I meet some troubles.
My code:
./SPrediXcan.py
--model_db_path 'data/dendritic_cell.db'
--covariance 'data/predict_db_Model_training_filtered.txt.gz'
--gwas_folder data/GWAS
--gwas_file_pattern ".*txt"
--snp_column VARID
--effect_allele_column A1
--non_effect_allele_column A2
--pvalue_column P
--or_column OR
--model_db_snp_key rsid
--keep_non_rsid
--output_file 'results/DC_result.csv'
My GWAS data:
SNP CHR BP N P OR MAF A1 A2 VARID
rs1000000 12 126890980 8 0.07901 0.9262 0.223658 G A chr12_126890980_A_G_b38
rs10000003 4 57561647 10 0.2496 1.045 0.297217 A G chr4_57561647_G_A_b38
rs10000005 4 85161558 8 0.002903 0.9166 NA NA NA chr4_85161558_NA_NA_b38
rs10000010 4 21618674 5 0.1527 0.9463 0.49006 T C chr4_21618674_C_T_b38
rs10000010 4 21618674 5 0.1527 0.9463 0.49006 T C chr4_21618674_C_T_b38
No matter I use SNP or VARID as snp_column, there is no results.
output:
WARNING - Missing --gwas_h2 and --gwas_N are required to calibrate the pvalue and zscore.
INFO - Processing GWAS command line parameters
INFO - Building beta for gwas_with_varid.txt and data/dendritic_cell.db
INFO - Reading input gwas with special handling: data/GWAS/gwas_with_varid.txt
INFO - Processing input gwas
INFO - Aligning GWAS to models
INFO - Trimming output
INFO - Successfully parsed input gwas in 9.162299307994545 seconds
INFO - Started metaxcan process
INFO - Loading model from: data/dendritic_cell.db
INFO - Loading covariance data from: data/predict_db_Model_training_filtered.txt.gz
INFO - Processing loaded gwas
INFO - Started metaxcan association
INFO - 0 % of model's snps used
I'm looking forward to your reply. Best wishes!@hakyimlab