Thèse diffusée initialement dans le cadre d'un projet pilote des Presses de l'Université de Montréal/Centre d'édition numérique UdeM (1997-2008) avec l'autorisation de l'auteur.

Cette thèse se trouve à l'intersection de deux sujets de recherches en économie qui ont connus des développements récents : d'un côté, la modélisation et l'étude des réseaux économiques, de l'autre, la théorie du transport et ses applications en économie. Les quatre chapitres de cette thèse déclinent plusieurs résultats théoriques en lien avec ces deux sujets, en insistant sur leurs connections. Le premier chapitre modélise des contrats bipartites sur un marché décentralisé comme un problème de flux. Ce travail prouve l’existence d'un équilibre compétitif et discute son efficacité. Nous présentons des interprétations de ce modèle dans le cas où les biens sont indivisibles. Ces spécifications sont utilisées pour proposer un modèle du marché du prêt interbancaire. Dans le second chapitre, co-écrit avec Alfred Galichon et Larry Samuelson, nous prouvons un théorème de statiques comparatives monotones. Nous appliquons ensuite ce résultat à plusieurs modèles économiques classiques (modèles d'appariement, modèles hédoniques et problèmes de flux de coût minimum). Le troisième chapitre, co-écrit avec Alfred Galichon et Arthur Charpentier, présente des outils pour résoudre les problèmes d'appariement sur des réseaux de grandes dimensions avant de les mettre en pratique. Enfin le quatrième chapitre, co-écrit avec Rakesh Vohra, montre comment un assureur monopolistique peut utiliser les externalités entre agents - modélisées comme un réseau - pour maximiser ses revenus. Bien qu’exerçant le niveau d'effort optimal d'un point de vue social, tous les agents préfèrent l'équilibre précédant l'introduction de l'assureur. / This dissertation lies at the intersection of two fields of research in economics that have recently substantially developed: on the one hand, the modeling and study of economic networks, and on the other hand, transport theory and its applications in economics. The four chapters of this dissertation present theoretical results in relation with these two topics and put stress on their connections. The first chapter models bipartite contracts on a decentralized market as an equilibrium flow problem. We prove the existence of a competitive equilibrium outcome and discuss its effciency. We interpret this equilibrium in the case of indivisible commodities. As an illustration, we build a model for the overnight interbank loan market with counterparty risk and collateralization costs. In the second one, written in collaboration with Alfred Galichon and Larry Samuelson, we prove a monotone comparative statics theorem that we then apply to several classical economic models (matching models, min-cost flow problems and hedonic models). The third chapter is a joint work with Alfred Galichon and Arthur Charpentier. It presents tools to solve for matching problems on large geographic networks before applying them to examples. Finally, the fourth chapter, written with Rakesh Vohra, shows how a monopolistic insurer can use externalities between agents - modeled as a network - to maximize his profit. We show that a monopolistic insurer decreases the welfare of all agents.

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Thèse diffusée initialement dans le cadre d'un projet pilote des Presses de l'Université de Montréal/Centre d'édition numérique UdeM (1997-2008) avec l'autorisation de l'auteur. / Dans cette thèse nous abordons la conception de nouveaux schémas de type volumes finis pour la résolution de systèmes hyperboliques non-linéaires pour la prédiction des écoulements compressibles instationnaires. Les nouveaux schémas présentés s'appuient tous sur les schémas proposés par Arminjon-Viallon et Arminjon-Stanescu-Viallon en 2 dimensions spatiales qui, eux, furent dérivés du schéma de Nessyahu-Tadmor en une dimension d'espace construit à partir du schéma décalé bien connu de Lax-Friedrichs. Ces schémas peuvent être considérés comme étant tous du type de Godunov et ont pour caractéristique principale d'éviter la résolution des problèmes de Riemann aux interfaces en utilisant 2 maillages différents pour, respectivement, les pas de temps pairs et impairs. Pour éviter la trop grande diffusion amenée par le schéma de Lax-Friedrichs, on a eu recours à l'utilisation d'une technique nommée MUSCL, originalement proposée par van Leer, consistant à reconstruire la solution constante par cellule en une solution linéaire par cellule tout en limitant les oscillations grâce à l'utilisation de fonctions non-linéaires. On obtient tout d'abord une extension en 3 dimensions spatiales sur des maillages cartésiens structurés. Ensuite, nous abordons le cas de maillages non-structurés composés de tétraèdres, et la formulation mathématique du schéma associé à ces cellules. Pour réduire les temps de calcul, un nouveau, schéma de type centré fondé sur celui de Nessyahu-Tadmor mais évitant l'utilisation d'un pas intermédiaire, et composé d'un nouveau flux est proposé en une et 2 dimensions spatiales pour des maillages structurés, puis en 3 dimensions sur des maillages non structurés composés de tétraèdres. Les résultats obtenus démontrent que les nouvelles méthodes sont moins sensibles aux maillages déformés et qu'elles sont plus simples à mettre en œuvre du fait que le problème de Riemann est évité et qu'aucune information sur la décomposition de la discontinuité en les différents champs caractéristiques du système n'est nécessaire.

