@@ -122,7 +122,7 @@ def rich_feeder_peripheral(x, sc, stat='median'):
122122 stat : str
123123 statistic to use over rich/feeder/peripheral links
124124 'mean' or 'median' (default)
125-
125+
126126 Returns
127127 -------
128128 rfp : (3, k) numpy.ndarray
@@ -139,7 +139,7 @@ def rich_feeder_peripheral(x, sc, stat='median'):
139139 This code was written by Justine Hansen who promises to fix and even
140140 optimize the code should any issues arise, provided you let her know.
141141 """
142-
142+
143143 nnodes = len (sc )
144144 mask = np .triu (np .ones (nnodes ), 1 ) > 0
145145 node_degree = degrees_und (sc )
@@ -177,11 +177,11 @@ def rich_feeder_peripheral(x, sc, stat='median'):
177177 rfp [2 , i ] = np .mean (x [rfp_label [:, i ] == 3 ]) # peripheral
178178
179179 # p-value (Welch's t-test)
180- _ , pvals [0 ,i ] = ttest_ind (x [rfp_label [:, i ] == 1 ], \
181- x [rfp_label [:, i ] != 1 ], equal_var = False )
182- _ , pvals [1 ,i ] = ttest_ind (x [rfp_label [:, i ] == 2 ], \
183- x [rfp_label [:, i ] == 3 ], equal_var = False )
184- _ , pvals [2 ,i ] = ttest_ind (x [rfp_label [:, i ] == 3 ], \
185- x [rfp_label [:, i ] == 2 ], equal_var = False )
180+ _ , pvals [0 , i ] = ttest_ind (x [rfp_label [:, i ] == 1 ],
181+ x [rfp_label [:, i ] != 1 ], equal_var = False )
182+ _ , pvals [1 , i ] = ttest_ind (x [rfp_label [:, i ] == 2 ],
183+ x [rfp_label [:, i ] == 3 ], equal_var = False )
184+ _ , pvals [2 , i ] = ttest_ind (x [rfp_label [:, i ] == 3 ],
185+ x [rfp_label [:, i ] == 2 ], equal_var = False )
186186
187187 return rfp , pvals
0 commit comments