-
Notifications
You must be signed in to change notification settings - Fork 62
Open
Description
Dear professor,
I came across the problem when runing mtag, and the err log file was shown as below.
I will bed truly appreciated fot your help.
Calling ./mtag.py \
--stream-stdout \
--n-min 0.0 \
--sumstats 1_OA2016_hm3samp_NEUR.txt,1_OA2016_hm3samp_SWB.txt \
--chunksize 10000 \
--out ./tutorial_results_1.1NS
Beginning MTAG analysis...
MTAG will use the Z column for analyses.
Read in Trait 1 summary statistics (185458 SNPs) from 1_OA2016_hm3samp_NEUR.txt ...
<><><<>><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Munging Trait 1 <><><<>><><><><><><><><><><><><><><><><><><><><><><><><><><
<><><<>><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Interpreting column names as follows:
snpid: Variant ID (e.g., rs number)
n: Sample size
a1: a1, interpreted as ref allele for signed sumstat.
pval: p-Value
a2: a2, interpreted as non-ref allele for signed sumstat.
z: Directional summary statistic as specified by --signed-sumstats.
Reading sumstats from provided DataFrame into memory 10000 SNPs at a time.
Read 185458 SNPs from --sumstats file.
Removed 0 SNPs with missing values.
Removed 0 SNPs with INFO <= None.
Removed 0 SNPs with MAF <= 0.01.
Removed 0 SNPs with SE <0 or NaN values.
Removed 0 SNPs with out-of-bounds p-values.
Removed 0 variants that were not SNPs. Note: strand ambiguous SNPs were not dropped.
185458 SNPs remain.
Removed 0 SNPs with duplicated rs numbers (185458 SNPs remain).
Removed 0 SNPs with N < 0.0 (185458 SNPs remain).
Segmentation fault (core dumped)
head -n 5 1_OA2016_hm3samp_NEUR.txt
snpid chr bpos a1 a2 freq z pval n
rs2736372 8 11106041 T C 0.4179 -7.71614161262 1.209e-14 111111.111111
rs2060465 8 11162609 T C 0.6194 7.69444599845 1.422e-14 62500.0
rs10096421 8 10831868 T G 0.4646 -7.561098219 3.989e-14 111111.111111
rs2409722 8 11039816 T G 0.4627 -7.38261601808 1.553e-13 111111.111111
head -n 5 1_OA2016_hm3samp_SWB.txt
snpid chr bpos a1 a2 freq z pval n
rs6648 20 47895242 A G 0.7761 5.56785113893 2.579e-08 62500.0
rs2295580 20 47626847 T C 0.6474 5.41782031863 6.033e-08 111111.111111
rs258681 7 81734211 A G 0.6213 -5.31715893243 1.054e-07 111111.111111
rs7294514 12 80821446 T C 0.09701 5.17914932223 2.229e-07 27777.7777778
Metadata
Metadata
Assignees
Labels
No labels