Dans ce mémoire, nous cherchons à modéliser des tables à deux entrées monotones en lignes et/ou en colonnes, pour une éventuelle application sur les tables de mortalité. Nous adoptons une approche bayésienne non paramétrique et représentons la forme fonctionnelle des données par splines bidimensionnelles. L’objectif consiste à condenser une table de mortalité, c’est-à-dire de réduire l’espace d’entreposage de la table en minimisant la perte d’information. De même, nous désirons étudier le temps nécessaire pour reconstituer la table. L’approximation doit conserver les mêmes propriétés que la table de référence, en particulier la monotonie des données. Nous travaillons avec une base de fonctions splines monotones afin d’imposer plus facilement la monotonie au modèle. En effet, la structure flexible des splines et leurs dérivées faciles à manipuler favorisent l’imposition de contraintes sur le modèle désiré. Après un rappel sur la modélisation unidimensionnelle de fonctions monotones, nous généralisons l’approche au cas bidimensionnel. Nous décrivons l’intégration des contraintes de monotonie dans le modèle a priori sous l’approche hiérarchique bayésienne. Ensuite, nous indiquons comment obtenir un estimateur a posteriori à l’aide des méthodes de Monte Carlo par chaînes de Markov. Finalement, nous étudions le comportement de notre estimateur en modélisant une table de la loi normale ainsi qu’une table t de distribution de Student. L’estimation de nos données d’intérêt, soit la table de mortalité, s’ensuit afin d’évaluer l’amélioration de leur accessibilité. / This master’s thesis is about the estimation of bivariate tables which are monotone within the rows and/or the columns, with a special interest in the approximation of life tables. This problem is approached through a nonparametric Bayesian regression model, in particular linear combinations of regression splines. By condensing a life table, our goal is to reduce its storage space without losing the entries’ accuracy. We will also study the reconstruction time of the table with our estimators. The properties of the reference table, specifically its monotonicity, must be preserved in the estimation. We are working with a monotone spline basis since splines are flexible and their derivatives can easily be manipulated. Those properties enable the imposition of constraints of monotonicity on our model. A brief review on univariate approximations of monotone functions is then extended to bivariate estimations. We use hierarchical Bayesian modeling to include the constraints in the prior distributions. We then explain the Markov chain Monte Carlo algorithm to obtain a posterior estimator. Finally, we study the estimator’s behaviour by applying our model on the Standard Normal table and the Student’s t table. We estimate our data of interest, the life table, to establish the improvement in data accessibility.

