Conservation of the Eurasian Lynx (Lynx lynx) in a fragmented landscape - habitat models, dispersal and potential distribution

Abstract: The expansion of a recovering population - whether re-introduced or spontaneously returning - is shaped by (i) biological (intrinsic) factors such as the land tenure system or dispersal, (ii) the distribution and availability of resources (e.g. prey), (iii) habitat and landscape features, and (iv) human attitudes and activities. In order to develop efficient conservation and recovery strategies, we need to understand all these factors and to predict the potential distribution and explore ways to reach it. An increased number of lynx in the north-western Swiss Alps in the nineties lead to a new controversy about the return of this cat. When the large carnivores were given legal protection in many European countries, most organizations and individuals promoting their protection did not foresee the consequences. Management plans describing how to handle conflicts with large predators are needed to find a balance between "overabundance" and extinction. Wildlife and conservation biologists need to evaluate the various threats confronting populations so that adequate management decisions can be taken. I developed a GIS probability model for the lynx, based on habitat information and radio-telemetry data from the Swiss Jura Mountains, in order to predict the potential distribution of the lynx in this mountain range, which is presently only partly occupied by lynx. Three of the 18 variables tested for each square kilometre describing land use, vegetation, and topography, qualified to predict the probability of lynx presence. The resulting map was evaluated with data from dispersing subadult lynx. Young lynx that were not able to establish home ranges in what was identified as good lynx habitat did not survive their first year of independence, whereas the only one that died in good lynx habitat was illegally killed. Radio-telemetry fixes are often used as input data to calibrate habitat models. Radio-telemetry is the only way to gather accurate and unbiased data on habitat use of elusive larger terrestrial mammals. However, it is time consuming and expensive, and can therefore only be applied in limited areas. Habitat models extrapolated over large areas can in turn be problematic, as habitat characteristics and availability may change from one area to the other. I analysed the predictive power of Ecological Niche Factor Analysis (ENFA) in Switzerland with the lynx as focal species. According to my results, the optimal sampling strategy to predict species distribution in an Alpine area lacking available data would be to pool presence cells from contrasted regions (Jura Mountains, Alps), whereas in regions with a low ecological variance (Jura Mountains), only local presence cells should be used for the calibration of the model. Dispersal influences the dynamics and persistence of populations, the distribution and abundance of species, and gives the communities and ecosystems their characteristic texture in space and time. Between 1988 and 2001, the spatio-temporal behaviour of subadult Eurasian lynx in two re-introduced populations in Switzerland was studied, based on 39 juvenile lynx of which 24 were radio-tagged to understand the factors influencing dispersal. Subadults become independent from their mothers at the age of 8-11 months. No sex bias neither in the dispersal rate nor in the distance moved was detected. Lynx are conservative dispersers, compared to bear and wolf, and settled within or close to known lynx occurrences. Dispersal distances reached in the high lynx density population - shorter than those reported in other Eurasian lynx studies - are limited by habitat restriction hindering connections with neighbouring metapopulations. I postulated that high lynx density would lead to an expansion of the population and validated my predictions with data from the north-western Swiss Alps where about 1995 a strong increase in lynx abundance took place. The general hypothesis that high population density will foster the expansion of the population was not confirmed. This has consequences for the re-introduction and recovery of carnivores in a fragmented landscape. To establish a strong source population in one place might not be an optimal strategy. Rather, population nuclei should be founded in several neighbouring patches. Exchange between established neighbouring subpopulations will later on take place, as adult lynx show a higher propensity to cross barriers than subadults. To estimate the potential population size of the lynx in the Jura Mountains and to assess possible corridors between this population and adjacent areas, I adapted a habitat probability model for lynx distribution in the Jura Mountains with new environmental data and extrapolated it over the entire mountain range. The model predicts a breeding population ranging from 74-101 individuals and from 51-79 individuals when continuous habitat patches < 50 km2 are disregarded. The Jura Mountains could once be part of a metapopulation, as potential corridors exist to the adjoining areas (Alps, Vosges Mountains, and Black Forest). Monitoring of the population size, spatial expansion, and the genetic surveillance in the Jura Mountains must be continued, as the status of the population is still critical. ENFA was used to predict the potential distribution of lynx in the Alps. The resulting model divided the Alps into 37 suitable habitat patches ranging from 50 to 18,711 km2, covering a total area of about 93,600 km2. When using the range of lynx densities found in field studies in Switzerland, the Alps could host a population of 961 to 1,827 residents. The results of the cost-distance analysis revealed that all patches were within the reach of dispersing lynx, as the connection costs were in the range of dispersal cost of radio-tagged subadult lynx moving through unfavorable habitat. Thus, the whole Alps could once be considered as a metapopulation. But experience suggests that only few disperser will cross unsuitable areas and barriers. This low migration rate may seldom allow the spontaneous foundation of new populations in unsettled areas. As an alternative to natural dispersal, artificial transfer of individuals across the barriers should be considered. Wildlife biologists can play a crucial role in developing adaptive management experiments to help managers learning by trial. The case of the lynx in Switzerland is a good example of a fruitful cooperation between wildlife biologists, managers, decision makers and politician in an adaptive management process. This cooperation resulted in a Lynx Management Plan which was implemented in 2000 and updated in 2004 to give the cantons directives on how to handle lynx-related problems. This plan was put into practice e.g. in regard to translocation of lynx into unsettled areas. Résumé: L'expansion d'une population en phase de recolonisation, qu'elle soit issue de réintroductions ou d'un retour naturel dépend 1) de facteurs biologiques tels que le système social et le mode de dispersion, 2) de la distribution et la disponibilité des ressources (proies), 3) de l'habitat et des éléments du paysage, 4) de l'acceptation de l'espèce par la population locale et des activités humaines. Afin de pouvoir développer des stratégies efficaces de conservation et de favoriser la recolonisation, chacun de ces facteurs doit être pris en compte. En plus, la distribution potentielle de l'espèce doit pouvoir être déterminée et enfin, toutes les possibilités pour atteindre les objectifs, examinées. La phase de haute densité que la population de lynx a connue dans les années nonante dans le nord-ouest des Alpes suisses a donné lieu à une controverse assez vive. La protection du lynx dans de nombreux pays européens, promue par différentes organisations, a entraîné des conséquences inattendues; ces dernières montrent que tout plan de gestion doit impérativement indiquer des pistes quant à la manière de gérer les conflits, tout en trouvant un équilibre entre l'extinction et la surabondance de l'espèce. Les biologistes de la conservation et de la faune sauvage doivent pour cela évaluer les différents risques encourus par les populations de lynx, afin de pouvoir rapidement prendre les meilleuresmdécisions de gestion. Un modèle d'habitat pour le lynx, basé sur des caractéristiques de l'habitat et des données radio télémétriques collectées dans la chaîne du Jura, a été élaboré afin de prédire la distribution potentielle dans cette région, qui n'est que partiellement occupée par l'espèce. Trois des 18 variables testées, décrivant pour chaque kilomètre carré l'utilisation du sol, la végétation ainsi que la topographie, ont été retenues pour déterminer la probabilité de présence du lynx. La carte qui en résulte a été comparée aux données télémétriques de lynx subadultes en phase de dispersion. Les jeunes qui n'ont pas pu établir leur domaine vital dans l'habitat favorable prédit par le modèle n'ont pas survécu leur première année d'indépendance alors que le seul individu qui est mort dans l'habitat favorable a été braconné. Les données radio-télémétriques sont souvent utilisées pour l'étalonnage de modèles d'habitat. C'est un des seuls moyens à disposition qui permette de récolter des données non biaisées et précises sur l'occupation de l'habitat par des mammifères terrestres aux moeurs discrètes. Mais ces méthodes de- mandent un important investissement en moyens financiers et en temps et peuvent, de ce fait, n'être appliquées qu'à des zones limitées. Les modèles d'habitat sont ainsi souvent extrapolés à de grandes surfaces malgré le risque d'imprécision, qui résulte des variations des caractéristiques et de la disponibilité de l'habitat d'une zone à l'autre. Le pouvoir de prédiction de l'Analyse Ecologique de la Niche (AEN) dans les zones où les données de présence n'ont pas été prises en compte dans le calibrage du modèle a été analysée dans le cas du lynx en Suisse. D'après les résultats obtenus, la meilleure mé- thode pour prédire la distribution du lynx dans une zone alpine dépourvue d'indices de présence est de combiner des données provenant de régions contrastées (Alpes, Jura). Par contre, seules les données sur la présence locale de l'espèce doivent être utilisées pour les zones présentant une faible variance écologique tel que le Jura. La dispersion influence la dynamique et la stabilité des populations, la distribution et l'abondance des espèces et détermine les caractéristiques spatiales et temporelles des communautés vivantes et des écosystèmes. Entre 1988 et 2001, le comportement spatio-temporel de lynx eurasiens subadultes de deux populations réintroduites en Suisse a été étudié, basé sur le suivi de 39 individus juvéniles dont 24 étaient munis d'un collier émetteur, afin de déterminer les facteurs qui influencent la dispersion. Les subadultes se sont séparés de leur mère à l'âge de 8 à 11 mois. Le sexe n'a pas eu d'influence sur le nombre d'individus ayant dispersés et la distance parcourue au cours de la dispersion. Comparé à l'ours et au loup, le lynx reste très modéré dans ses mouvements de dispersion. Tous les individus ayant dispersés se sont établis à proximité ou dans des zones déjà occupées par des lynx. Les distances parcourues lors de la dispersion ont été plus courtes pour la population en phase de haute densité que celles relevées par les autres études de dispersion du lynx eurasien. Les zones d'habitat peu favorables et les barrières qui interrompent la connectivité entre les populations sont les principales entraves aux déplacements, lors de la dispersion. Dans un premier temps, nous avons fait l'hypothèse que les phases de haute densité favorisaient l'expansion des populations. Mais cette hypothèse a été infirmée par les résultats issus du suivi des lynx réalisé dans le nord-ouest des Alpes, où la population connaissait une phase de haute densité depuis 1995. Ce constat est important pour la conservation d'une population de carnivores dans un habitat fragmenté. Ainsi, instaurer une forte population source à un seul endroit n'est pas forcément la stratégie la plus judicieuse. Il est préférable d'établir des noyaux de populations dans des régions voisines où l'habitat est favorable. Des échanges entre des populations avoisinantes pourront avoir lieu par la suite car les lynx adultes sont plus enclins à franchir les barrières qui entravent leurs déplacements que les individus subadultes. Afin d'estimer la taille de la population de lynx dans le Jura et de déterminer les corridors potentiels entre cette région et les zones avoisinantes, un modèle d'habitat a été utilisé, basé sur un nouveau jeu de variables environnementales et extrapolé à l'ensemble du Jura. Le modèle prédit une population reproductrice de 74 à 101 individus et de 51 à 79 individus lorsque les surfaces d'habitat d'un seul tenant de moins de 50 km2 sont soustraites. Comme des corridors potentiels existent effectivement entre le Jura et les régions avoisinantes (Alpes, Vosges, et Forêt Noire), le Jura pourrait faire partie à l'avenir d'une métapopulation, lorsque les zones avoisinantes seront colonisées par l'espèce. La surveillance de la taille de la population, de son expansion spatiale et de sa structure génétique doit être maintenue car le statut de cette population est encore critique. L'AEN a également été utilisée pour prédire l'habitat favorable du lynx dans les Alpes. Le modèle qui en résulte divise les Alpes en 37 sous-unités d'habitat favorable dont la surface varie de 50 à 18'711 km2, pour une superficie totale de 93'600 km2. En utilisant le spectre des densités observées dans les études radio-télémétriques effectuées en Suisse, les Alpes pourraient accueillir une population de lynx résidents variant de 961 à 1'827 individus. Les résultats des analyses de connectivité montrent que les sous-unités d'habitat favorable se situent à des distances telles que le coût de la dispersion pour l'espèce est admissible. L'ensemble des Alpes pourrait donc un jour former une métapopulation. Mais l'expérience montre que très peu d'individus traverseront des habitats peu favorables et des barrières au cours de leur dispersion. Ce faible taux de migration rendra difficile toute nouvelle implantation de populations dans des zones inoccupées. Une solution alternative existe cependant : transférer artificiellement des individus d'une zone à l'autre. Les biologistes spécialistes de la faune sauvage peuvent jouer un rôle important et complémentaire pour les gestionnaires de la faune, en les aidant à mener des expériences de gestion par essai. Le cas du lynx en Suisse est un bel exemple d'une collaboration fructueuse entre biologistes de la faune sauvage, gestionnaires, organes décisionnaires et politiciens. Cette coopération a permis l'élaboration du Concept Lynx Suisse qui est entré en vigueur en 2000 et remis à jour en 2004. Ce plan donne des directives aux cantons pour appréhender la problématique du lynx. Il y a déjà eu des applications concrètes sur le terrain, notamment par des translocations d'individus dans des zones encore inoccupées.

