Skip to content
This repository was archived by the owner on Dec 11, 2022. It is now read-only.
This repository was archived by the owner on Dec 11, 2022. It is now read-only.

Failed to run haps_filters through multipop_selection_pipeline #28

@AhmedElfarargi

Description

@AhmedElfarargi

I am testing the multipop_selection_pipeline with the test lactase data. After running the following command, I noticed that there's something went wrong with the selection_pipeline

multipop_selection_pipeline -p CEU_ids.txt -p YRI_ids.txt -i CEU_YRI_lactase.vcf --config-file defaults.cfg --a "--phased-vcf" -c 2

multi_population.log
selection_pipeline command = selection_pipeline -c 2 -i /SelectionPipelineTestData/lactase/CEU.vcf --population CEU --config-file defaults.cfg --phased-vcf
selection_pipeline failed to run selection_pipeline -c 2 -i /SelectionPipelineTestData/lactase/CEU.vcf --population CEU --config-file defaults.cfg --phased-vcf

CEU2_selection_pipeline.log
haps filter command = haps_filters --maf 0.01 --hwe 0.0001 --chi-sq --missing 0.99 --output CEU2_filtered.haps --haps CEU2vcf_to_haps.haps
haps filter failed to run haps_filters --maf 0.01 --hwe 0.0001 --chi-sq --missing 0.99 --output CEU2_filtered.haps --haps CEU2vcf_to_haps.haps
haps filter STDERR: Traceback (most recent call last):
haps filter STDERR: File "/bin/haps_filters", line 11, in
haps filter STDERR: load_entry_point('selectionTools==1.1.1', 'console_scripts', 'haps_filters')()
haps filter STDERR: File "/lib/python2.7/site-packages/selectionTools-1.1.1-py2.7.egg/selection_pipeline/haps_filters.py", line 123
haps filter STDERR: filter_haps_file(args)
haps filter STDERR: File "/lib/python2.7/site-packages/selectionTools-1.1.1-py2.7.egg/selection_pipeline/haps_filters.py", line 83,
haps filter STDERR: hardy_weinberg_asymptotic(countAB, countAA, countBB)
haps filter STDERR: File "/lib/python2.7/site-packages/selectionTools-1.1.1-py2.7.egg/selection_pipeline/haps_filters.py", line 37,
haps filter STDERR: return stats.chisqprob(chi_sq_total, 1)
haps filter STDERR: AttributeError: 'module' object has no attribute 'chisqprob'

should I modify something within haps_filters.py file?

Thanks

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions