Update compute_anova.py

give list of anova selective as output

src/numan_plus/compute_anova.py

@@ -224,7 +224,7 @@ def compute_anova_neurons(Q, C, Hf, alpha_level, n_permutations, filtered_idx):#
     #anova_df.to_csv(f'./caiman_final_datasets/numerosityCells_{brain_region_tag}.csv')
     #print('\033[1m\nYour number units are calculated.\033[0m\nYou can find them in ./processed/caiman_final_dataset')
 
-    return R, chance_lev, anova_df, pref_num, excitatory_or_inhibitory#, save_df
+    return anova_cells, R, chance_lev, anova_df, pref_num, excitatory_or_inhibitory#, save_df
 
 def compute_shuffled_anova_neurons(Q, C, Hf, alpha_level, n_permutations):#, save_df):
 
@@ -249,7 +249,7 @@ def compute_shuffled_anova_neurons(Q, C, Hf, alpha_level, n_permutations):#, sav
     #for n in range(6):
     #    save_df[f'Preferring_{n}'].append(sum(pref_num_shuffled==n))
 
-    return Q_S, C_S, R_S#, save_df
+    return anova_cells_shuffled, Q_S, C_S, R_S#, save_df
 
 
 def replot_tuning_curves(output_real, output_shuffled):