Index de Brillouin (langage R)

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 R 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)).

Partagez : These icons link to social bookmarking sites where readers can share and discover new web pages.
  • Digg
  • del.icio.us
  • Facebook
  • Google
  • blogmarks
  • BlogMemes Fr
  • StumbleUpon
  • Technorati
  • TwitThis
  • connotea
  • Live
  • Reddit
  • Scoopeo
  • Slashdot

2 réponses pour le moment

2 réponses à “Index de Brillouin (langage R)”

  1. 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. 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 !

Sitemap