a-slide / pycoMeth

DNA methylation analysis downstream to Nanopolish for Oxford Nanopore DNA sequencing datasets
GNU General Public License v3.0
34 stars 2 forks source link

Documentation of output files #19

Closed BenjaminSchwessinger closed 4 years ago

BenjaminSchwessinger commented 4 years ago

Hi, Nice tool and I think I wrapped my head round it. Is there a chance to describe the outputfiles a bit more e.g. llr and such. I can guess what it all means and some clear explanation and docs will go a long way towards knowing.

Thanks!

Benjamin

a-slide commented 4 years ago

Hi @BenjaminSchwessinger, I have a section describing the header of each tool in he documentation:

Did you mean you'd like more than that ?

Cheers

BenjaminSchwessinger commented 4 years ago

Hi,

I missed this page. Apologies for this. Looks pretty comprehensive.

Now I still have some questions:

https://a-slide.github.io/pycoMeth/CpG_Aggregate/usage/

In the TSV section:

median_llr: Median of log likelihood ratios for all read mapped <- maybe explain your cut-offs for methy/unmeth like below for beds.

https://a-slide.github.io/pycoMeth/Meth_Comp/usage/

In TSV section:

neg_med / pos_med / ambiguous_med: Number of samples with a median below the negative llr threshold / above the positive llr threshold or with and ambiguous median between the 2 thresholds <- might be worth spelling out the thresholds again (>2 and <-2).

Thanks!

Great tool.

Benjamin

On 31/3/20 7:45 pm, Adrien Leger wrote:

Hi @BenjaminSchwessingerhttps://github.com/BenjaminSchwessinger, I have a section describing the header of each tool in he documentation:

Did you mean you'd like more than that ?

Cheers

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://github.com/a-slide/pycoMeth/issues/19#issuecomment-606485862, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABRMZBVCBD7XXTPPOQLNKP3RKGURDANCNFSM4LXC4LYA.

--

Benjamin Schwessinger PhD. ARC Future Fellow 2018 Senior Lecturer Division of Plant Science Research School of Biology College of Science Linnaeus Building (134), Linnaeus Way The Australian National University Canberra ACT 0200 Australia

P: +61 2 6125 7794 (Office, 3.085) benjamin.schwessinger@anu.edu.aumailto:benjamin.schwessinger@anu.edu.au lab webpagehttp://teamschwessinger.com/, @schwessingerhttp://twitter.com/schwessinger, bloghttp://blushgreengrassatafridayafternoon.wordpress.com/, google scholarhttp://scholar.google.com.au/citations?user=lEhYW3QAAAAJ&hl=en

Ambassador of ASAPbio<:"https://asapbio.org//"> Member of eLife's Early Career Advisory Group<:"https://elifesciences.org/about/people/early-career"> Co-founder of Reproducibility 4 Everyone<:"https://www.repro4everyone.org/">

CRICOS Provider # 00120C

a-slide commented 4 years ago

Hi @BenjaminSchwessinger,

I slightly updated the doc in the last version of PycoMeth

Regarding CpG_Aggregate usage, there is no cuttoff for the median, it is just the median of all observed llr for the position As for Meth_Comp usage, the default threshold is indeed 2 (like in Nanopolish) but it can actually be redefined by users with min_diff_llr.

The last version can generate HTML reports of the significant candidates.

Cheers

Ad