Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Peak Scores when using LanceOtron #24

Open
zhanglongnation opened this issue Jun 5, 2024 · 2 comments
Open

Peak Scores when using LanceOtron #24

zhanglongnation opened this issue Jun 5, 2024 · 2 comments

Comments

@zhanglongnation
Copy link

Hi,
I have an H3K27ac CUT&Tag sequecing data for a non-model animal. A bw file was generated accoding to the protocol. After upload to the web, the process fail. Do you have any idea about the problem?
I install the LanceOtron on my server, and it runs successfully. A bed file was generated, and the header contains many contents. Can you give a detailed explanation about the contents meaning? I noticed that overall_peak_score was in the generated bed file, is it resonable to exclued some peaks under a threashold?
Many thanks

@LHentges
Copy link
Owner

LHentges commented Jun 7, 2024

My co-creator is looking into your issue using the webtool.

The overall_peak_score is the metric we published on in our paper, and is a combined assessment of both peak shape and multi-way enrichment. To answer your question, I think this is the best single score to filter on. In the paper we used a score of 0.5, which is equivalent to the model believing there is at least a 50% chance that the enriched region is derived from a biological signal rather than noise.

@zhanglongnation
Copy link
Author

OK, I received the mail and it works now. Thanks a lot. Another question, peak calling using macs2 and Lanceotron may give different enrichment regions for the same region, do you have any suggestions for distinguishing which is more correct called peak? like the peak called in the attached file. LanceOtron give 4 peaks in the region, but macs2 gives almost 8 peaks, although it is more scatter.
Peak_call_mcas2&LanceOtron.docx

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants