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

Error running Tutorial for VcfHunter Genetic linkage analysis #4

Open
jhtrujillo opened this issue Jun 7, 2022 · 2 comments
Open

Comments

@jhtrujillo
Copy link

I got the following error, when I am running the tutorial:

raceback (most recent call last):
File "/biodata5/mapa-genetico-jhtrujillo/tools/VcfHunter/TestTools/../bin/vcf2popNew.1.0.py", line 881, in
if name == "main": main()
File "/biodata5/mapa-genetico-jhtrujillo/tools/VcfHunter/TestTools/../bin/vcf2popNew.1.0.py", line 804, in main
Khi_stat = KhiSquare(liste[0], Accession_header, dico_acc, DicoSegregation)
File "/biodata5/mapa-genetico-jhtrujillo/tools/VcfHunter/TestTools/../bin/vcf2popNew.1.0.py", line 326, in KhiSquare
chi2, p = stats.chisquare(FreqObs, f_exp=FreqTheor)
File "/home/estuvar4/.local/lib/python3.9/site-packages/scipy/stats/_stats_py.py", line 6911, in chisquare
return power_divergence(f_obs, f_exp=f_exp, ddof=ddof, axis=axis,
File "/home/estuvar4/.local/lib/python3.9/site-packages/scipy/stats/_stats_py.py", line 6753, in power_divergence
raise ValueError(msg)
ValueError: For each axis slice, the sum of the observed frequencies must agree with the sum of the expected frequencies to a relative tolerance of 1e-08, but the percent differences are

@guiguimartin
Copy link
Member

I assume you are running the tutorial_GeneticLinkage.md? I see that you are using python3.9. Maybe your version of python is too up to date for those programs that were developed on python3.4.10... I have observed that some functions have changed when switching to a newer version but I have not yet had time to update the scripts to these new versions of python. You can try to use them with python 3.4.10, it should work.

@sanderdebacker
Copy link

I would advise creating a virtual environment based on Python3.4.10.
Working in this environment will not mess with your default Python3 installation.

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

3 participants