R : index de Brillouin

Posté par Timothée le 4 May 2007 | , , ,

Je suis en pleine séance de stats, et j’ai besoin de calculer l’index de diversité de Brillouin, qui me donne la diversité des parasites présents sur chaque individu. Excel ayant quelques difficultés, j’ai écrit une petite fonction en qui me fait ça tout à fait bien.

Index de Brillouin, ou N est le nombre total de parasites, et ni le nombre de parasites pour chaque espèce.Dans cet index, N représente le nombre total d’individus parasites sur chaque individu hôte, et ni est le nombre d’individus d’une espèce i. Je passe sur les raisons de son utilisation, pour donner la formule:

brillouin <- function(pop) {(log(factorial(sum(pop)))-sum(log(factorial(pop))))/sum(pop)}

ou pop est un array qui contient l’ensemble des valeurs ni. La fonction s’utilise en faisant brillouin(variable), ou directement brillouin(c(n1,n2,…nk)).

Billets similaires

2 réponses pour le moment

2 réponses à “R : index de Brillouin”

  1. # Tom Roud

    Remarque peut-être sans intérêt, mais à mon avis ce serait beaucoup plus rapide et précis d’utiliser le fait que log(factorial(N))= sum(log i) pour i de 1 à N.

    04 May 2007 à 3:33 pm

  2. # Benjamin

    Comme il n’y a pas d’algo dedans, ce serait peut être plus pratique avec une cellule excel dans le tableau où tu rentres tes chiffres, non?

    10 May 2007 à 5:51 pm

Trackback URI | Comments RSS

Commentez !