Rockfall susceptibility assessment and remote geological mapping with LiDAR point clouds

Characterizing the geological features and structures in three dimensions over inaccessible rock cliffs is needed to assess natural hazards such as rockfalls and rockslides and also to perform investigations aimed at mapping geological contacts and building stratigraphy and fold models. Indeed, the detailed 3D data, such as LiDAR point clouds, allow to study accurately the hazard processes and the structure of geologic features, in particular in vertical and overhanging rock slopes. Thus, 3D geological models have a great potential of being applied to a wide range of geological investigations both in research and applied geology projects, such as mines, tunnels and reservoirs. Recent development of ground-based remote sensing techniques (LiDAR, photogrammetry and multispectral / hyperspectral images) are revolutionizing the acquisition of morphological and geological information. As a consequence, there is a great potential for improving the modeling of geological bodies as well as failure mechanisms and stability conditions by integrating detailed remote data. During the past ten years several large rockfall events occurred along important transportation corridors where millions of people travel every year (Switzerland: Gotthard motorway and railway; Canada: Sea to sky highway between Vancouver and Whistler). These events show that there is still a lack of knowledge concerning the detection of potential rockfalls, making mountain residential settlements and roads highly risky. It is necessary to understand the main factors that destabilize rocky outcrops even if inventories are lacking and if no clear morphological evidences of rockfall activity are observed. In order to increase the possibilities of forecasting potential future landslides, it is crucial to understand the evolution of rock slope stability. Defining the areas theoretically most prone to rockfalls can be particularly useful to simulate trajectory profiles and to generate hazard maps, which are the basis for land use planning in mountainous regions. The most important questions to address in order to assess rockfall hazard are: Where are the most probable sources for future rockfalls located? What are the frequencies of occurrence of these rockfalls? I characterized the fracturing patterns in the field and with LiDAR point clouds. Afterwards, I developed a model to compute the failure mechanisms on terrestrial point clouds in order to assess the susceptibility to rockfalls at the cliff scale. Similar procedures were already available to evaluate the susceptibility to rockfalls based on aerial digital elevation models. This new model gives the possibility to detect the most susceptible rockfall sources with unprecented detail in the vertical and overhanging areas. The results of the computation of the most probable rockfall source areas in granitic cliffs of Yosemite Valley and Mont-Blanc massif were then compared to the inventoried rockfall events to validate the calculation methods. Yosemite Valley was chosen as a test area because it has a particularly strong rockfall activity (about one rockfall every week) which leads to a high rockfall hazard. The west face of the Dru was also chosen for the relevant rockfall activity and especially because it was affected by some of the largest rockfalls that occurred in the Alps during the last 10 years. Moreover, both areas were suitable because of their huge vertical and overhanging cliffs that are difficult to study with classical methods. Limit equilibrium models have been applied to several case studies to evaluate the effects of different parameters on the stability of rockslope areas. The impact of the degradation of rockbridges on the stability of large compartments in the west face of the Dru was assessed using finite element modeling. In particular I conducted a back-analysis of the large rockfall event of 2005 (265'000 m3) by integrating field observations of joint conditions, characteristics of fracturing pattern and results of geomechanical tests on the intact rock. These analyses improved our understanding of the factors that influence the stability of rock compartments and were used to define the most probable future rockfall volumes at the Dru. Terrestrial laser scanning point clouds were also successfully employed to perform geological mapping in 3D, using the intensity of the backscattered signal. Another technique to obtain vertical geological maps is combining triangulated TLS mesh with 2D geological maps. At El Capitan (Yosemite Valley) we built a georeferenced vertical map of the main plutonio rocks that was used to investigate the reasons for preferential rockwall retreat rate. Additional efforts to characterize the erosion rate were made at Monte Generoso (Ticino, southern Switzerland) where I attempted to improve the estimation of long term erosion by taking into account also the volumes of the unstable rock compartments. Eventually, the following points summarize the main out puts of my research: The new model to compute the failure mechanisms and the rockfall susceptibility with 3D point clouds allows to define accurately the most probable rockfall source areas at the cliff scale. The analysis of the rockbridges at the Dru shows the potential of integrating detailed measurements of the fractures in geomechanical models of rockmass stability. The correction of the LiDAR intensity signal gives the possibility to classify a point cloud according to the rock type and then use this information to model complex geologic structures. The integration of these results, on rockmass fracturing and composition, with existing methods can improve rockfall hazard assessments and enhance the interpretation of the evolution of steep rockslopes. -- La caractérisation de la géologie en 3D pour des parois rocheuses inaccessibles est une étape nécessaire pour évaluer les dangers naturels tels que chutes de blocs et glissements rocheux, mais aussi pour réaliser des modèles stratigraphiques ou de structures plissées. Les modèles géologiques 3D ont un grand potentiel pour être appliqués dans une vaste gamme de travaux géologiques dans le domaine de la recherche, mais aussi dans des projets appliqués comme les mines, les tunnels ou les réservoirs. Les développements récents des outils de télédétection terrestre (LiDAR, photogrammétrie et imagerie multispectrale / hyperspectrale) sont en train de révolutionner l'acquisition d'informations géomorphologiques et géologiques. Par conséquence, il y a un grand potentiel d'amélioration pour la modélisation d'objets géologiques, ainsi que des mécanismes de rupture et des conditions de stabilité, en intégrant des données détaillées acquises à distance. Pour augmenter les possibilités de prévoir les éboulements futurs, il est fondamental de comprendre l'évolution actuelle de la stabilité des parois rocheuses. Définir les zones qui sont théoriquement plus propices aux chutes de blocs peut être très utile pour simuler les trajectoires de propagation des blocs et pour réaliser des cartes de danger, qui constituent la base de l'aménagement du territoire dans les régions de montagne. Les questions plus importantes à résoudre pour estimer le danger de chutes de blocs sont : Où se situent les sources plus probables pour les chutes de blocs et éboulement futurs ? Avec quelle fréquence vont se produire ces événements ? Donc, j'ai caractérisé les réseaux de fractures sur le terrain et avec des nuages de points LiDAR. Ensuite, j'ai développé un modèle pour calculer les mécanismes de rupture directement sur les nuages de points pour pouvoir évaluer la susceptibilité au déclenchement de chutes de blocs à l'échelle de la paroi. Les zones sources de chutes de blocs les plus probables dans les parois granitiques de la vallée de Yosemite et du massif du Mont-Blanc ont été calculées et ensuite comparés aux inventaires des événements pour vérifier les méthodes. Des modèles d'équilibre limite ont été appliqués à plusieurs cas d'études pour évaluer les effets de différents paramètres sur la stabilité des parois. L'impact de la dégradation des ponts rocheux sur la stabilité de grands compartiments de roche dans la paroi ouest du Petit Dru a été évalué en utilisant la modélisation par éléments finis. En particulier j'ai analysé le grand éboulement de 2005 (265'000 m3), qui a emporté l'entier du pilier sud-ouest. Dans le modèle j'ai intégré des observations des conditions des joints, les caractéristiques du réseau de fractures et les résultats de tests géoméchaniques sur la roche intacte. Ces analyses ont amélioré l'estimation des paramètres qui influencent la stabilité des compartiments rocheux et ont servi pour définir des volumes probables pour des éboulements futurs. Les nuages de points obtenus avec le scanner laser terrestre ont été utilisés avec succès aussi pour produire des cartes géologiques en 3D, en utilisant l'intensité du signal réfléchi. Une autre technique pour obtenir des cartes géologiques des zones verticales consiste à combiner un maillage LiDAR avec une carte géologique en 2D. A El Capitan (Yosemite Valley) nous avons pu géoréferencer une carte verticale des principales roches plutoniques que j'ai utilisé ensuite pour étudier les raisons d'une érosion préférentielle de certaines zones de la paroi. D'autres efforts pour quantifier le taux d'érosion ont été effectués au Monte Generoso (Ticino, Suisse) où j'ai essayé d'améliorer l'estimation de l'érosion au long terme en prenant en compte les volumes des compartiments rocheux instables. L'intégration de ces résultats, sur la fracturation et la composition de l'amas rocheux, avec les méthodes existantes permet d'améliorer la prise en compte de l'aléa chute de pierres et éboulements et augmente les possibilités d'interprétation de l'évolution des parois rocheuses.

Mélanges historiques.

Contient : Traité sur la charge d'ambassadeur ; « De la prélation et retenue feudale, pour messire Henry de La Tour, duc de Bouillon, premier mareschal de France, vicomte de Turenne, demandeur, contre le sieur de Noailles, deffendeur » [par Nicolas Rigault], et autres pièces, inventaire de productions, etc., concernant le procès Bouillon-Noailles ; Pièces concernant les biens et dépendances de l'abbaye de Saint-Wandrille ; « Très humble et très importante remonstrance au Roy [Louis XIV] sur la remise des places maritimes de Flandres entre les mains des Anglois », au temps de Cromwell ; « Journal de ce qui s'est passé à la conférence pour la paix..., à Munster..., l'an 1644 », et autres pièces concernant les négociations de Münster, etc., 1639-1644 ; Pièces diplomatiques, concernant diverses questions de préséance ; Extrait d'une lettre de Charles Brulart, prieur de Léon, à Villeroy, Venise, 1611 ; « Estat par estimation des train et suitte, livrées, équipage et ameublement d'ambassadeur à Rome », etc

Recueil de lettres et de pièces originales, et de copies de pièces indiquées comme telles dans le dépouillement qui suit.

Contient : 1 Lettre de « HENRY [duc D'ANJOU]... à mon cousin monseigneur le duc de Nemoux,... Escript au camp à St Mesmyn, ce VIIe jour de fevrier » ; 2 Lettre de « HENRY [duc D'ANJOU]... à mon cousin monseigneur de Nemoux,... Escript au camp à Pons sur Senne, ce VIIe jour de fevrier » ; 3 Lettre de « HENRY [duc D'ANJOU]... à mon cousin monseigneur le duc de Nemoux,... Escript au camp à Nogen sur Senne, ce Xme jour de fevrier » ; 4 Lettre de « HENRY [duc D'ANJOU]... à mon cousin monseigneur le duc de Nemoux,... Escript au camp à Nogen sur Senne, ce VIIIme jour de fevrier » ; 5 Lettre de « JEHANNE » D'ALBRET, reine de Navarre ; 6 Lettre de « MARGUERITE [reine de Navarre]... à mon cousin monsieur le duc de Nevers,... Du 16 may 1573 » ; 7 Lettre d'ARMAND DE GONTAUT, baron DE « BIRON,... à monsieur... de Cossé, conte de Segondiny et mareschal de France... De La Rochelle, se dimanche matin » ; 8 Lettre du roi « CHARLES [IX]... à mon cousin le duc de Montpansier, gouverneur... en Bretaigne... Escript à St Leger, le IIIIe jour de mars 1573 » ; 9 Lettre de « CATERINE [DE MEDICIS]... à mon cousin monseigneur le duc de Montpensier,... gouverneur et lieutenant general du roy... en son pays et duché de Bretaigne... Escript à Paris, le XIIe jour de fevrier 1573 » ; 10 Lettre de « HENRY [duc D'ANJOU]... à mon cousin monseigneur le duc de Montpensier,... gouverneur... en Bretaigne... Escrit au camp de Nyeul, ce VIIIe may 1573 » ; 11 Lettre du roi « CHARLES [IX]... à mon cousin le duc de Montpensier,... gouverneur... en Bretaigne... Escript à St Liger, le IIe jour de mars 1573 » ; 12 Lettre du roi « CHARLES [IX]... à mon cousin le duc de Montpensier,... gouverneur... en mon pays et duché de Bretaigne... Escript à Fontainebleau, ce XXVe jour d'avril 1573 » ; 13 Lettre de « CATERINE [DE MEDICIS]... à mon cousin monseigneur le duc de Montpensier,... gouverneur... en Bretaigne... Escript à Fontainebleau, le dernier jour de may 1573 » ; 14 Lettre du roi « CHARLES [IX]... à mon cousin le duc de Montpensier,... gouverneur... en Bretaigne... Escrit à Fontainebleau,... le IIe jour de may 1573 » ; 15 Lettre de « CATERINE [DE MEDICIS]... à mon cousin le duc de Montpensier,... gouverneur... en Bretaigne... Escrit à Paris, le XVIe jour de fevrier 1573 » ; 16 Lettre de « CATERINE [DE MEDICIS]... à mon cousin monseigneur le duc de Montpensier,... gouverneur... en Bretaigne... Escript à St Liger, le IIe jour de mars 1573 » ; 17 Lettre du roi « CHARLES [IX]... à mon cousin le duc de Montpensier,... gouverneur... en mon pays et duché de Bretaigne... Escript à Paris, ce XVIIe jour de janvier 1573 » ; 18 Lettre de « HENRY [duc D'ANJOU]... à mon cousin le duc de Montpensier,... gouverneur et lieutenant general du roy monseigneur et frere en son pays et duché de Bretaigne... Escript au camp devant La Rochelle, le VIIe jour de juing 1573 » ; 19 Lettre de « HENRY [duc D'ANJOU]... à mon cousin monseigneur le duc de Montpensier,... gouverneur... en Bretaigne... Escript au camp devant La Rochelle, le XXIXme jour de may 1573 » ; 20 Lettre de GASPARD DE COLIGNY, Sr DE « CHASTILLON,... à madame... la duchesse de Ferrare,... De Chastillon, ce XIIIme de septembre 1566 » ; 21 Lettre de GASPARD DE COLIGNY, Sr DE « CHASTILLON », et de FRANÇOIS DE COLIGNY, Sr D'« ANDELOT,... à madame... la duchesse de Ferrare,... De Chastillon sur Loing, ce XXVIIIe de may 1563 » ; 22 Lettre d'ODET DE COLIGNY, « cardinal DE CHASTILLON,... à madame... la duchesse de Ferrare,... De Moulins, ce XXIIIIe jour de janvier 1566 » ; 23 « Memoriale ad S. R. Majestatem per dominum episcopum Valentinum ». En latin. Copie ; 24 « Comitive millitum in terris Podoliae equestres et pedestres quibus ex senatusconsulto Varschoviae in ellectione habito conscriptio designata erat per magnificum Do. Georgium Iazlowieczky de Buczacz, pallatinum Russiaeet exercituum regni capitaneum, a die 4 septembris 1573... ad diem 19 february 1574 », etc. En latin ; 25 Promesse du roi de Pologne de ne pas établir de nouvelles constitutions sans l'assentiment du pays. En latin. Copie ; 26 Mémoire en latin sur le couronnement du roi de Pologne et sur l'assemblée générale tenue à cette occasion. Copie ; 27 « Roole des soldats à cheval entretenus en la Podolie vers les Tartares ; aussy somaire d'aulcuns soldats tant à cheval que à pied entretenus aux frontieres de Silesie, Ongarie et duché de Cracovie et Lictuanye ». En latin ; 28 Discours de JACQUES LASKOWSKI sur les différentes parties de l'administration du royaume de Pologne. En latin ; 29 Copie d'un decret des ordres du royaume de Pologne en faveur de Jean Saborowski. « Datum in comitiis generalibus electionis novi regis sub Warschovia, ad villam Kamieci, die 14 maii, anno 1573 ». En latin ; 30 « Coppie de l'acte delivré de la protestation faicte par l'evesque de Posnanie » contre la formule du serment prêté à Paris par Henri, élu roi de Pologne. « Actum et datum Parisiis, die Jovis decima mensis septembris, anno... millesimo quingentesimo septuagesimo tertio ». En latin ; 31 Ordonnance faicte par le duc DE NEVERS pour la discipline des gens de guerre ; 32 Copie du serment d'obéissance prêté par les soldats suisses au roi Henri II et à son capitaine général ; 33 Requête adressée au roi par GILBERT DE LEVIS, duc DE « VANTADOUR » ; 34 « Estimation des chesnes faictes par Jacques Fournier, Me orphevre, à Paris », données par le roy Charles IX au « roy de Pologne », son frère; suivie d'une décharge donnée par « René de Villequier,... grant chambellan... du roy de Poloigne... à Evry, le VIIIme jour de novembre M.V.C. soixante treize » ; 35 « Registro di tutti li castelli delli confini di quanti soldati li sono hora et quanti li era inanti la tregua fatta... 1573 ». En italien ; 36 « Registro di tutti li castelli di Russia et quanti soldati gli e per ogni castello et quanti fiorini li va ogni doi quartalli ». En italien ; 37 « Registro di vile che tengono molti de le quale S. M. non ne a utile alcuno, e pol contentar tutti li capitani di Lithuania ». En italien ; 38 Mémoire, en italien, relatif à une entreprise contre les Turcs ; suivi d'une énumération des forces militaires de l'expédition ; 39 Lettre d'« ADAMO CONARSKI, vescovo di Posnania... a l'excellentissimo signor duca di Nivers,... Di Pszczcw, a li 25 di gennaro 1574 ». En italien ; 40 Lettre d' « il vescovo di Ploscha... all... signor duca di Nevers,... Da Pultovia, alli 29 di settembre 1574 ». En italien ; 41 Noms des palatinats de Pologne et de Lithuanie. En latin ; 42 « Certificat de mons. de Valance [JEAN DE MONTLUC, envoyé des rois de France et de Pologne], pour Jo. Zaborouski [prévôt de Sandomir]... Datum Warssoviae, 26 mensis maii, anno Domini M.D.LXXIII ». En latin ; 43 Mémoire sur le sel de Pologne. En latin ; 44 Observations sur un mémoire relatif à la conduite du roi de Pologne dans son gouvernement. En latin ; 45 Mémoire adressé au roi de Pologne sur la conduite à tenir dans son gouvernement. En latin ; 46 Formulaire pour les lettres à envoyer à l'archevêque de Gnesne et autres personnages de Pologne. En latin ; 47 Lettre du duc DE NEVERS au comte Tenchini. « Datum Lutetiae, 27a die martii 1575 ». En latin ; 48 Lettre de « HIERAUSME, conte DE ROSDRAZOW, prowost de Wratislawie... à monsieur... le duc de Nevers,... Donné le 17 de febvrier 1574 » ; 49 « Brevis descriptio status et dispositionis... Nyvernensium ducis ab anno 1573 ad 87 ». En latin ; 50 Copie inachevée d'un mémoire justificatif de la conduite de monsieur DE NEVERS au retour de Henri, roi de Pologne ; 51 « Machines, ustilz et choses necessaires pour les entreprises » de guerre ; 52 « Memoire pour les oultilz qu'il faut pour les mines et sappes » ; 53 Mémoire sur les « estoffes pour fondre le canon » ; 54 « Avis du duc DE ZEBARAZ » sur l'« adsistance des Polonoys et des Kosaques » dans la guerre contre les Turcs ; 55 Lettre d'ADAMO CONARSKI, « vescovo di Posnania... al... signor Ludovico Gonzaga, principe di Mantoa e ducca di Nivers,... Di Isbiza, XVIII settembre 1574 ». En italien ; 56 Lettre d'ADAMO CONARSKI, « vescovo di Posnania... al... signor Ludovico Gonzaga, principe di Mantoa, duca di Nivers,... Di Isbiza, XXIIII settembre 1574 ». En italien ; 57 Mémoire sur les comices de Varsovie. En latin ; 58 Lettre de « JOANNES BORNKOWSKI » et « JOANNES DEMETRIUS SOLIKOWSKI,... domino Ludovico Gonzaga, Nivernensium duci, Manthuae et Registensi principi... Cracovie, 21 aprilis 1574 ». En latin ; 59 Lettre de « JOANNES BORNKOWSKI » et « JOANNES DEMETRIUS SOLIKOWSKI,... domino Martino Masparaultio,... Cracovie, ex arce, 17 maii 1574 ». En latin ; 60 Lettre de « JOANNES ZABOROWSKI,... domino Ludowico, duci Nivernensium... Datum Cracoviae, octava maii, anno 1574 ». En latin ; 61 Lettre d'« OSVALDUS REDLING, rector... Datum Fuldae, quinto calendas januarii anno 1573 ». En latin ; 62 Lettre d'« HIERONIMUS ROZORAZOROSKII,... praepositus Plocensis... duci Nivernensi... Divioni, 8 februarii 1575 ». En latin ; 63 Lettre de « STANISLAUS CARXCOVIUS, episcopus Cujavie et Pomeranie... domino duci Ninvernensi,... Datum Volborio, VI die mensis decembris anno 1574 ». En latin ; 64 Lettre de « JOANNES BORNKOWSKI » et « JOANNES DEMETRIUS SOLIKOWSKI,... domino Ludovico Gonzaga, Nivernensium duci, Mantuae et Registensi principi... Datum Cracoviae, 18a maii 1574 ». En latin ; 65 Lettre d'ADAMO CONARSKI, « vescovo di Posnania ». En italien ; 66 Extrait des constitutions du palatinat de Lublin, de l'an 1569, sur la tenue des comices généraux du royaume de Pologne. En latin ; 67 « Copie de la negotiation faitte par LE BOURG en Constantinople pres le Grand Seigneur, l'an passé 1569 » ; 68 Capitulation proposée par le « sultam SELIN, filz de sultam SOLIMAN, roy », au roi Charles IX. « Escript en la ville et cité de Constantinople... l'an neuf cens soixante dix sept et de Christ mil V.C. soixante neuf au mois d'octobre. Traduction faicte à l'original... en langue arabique... par DOMINICO OLIVIERI,... truchement et interprette du roy en ladicte langue ». Copie ; 69 Lettre de « CARLO SODERINI,... al... signor duca di Nevers,... Di Lione, l'ultimo del anno [15]74 ». En italien

Computation and Analysis of Multiple Structural-Change Models

In a recent paper, Bai and Perron (1998) considered theoretical issues related to the limiting distribution of estimators and test statistics in the linear model with multiple structural changes. In this companion paper, we consider practical issues for the empirical applications of the procedures. We first address the problem of estimation of the break dates and present an efficient algorithm to obtain global minimizers of the sum of squared residuals. This algorithm is based on the principle of dynamic programming and requires at most least-squares operations of order O(T 2) for any number of breaks. Our method can be applied to both pure and partial structural-change models. Secondly, we consider the problem of forming confidence intervals for the break dates under various hypotheses about the structure of the data and the errors across segments. Third, we address the issue of testing for structural changes under very general conditions on the data and the errors. Fourth, we address the issue of estimating the number of breaks. We present simulation results pertaining to the behavior of the estimators and tests in finite samples. Finally, a few empirical applications are presented to illustrate the usefulness of the procedures. All methods discussed are implemented in a GAUSS program available upon request for non-profit academic use.

Methods to Estimate Dynamic Stochastic General Equilibrium Models

This paper employs the one-sector Real Business Cycle model as a testing ground for four different procedures to estimate Dynamic Stochastic General Equilibrium (DSGE) models. The procedures are: 1 ) Maximum Likelihood, with and without measurement errors and incorporating Bayesian priors, 2) Generalized Method of Moments, 3) Simulated Method of Moments, and 4) Indirect Inference. Monte Carlo analysis indicates that all procedures deliver reasonably good estimates under the null hypothesis. However, there are substantial differences in statistical and computational efficiency in the small samples currently available to estimate DSGE models. GMM and SMM appear to be more robust to misspecification than the alternative procedures. The implications of the stochastic singularity of DSGE models for each estimation method are fully discussed.

Internet et le paradigme juridictionnel

"Mémoire présenté à la faculté des études supérieures en vue de l'obtention du grade de maîtrise en droit (LL.M.)"

Developmental trajectories of body mass index in early childhood : an 8-year longitudinal study

Trajectoires développementales de l’IMC durant l’enfance: Une étude longitudinale sur 8 ans. Introduction : L’obésité infantile, origine de nombreux problèmes de santé, représente un grand défi en santé publique. Récemment, l’importance d’étudier l’évolution du surpoids durant l’enfance ainsi que les facteurs de risques précoces pour l’obésité a été reconnue. Les trajectoires développementales d’indice de masse corporelle (IMC) chez les jeunes représentent une approche innovatrice qui nous permet de mieux comprendre cette problématique importante. Objectifs: 1) Identifier des trajectoires développementales distinctes de groupes d’enfants selon leur IMC durant l’enfance, et 2) Explorer les facteurs de risques précoces qui prédisent l’appartenance de l’enfant à la trajectoire d’IMC le plus élevé Hypothèses: 1) On s’attend à retrouver un groupe d’enfants qui suit une trajectoire d’IMC élevée durant l’enfance. 2) On s’attend à ce que certaines caractéristiques de la mère (ex : tabac pendant la grossesse et IMC élevé), soient associées à l’appartenance de l’enfant au groupe ayant la trajectoire «IMC élevé ». Méthodes: Estimation des trajectoires développementales d’IMC d’enfants, dans un échantillon populationnel (n=1957) au Québec (ELDEQ). Les IMC ont été calculés à partir de données fournies par les mères des enfants et recueillis chaque année sur une durée de 8 ans. Des données propres à l’enfant sa mère, ainsi que socioéconomiques, ont étés recueillies. Une régression logistique multinomiale a été utilisée pour distinguer les enfants avec un IMC élevé des autres enfants, selon les facteurs de risques précoces. Les programmes PROC TRAJ (extension de SAS), SPSS (version 16), et SAS (version 9.1.3) ont été utilisés pour ces analyses. Résultats: Trois trajectoires d’IMC ont étés identifiées : IMC « bas-stable » (54,5%), IMC « modéré » (41,0%) et IMC « élevé et en hausse » (4,5%). Le groupe « élevé et en hausse » incluait des enfants pour qui l’IMC à 8 ans dépassait la valeur limite pour l’obésité. Les analyses de régression logistique ont révélé que deux facteurs de risques maternels étaient significativement associés avec la trajectoire “en hausse” par rapport aux deux autres groupes : le tabac durant la grossesse et le surpoids maternel. Conclusions: Des risques d’obésité infantile peuvent êtres identifiés dès la grossesse. Des études d’intervention sont requises pour identifier la possibilité de réduire le risque d’obésité chez l’enfant en ciblant le tabac et le surpoids maternelle durant la grossesse. Mots clés: Indice de masse corporelle (IMC), obésité infantile, trajectoires développementales de groupe, facteurs de risque précoce, étude populationnelle, tabac pendant la grossesse, obésité maternelle.

Évaluation et contrôle de l'irrégularité de la prise médicamenteuse : proposition et développement de stratégies rationnelles fondées sur une démarche de modélisations pharmacocinétiques et pharmacodynamiques

L'hétérogénéité de réponses dans un groupe de patients soumis à un même régime thérapeutique doit être réduite au cours d'un traitement ou d'un essai clinique. Deux approches sont habituellement utilisées pour atteindre cet objectif. L'une vise essentiellement à construire une observance active. Cette approche se veut interactive et fondée sur l'échange ``médecin-patient '', ``pharmacien-patient'' ou ``vétérinaire-éleveurs''. L'autre plutôt passive et basée sur les caractéristiques du médicament, vise à contrôler en amont cette irrégularité. L'objectif principal de cette thèse était de développer de nouvelles stratégies d'évaluation et de contrôle de l'impact de l'irrégularité de la prise du médicament sur l'issue thérapeutique. Plus spécifiquement, le premier volet de cette recherche consistait à proposer des algorithmes mathématiques permettant d'estimer efficacement l'effet des médicaments dans un contexte de variabilité interindividuelle de profils pharmacocinétiques (PK). Cette nouvelle méthode est fondée sur l'utilisation concommitante de données \textit{in vitro} et \textit{in vivo}. Il s'agit de quantifier l'efficience ( c-à-dire efficacité plus fluctuation de concentrations \textit{in vivo}) de chaque profil PK en incorporant dans les modèles actuels d'estimation de l'efficacité \textit{in vivo}, la fonction qui relie la concentration du médicament de façon \textit{in vitro} à l'effet pharmacodynamique. Comparativement aux approches traditionnelles, cette combinaison de fonction capte de manière explicite la fluctuation des concentrations plasmatiques \textit{in vivo} due à la fonction dynamique de prise médicamenteuse. De plus, elle soulève, à travers quelques exemples, des questions sur la pertinence de l'utilisation des indices statiques traditionnels ($C_{max}$, $AUC$, etc.) d'efficacité comme outil de contrôle de l'antibiorésistance. Le deuxième volet de ce travail de doctorat était d'estimer les meilleurs temps d'échantillonnage sanguin dans une thérapie collective initiée chez les porcs. Pour ce faire, nous avons développé un modèle du comportement alimentaire collectif qui a été par la suite couplé à un modèle classique PK. À l'aide de ce modèle combiné, il a été possible de générer un profil PK typique à chaque stratégie alimentaire particulière. Les données ainsi générées, ont été utilisées pour estimer les temps d'échantillonnage appropriés afin de réduire les incertitudes dues à l'irrégularité de la prise médicamenteuse dans l'estimation des paramètres PK et PD . Parmi les algorithmes proposés à cet effet, la méthode des médianes semble donner des temps d'échantillonnage convenables à la fois pour l'employé et pour les animaux. Enfin, le dernier volet du projet de recherche a consisté à proposer une approche rationnelle de caractérisation et de classification des médicaments selon leur capacité à tolérer des oublis sporadiques. Méthodologiquement, nous avons, à travers une analyse globale de sensibilité, quantifié la corrélation entre les paramètres PK/PD d'un médicament et l'effet d'irrégularité de la prise médicamenteuse. Cette approche a consisté à évaluer de façon concomitante l'influence de tous les paramètres PK/PD et à prendre en compte, par la même occasion, les relations complexes pouvant exister entre ces différents paramètres. Cette étude a été réalisée pour les inhibiteurs calciques qui sont des antihypertenseurs agissant selon un modèle indirect d'effet. En prenant en compte les valeurs des corrélations ainsi calculées, nous avons estimé et proposé un indice comparatif propre à chaque médicament. Cet indice est apte à caractériser et à classer les médicaments agissant par un même mécanisme pharmacodynamique en terme d'indulgence à des oublis de prises médicamenteuses. Il a été appliqué à quatre inhibiteurs calciques. Les résultats obtenus étaient en accord avec les données expérimentales, traduisant ainsi la pertinence et la robustesse de cette nouvelle approche. Les stratégies développées dans ce projet de doctorat sont essentiellement fondées sur l'analyse des relations complexes entre l'histoire de la prise médicamenteuse, la pharmacocinétique et la pharmacodynamique. De cette analyse, elles sont capables d'évaluer et de contrôler l'impact de l'irrégularité de la prise médicamenteuse avec une précision acceptable. De façon générale, les algorithmes qui sous-tendent ces démarches constitueront sans aucun doute, des outils efficients dans le suivi et le traitement des patients. En outre, ils contribueront à contrôler les effets néfastes de la non-observance au traitement par la mise au point de médicaments indulgents aux oublis

