novoalab / EpiNano

Detection of RNA modifications from Oxford Nanopore direct RNA sequencing reads (Liu*, Begik* et al., Nature Comm 2019)
GNU General Public License v2.0
108 stars 31 forks source link

Unable to plot results #69

Closed Wardale24 closed 3 years ago

Wardale24 commented 3 years ago

After running Epinano_Predict.py with the provided rrach.q3.mis3.del3.linear.dump file, I attempted to plot my results after filtering the rrach motifs.

First of all I wasn't sure whether to use EPINANO_HOME/Epinano_Plot.R or EPINANO_HOME/misc/Epinano_Plot_Prob.R , since in the wiki you mention "Rscript $EPINANO_HOME/misc/Epinano_Plot.R SVM_prediction.csv" and after "Rscript $EPINANO_HOME/Epinano_Plot.R diffErr.delta-mis.prediction.csv". Although I'm pretty sure now it is EPINANO_HOME/Epinano_Plot.R after looking through the script, I would still appreciate confirmation this is correct.

However, when running "Rscript EPINANO_HOME/Epinano_Plot.R my_rrach_modification_prediction.csv" I get no output, no errors, no empty pdf.

I am using Epinano1.2 since the plotting scripts were only available in this update. And the head of my rrach.csv file is:

GAACA,6360-6364,chr1,-,1:1:1:1:1,7.0,7.0,7.0,21.0,10.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,mod,-2.9078354219478157,0.9835656588992987,0.016434341100701327 GAACA,6372-6376,chr1,-,1:1:1:1:1,3.0,4.0,8.0,19.0,7.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,mod,-2.5956296188947974,0.9741645009292426,0.02583549907075754 AGACT,6389-6393,chr1,-,1:1:1:1:1,2.0,2.0,0.0,8.0,10.5,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,1.0,0.0,0.0,mod,-5.150868401500868,0.9993704871974265,0.000629512802573528 AAACC,6412-6416,chr1,-,2:2:2:2:2,12.0,10.0,16.5,20.0,23.0,0.0,0.0,0.0,0.0,0.0,0.5,0.5,0.0,0.0,0.0,0.0,0.0,0.0,0.5,0.5,unm,0.05811970706343317,0.44806966643190144,0.5519303335680985 GAACA,6441-6445,chr1,-,2:2:2:2:2,12.0,11.0,19.5,27.0,26.5,0.5,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,unm,0.9947371162251883,0.17200618172401694,0.8279938182759831

Wardale24 commented 3 years ago

When filtering for the drach motifs I did not include the header line, my mistake.