Dans ce mémoire, nous cherchons à modéliser des tables à deux entrées monotones en lignes et/ou en colonnes, pour une éventuelle application sur les tables de mortalité. Nous adoptons une approche bayésienne non paramétrique et représentons la forme fonctionnelle des données par splines bidimensionnelles. L’objectif consiste à condenser une table de mortalité, c’est-à-dire de réduire l’espace d’entreposage de la table en minimisant la perte d’information. De même, nous désirons étudier le temps nécessaire pour reconstituer la table. L’approximation doit conserver les mêmes propriétés que la table de référence, en particulier la monotonie des données. Nous travaillons avec une base de fonctions splines monotones afin d’imposer plus facilement la monotonie au modèle. En effet, la structure flexible des splines et leurs dérivées faciles à manipuler favorisent l’imposition de contraintes sur le modèle désiré. Après un rappel sur la modélisation unidimensionnelle de fonctions monotones, nous généralisons l’approche au cas bidimensionnel. Nous décrivons l’intégration des contraintes de monotonie dans le modèle a priori sous l’approche hiérarchique bayésienne. Ensuite, nous indiquons comment obtenir un estimateur a posteriori à l’aide des méthodes de Monte Carlo par chaînes de Markov. Finalement, nous étudions le comportement de notre estimateur en modélisant une table de la loi normale ainsi qu’une table t de distribution de Student. L’estimation de nos données d’intérêt, soit la table de mortalité, s’ensuit afin d’évaluer l’amélioration de leur accessibilité. / This master’s thesis is about the estimation of bivariate tables which are monotone within the rows and/or the columns, with a special interest in the approximation of life tables. This problem is approached through a nonparametric Bayesian regression model, in particular linear combinations of regression splines. By condensing a life table, our goal is to reduce its storage space without losing the entries’ accuracy. We will also study the reconstruction time of the table with our estimators. The properties of the reference table, specifically its monotonicity, must be preserved in the estimation. We are working with a monotone spline basis since splines are flexible and their derivatives can easily be manipulated. Those properties enable the imposition of constraints of monotonicity on our model. A brief review on univariate approximations of monotone functions is then extended to bivariate estimations. We use hierarchical Bayesian modeling to include the constraints in the prior distributions. We then explain the Markov chain Monte Carlo algorithm to obtain a posterior estimator. Finally, we study the estimator’s behaviour by applying our model on the Standard Normal table and the Student’s t table. We estimate our data of interest, the life table, to establish the improvement in data accessibility.

.

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