Approximation de la distribution a posteriori d'un modèle Gamma-Poisson hiérarchique à effets mixtes

La méthode que nous présentons pour modéliser des données dites de "comptage" ou données de Poisson est basée sur la procédure nommée Modélisation multi-niveau et interactive de la régression de Poisson (PRIMM) développée par Christiansen et Morris (1997). Dans la méthode PRIMM, la régression de Poisson ne comprend que des effets fixes tandis que notre modèle intègre en plus des effets aléatoires. De même que Christiansen et Morris (1997), le modèle étudié consiste à faire de l'inférence basée sur des approximations analytiques des distributions a posteriori des paramètres, évitant ainsi d'utiliser des méthodes computationnelles comme les méthodes de Monte Carlo par chaînes de Markov (MCMC). Les approximations sont basées sur la méthode de Laplace et la théorie asymptotique liée à l'approximation normale pour les lois a posteriori. L'estimation des paramètres de la régression de Poisson est faite par la maximisation de leur densité a posteriori via l'algorithme de Newton-Raphson. Cette étude détermine également les deux premiers moments a posteriori des paramètres de la loi de Poisson dont la distribution a posteriori de chacun d'eux est approximativement une loi gamma. Des applications sur deux exemples de données ont permis de vérifier que ce modèle peut être considéré dans une certaine mesure comme une généralisation de la méthode PRIMM. En effet, le modèle s'applique aussi bien aux données de Poisson non stratifiées qu'aux données stratifiées; et dans ce dernier cas, il comporte non seulement des effets fixes mais aussi des effets aléatoires liés aux strates. Enfin, le modèle est appliqué aux données relatives à plusieurs types d'effets indésirables observés chez les participants d'un essai clinique impliquant un vaccin quadrivalent contre la rougeole, les oreillons, la rub\'eole et la varicelle. La régression de Poisson comprend l'effet fixe correspondant à la variable traitement/contrôle, ainsi que des effets aléatoires liés aux systèmes biologiques du corps humain auxquels sont attribués les effets indésirables considérés.

