@@ -479,26 +479,22 @@ def phys2cvr(
479479 if lag_map :
480480 lag , _ , _ = io .load_nifti_get_mask (lag_map )
481481 if func .shape [:3 ] != lag .shape :
482- raise ValueError (
483- f'{ lag_map } and { fname_func } have different sizes!'
484- )
482+ raise ValueError (f'{ lag_map } and { fname_func } have different sizes!' )
485483
486484 # Read lag_step and lag_max from file (or try to)
487485 lag = lag * mask
488-
486+
489487 lag_list = np .unique (lag [mask > 0 ])
490488
491489 if lag_step is None :
492- lag_step = np .unique (np .round (lag_list [1 :] - lag_list [:- 1 ],3 ))
490+ lag_step = np .unique (np .round (lag_list [1 :] - lag_list [:- 1 ], 3 ))
493491 if lag_step .size > 1 :
494492 raise ValueError (
495493 f'phys2cvr found different delta lags in { lag_map } '
496494 )
497495 else :
498- lag_step = lag_step [0 ]
499- LGR .warning (
500- f'phys2cvr detected a delta lag of { lag_step } seconds'
501- )
496+ lag_step = lag_step [0 ]
497+ LGR .warning (f'phys2cvr detected a delta lag of { lag_step } seconds' )
502498 else :
503499 LGR .warning (f'Forcing delta lag to be { lag_step } ' )
504500
@@ -634,20 +630,19 @@ def phys2cvr(
634630
635631 # If user specified a lag map, use that one to regress things
636632 if lag_map :
637-
638633 lag_idx = np .round ((lag + lag_max ) * freq / step ).astype (int )
639634 lag_idx_list = np .unique (lag_idx )
640-
635+
641636 # Prepare empty matrices
642637 beta = np .empty_like (lag , dtype = 'float32' )
643638 tstat = np .empty_like (lag , dtype = 'float32' )
644-
639+
645640 for index , i in enumerate (lag_idx_list ):
646641 LGR .info (f'Perform L-GLM number { index + 1 } of { len (lag_idx_list )} ' )
647642 regr = regr_shifts [(i * step ), :, np .newaxis ]
648643
649644 x1D = os .path .join (outdir , 'mat' , f'mat_{ i :04g} .1D' )
650-
645+
651646 (beta [lag_idx == i ], tstat [lag_idx == i ], _ ) = stats .regression (
652647 func [lag_idx == i ],
653648 regr ,
0 commit comments