Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

About hicPCA error #898

Open
HW-Superior opened this issue Mar 31, 2024 · 1 comment
Open

About hicPCA error #898

HW-Superior opened this issue Mar 31, 2024 · 1 comment

Comments

@HW-Superior
Copy link

HW-Superior commented Mar 31, 2024

Hi ,

I'm experiencing an issue with hicPCA, where it gives a ValueError. I'm using python version 3.8.19 and hicexplorer version 3.7.2. Here is the command I run:
hicPCA --matrix CX6_200000.h5 --whichEigenvectors 1 --outputFileName CX6.200k.hicpca.bedgraph --format bedgraph --extraTrack gene.bed --pearsonMatrix CX6.200k.hicpca.pm --obsexpMatrix CX6.200k.hicpca.oem

The error as Traceback (most recent call last):
File "/home/wanghao/miniconda3/envs/hicexplorer/bin/hicPCA", line 7, in
main()
File "/home/wanghao/miniconda3/envs/hicexplorer/lib/python3.8/site-packages/hicexplorer/hicPCA.py", line 355, in main
vecs_list = correlateEigenvectorWithGeneTrack(ma, vecs_list, args.extraTrack)
File "/home/wanghao/miniconda3/envs/hicexplorer/lib/python3.8/site-packages/hicexplorer/hicPCA.py", line 190, in correlateEigenvectorWithGeneTrack
_correlation = pearsonr(np.array(eigenvector[bin_id[0]:bin_id[1]]).real,
File "/home/wanghao/miniconda3/envs/hicexplorer/lib/python3.8/site-packages/scipy/stats/_stats_py.py", line 4411, in pearsonr
raise ValueError('x and y must have length at least 2.')
ValueError: x and y must have length at least 2.
The picture below shows the gene.bed file. When I remove --extraTrack gene.bed parameter, the output will be normal. How should I solve this problem?
Thanks
image

@shRDJr
Copy link

shRDJr commented Aug 12, 2024

Hello, I also encountered the same mistake, how did you solve it later?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants