hardingnj / xpclr

Code to compute the XP-CLR statistic to infer natural selection
MIT License
85 stars 26 forks source link

ValueError: operands could not be broadcast together with shapes #95

Open maizecxwu opened 1 year ago

maizecxwu commented 1 year ago

My Code:

xpclr \ --format txt \ --map PANZ_SNP_Flt_AMP.chr2.xpclr.space.inf \ --popA PANZ_SNP_Flt_AMP.rmfmt.TEO_ID.chr2.vcf.gz_xpclr.V2.geno \ --popB PANZ_SNP_Flt_AMP.rmfmt.Zmay_ID.chr2.vcf.gz_xpclr.V2.geno \ --chr 1 \ --maxsnps 100 --size 1000 --step 100 --out \ /public/home/cxwu/yhdu/chr2_TEO_Zmay.xpclr

My Error: 2023-04-26 19:25:53 : INFO : running xpclr v1.1.1 2023-04-26 19:25:53 : INFO : Loading TXT 2023-04-26 19:39:08 : INFO : TXT loading complete 2023-04-26 19:39:12 : INFO : 183 SNPs in total are in the provided input files Traceback (most recent call last): File "/public/home/software/opt/bio/software/Python/3.6.5//bin/xpclr", line 4, in import('pkg_resources').run_script('xpclr==1.1.1', 'xpclr') File "/public/home/software/opt/bio/software/Python/3.6.5/lib/python3.6/site-packages/pkg_resources/init.py", line 667, in run_script self.require(requires)[0].run_script(script_name, ns) File "/public/home/software/opt/bio/software/Python/3.6.5/lib/python3.6/site-packages/pkg_resources/init.py", line 1471, in run_script exec(script_code, namespace, namespace) File "/public/home/software/opt/bio/software/xpclr/v1.1.1/lib/python3.6/site-packages/xpclr-1.1.1-py3.6.egg/EGG-INFO/scripts/xpclr", line 195, in File "/public/home/software/opt/bio/software/xpclr/v1.1.1/lib/python3.6/site-packages/xpclr-1.1.1-py3.6.egg/EGG-INFO/scripts/xpclr", line 143, in main ValueError: operands could not be broadcast together with shapes (183,) (506,) 2023-04-26 19:46:25 : INFO : running xpclr v1.1.2 2023-04-26 19:46:25 : INFO : Loading TXT 2023-04-26 19:57:55 : INFO : TXT loading complete 2023-04-26 19:57:59 : INFO : 183 SNPs in total are in the provided input files Traceback (most recent call last): File "/public/home/cxwu/~/python_package/bin/xpclr", line 4, in import('pkg_resources').run_script('xpclr==1.1.2', 'xpclr') File "/public/home/software/opt/bio/software/Python/3.6.9/lib/python3.6/site-packages/pkg_resources/init.py", line 664, in run_script self.require(requires)[0].run_script(script_name, ns) File "/public/home/software/opt/bio/software/Python/3.6.9/lib/python3.6/site-packages/pkg_resources/init.py", line 1451, in run_script exec(script_code, namespace, namespace) File "/public/home/cxwu/~/python_package/lib/python3.6/site-packages/xpclr-1.1.2-py3.6.egg/EGG-INFO/scripts/xpclr", line 196, in File "/public/home/cxwu/~/python_package/lib/python3.6/site-packages/xpclr-1.1.2-py3.6.egg/EGG-INFO/scripts/xpclr", line 143, in main ValueError: operands could not be broadcast together with shapes (183,) (506,)

My guess : maybe .geno file cannot have 9 (meaning NA) ?

Jiangjiangzhang6 commented 10 months ago

have you solve this bug, I met the same error

maizecxwu commented 10 months ago

no,i use the xpclr(python version)to calculate

发自我的iPhone

------------------ Original ------------------ From: Jiangjiangzhang6 @.> Date: Sat,Sep 16,2023 4:15 PM To: hardingnj/xpclr @.> Cc: maizecxwu @.>, Author @.> Subject: Re: [hardingnj/xpclr] ValueError: operands could not be broadcasttogether with shapes (Issue #95)

have you solve this bug, I met the same error

— Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you authored the thread.Message ID: @.***>

Jiangjiangzhang6 commented 10 months ago

this xpclr I use was the python version, I used the orignal version ,but now it wrong with bug too

| ---- Replied Message ---- | From | @.> | | Date | 9/16/2023 16:34 | | To | @.> | | Cc | @.> , @.> | | Subject | Re: [hardingnj/xpclr] ValueError: operands could not be broadcast together with shapes (Issue #95) |

no,i use the xpclr(python version)to calculate

发自我的iPhone

------------------ Original ------------------ From: Jiangjiangzhang6 @.> Date: Sat,Sep 16,2023 4:15 PM To: hardingnj/xpclr @.> Cc: maizecxwu @.>, Author @.> Subject: Re: [hardingnj/xpclr] ValueError: operands could not be broadcasttogether with shapes (Issue #95)

have you solve this bug, I met the same error

— Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you authored the thread.Message ID: @.***>

— Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you commented.Message ID: @.***>