Skip to content
Snippets Groups Projects
Commit 2e2dcd54 authored by Julien  GUGLIELMINI's avatar Julien GUGLIELMINI
Browse files

README update

parent 528c7733
No related branches found
No related tags found
No related merge requests found
...@@ -7,11 +7,12 @@ wGRR is a set of scripts to calculate the weighted Gene Repertoire Relatedness ( ...@@ -7,11 +7,12 @@ wGRR is a set of scripts to calculate the weighted Gene Repertoire Relatedness (
wGRR(A,B) = \frac{\sum_{i}{id(A_i,B_i)}}{min(P_A,P_B)} wGRR(A,B) = \frac{\sum_{i}{id(A_i,B_i)}}{min(P_A,P_B)}
``` ```
where where
* $`id(A_i,B_i)`$ is the identity score for each BBH pair between element $`A`$ and element $`B`$ * $`i`$ is the number of BBH between elements $`A`$ and $`B`$
* $`min(P_A,P_B)`$ is the number of proteins in the smallest of $`A`$ and $`B`$ elements. * $`id(A_i,B_i)`$ is the identity score of BBH $`i`$
* $`P_A`$ and $`P_B`$ are the numbers of proteins of $`A`$ and $`B`$ respectively
**NOTE** **NOTE**
`wGRR` calculates 3 versions of the wGRR depending on which BBH pairs are considered (numerator) and what proteins should be counted (denominator). See the `Output` section of this manual for more explanations. `wGRR` calculates 3 versions of the wGRR score depending on which BBH are considered (numerator) and what proteins should be counted (denominator). See the `Output` section of this manual for more explanations.
### Dependencies ### Dependencies
BBH are defined by all versus all protein comparisons using [MMseqs2][1]. BBH are defined by all versus all protein comparisons using [MMseqs2][1].
...@@ -105,6 +106,6 @@ wGRR3 also uses the protein families. But this time, if two BBH pairs are found ...@@ -105,6 +106,6 @@ wGRR3 also uses the protein families. But this time, if two BBH pairs are found
Common3 is the number of protein families containing at least one BBH divided by the mean nProt3 of two elements. Common3 is the number of protein families containing at least one BBH divided by the mean nProt3 of two elements.
![wgrr.png](https://gitlab.pasteur.fr/jgugliel/wgrr/-/raw/main/wgrr.png) ![wgrr.png](https://uc19a2ea477fef500ecb8a39fd57.dl.dropboxusercontent.com/cd/0/inline/BnUpL9u4vppyAsK9XPDgE2INoPoLdGWQCmgFf8ztrB4q2Y10nP1ElHmdoNujNc8pwMB9j6iGgm_xme45W12ivgAe7ffaK8HZ88-PDEWkvdmUCyBcK6EifEeCCYEEdM3WVksE49oSD1nKzzjF9WJAwIeZf6Idm1gdlDye_9bvb4uF9A/file#)
[1]: https://mmseqs.com/ [1]: https://mmseqs.com/
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment