@@ -167,17 +167,20 @@ def rich_feeder_peripheral(x, sc, stat='median'):
167167 redfunc = np .median if stat == 'median' else np .mean
168168 for linktype in range (3 ):
169169 rfp [linktype , degthresh ] = redfunc (x [rfp_label [:, degthresh ] ==
170- linktype + 1 ])
170+ linktype + 1 ])
171171
172172 # p-value (one-sided Welch's t-test)
173173 _ , pvals [0 , degthresh ] = ttest_ind (x [rfp_label [:, degthresh ] == 1 ],
174- x [rfp_label [:, degthresh ] != 1 ],
175- equal_var = False , alternative = 'greater' )
174+ x [rfp_label [:, degthresh ] != 1 ],
175+ equal_var = False ,
176+ alternative = 'greater' )
176177 _ , pvals [1 , degthresh ] = ttest_ind (x [rfp_label [:, degthresh ] == 2 ],
177- x [rfp_label [:, degthresh ] == 3 ],
178- equal_var = False , alternative = 'greater' )
178+ x [rfp_label [:, degthresh ] == 3 ],
179+ equal_var = False ,
180+ alternative = 'greater' )
179181 _ , pvals [2 , degthresh ] = ttest_ind (x [rfp_label [:, degthresh ] == 3 ],
180- x [rfp_label [:, degthresh ] == 2 ],
181- equal_var = False , alternative = 'greater' )
182+ x [rfp_label [:, degthresh ] == 2 ],
183+ equal_var = False ,
184+ alternative = 'greater' )
182185
183186 return rfp , pvals
0 commit comments