International audienceLes communautés de macrophytes (algues et phanérogames), dont certaines espèces sont protégées, constituent un des indicateurs clés de l'état écologique des lagunes, tel que défini par la Directive Européenne Cadre sur l'Eau (DCE). Par ailleurs, la restauration de ces communautés de macrophytes, et des herbiers de Zostera en particulier, est l'un des objectifs principaux de la réhabilitation de l'Etang de Berre (Bouches-du-Rhône) portée par le Gipreb (Groupement d'intérêt public pour la réhabilitation de l'Etang de Berre). C'est pourquoi celui-ci met en oeuvre, depuis 1996, le suivi des principales espèces de macrophytes encore présentes dans l'étang. La densité de chaque espèce ou groupe d'espèces a été évaluée visuellement par des plongeurs le long de 31 transects (composés chacun de 20 segments de même longueur) répartis sur le pourtour de l'étang, et ce pendant 11 années. L'évaluation consiste en un codage en 6 modalités du taux de couverture du fonds par l'espèce. L'état annuel de la population de chaque espèce est finalement décrit par un tableau de type 6x20, croisant la densité avec la position : deux variables ordinales. La suite des 11 tableaux sera donc ici décrite via l'Analyse des Correspondances Ordinales de Beh (1997). Les résultats obtenus seront comparés avec ceux de l'AFC. Bibliographie : E. J. Beh, Simple Correspondence Analysis of ordinal cross-classifications using orthogonal polynomials, Biometrical Journal, 39, 5, 589-613 (1997)

Bernard, Guillaume

Bonhomme, Patrick

Manté, Claude

HAL AMU