Le Mercosur et la protection de la propriété intellectuelle : formation, évolution historique et importance dans ses relations internationales

Depuis le début de son histoire, l’être humain cherche, pour de nombreuses raisons à s’associer et développer le milieu où il vit. Ces caractéristiques culminèrent dans les phénomènes de mondialisation et de révolution technologique responsables de la diminution des distances physiques et temporelles du monde contemporain. Ces facteurs ont stimulé les processus d’intégration régionaux, lesquels ont pour but le renforcement de ses adhérents face à la nouvelle réalité mondiale. Parmi ces processus, on note le Mercosur qui vise à la création d’un marché commun entre l’Argentine, le Brésil, le Paraguay et l’Uruguay. Fondé sur l’intergouvernabilité et avec une structure organisationnelle définie, ce bloc économique possède un cadre normatif en évolution, cherchant constamment la sécurité juridique dans divers secteurs. Quant à la protection des droits de propriété intellectuelle, on constate une ancienne préoccupation de la région sur le thème. Toutefois, le peu de normes existantes au sein du bloc confèrent des droits subjectifs minimaux et la plupart d’entre eux ne sont pas intégrés dans les législations nationales, lesquelles sont alignées sur les traités internationaux. Dans ses relations internationales, le Mercosur a une longue histoire de négociations avec l’Union Européenne et une histoire plus modeste avec le Canada, dans lesquelles l’attention portée aux droits de propriété intellectuelle n’est pas uniforme. De plus, sa participation dans les organismes internationaux responsables de ces droits est faible. Dans ce contexte, même avec un futur incertain le Mercosur présente de grandes expectatives. Toutefois, la préoccupation à protéger les droits de propriété intellectuelle est présente et réelle et commence à prendre de l’ampleur dans ses relations internes et externes.

