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

Changing marker size for bonferoni correction #54

Closed
mtross2 opened this issue Nov 16, 2020 · 2 comments
Closed

Changing marker size for bonferoni correction #54

mtross2 opened this issue Nov 16, 2020 · 2 comments

Comments

@mtross2
Copy link

mtross2 commented Nov 16, 2020

Hi,

I'm calculated the number of effective SNPs and was wondering how I would be able to change marker size to this value so that I can have a corrected Bonferroni threshold according to the number of effective SNPs. How do I change that in rMVP? I tried doing threshold = 0.05/effective size but that seems to inflate the threshold value even more. Any suggestions would be appreciated!

@hyacz
Copy link
Collaborator

hyacz commented Nov 17, 2020

Thank you for using rMVP. Actually, the threshold is equal to -log10(threshold/n_marker). The number of markers will be automatically counted and used for Bonferroni correction. If you really want to modify it, you can set threshold = 0.05 * n_marker_old / n_marker_new.

@hyacz hyacz closed this as completed Mar 31, 2023
@ramendra-sarma
Copy link

Thank you for using rMVP. Actually, the threshold is equal to -log10(threshold/n_marker). The number of markers will be automatically counted and used for Bonferroni correction. If you really want to modify it, you can set threshold = 0.05 * n_marker_old / n_marker_new.

What code shouls I use as " threshold = 0.05 * n_marker_old / n_marker_new" generate error

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

3 participants