diff --git a/RNA_Seq_Cecere/RNA-seq.snakefile b/RNA_Seq_Cecere/RNA-seq.snakefile index 1bb65eeed75ac62ac2c3714917eb35d41a57a6bb..51607f95d10a96cabcb4dc81951e9ab0c4323663 100644 --- a/RNA_Seq_Cecere/RNA-seq.snakefile +++ b/RNA_Seq_Cecere/RNA-seq.snakefile @@ -1276,8 +1276,14 @@ rule compare_replicates: (r_val, p_val) = pearsonr(data[rep_a], data[rep_b]) corr_file.write(f"{rep_a}{warn_a}_vs_{rep_b}{warn_b}\t{r_val} ({p_val})\n") filtered = log_data[[rep_a, rep_b]].dropna() - (log_r_val, log_p_val) = pearsonr(filtered[rep_a], filtered[rep_b]) - corr_file.write(f"log10_{rep_a}{warn_a}_vs_log10_{rep_b}{warn_b}\t{log_r_val} ({log_p_val})\n") + try: + (log_r_val, log_p_val) = pearsonr(filtered[rep_a], filtered[rep_b]) + corr_file.write(f"log10_{rep_a}{warn_a}_vs_log10_{rep_b}{warn_b}\t{log_r_val} ({log_p_val})\n") + except ValueError as err: + if str(err) == "x and y must have length at least 2.": + corr_file.write(f"Only {len(filtered)} points in log10 data.\n") + else: + raise corr_file.write("(*): not normal\n(**): log10 not normal\n")