Script for generating the "sharkfin plot of the p-values compared with Log Odds Ratio"

[sbridgett]

Script for generating the "Sharkfin plot (Figure 4A) of the p-values compared with Log Odds Ratio"

Thank you very much for developing and publishing the nanocompore pipeline and its python wrapper module.

I've run the steps given in the Data preparation and usage, to generate the database and the plots for individual genes and regions.

I would like to generate a plot, or plots, representing the modifications found overall, and identify which genes are most modified. I think that I can extract that information from the report file produced by the "SampCompDB.save_report()" command, then plot it. I noticed the "Sharkfin plot" in Figure 4A of your paper and was wondering how to generate it. I noticed in your paper that the sharkfin plot was included in the python wrapper, but I haven't found that function. I may have overlooked it or misunderstood:

"Finally, we provide a convenient python wrapper over the GDBM database, allowing users to interactively access simple high level functions to plot and export the results (https://nanocompore.rna.rocks/demo/SampCompDB_usage/). ........ At the time of publication the wrapper allows to generate 6 different types of publication ready plots for a given transcript including ............ and (6) the sharkfin plot of the p-values compared with Log Odds Ratio (for the GMM method).

Thank you for your help with this.