L'objectif de cette thèse est de modéliser la transmission du paludisme dans la région montagneuse de l'ouest du Kenya, en se servant des outils de systèmes dynamiques. Nous considérons deux modèles mathématiques. Le premier prend en compte une susceptibilité et une infectivité différentielle dans les métapopulations, et le second un taux de saturation des repas sanguins dans la population des moustiques. Dans le premier modèle, nous considérons plusieurs écosystèmes identifiés comme zones sensibles dans la région montagneuse de l'ouest du Kenya. Dans ce modèle, ces zones sensibles sont considérées comme nos différents patchs. Les populations de chaque patch sont divisées en deux : les enfants et les adultes. Le modèle nous permet d'évaluer le rôle de l'hétérogénéité de l'écosystème et la persistance de l'épidémie dans la région, due à la structuration d'âge. Nous prenons en compte la susceptibilité et l'infectivité différentielle afin d'étendre le modèle d'un patch en un modèle à plusieurs patchs. Après avoir subdivisé la région en n zones sensibles, nous faisons une analyse mathématique du modèle obtenu. Pour effectuer cette analyse, nous utilisons la théorie des systèmes triangulaires, des systèmes dynamiques monotones, des systèmes dynamiques non linéaires anti-monotones et le principe d'invariance de LaSalle. Un des éléments très utilisés dans notre analyse qui est un concept clé en épidémiologie, est le taux de reproduction de base, très souvent noté Ro. Cette quantité, sans dimension, est le nombre moyen de cas secondaires, engendré par un individu infectieux typique durant sa période d'infectiosité, quand il est introduit dans une population constituée entièrement de susceptibles. L'existence et la stabilité du point d'équilibre sans maladie (DFE) sont établies et nous prouvons que le DFE est globalement asymptotiquement stable lorsque Ro<1. Lorsque Ro>1, le modèle admet un point d'équilibre endémique qui est globalement asymptotiquement stable. L'analyse de notre modèle montre que la structuration d'âge réduit l'ampleur de l'infection. En utilisant les données relevées, nous faisons quelques simulations numériques afin de montrer l'impact de la métapopulation et de la structuration d'âge sur le taux de reproduction de base. Dans la seconde partie, nous formulons un modèle de paludisme avec saturation du taux d'alimentation des moustiques qui nous conduit à une incidence non linéaire. Nous démontrons que DFE est globalement asymptotiquement stable si Ro<1. Lorsque Ro>1, il existe un unique point d'équilibre endémique qui est globalement asymptotiquement stable. Des simulations numériques sont faites afin d'illustrer l'impact de la saturation d'alimentation sur le taux de reproduction de base / The objective of this thesis is to model highland malaria in western Kenya using dynamical systems. Two mathematical models are formulated ; one, on differentiated susceptibility and differentiated infectivity in a metapopulation setting with age structure, the other, a saturated vector feeding rate model with disease induced deaths and varying host and vector populations. In the first model, we consider the different ecosystems identified as malaria hotspots in the western Kenya highlands and consider the ecosystems as different patches. The population in each patch is classified as, either child or, adult. The model will aid in examining the role of ecosystem heterogeneity and age structure to the persistent malaria epidemics in the highlands. We formulate the differentiated susceptibility and infectivity model that extend to multiple patches the well known epidemiological models in one patch. Classifying the hot spots as n patches, we give its mathematical analysis using the theory of triangular system, monotone non-linear dynamical systems, and Lyapunov-Lasalle invariance principle techniques. Key to our analysis is the definition of a reproductive number, Ro, the number of new infections caused by one individual in an otherwise fully susceptible population throughout the duration of the infectious period. The existence and stability of disease-free and endemic equilibrium is established. We prove that the disease free state of the systems is globally asymptotically stable when the basic reproduction number Ro<1, and when Ro>1 an endemic equilibrium is established which is locally and globally asymptotically stable. The model shows that the age structuring reduces the magnitude of infection. Using relevant data we did some simulation, to demonstrate the role played by metapopulation and age structuring on the incidence and Ro. In the second part we formulate a model for malaria with saturation on the vector feeding rates that lead to a nonlinear function in the infection term. The vector feeding rate is assumed, as in the predator prey models, to rise linearly as a function of the host-vector ratio until it reaches a threshold Qv, after which the vector feeds freely at its desired rate. The two populations are variable and drive malaria transmission, such that when the vectors are fewer than hosts, the rate of feeding is determined by the vectors feeding desire, whereas, when the hosts are more than the vectors, the feeding rate is limited by host availability and other feeding sources may have to be sought by the vector. Malaria induced deaths are introduced in the host population, while the vector is assumed to survive with the parasite till its death. We prove that the Disease Free Equilibrium is locally and globally asymptotically stable if Ro<1 and when Ro>1, an endemic equilibrium emerges, which is unique, locally and globally asymptotically stable. The role of the saturated mosquito feeding rate is explored with simulation showing the crucial role it plays especially on the basic reproduction number

On s'intéresse dans cette thèse à l'étude théorique et numérique des équations différentielles stochastiques multivoques et leurs applications à la modélisation de structures mécaniques sous sollicitations aléatoires. Les équations différentielles stochastiques considérées comportent dans le terme de dérive un opérateur multivoque maximal monotone pour lesquelles l'existence et l'unicité de solutions ont déjà été obtenues dans un cadre euclidien. Pour de telles équations on montre la convergence d'un schéma numérique, faisant intervenir grâce à la maximalité et à la monotonie des opérateurs considérés, des applications exclusivement univoques, rendant son implémentation aisée. Un ordre de convergence est de plus obtenu sous certaines conditions sur le coefficient de diffusion. Pour enrichir la modélisation, on envisage des équations différentielles stochastiques multivoques d'ordre 2 évoluant sur une variété riemanienne pour lesquelles ont été obtenues l'existence et l'unicité d'une solution. Des simulations numériques sur des modèles d'association en série ou en parallèle de ressorts, amortisseurs et patins (ou éléments de Saint-Venant), dont la formalisation mathématique fait intervenir des équations différentielles stochastiques multivoques, ont permis de valider des méthodes d'identification de paramètres à partir de cycles d'hystérésis.