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

Add warning for negative distances during rate-adjustment #49

Merged
merged 2 commits into from
May 10, 2023

Conversation

Cecilia-Sensalari
Copy link
Collaborator

@Cecilia-Sensalari Cecilia-Sensalari commented May 10, 2023

During rate-adjustment, ksrates prints a warning if there are negative values for Ks distances, stating that they are artefacts.
More details below.

The rate-adjusted ortholog peak (K_AB) is derived through the following mathematical formulas, using three estimated ortholog distribution peaks (A: focal species, B: sister species, C: outgroup):
image
image

The numerator in (2) should always in principle be a positive number, as it is a measure of genetic distance. However, poor peak estimate for the three AC, AB and BC ortholog distributions might produce a negative sum, and therefore a negative K_OA value. The negative sign is then maintained also in (4).
Poor peak estimate is often due to the very old divergence age between the involved species or by bad quality sequence data.

@Cecilia-Sensalari Cecilia-Sensalari merged commit 46297d7 into master May 10, 2023
@Cecilia-Sensalari Cecilia-Sensalari deleted the warn_neg_rate branch May 10, 2023 14:17
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

Successfully merging this pull request may close these issues.

1 participant