@@ -226,18 +226,18 @@ def get_test_data(matrix_text, par):
226
226
x = convert_to_numpy (matrix_text )
227
227
if par == "switch" :
228
228
# Switch probability, one per site
229
- return x [:, 2 ].reshape ((4 , 4 ))[:, 2 ]
229
+ return x [:, 2 ].reshape ((4 , 4 ))[:, 0 ]
230
230
elif par == "mismatch" :
231
- # Mismatch probability
232
- return x [:, 2 ].reshape ((4 , 4 ))[:, 4 ]
231
+ # Mismatch probability, one per site
232
+ return x [:, 4 ].reshape ((4 , 4 ))[:, 0 ]
233
233
elif par == "ref_hap_allele" :
234
234
# Allele in haplotype in reference panel
235
235
# 0 = ref allele, 1 = alt allele
236
- return x [:, 2 ].reshape ((4 , 4 ))[:, 6 ]
236
+ return x [:, 6 ].reshape ((4 , 4 ))
237
237
elif par == "query_hap_allele" :
238
238
# Allele in haplotype in query
239
239
# 0 = ref allele, 1 = alt allele
240
- return x [:, 2 ].reshape ((4 , 4 ))[:, 7 ]
240
+ return x [:, 7 ].reshape ((4 , 4 ))[:, 0 ]
241
241
elif par == "shift" :
242
242
# Shift factor
243
243
# TODO
@@ -248,6 +248,6 @@ def get_test_data(matrix_text, par):
248
248
pass
249
249
elif par == "sum" :
250
250
# Sum of values over haplotypes
251
- return x [:, 2 ].reshape ((4 , 4 ))[:, 10 ]
251
+ return x [:, 10 ].reshape ((4 , 4 ))[:, 0 ]
252
252
else :
253
253
raise ValueError (f"Unknown parameter: { par } " )
0 commit comments