Application de l’Analyse des Correspondances Ordinalesau suivi d’espe`ces ve´ge´tales aquatiques.Claude Mante´, Guillaume Bernard, Patrick BonhommeTo cite this version:Claude Mante´, Guillaume Bernard, Patrick Bonhomme. Application de l’Analyse des Corre-spondances Ordinales au suivi d’espe`ces ve´ge´tales aquatiques.. 42e`mes Journe´es de Statistique,2010, Marseille, France, France. 2010. <inria-00494699>HAL Id: inria-00494699https://hal.inria.fr/inria-00494699Submitted on 24 Jun 2010HAL is a multi-disciplinary open accessarchive for the deposit and dissemination of sci-entific research documents, whether they are pub-lished or not. The documents may come fromteaching and research institutions in France orabroad, or from public or private research centers.L’archive ouverte pluridisciplinaire HAL, estdestine´e au de´poˆt et a` la diffusion de documentsscientifiques de niveau recherche, publie´s ou non,e´manant des e´tablissements d’enseignement et derecherche franc¸ais ou e´trangers, des laboratoirespublics ou prive´s.APPLICATION DE L’ANALYSE DES CORRESPONDANCES ORDINALES AU SUIVI D’ESPÈCES VÉGÉTALES AQUATIQUES. Claude Manté1, Guillaume Bernard2 et Patrick Bonhomme3  1LMGEM, UMR CNRS 6117, Campus de Luminy - 13009 Marseille 2GIPREB, Cours Mirabeau - 13130 Berre-l'étang 3GIS Posidonie, Campus de Luminy - 13009 Marseille  Résumé Les communautés de macrophytes (algues et phanérogames), dont certaines espèces sont protégées, constituent un des indicateurs clés de l’état écologique des lagunes, tel que défini par la Directive Européenne Cadre sur l’Eau (DCE). Par ailleurs, la restauration de ces communautés de macrophytes, et des herbiers de Zostera en particulier, est l’un des objectifs principaux du plan de réhabilitation de l’Etang de Berre (Bouches-du-Rhône) porté par le Gipreb (Groupement d’intérêt public pour la réhabilitation de l’Etang de Berre). C’est pourquoi celui-ci met en œuvre, depuis 1996, le suivi annuel des principales espèces de macrophytes encore présentes dans l’étang, afin d’évaluer l’impact sur ce compartiment de mesures de régulation des rejets d’eau douce d’une centrale hydroélectrique, principale source de perturbation du système. Une suite de onze tableaux annuels relatifs à la densité de l’espèce Zostera noltii est ici étudiée via l’Analyse des Correspondances Ordinales de Beh (1997, 2004), dont les principes sont rappelés. Les résultats obtenus seront comparés avec ceux de l’AFC classique. Mots clés : Analyse des données, Environnement Abstract Macrophytes communities, (seaweeds and seagrass), some of them are protected, consist in a key indicator of the ecological health of lagoons, as stated by the European Water Framework Directive. Furthermore, the restoration of these communities, especially the Zostera meadows is one of the main objectives of the Berre lagoon restoration plan, leaded by the Gipreb (Groupement d’intérêt public pour la réhabilitation de l’Etang de Berre). Consequently, this organism is monitoring since 1996 the main macrophytes species still present in the lagoon, in order to evaluate the impact of freshwater releases from an electric power plant on this ecosystem.  A sequence of eleven annual tables of density of Zostera noltii is studied through the Ordinal Correspondence Analysis of Beh (1997, 2004), whose principles are first recalled. The obtained results are compared with those of classical Correspondence Analysis. Keywords : Data analysis, Environment Introduction L’étang de Berre, une des plus grandes lagunes méditerranénnes, est soumise depuis 1966 à des apports massifs d’eau douce et de limons par la centrale hydroélectrique de St-Chamas. Il en a résulté des modifications brutales de la salinité et de la transparence des eaux de surface, ce qui a conduit à la régression des communautés de macrophytes de l’étang. Ces peuplements sont suivis depuis 1996 par le Gipreb: la localisation et l’abondance de 8 espèces ou groupes d’espèces de macrophytes sont ainsi annuellement relevés (en période estivale, période de plus forte végétation) le long du littoral de l’étang. La densité de chaque espèce ou groupe d’espèces a été évaluée visuellement par des plongeurs le long de 31 transects (composés chacun de 20 segments de 5 m) répartis régulièrement sur l’ensemble du pourtour de l’étang, et ce pendant 11 années (1996 à 2008). L’évaluation de la densité consiste en un codage en 6 modalités (cf. figure 4) du taux de couverture du fonds par l’espèce. L’état annuel de la population de chaque espèce est donc décrit par un tableau de type 6x20, croisant la densité avec la position : deux variables ordinales. Il s’agit pour nous de tester aussi rigoureusement que possible l’existence d’effets spatio-temporels pour chaque macrophyte. Nous tenterons donc, pour l’espèce Zostera noltii (espèce reconnue d’importance patrimoniale, protégée et qui par sa présence indique un bon fonctionnement écologique du système lagunaire), d’apporter une réponse aux questions suivantes : 1. la densité de l’espèce dépend-elle de la distance à la côte ? 2. si oui, quelles sont les zones du transect importantes pour sa dynamique? 3. y a-t-il eu un changement à partir de 2005 (année n°8 et suivantes), année à compter de laquelle des nouvelles régulations des apports d’eau douce par la centrale hydroélectrique de Saint-Chamas ont été appliquées ? Une manière classique d’analyser de ce genre de données est de procéder à l’AFC du tableau T de type 206X  associé à chaque espèce, obtenu en sommant les onze tableaux annuels. Cette méthode (Benzécri, 1967 ; Beh, 2004) permet d’analyser un tableau de contingence T à I lignes et J colonnes, d’effectif total n  (ici, I=6, J=20, et n=11x31x20=6820). On note nTP /:=  la probabilité associée, et ∑==Jj jii PP 1 ,:  (resp. ∑ ==Ii jij PP 1 ,:  ) le poids de la ième ligne (resp. de la jème colonne). L’AFC met en évidence la manière dont la relation entre lignes et colonnes de T  s’écarte de l’indépendance, en décomposant le tableau Θ  de terme généraljijiji PPP×=,,:θ . Dans la pratique, ( )2,,2 1: ∑ −×=Χ ji jiji PPn θ  suit (sous hypothèse d’indépendance) une loi ( ))1()1(2 −− JIχ , et l’AFC consiste essentiellement à déterminer des couples d’axes factoriels ( ) MmBA mm ≤≤1:, , avec ]1,1[: −−= JIMinM , tels que l’on ait : jmmMmimji BA ,0,,λθ ∑== . On a alors les relations fondamentales : ∑==ΧMmmn122 λ , et : ∑ ×∈= JIji jpjiimpmm BPA),( ,,,δλ . La méthode de Beh pour les tableaux ordonnés Mais les facteurs d’une AFC ne sont pas toujours faciles à interpréter, et surtout leur construction ne dépend pas de l’ordre des modalités. C’est pourquoi Beh (1997, 2004) a proposé une méthode alternative de décomposition de 2Χ , adaptée au cas ou les lignes et/ou les colonnes du tableau traité sont ordonnées. 2Χ  n’est plus alors décomposée en fonction d’axes factoriels dépendant de l’ensemble du tableau analysé, mais en fonction de polynômes dépendant des variables ordinales structurant le tableau au travers de scores, et du tableau au travers de ses marges { }IiPi ≤≤1:  et { }JjPj ≤≤1:  uniquement. Ces polynômes, ici notés mA  pour les lignes (densité) et pB  pour les colonnes (segments), jouent le rôle de facteurs. Comme des axes factoriels, ils sont centrés et orthonormés. Par contre, ce ne sont pas des vecteurs singuliers de Θ , ce qui fait qu’il est impossible de parler de paires d’axes associés à une part de variance commune, comme c’est le cas en AFC. Pour la même raison, des polynômes de degrés différents peuvent être corrélés. Si l’on note µ  (resp.η ) les parts de variance associés aux polynômes sur la position (resp. la densité), nous savons seulement que : ∑∑−=−===Χ 11112Jp pIm mnµη . On a donc ici deux autres manières de décomposer 2Χ , distinctes de celle de l’AFC.  Les tests Il est bien connu qu’il n’existe pas de test rigoureux de significativité des axes factoriels en AFC, en raison de la complexité de la loi des valeurs propres associées. A notre connaissance, il n’existe guère que la procédure classique basée sur le cumul des valeurs propres proposée par Benzécri (1976). En revanche, dans le cas de la méthode de Beh, il est possible de tester l’effet des polynômes et de leurs interactions, quantifiées par des moments bivariés. Le moment bivarié d’ordre (m,p) est : (Beh, 2004) ∑ ×∈= JIji jpjiimpm BPA),( ,,,, :λ . Remarquons que cette formule est formellement identique à celle donnant la corrélation entre les axes factoriels de l’AFC. Mais ici les axes (polynômes sur les deux ensembles) ne sont plus décorrélés : toutes les interactions sont possibles. Rayner et Best (1996) ont proposé des tests pour ces moments, en démontrant qu’ils sont asymptotiquement indépendants et de loi )1,0(N , sous hypothèse d’indépendance des lignes et des colonnes. On en déduit une statistique permettant sous cette hypothèse de tester la significativité du polynôme 11, −≤≤∀ ImAm  : ∑−=−≈=112)1(2,Jp Jpmmχλη . On a de même pour les colonnes la statistique : ∑−=−≈=−≤≤∀ 112)1(2,,11 Im IpmpJp χλµ . Un exemple : l’analyse des transects de Zostères Zostera noltii est une espèce qualifiée dans l’étang de Berre de “pseudo-annuelle”, les herbiers qu’elle constitue présentent dont un caractére pérenne d’une année à l’autre, en terme au moins de localisation géographique et d’extension. Le recouvrement sur le fond des faisceaux de feuilles, en revanche, connait des fluctuations saisonnières et interannuelles, qui témoignent de la vitalité et de la dynamique de l’espèce d’une année à l’autre. Dans l’étang de Berre, son extension vers le large semble contrainte par la profondeur, directement en lien avec la lumière disponible sur le fond (Bernard et al., 2007). Les zostères étant totalement absentes des 10 derniers segments des transects, nous les avons supprimés de cette analyse, et donc exceptionnellement J=10. Le test du Chi-deux global ( )( )422 )45( 10.312.84 −≈=> XP χ  montre que la densité et la distance au rivage sont en relation forte, et le test sur la trace ne donne pour significatif que le 1er facteur de l’AFC (72.6% de l’inertie), le 2ème correspondant à un taux d’inertie de 12.%.  0-505-1010-1515-2020-2525-3030-3535-4040-4545-50123456-0.5 0.0 0.5 1.0 1.5 2.0-1.5-1.0-0.50.00.51.01.5Facteur 1Facteur 2 Figure 1 : 1er plan de l'Analyse des Correspondances classique  Il apparaît clairement sur ce plan que la densité « 1 » (nulle) est très fréquente (proche de la moyenne), alors que les fortes densités de Zostères sont associées aux segment 5-10m et surtout 10-15m. Les scores et les polynômes orthonormés On remarque sur la figure 2 ci-dessous que le poids des segments, contrairement à celui des niveaux, est constant. En effet, chaque colonne de chaque tableau à toujours pour effectif le nombre de stations, 31. ææ æ æ æ æ2 3 4 5 60.20.40.60.8æ æ æ æ æ æ æ æ æ æ2 4 6 8 100.050.100.150.20 Figure 2 : poids des niveaux de densité et des segments  Beh (1998) quantifie la différence entre les modalités consécutive d’une des variable ordinales au moyen d’un score S. Le score « naturel » est jjS =:)( , où la distance entre modalités est constante. C’est celui que nous avons choisi pour coder la position le long du transect, car les segments considérés sont de même largeur. Les polynômes orthonormés sur les transects, { }91 ,, BB L  ont donc une allure parfaitement symétrique : ce sont (en gros) des polynômes de Legendre discrétisés et changés d’échelles. La situation est différente pour la variable “densité”. D’abord, notons (voir ci-dessus) que le poids de la première modalité est proche de 1 (0.96) : très fréquemment, les Zostères sont absentes. Pour ce qui est du score, nous avons associé à chaque classe le taux de recouvrement central évalué par le plongeur (voir Bonhomme et al. (2008)), ce qui donne la figure 3 ci-dessous. La non-linéarité du score, et le poids considérable de la modalité « 1 », font que les polynôme relatifs à cette variable ne sont plus du tout symétriques.  ææææææ1 2 3 4 5 620406080 Figure 3 : score associé au taux de recouvrement  Remarquons maintenant que sur la figure 4 les taux d’inertie de cette analyse ne sont pas systématiquement décroissants, contrairement à ce qui se produit en AFC. Ces graphiques montrent que les variations de structure des transects sont expliquées le mieux possible par{ }51 , AA , alors que les variations de densité sont expliquées par les polynômes{ }431 ,, BBB . CCC CCCCC C2 4 6 8 105101520Contribution des Bp à l'inertie des lignesLL L LL2 4 6 8 101020304050Contribution des Am à l'inertie des colonnesTests du c2 sur les composantes de Beh HRouge: 0.95, Vert: médianeL Figure 4 : épreuve de validité des composantes de Beh  Les tests portant sur les interactions permettent davantage de précision. Ils sont résumés par la figure 5 ci-dessous, représentant la matrice Λ  des moments bivariés. La case de coordonnées (m,p) est blanche lorsque l’hypothèse : )1,0(,Npm ≈λ  n’est pas rejetée au seuil 0.05 (test bilatère) ; sinon elle est d’autant plus rouge que l’interaction est statistiquement significative. 1 2 3 4 5 6 7 8 9123451 2 3 4 5 6 7 8 912345Degré&DistanceDegré&Densité Figure 5 :  test de significativité des éléments de la matrice Λ .  D’après la figure 4, le plan factoriel le plus pertinent (74% de la variance) pour projeter les segments est le plan (1,5), représenté sur la figure 6. Sur ce plan, les segments sont représentés en abscisse (resp. en ordonnée) par la composante de Beh (1997) )1(*G  (resp. )5(*G ) dont la jième coordonnée est la projection du jième segment sur 1A  (resp. 5A ).   0-505-1010-1515-2020-2525-3030-3535-4040-4545-50-0.2 -0.1 0.0 0.1-0.10-0.050.000.05A 1A5 Figure 6 : plan (1,5) de l'analyse ordinale de Beh.  Enfin, il est aussi possible d’introduire le temps dans l’analyse, en remarquant que ∑==111aaa PP α , où aP  désigne la probabilité associée au tableau partiel des mesures de l’année « a », et 11/1=aα  (les tableaux partiels ont ici tous la même masse). Il s’ensuit que la position de chaque moyenne temporelle est au barycentre des points annuels correspondant à ce segment. Cela permet de tracer une ultime figure (non représentée, faute de place) mettant en évidence le(s) segment(s) où la dynamique temporelle est la plus agitée, et une éventuelle rupture en 2005 (année 8 et suivantes). +Conclusions : L’analyse en correspondance ordinale s’avère être un outil intéressant pour le suivi à long terme des communautés littorales (ici les macrophytes). Appliquée aux herbiers de zostères (Z. Noltii), elle permet de mettre clairement en évidence la dépendance de l’abondance des Zostères en un point avec sa distance au rivage. Cela confirme la contrainte liée à la profondeur (et donc à l’éclairement) pour leur développement et leur extension vers le large. En première approximation, leur densité décroît le long du transect (car 71.31,1 −=λ ). Cependant, comme on peut le voir sur la figure 5, il existe aussi des interactions de degrés 3 et 4 significatives. Les segments où la variabilité temporelle des Zostères est maximale sont ceux correspondant à une distance au rivage de 0-25m, surtout le segment 10-15m (très proches de la côte). En terme d’évolution temporelle, d’après l’interprétation de la 1ère composante, la régression de l’espèce observée depuis 1996 se poursuit en 2005, et après. Toutefois, on ne peut constater de changement clair que sur le segment 10-15m.  Bibliographie   [1] Beh, E. J. (1997) Simple Correspondence Analysis of ordinal cross-classifications using orthogonal polynomials, Biometrical Journal, 39, 5, 589-613. [2] Beh, E. J. (1998) A comparative study of Scores for Correspondence Analysis with ordered categories, Biometrical Journal, 40, 4, 413-429. [3] Beh, E. J. (2004) Simple Correspondence Analysis : a bibliographic review, International Statistical Review, 72, 2, 257-284. [4] Benzécri, J. P. (1976) L’Analyse des Données, tome 2, 2ème ed., Dunod. [5] Bernard G., Boudouresque CF., Picon P. (2007) Long term changes in Zostera meadows in the Berre lagoon (Provence, Mediterranean Sea). Estuarine, Coastal and Shelf Science, 73:617-629. [5] Bonhomme, P., Gravel, R. Ruitton, S., Bonhomme, Boudouresque, C.-F. (2008) Système de surveillance des principales espèces de macrophytes dans les étangs de Berre et de Vaïne : Suivi 2007, Contrat GIS Posidonie/GIPREB. GIS Posidonie publ., Marseille, Fr. : 1- 47+ Annexes 1-3. [6] Rayner, J.C.W. and Best, D.J. (1996) Smooth extensions of Pearson’s product moment correlation and Spearman’s Rho, Statistics & Probability Letters, 30, 171-177.  

Application de l'Analyse des Correspondances Ordinales au suivi d'espèces végétales aquatiques.

Abstract

Similar works

Full text

Available Versions

HAL AMU