-
Notifications
You must be signed in to change notification settings - Fork 10
Description
I was in trouble in step 1, seeing below:
(HESS) [lijunpeng@login05 step2]$ DIR1="/public/home/lijunpeng/HESS/Cleandata/scz/covariance_correlation/step1/"
(HESS) [lijunpeng@login05 step2]$ OUTPUT="/public/home/lijunpeng/HESS/Cleandata/scz/covariance_correlation/step2/"
(HESS) [lijunpeng@login05 step2]$ python /public/home/lijunpeng/HESS/hess-0.5.3-beta/hess.py --prefix ${DIR1}SCZ_PUD_step1_trait1 --out ${OUTPUT}SCZ_PUD_step2_trait1
[INFO] @----------------------------------------------------------@
| HESS | v0.5 | 9/October/2017 |
|----------------------------------------------------------|
| (C) 2017 Huwenbo Shi, GNU General Public License, v3 |
|----------------------------------------------------------|
| For documentation, citation & bug-report instructions: |
| http://bogdan.bioinformatics.ucla.edu/software/hess/ |
@----------------------------------------------------------@
[INFO] Command started at: Tue, 25 Jun 2024 17:15:37
[INFO] Command issued:
/public/home/lijunpeng/HESS/hess-0.5.3-beta/hess.py
--prefix /public/home/lijunpeng/HESS/Cleandata/scz/covariance_correlation/step1/SCZ_PUD_step1_trait1
--out /public/home/lijunpeng/HESS/Cleandata/scz/covariance_correlation/step2/SCZ_PUD_step2_trait1
/public/home/lijunpeng/HESS/hess-0.5.3-beta/src/estimation.py:126: FutureWarning: read_table is deprecated, use read_csv instead, passing sep='\t'.
compression='gzip', names=['start', 'stop', 'nsnp', 'rank', 'N'])
[INFO] Loaded results for 1703 loci from step 1
[INFO] Using 4386216 SNPs with average sample size 58740.1388381
[INFO] Re-inflate the summary statistics with lambda_gc: 1
Traceback (most recent call last):
File "/public/home/lijunpeng/HESS/hess-0.5.3-beta/hess.py", line 217, in
main()
File "/public/home/lijunpeng/HESS/hess-0.5.3-beta/hess.py", line 40, in main
argmap['gwse-thres'], argmap['tot-hsqg'], argmap['out'])
File "/public/home/lijunpeng/HESS/hess-0.5.3-beta/src/estimation.py", line 174, in local_hsqg_step2
min_eigval, reinflate, gwse_thres)
File "/public/home/lijunpeng/HESS/hess-0.5.3-beta/src/estimation.py", line 281, in local_hsqg_step2_helper
local_hsqg_est = np.dot(np.linalg.pinv(A), b)
File "/public/home/lijunpeng/miniforge3/envs/HESS/lib/python2.7/site-packages/numpy/linalg/linalg.py", line 1938, in pinv
u, s, vt = svd(a, full_matrices=False)
File "/public/home/lijunpeng/miniforge3/envs/HESS/lib/python2.7/site-packages/numpy/linalg/linalg.py", line 1612, in svd
u, s, vh = gufunc(a, signature=signature, extobj=extobj)
File "/public/home/lijunpeng/miniforge3/envs/HESS/lib/python2.7/site-packages/numpy/linalg/linalg.py", line 106, in _raise_linalgerror_svd_nonconvergence
raise LinAlgError("SVD did not converge")
numpy.linalg.LinAlgError: SVD did not converge
I do not know what was going wrong? Could you do me a favor?