diff --git a/PRO-seq/PRO-seq.snakefile b/PRO-seq/PRO-seq.snakefile index c992b56da37e453d9025878de14ac176f3c7b813..0f3ffef1e86e02933c818ca8b05e13988582d285 100644 --- a/PRO-seq/PRO-seq.snakefile +++ b/PRO-seq/PRO-seq.snakefile @@ -1198,7 +1198,7 @@ rule merge_bigwig_reps: from rpy2.robjects import Formula, StrVector -from rpy2.rinterface import RRuntimeError +#from rpy2.rinterface import RRuntimeError rule differential_expression: input: counts_table = source_counts, @@ -1227,7 +1227,8 @@ rule differential_expression: contrast = StrVector(["lib", cond, ref]) formula = Formula("~ lib") res, size_factors = do_deseq2(COND_NAMES, CONDITIONS, counts_data, formula=formula, contrast=contrast) - except RRuntimeError as e: + #except RRuntimeError as e: + except RuntimeError as e: warnings.warn( "Probably not enough usable data points to perform DESeq2 analyses:\n%s\nSkipping %s_%s_%s" % ( str(e), wildcards.contrast, wildcards.orientation, wildcards.biotype)) diff --git a/RNA_Seq_Cecere/RNA-seq.snakefile b/RNA_Seq_Cecere/RNA-seq.snakefile index 9d94427fda04677dac8842934321a21a73155f97..2cdada560f90ec7110fcf1c80fd8039bea1832dc 100644 --- a/RNA_Seq_Cecere/RNA-seq.snakefile +++ b/RNA_Seq_Cecere/RNA-seq.snakefile @@ -1733,7 +1733,7 @@ rule compute_RPM_folds: from rpy2.robjects import Formula, StrVector -from rpy2.rinterface import RRuntimeError +#from rpy2.rinterface import RRuntimeError rule differential_expression: input: counts_table = source_counts, @@ -1769,7 +1769,8 @@ rule differential_expression: contrast = StrVector(["lib", cond, ref]) try: res, size_factors = do_deseq2(COND_NAMES, CONDITIONS, counts_data, formula=formula, contrast=contrast) - except RRuntimeError as e: + #except RRuntimeError as e: + except RuntimeError as e: warnings.warn( "Probably not enough usable data points to perform DESeq2 analyses:\n%s\nSkipping %s_%s_%s" % ( str(e), wildcards.contrast, wildcards.orientation, wildcards.biotype)) diff --git a/Ribo-seq/Ribo-seq.snakefile b/Ribo-seq/Ribo-seq.snakefile index 6ccfee8e4f4e8dd8d487b14789576d3a24ea34c9..f9391b595ba202d88373dc1fb3b46e432570b4ab 100644 --- a/Ribo-seq/Ribo-seq.snakefile +++ b/Ribo-seq/Ribo-seq.snakefile @@ -1052,7 +1052,7 @@ rule gather_read_counts_summaries: from rpy2.robjects import Formula, StrVector -from rpy2.rinterface import RRuntimeError +#from rpy2.rinterface import RRuntimeError rule differential_expression: input: counts_table = source_counts, @@ -1092,7 +1092,8 @@ rule differential_expression: contrast = StrVector(["lib", cond, ref]) try: res, size_factors = do_deseq2(COND_NAMES, CONDITIONS, counts_data, formula=formula, contrast=contrast) - except RRuntimeError as e: + #except RRuntimeError as e: + except RuntimeError as e: warnings.warn( "Probably not enough usable data points to perform DESeq2 analyses:\n%s\nSkipping %s_%s_%s_%s" % ( str(e), wildcards.read_type, wildcards.contrast, wildcards.orientation, wildcards.biotype)) @@ -2575,7 +2576,7 @@ def plot_counts_scatters(counts_and_res, cols): ################################################### from rpy2.robjects import Formula, StrVector -from rpy2.rinterface import RRuntimeError +#from rpy2.rinterface import RRuntimeError rule small_RNA_differential_expression: input: counts_table = source_counts, @@ -2611,7 +2612,8 @@ rule small_RNA_differential_expression: # raise NotImplementedError("Unknown contrast: %s" % wildcards.contrast) try: res, size_factors = do_deseq2(COND_NAMES, CONDITIONS, counts_data, formula=formula, contrast=contrast) - except RRuntimeError as e: + #except RRuntimeError as e: + except RuntimeError as e: warn( "Probably not enough usable data points to perform DESeq2 analyses:\n%s\nSkipping %s_%s_%s" % ( str(e), wildcards.contrast, wildcards.orientation, wildcards.small_type)) diff --git a/small_RNA-seq/small_RNA-seq.snakefile b/small_RNA-seq/small_RNA-seq.snakefile index 8c516108b200ccb36029a6f17f201a6354d8e286..547017dbc99df73717ecd2395b7627fef0d36f61 100644 --- a/small_RNA-seq/small_RNA-seq.snakefile +++ b/small_RNA-seq/small_RNA-seq.snakefile @@ -3636,7 +3636,7 @@ def plot_counts_scatters(counts_and_res, cols): ################################################### from rpy2.robjects import Formula, StrVector -from rpy2.rinterface import RRuntimeError +#from rpy2.rinterface import RRuntimeError rule small_RNA_differential_expression: input: counts_table = source_small_RNA_counts, @@ -3668,7 +3668,8 @@ rule small_RNA_differential_expression: contrast = StrVector(["lib", cond, ref]) try: res, size_factors = do_deseq2(COND_NAMES, CONDITIONS, counts_data, formula=formula, contrast=contrast) - except RRuntimeError as e: + #except RRuntimeError as e: + except RuntimeError as e: warn( "Probably not enough usable data points to perform DESeq2 analyses:\n%s\nSkipping %s_%s_%s" % ( str(e), wildcards.contrast, wildcards.orientation, wildcards.small_type))