ValueError: operands could not be broadcast together with shapes (30,30) (2,2)

[mateushg1]

The number of chromosomes (contigs in .vcf HEADER): 22

Two Pulse Model:
Group_1 - admixed population! Sample size: 2983
north - admixed one time! Sample size: 20
south - admixed two times! Sample size: 28

estimating weighted LD...
CHR: 1
variants count (.vcf): 10157
10157
0:00:12.928536 - data read.
0:00:00.072327 - estimation.
1/22
CHR: 2
variants count (.vcf): 11004
11004
0:00:14.016453 - data read.
0:00:00.076937 - estimation.
2/22
CHR: 3
variants count (.vcf): 9235
9235
0:00:11.930279 - data read.
0:00:00.068707 - estimation.
3/22
CHR: 4
variants count (.vcf): 7692
7692
0:00:09.927802 - data read.
0:00:00.058400 - estimation.
4/22
CHR: 5
variants count (.vcf): 8169
8169
0:00:10.525183 - data read.
0:00:00.060438 - estimation.
5/22
CHR: 6
variants count (.vcf): 8475
8475
0:00:10.902056 - data read.
0:00:00.060794 - estimation.
6/22
CHR: 7
variants count (.vcf): 6927
6927
0:00:09.029840 - data read.
0:00:00.052100 - estimation.
7/22
CHR: 8
variants count (.vcf): 7376
7376
0:00:09.598808 - data read.
0:00:00.054676 - estimation.
8/22
CHR: 9
variants count (.vcf): 6292
6292
0:00:08.288601 - data read.
0:00:00.048536 - estimation.
9/22
CHR: 10
variants count (.vcf): 7151
7151
0:00:09.254934 - data read.
0:00:00.054075 - estimation.
10/22
CHR: 11
variants count (.vcf): 6587
6587
0:00:08.577645 - data read.
0:00:00.049995 - estimation.
11/22
CHR: 12
variants count (.vcf): 6587
6587
0:00:08.592770 - data read.
0:00:00.050457 - estimation.
12/22
CHR: 13
variants count (.vcf): 5123
5123
0:00:06.776116 - data read.
0:00:00.042195 - estimation.
13/22
CHR: 14
variants count (.vcf): 4443
4443
0:00:05.895077 - data read.
0:00:00.038723 - estimation.
14/22
CHR: 15
variants count (.vcf): 4068
4068
0:00:05.491133 - data read.
0:00:00.036777 - estimation.
15/22
CHR: 16
variants count (.vcf): 4136
4136
0:00:05.561663 - data read.
0:00:00.037018 - estimation.
16/22
CHR: 17
variants count (.vcf): 2974
2974
0:00:04.143150 - data read.
0:00:00.030630 - estimation.
17/22
CHR: 18
variants count (.vcf): 4037
4037
0:00:05.458612 - data read.
0:00:00.036575 - estimation.
18/22
CHR: 19
variants count (.vcf): 1469
1469
0:00:02.278955 - data read.
0:00:00.022641 - estimation.
19/22
CHR: 20
variants count (.vcf): 3515
3515
0:00:04.787149 - data read.
0:00:00.034019 - estimation.
20/22
CHR: 21
variants count (.vcf): 1841
1841
0:00:02.747813 - data read.
0:00:00.023686 - estimation.
21/22
CHR: 22
variants count (.vcf): 1658
1658
0:00:02.498437 - data read.
0:00:00.022554 - estimation.
22/22

estimating parameters...
cm min: 0.5
cm max: 30

free variables: t1 t2 m1 m2 mt
Traceback (most recent call last):
File "/gpfs/gsfs11/users/gouveiamh/softwares/LaNeta/LaNeta/laneta.py", line 130, in
parameters = exp.estimateTwoPulse(cm_max=clargs.max,
File "/gpfs/gsfs11/users/gouveiamh/softwares/LaNeta/LaNeta/LaNeta/ThLd.py", line 835, in estimateTwoPulse
parameters = self.estimate_parameters(T1=T1, T2=T2, M1=M1, M2=M2,
File "/gpfs/gsfs11/users/gouveiamh/softwares/LaNeta/LaNeta/LaNeta/ThLd.py", line 981, in estimate_parameters
res = least_squares(metr, x0, bounds=(bounds, bounds))
File "/usr/local/Anaconda/envs/py3.9/lib/python3.9/site-packages/scipy/optimize/_lsq/least_squares.py", line 830, in least_squares
f0 = fun_wrapped(x0)
File "/usr/local/Anaconda/envs/py3.9/lib/python3.9/site-packages/scipy/optimize/_lsq/least_squares.py", line 825, in fun_wrapped
return np.atleast_1d(fun(x, *args, **kwargs))
File "/gpfs/gsfs11/users/gouveiamh/softwares/LaNeta/LaNeta/LaNeta/ThLd.py", line 957, in metr
return (th_E[P:P, P:P] - self.WLD.ld[P:P, P:P]).reshape(-1)
ValueError: operands could not be broadcast together with shapes (30,30) (2,2)