@@ -202,7 +202,7 @@ def _genefusion_fuse(self,modules,chrs=None):
202202
203203 return COMPUTE_SET , FUSION_SET
204204
205- def _genefusion (self ,modules ,method = 'auto' ,mode = 'auto' ,reqacc = 1e-100 ,parallel = 1 ,nobar = False ,chrs = None ):
205+ def _genefusion (self ,modules ,method = 'auto' ,mode = 'auto' ,reqacc = 1e-100 ,parallel = 1 ,nobar = False ,chrs = None , autorescore = False ):
206206
207207 COMPUTE_SET , FUSION_SET = self ._genefusion_fuse (modules ,chrs )
208208
@@ -214,7 +214,7 @@ def _genefusion(self,modules,method='auto',mode='auto',reqacc=1e-100,parallel=1,
214214 print ("Scoring" ,len (SET ),"missing (meta)-genes" )
215215
216216 # Compute missing (meta)-genes
217- R = self ._genescorer .score (SET ,method = method ,mode = mode ,reqacc = reqacc ,parallel = parallel ,nobar = nobar ,autorescore = True )
217+ R = self ._genescorer .score (SET ,method = method ,mode = mode ,reqacc = reqacc ,parallel = parallel ,nobar = nobar ,autorescore = autorescore )
218218
219219 #print(R)
220220 #print(FUSION_SET)
@@ -263,7 +263,7 @@ class chi2rank(pathwayscorer):
263263
264264 """
265265
266- def score (self ,modules ,method = 'saddle' ,mode = 'auto' ,reqacc = 1e-100 ,parallel = 1 ,nobar = False ,genes_only = False ,chrs_only = None ):
266+ def score (self ,modules ,method = 'saddle' ,mode = 'auto' ,reqacc = 1e-100 ,parallel = 1 ,nobar = False ,genes_only = False ,chrs_only = None , autorescore = True ):
267267 """
268268 Scores a set of pathways/modules
269269
@@ -274,15 +274,17 @@ def score(self,modules,method='saddle',mode='auto',reqacc=1e-100,parallel=1,noba
274274 method(string): Method to use to evaluate tail probability ('auto','davies','ruben','satterthwaite','pearson','saddle')
275275 mode(string): Precision mode to use ('','128b','100d','auto')
276276 reqacc(float): requested accuracy
277+ autorescore(bool): Automatically try to re-score failed genes
277278 nobar(bool): Show progress bar
278279 genes_only(bool): Compute only (fused)-genescores (accessible via genescorer method)
279- chrs_only(list): Only consider genes on listed chromosomes. None for all.
280+ chrs_only(list): Only consider genes on listed chromosomes. None for all
280281 """
281282 tic = time .time ()
282283
284+
283285 # Compute fusion sets
284286 if self ._fuse :
285- COMPUTE_SET ,FUSION_SET ,R = self ._genefusion (modules ,method = method ,mode = mode ,reqacc = reqacc ,parallel = parallel ,nobar = nobar ,chrs = chrs_only )
287+ COMPUTE_SET ,FUSION_SET ,R = self ._genefusion (modules ,method = method ,mode = mode ,reqacc = reqacc ,parallel = parallel ,nobar = nobar ,chrs = chrs_only , autorescore = autorescore )
286288 else :
287289 COMPUTE_SET ,FUSION_SET ,R = self ._nogenefusion (modules )
288290
@@ -400,7 +402,7 @@ class chi2perm(pathwayscorer):
400402 Genes in the background gene sets are NOT fused.
401403 """
402404
403- def score (self ,modules ,samples = 100000 ,method = 'saddle' ,mode = 'auto' ,reqacc = 1e-100 ,parallel = 1 ,nobar = False ):
405+ def score (self ,modules ,samples = 100000 ,method = 'saddle' ,mode = 'auto' ,reqacc = 1e-100 ,parallel = 1 ,nobar = False , autorescore = True ):
404406 """
405407 Scores a set of pathways/modules
406408
@@ -412,11 +414,12 @@ def score(self,modules,samples=100000,method='saddle',mode='auto',reqacc=1e-100,
412414 mode(string): Precision mode to use ('','128b','100d','auto')
413415 reqacc(float): requested accuracy
414416 nobar(bool): Show progress bar
415-
417+ autorescore(bool): Automatically try to re-score failed genes
418+
416419 """
417420 # Compute fusion sets
418421 if self ._fuse :
419- COMPUTE_SET ,FUSION_SET ,R = self ._genefusion (modules ,method = method ,mode = mode ,reqacc = reqacc ,parallel = parallel ,nobar = nobar )
422+ COMPUTE_SET ,FUSION_SET ,R = self ._genefusion (modules ,method = method ,mode = mode ,reqacc = reqacc ,parallel = parallel ,nobar = nobar , autorescore = autorescore )
420423 else :
421424 COMPUTE_SET ,FUSION_SET ,R = self ._nogenefusion (modules )
422425
0 commit comments