Post orthodontic effects of SARPE on sleep-disordered breathing in young adults as observed in a sleep laboratory

Introduction: L’expansion palatine du maxillaire a beaucoup d’effets positifs sur la respiration et la qualité du sommeil, mais peu d'études ont examiné ces données sur des adultes ayant dépassé l’âge permettant de bénéficier d'une expansion palatine conventionnelle. Le but de cette recherche est d’évaluer la stabilité de l’EPRAC (expansion palatine rapide assistée chirurgicalement) et son effet sur les troubles respiratoires après l’ablation des appareils orthodontiques. Méthodes: Neuf patients (Âge moyen 21, entre 16-39 ans) nécessitant une EPRAC ont passé des nuits dans un laboratoire de sommeil, et ce avant l’EPRAC, après l’EPRAC, et après l’ablation des appareils fixes. Les radiographies céphalométriques postéroantérieures ainsi que les modèles d’étude ont été pris pendant ces trois périodes de temps. Résultats: L’analyse des modèles d’étude a démontré une récidive significative au niveau des distances inter-molaires et inter-canines au niveau du maxillaire seulement. Les analyses céphalométriques ont démontré une récidive au niveau de la largeur maxillaire. Aucun changement important n'a été observé dans les stades de sommeil, mais une réduction importante dans l’index de ronflement a été notée. De plus, il y avait moins de changements entre les stades de sommeil. Conclusions: La récidive squelettique est minime et cliniquement non significative. Par contre, les changements dans les distances intermolaires et intercanines sont cliniquement importants. Il semble également qu'une EPRAC ait un effet positif sur la qualité de sommeil par la réduction de l’indice de ronflement ainsi que sur la diminution des changements entre les stades de sommeil.

