Hi guys,
I'm performing Multi-Trait Analysis of GWAS (mtag) (see the following link:
github.com/JonJala/mtag

After running the command:

python ./mtag.py  --sumstats Sumstats/PGC_UKB_depression_genome-wide_mtag.txt,Sumstats/CRP_UKBB.imputed_v3.both_sexes_mtag.txt --out ./MDD_vs_CRP --n_min 0.0 --snp_name SNP --use_beta_se --beta_name BETA --se_name SE --force --maf_min 0.01 --stream_stdout --fdr  --stream_stdout & 

I found the following error message:

2021/06/03/10:48:00 AM attempt to get argmax of an empty sequence
Traceback (most recent call last):
  File "./mtag.py", line 1567, in <module>
    mtag(args)
  File "./mtag.py", line 1455, in mtag
    fdr(args, Ns, Zs)
  File "./mtag.py", line 1269, in fdr
    ind_max = np.argmax(FDR, axis=0)
  File "/home/platone1/.local/lib/python2.7/site-packages/numpy/core/fromnumeric.py", line 1103, in argmax
    return _wrapfunc(a, 'argmax', axis=axis, out=out)
  File "/home/platone1/.local/lib/python2.7/site-packages/numpy/core/fromnumeric.py", line 56, in _wrapfunc
    return getattr(obj, method)(*args, **kwds)
ValueError: attempt to get argmax of an empty sequence
2021/06/03/10:48:00 AM Analysis terminated from error at Thu Jun  3 10:48:00 2021
2021/06/03/10:48:00 AM Total time elapsed: 10.0m:9.81s

How can I overcome this error? May it be linked to python ?



Source link