maenas (intertidal zone) and Segonzacia mesatlantica - Station ...

232 ANNEXE B. DISTRIBUTIONS DE MASSE DES HEXAMÈRES
1 Principe
La distribution des masses des hexamères d’Hc à partir des données de spectrométrie de masse
peut être établie avec le modèle suivant : un hexamère est formé de l’assemblage de 6 sous-unités
prises au hasard dans un ensemble de 4 types de sous-unités principales, dont les masses sont 73922,
74043, 75073 et 75187 Da. On peut déterminer l’ensemble des masses possibles des complexes en
établissant la somme des masses de toutes les combinaisons possibles de sous-unités ; si on considère
que toutes les sous-unités ont la même probabilité d’être intégrées dans un hexamère, la probabilité
d’avoir un hexamère d’une masse donnée est le rapport entre le nombre de combinaisons de sousunités
donnant cette masse et le nombre total de combinaisons possibles. On peut établir des modèles
différents en imposant par exemple la présence d’une ou plusieurs sous-unités dans l’hexamère (imposer
la présence de sous-unités lourdes ou légères par exemple).
2 Implémentation dans le langage R
La distribution des masses selon les différents modèles a été déterminée en utilisant le langage R
pour calculer les masses de chaque combinaison.
Modèle 1 : 6 sous-unités aléatoires
masses_hexa=function() { ### 6 su aléatoires
liste_masses_hexa=NULL
masses_su=c(73922,74043,75073,75184)
for (i in 1:4) {
for (j in 1:4) {
for (k in 1:4) {
for (r in 1:4) {
for (s in 1:4) {
for (q in 1:4) {
liste_masses_hexa=append(liste_masses_hexa,masses_su[i]
+masses_su[j]+masses_su[k]+masses_su[r]
+masses_su[s]+masses_su[q])}}}}}}
return(liste_masses_hexa)}
masses_hexa1=sort(masses_hexa())