Algorithmes pour la reconstruction de génomes ancestraux

L’inférence de génomes ancestraux est une étape essentielle pour l’étude de l’évolution des génomes. Connaissant les génomes d’espèces éteintes, on peut proposer des mécanismes biologiques expliquant les divergences entre les génomes des espèces modernes. Diverses méthodes visant à résoudre ce problème existent, se classant parmis deux grandes catégories : les méthodes de distance et les méthodes de synténie. L’état de l’art des distances génomiques ne permettant qu’un certain répertoire de réarrangements pour le moment, les méthodes de synténie sont donc plus appropriées en pratique. Nous proposons une méthode de synténie pour la reconstruction de génomes ancestraux basée sur une définition relaxée d’adjacences de gènes, permettant un contenu en gène inégal dans les génomes modernes causé par des pertes de gènes de même que des duplications de génomes entiers (DGE). Des simulations sont effectuées, démontrant une capacité de former une solution assemblée en un nombre réduit de régions ancestrales contigües par rapport à d’autres méthodes tout en gardant une bonne fiabilité. Des applications sur des données de levures et de plantes céréalières montrent des résultats en accord avec d’autres publications, notamment la présence de fusion imbriquée de chromosomes pendant l’évolution des céréales.

Essays on numerically efficient inference in nonlinear and non-Gaussian state space models, and commodity market analysis.

The first two articles build procedures to simulate vector of univariate states and estimate parameters in nonlinear and non Gaussian state space models. We propose state space speci fications that offer more flexibility in modeling dynamic relationship with latent variables. Our procedures are extension of the HESSIAN method of McCausland[2012]. Thus, they use approximation of the posterior density of the vector of states that allow to : simulate directly from the state vector posterior distribution, to simulate the states vector in one bloc and jointly with the vector of parameters, and to not allow data augmentation. These properties allow to build posterior simulators with very high relative numerical efficiency. Generic, they open a new path in nonlinear and non Gaussian state space analysis with limited contribution of the modeler. The third article is an essay in commodity market analysis. Private firms coexist with farmers' cooperatives in commodity markets in subsaharan african countries. The private firms have the biggest market share while some theoretical models predict they disappearance once confronted to farmers cooperatives. Elsewhere, some empirical studies and observations link cooperative incidence in a region with interpersonal trust, and thus to farmers trust toward cooperatives. We propose a model that sustain these empirical facts. A model where the cooperative reputation is a leading factor determining the market equilibrium of a price competition between a cooperative and a private firm

Les modèles vectoriels et multiplicatifs avec erreurs non-négatives de séries chronologiques.

L'objectif du présent mémoire vise à présenter des modèles de séries chronologiques multivariés impliquant des vecteurs aléatoires dont chaque composante est non-négative. Nous considérons les modèles vMEM (modèles vectoriels et multiplicatifs avec erreurs non-négatives) présentés par Cipollini, Engle et Gallo (2006) et Cipollini et Gallo (2010). Ces modèles représentent une généralisation au cas multivarié des modèles MEM introduits par Engle (2002). Ces modèles trouvent notamment des applications avec les séries chronologiques financières. Les modèles vMEM permettent de modéliser des séries chronologiques impliquant des volumes d'actif, des durées, des variances conditionnelles, pour ne citer que ces applications. Il est également possible de faire une modélisation conjointe et d'étudier les dynamiques présentes entre les séries chronologiques formant le système étudié. Afin de modéliser des séries chronologiques multivariées à composantes non-négatives, plusieurs spécifications du terme d'erreur vectoriel ont été proposées dans la littérature. Une première approche consiste à considérer l'utilisation de vecteurs aléatoires dont la distribution du terme d'erreur est telle que chaque composante est non-négative. Cependant, trouver une distribution multivariée suffisamment souple définie sur le support positif est plutôt difficile, au moins avec les applications citées précédemment. Comme indiqué par Cipollini, Engle et Gallo (2006), un candidat possible est une distribution gamma multivariée, qui impose cependant des restrictions sévères sur les corrélations contemporaines entre les variables. Compte tenu que les possibilités sont limitées, une approche possible est d'utiliser la théorie des copules. Ainsi, selon cette approche, des distributions marginales (ou marges) peuvent être spécifiées, dont les distributions en cause ont des supports non-négatifs, et une fonction de copule permet de tenir compte de la dépendance entre les composantes. Une technique d'estimation possible est la méthode du maximum de vraisemblance. Une approche alternative est la méthode des moments généralisés (GMM). Cette dernière méthode présente l'avantage d'être semi-paramétrique dans le sens que contrairement à l'approche imposant une loi multivariée, il n'est pas nécessaire de spécifier une distribution multivariée pour le terme d'erreur. De manière générale, l'estimation des modèles vMEM est compliquée. Les algorithmes existants doivent tenir compte du grand nombre de paramètres et de la nature élaborée de la fonction de vraisemblance. Dans le cas de l'estimation par la méthode GMM, le système à résoudre nécessite également l'utilisation de solveurs pour systèmes non-linéaires. Dans ce mémoire, beaucoup d'énergies ont été consacrées à l'élaboration de code informatique (dans le langage R) pour estimer les différents paramètres du modèle. Dans le premier chapitre, nous définissons les processus stationnaires, les processus autorégressifs, les processus autorégressifs conditionnellement hétéroscédastiques (ARCH) et les processus ARCH généralisés (GARCH). Nous présentons aussi les modèles de durées ACD et les modèles MEM. Dans le deuxième chapitre, nous présentons la théorie des copules nécessaire pour notre travail, dans le cadre des modèles vectoriels et multiplicatifs avec erreurs non-négatives vMEM. Nous discutons également des méthodes possibles d'estimation. Dans le troisième chapitre, nous discutons les résultats des simulations pour plusieurs méthodes d'estimation. Dans le dernier chapitre, des applications sur des séries financières sont présentées. Le code R est fourni dans une annexe. Une conclusion complète ce mémoire.