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

Confusion about expected number of CNV calls #58

Open
tsin-lau opened this issue Sep 21, 2023 · 0 comments
Open

Confusion about expected number of CNV calls #58

tsin-lau opened this issue Sep 21, 2023 · 0 comments

Comments

@tsin-lau
Copy link

Thank you for developing this excellent tool for CNV calling and it works well with my WES data. However, I'm still confused about the number of CNVs of each exome.
As mentioned in ExomeDepth-vignette, you usually obtained 150-280 CNV calls per exome sample. But according to the reference manual, the argument, transition probability, in the CallCNVs function is described as follows —— Transition probability of the hidden Markov Chain from the normal copy number state to either a deletion or a duplication. The default (0.0001) expect approximately 20 CNVs genome-wide.
I wonder where the discrepancy comes from and how the transition probability argument can be modified.
THANKS!

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

1 participant