@@ -72,11 +72,7 @@ def calc_loccenter(x, lin_log_flag):
7272 w = x [i ,ind ]/ sum (x [i ,ind ], 0 )
7373 else :
7474 w = (2 ** x [i ,ind ])/ sum (2 ** x [i ,ind ], 0 )
75- < << << << HEAD
7675 loc_center [i ] = sum (w * ind , 0 )
77- == == == =
78- loc_center [i ] = sum (w * ind , 0 )
79- >> >> >> > master
8076 else :
8177 loc_center [i ] = 0
8278
@@ -490,19 +486,12 @@ def _divide_to_2and_resort(sorted_data, wid, iters_spin=8, stop_const = 1.15, lo
490486 gr1 = arange (N )[:breakp1 ]
491487 gr2 = arange (N )[breakp1 :]
492488 # and assign the genes into the two groups
493- < << << << HEAD
494489 sorted_gr1 = sorted_data [:,gr1 ]
495490 sorted_gr2 = sorted_data [:,gr2 ]
496491 mean_gr1 = sorted_gr1 .mean (1 )
497492 mean_gr2 = sorted_gr2 .mean (1 )
498493 concat_loccenter_gr1 = c_ [ calc_loccenter (sorted_gr1 , 2 ), calc_loccenter (sorted_gr1 [...,::- 1 ], 2 ) ]
499494 concat_loccenter_gr2 = c_ [ calc_loccenter (sorted_gr2 , 2 ), calc_loccenter (sorted_gr2 [...,::- 1 ], 2 ) ]
500- == == == =
501- mean_gr1 = sorted_data [:,gr1 ].mean (1 )
502- mean_gr2 = sorted_data [:,gr2 ].mean (1 )
503- concat_loccenter_gr1 = c_ [ calc_loccenter (sorted_data [:,gr1 ], 2 ), calc_loccenter (sorted_data [:,gr1 ][...,::- 1 ], 2 ) ]
504- concat_loccenter_gr2 = c_ [ calc_loccenter (sorted_data [:,gr2 ], 2 ), calc_loccenter (sorted_data [:,gr2 ][...,::- 1 ], 2 ) ]
505- >> >> >> > master
506495 center_gr1 , flip_flag1 = concat_loccenter_gr1 .min (1 ), concat_loccenter_gr1 .argmin (1 )
507496 center_gr2 , flip_flag2 = concat_loccenter_gr2 .max (1 ), concat_loccenter_gr2 .argmax (1 )
508497 sorted_data_tmp = array ( sorted_data )
@@ -955,7 +944,7 @@ def usage():
955944
956945== == == =
957946
958-
947+
959948>> >> >> > master
960949
961950 output_cef .writeCEF ( outfiles_path )
0 commit comments