Modélisation bayésienne des changements aux niches écologiques causés par le réchauffement climatique

Cette thèse présente des méthodes de traitement de données de comptage en particulier et des données discrètes en général. Il s'inscrit dans le cadre d'un projet stratégique du CRNSG, nommé CC-Bio, dont l'objectif est d'évaluer l'impact des changements climatiques sur la répartition des espèces animales et végétales. Après une brève introduction aux notions de biogéographie et aux modèles linéaires mixtes généralisés aux chapitres 1 et 2 respectivement, ma thèse s'articulera autour de trois idées majeures. Premièrement, nous introduisons au chapitre 3 une nouvelle forme de distribution dont les composantes ont pour distributions marginales des lois de Poisson ou des lois de Skellam. Cette nouvelle spécification permet d'incorporer de l'information pertinente sur la nature des corrélations entre toutes les composantes. De plus, nous présentons certaines propriétés de ladite distribution. Contrairement à la distribution multidimensionnelle de Poisson qu'elle généralise, celle-ci permet de traiter les variables avec des corrélations positives et/ou négatives. Une simulation permet d'illustrer les méthodes d'estimation dans le cas bidimensionnel. Les résultats obtenus par les méthodes bayésiennes par les chaînes de Markov par Monte Carlo (CMMC) indiquent un biais relatif assez faible de moins de 5% pour les coefficients de régression des moyennes contrairement à ceux du terme de covariance qui semblent un peu plus volatils. Deuxièmement, le chapitre 4 présente une extension de la régression multidimensionnelle de Poisson avec des effets aléatoires ayant une densité gamma. En effet, conscients du fait que les données d'abondance des espèces présentent une forte dispersion, ce qui rendrait fallacieux les estimateurs et écarts types obtenus, nous privilégions une approche basée sur l'intégration par Monte Carlo grâce à l'échantillonnage préférentiel. L'approche demeure la même qu'au chapitre précédent, c'est-à-dire que l'idée est de simuler des variables latentes indépendantes et de se retrouver dans le cadre d'un modèle linéaire mixte généralisé (GLMM) conventionnel avec des effets aléatoires de densité gamma. Même si l'hypothèse d'une connaissance a priori des paramètres de dispersion semble trop forte, une analyse de sensibilité basée sur la qualité de l'ajustement permet de démontrer la robustesse de notre méthode. Troisièmement, dans le dernier chapitre, nous nous intéressons à la définition et à la construction d'une mesure de concordance donc de corrélation pour les données augmentées en zéro par la modélisation de copules gaussiennes. Contrairement au tau de Kendall dont les valeurs se situent dans un intervalle dont les bornes varient selon la fréquence d'observations d'égalité entre les paires, cette mesure a pour avantage de prendre ses valeurs sur (-1;1). Initialement introduite pour modéliser les corrélations entre des variables continues, son extension au cas discret implique certaines restrictions. En effet, la nouvelle mesure pourrait être interprétée comme la corrélation entre les variables aléatoires continues dont la discrétisation constitue nos observations discrètes non négatives. Deux méthodes d'estimation des modèles augmentés en zéro seront présentées dans les contextes fréquentiste et bayésien basées respectivement sur le maximum de vraisemblance et l'intégration de Gauss-Hermite. Enfin, une étude de simulation permet de montrer la robustesse et les limites de notre approche.

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.

Sur les prolongements de sous-copules

L’objet du travail est d’étudier les prolongements de sous-copules. Un cas important de l’utilisation de tels prolongements est l’estimation non paramétrique d’une copule par le lissage d’une sous-copule (la copule empirique). Lorsque l’estimateur obtenu est une copule, cet estimateur est un prolongement de la souscopule. La thèse présente au chapitre 2 la construction et la convergence uniforme d’un estimateur bona fide d’une copule ou d’une densité de copule. Cet estimateur est un prolongement de type copule empirique basé sur le lissage par le produit tensoriel de fonctions de répartition splines. Le chapitre 3 donne la caractérisation de l’ensemble des prolongements possibles d’une sous-copule. Ce sujet a été traité par le passé; mais les constructions proposées ne s’appliquent pas à la dépendance dans des espaces très généraux. Le chapitre 4 s’attèle à résoudre le problème suivant posé par [Carley, 2002]. Il s’agit de trouver la borne supérieure des prolongements en dimension 3 d’une sous-copule de domaine fini.

Implementación de cópulas para la estimación del valor en riesgo

La dependencia entre las series financieras, es un parámetro fundamental para la estimación de modelos de Riesgo. El Valor en Riesgo (VaR) es una de las medidas más importantes utilizadas para la administración y gestión de Riesgos Financieros, en la actualidad existen diferentes métodos para su estimación, como el método por simulación histórica, el cual no asume ninguna distribución sobre los retornos de los factores de riesgo o activos, o los métodos paramétricos que asumen normalidad sobre las distribuciones. En este documento se introduce la teoría de cópulas, como medida de dependencia entre las series, se estima un modelo ARMA-GARCH-Cópula para el cálculo del Valor en Riesgo de un portafolio compuesto por dos series financiera, la tasa de cambio Dólar-Peso y Euro-Peso. Los resultados obtenidos muestran que la estimación del VaR por medio de copulas es más preciso en relación a los métodos tradicionales.

Diversification when it hurts? The joint distributions of real estate and equity markets

Much of the literature on the construction of mixed asset portfolios and the case for property as a risk diversifier rests on correlations measured over the whole of a given time series. Recent developments in finance, however, focuses on dependence in the tails of the distribution. Does property offer diversification from equity markets when it is most needed - when equity returns are poor. The paper uses an empirical copula approach to test tail dependence between property and equity for the UK and for a global portfolio. Results show strong tail dependence: in the UK, the dependence in the lower tail is stronger than in the upper tail, casting doubt on the defensive properties of real estate stocks.

The heliospheric magnetic field

The heliospheric magnetic field (HMF) is the extension of the coronal magnetic field carried out into the solar system by the solar wind. It is the means by which the Sun interacts with planetary magnetospheres and channels charged particles propagating through the heliosphere. As the HMF remains rooted at the solar photosphere as the Sun rotates, the large-scale HMF traces out an Archimedean spiral. This pattern is distorted by the interaction of fast and slow solar wind streams, as well as the interplanetary manifestations of transient solar eruptions called coronal mass ejections. On the smaller scale, the HMF exhibits an array of waves, discontinuities, and turbulence, which give hints to the solar wind formation process. This review aims to summarise observations and theory of the small- and large-scale structure of the HMF. Solar-cycle and cycle-to-cycle evolution of the HMF is discussed in terms of recent spacecraft observations and pre-spaceage proxies for the HMF in geomagnetic and galactic cosmic ray records.

Morphology and phylogenetic relationships of a remarkable new genus and two new species of Neotropical freshwater stingrays from the Amazon basin (Chondrichthyes: Potamotrygonidae)

The morphology and phylogenetic relationships of a new genus and two new species of Neotropical freshwater stingrays, family Potamotrygonidae, are investigated and described in detail. The new genus, Heliotrygon, n. gen., and its two new species, Heliotrygon gomesi, n. sp. (type-species) and Heliotrygon rosai, n. sp., are compared to all genera and species of potamotrygonids, based on revisions in progress. Some of the derived features of Heliotrygon include its unique disc proportions (disc highly circular, convex anteriorly at snout region, its width and length very similar), extreme subdivision of suborbital canal (forming a complex honeycomb-like pattern anterolaterally on disc), stout and triangular pelvic girdle, extremely reduced caudal sting, basibranchial copula with very slender and acute anterior extension, and precerebral and frontoparietal fontanellae of about equal width, tapering very little posteriorly. Both new species can be distinguished by their unique color patterns: Heliotrygon gomesi is uniform gray to light tan or brownish dorsally, without distinct patterns, whereas Heliotrygon rosai is characterized by numerous white to creamy-white vermiculate markings over a light brown, tan or gray background color. Additional proportional characters that may further distinguish both species are also discussed. Morphological descriptions are provided for dermal denticles, ventral lateral-line canals, skeleton, and cranial, hyoid and mandibular muscles of Heliotrygon, which clearly corroborate it as the sister group of Paratrygon. Both genera share numerous derived features of the ventral lateral-line canals, neurocranium, scapulocoracoid, pectoral basals, clasper morphology, and specific patterns of the adductor mandibulae and spiracularis medialis muscles. Potamotrygon and Plesiotrygon are demonstrated to share derived characters of their ventral lateral-line canals, in addition to the presence of angular cartilages. Our morphological phylogeny is further corroborated by a molecular phylogenetic analysis of cytochrome b based on four sequences (637 base pairs in length), representing two distinct haplotypes for Heliotrygon gomesi. Parsimony analysis produced a single most parsimonious tree revealing Heliotrygon and Paratrygon as sister taxa (boot-strap proportion of 70%), which together are the sister group to a clade including Plesiotrygon and species of Potamotrygon. These unusual stingrays highlight that potamotrygonid diversity, both in terms of species composition and undetected morphological and molecular patterns, is still poorly known.

Automatic spiral analysis for objective assessment of motor symptoms in Parkinson's disease

Objective: To develop a method for objective quantification of PD motor symptoms related to Off episodes and peak dose dyskinesias, using spiral data gathered by using a touch screen telemetry device. The aim was to objectively characterize predominant motor phenotypes (bradykinesia and dyskinesia), to help in automating the process of visual interpretation of movement anomalies in spirals as rated by movement disorder specialists. Background: A retrospective analysis was conducted on recordings from 65 patients with advanced idiopathic PD from nine different clinics in Sweden, recruited from January 2006 until August 2010. In addition to the patient group, 10 healthy elderly subjects were recruited. Upper limb movement data were collected using a touch screen telemetry device from home environments of the subjects. Measurements with the device were performed four times per day during week-long test periods. On each test occasion, the subjects were asked to trace pre-drawn Archimedean spirals, using the dominant hand. The pre-drawn spiral was shown on the screen of the device. The spiral test was repeated three times per test occasion and they were instructed to complete it within 10 seconds. The device had a sampling rate of 10Hz and measured both position and time-stamps (in milliseconds) of the pen tip. Methods: Four independent raters (FB, DH, AJ and DN) used a web interface that animated the spiral drawings and allowed them to observe different kinematic features during the drawing process and to rate task performance. Initially, a number of kinematic features were assessed including ‘impairment’, ‘speed’, ‘irregularity’ and ‘hesitation’ followed by marking the predominant motor phenotype on a 3-category scale: tremor, bradykinesia and/or choreatic dyskinesia. There were only 2 test occasions for which all the four raters either classified them as tremor or could not identify the motor phenotype. Therefore, the two main motor phenotype categories were bradykinesia and dyskinesia. ‘Impairment’ was rated on a scale from 0 (no impairment) to 10 (extremely severe) whereas ‘speed’, ‘irregularity’ and ‘hesitation’ were rated on a scale from 0 (normal) to 4 (extremely severe). The proposed data-driven method consisted of the following steps. Initially, 28 spatiotemporal features were extracted from the time series signals before being presented to a Multilayer Perceptron (MLP) classifier. The features were based on different kinematic quantities of spirals including radius, angle, speed and velocity with the aim of measuring the severity of involuntary symptoms and discriminate between PD-specific (bradykinesia) and/or treatment-induced symptoms (dyskinesia). A Principal Component Analysis was applied on the features to reduce their dimensions where 4 relevant principal components (PCs) were retained and used as inputs to the MLP classifier. Finally, the MLP classifier mapped these components to the corresponding visually assessed motor phenotype scores for automating the process of scoring the bradykinesia and dyskinesia in PD patients whilst they draw spirals using the touch screen device. For motor phenotype (bradykinesia vs. dyskinesia) classification, the stratified 10-fold cross validation technique was employed. Results: There were good agreements between the four raters when rating the individual kinematic features with intra-class correlation coefficient (ICC) of 0.88 for ‘impairment’, 0.74 for ‘speed’, 0.70 for ‘irregularity’, and moderate agreements when rating ‘hesitation’ with an ICC of 0.49. When assessing the two main motor phenotype categories (bradykinesia or dyskinesia) in animated spirals the agreements between the four raters ranged from fair to moderate. There were good correlations between mean ratings of the four raters on individual kinematic features and computed scores. The MLP classifier classified the motor phenotype that is bradykinesia or dyskinesia with an accuracy of 85% in relation to visual classifications of the four movement disorder specialists. The test-retest reliability of the four PCs across the three spiral test trials was good with Cronbach’s Alpha coefficients of 0.80, 0.82, 0.54 and 0.49, respectively. These results indicate that the computed scores are stable and consistent over time. Significant differences were found between the two groups (patients and healthy elderly subjects) in all the PCs, except for the PC3. Conclusions: The proposed method automatically assessed the severity of unwanted symptoms and could reasonably well discriminate between PD-specific and/or treatment-induced motor symptoms, in relation to visual assessments of movement disorder specialists. The objective assessments could provide a time-effect summary score that could be useful for improving decision-making during symptom evaluation of individualized treatment when the goal is to maximize functional On time for patients while minimizing their Off episodes and troublesome dyskinesias.

Modelando contágio financeiro através de cópulas

This article aims to test the hypothesis of contagion between the indices of nancial markets from the United States to Brazil, Japan and England for the period 2000 to 2009. Time varying copulas were used to capture the impact of Sub-prime crisis in the dependence between markets. The implemented model was a ARMA(1,0) st-ARCH(1,2) to the marginal distributions and Normal and Joe Clayton (SJC) copulas for the joint distribution. The results obtained allow to conclude that both for the gaussiana copula and for the SJC copula there is evidence of contagion between the American market and the Brazilian market. For the other two markets Londoner and Japanese, the evidence of the presence of contagion between these markets and the American has not been suf ciently clear in both copula

Essays on infrastructure in Brazil

A presente tese engloba cinco trabalhos sobre dois setores brasileiros de infraestrutura: setor de energia e de transportes. No primeiro trabalho, modela-se o lance dos leilões de linhas de transmissão de energia, buscando compreender porque os lances ganhadores têm se situado em níveis abaixo do que se espera. No segundo paper, estima-se a demanda de energia do Brasil e as elasticidades preço e renda que os consumidores apresentem em relação a essa demanda, usando uma técnica ainda não aplicada na literatura e incluindo o período da Crise do Racionamento – a qual pode ter mudado o padrão de consumo dos agentes. No terceiro trabalho, estuda-se a evolução do mercado de gás natural liquefeito (GNL) ao redor do mundo. O GNL pode ser o link que faltava entre os mercados de gás natural e essa hipótese é testada por meio de análises de séries de tempo e de cópulas. O quarto paper discute a entrada do GNL no sistema energético brasileiro e suas eventuais consequências. Por fim, o quinto paper analisa a experiência de leilões de rodovias no Brasil. O trabalho levanta relevantes insights a respeito do modelo de leilão utilizado e dos detalhes contratuais dessas licitações.

Statistical modelling of extreme rainfall in Madeira Island

Extreme rainfall events have triggered a significant number of flash floods in Madeira Island along its past and recent history. Madeira is a volcanic island where the spatial rainfall distribution is strongly affected by its rugged topography. In this thesis, annual maximum of daily rainfall data from 25 rain gauge stations located in Madeira Island were modelled by the generalised extreme value distribution. Also, the hypothesis of a Gumbel distribution was tested by two methods and the existence of a linear trend in both distributions parameters was analysed. Estimates for the 50– and 100–year return levels were also obtained. Still in an univariate context, the assumption that a distribution function belongs to the domain of attraction of an extreme value distribution for monthly maximum rainfall data was tested for the rainy season. The available data was then analysed in order to find the most suitable domain of attraction for the sampled distribution. In a different approach, a search for thresholds was also performed for daily rainfall values through a graphical analysis. In a multivariate context, a study was made on the dependence between extreme rainfall values from the considered stations based on Kendall’s τ measure. This study suggests the influence of factors such as altitude, slope orientation, distance between stations and their proximity of the sea on the spatial distribution of extreme rainfall. Groups of three pairwise associated stations were also obtained and an adjustment was made to a family of extreme value copulas involving the Marshall–Olkin family, whose parameters can be written as a function of Kendall’s τ association measures of the obtained pairs.

As sentenças com é ruim que no português brasileiro

Based on studies on the structure of copulates sentences in Brazilian Portuguese (PB), according to grammar theory, it s analyzed on this dissertation, the complex copulate sentences what, in the PB, initiated by the formation é ruim que in the PB, bringing the hypothesis that they can have one predicacional reading (PRED) and other especificacional reading (ESP). It s revealed what, though they are similar slightly, that makes the different that are triggered is structural configuration of sentence and the form as come emerged the its constituents: that of PRED interpretation, that we defined as Common Copulate Sentence (or SCC), it s checked that the constituent ruim, and only it, is predicate of a Small Clause, and appears in the structure on even position wherein is pronounced, not passing by the method of movement, from where introduce wide predication on all subject, that is the CP built-in; those sentences é ruim que understanding ESP, that denominate of Negative Copulate Sentence (or SCN), whereas the fixed expression contradict or right what is asserted in built-in, it s perceived what the ruim is a part of this fixed expression, and surely crystallized, which appears in structure per movement, germinated as adjunct of Inflexional Phrase (IP) and being elevation to specifier to FocP (or SpecFocP), where acquires discursive interpretation of focus. Besides the mode as the ruim or the fixed expression emerge in that sentences, the relation between the copula and main verb likewise contribute to distinguish the sentences: 1) regarding flexion of mode, when is a SCC, the copula must be in the mode indicative, and the main verb, in the subjunctive; when is a SCN, copula and main verb must be always on indicative; 2) regarding flexion of temp, both the copula as main verb of the SCC can be flexed, whereas at the SCN, the main verb can vary in present, past and future, but the copula should appear, necessarily, in third persona of present of indicative, what confirms our hypothesis that there is an fixed expression at the copulate sentences with é ruim que ESP. Other two evidences are pointed as characteristics that distinguish the sentence with é ruim que PRED of sentence with é ruim que ESP: 1) in semantics, the constituent ruim equals not good, that has appreciation, when the reading will PRED; already the pair é ruim equals not, at the ESP; 2) in prosody, there discrete sound elevation the ruim on the other constituents of SCC, PRED, while there is accentuated acoustic elevation on the ruim of SCN, ESP. Our search it is grounded in authors as Zanfeliz (2000), Modesto (2001), Mioto (2003), Kato & Ribeiro (2006), Lobo (2006), Quarezemin (2006, 2009, 2011, 2012) e Resenes (2009), researchers that devoted their attention on studies the formation and organization of the constituents of cleaved sentences and focalization of constituent, basing itself in approach generative of linguistic

Evaluation of the courtship and of the hybrid male sterility among Drosophila buzzatii cluster species (Diptera, Drosophilidae)

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Ordenação de constituintes em construções categorial, tética e apresentativa

A ordenação de constituintes oracionais nas variedades portuguesas é aqui tratada considerando-se três moldes de conteúdo, Tético, Apresentativo e Categorial. A linearização das estruturas oracionais refletem decisões assumidas na formulação do Nível Interpessoal e configuram predominantemente Molde de Conteúdo Categorial Tópico-orientado, com o constituinte Tópico ocupando a posição P I, a palavra verbal, a posição P M, e os outros constituintes, as posições à direita e à esquerda de P M. em Construções Téticas, a oração toda, por ser focal, ocupa o domínio de P F, já em Construções Apresentativas a cópula vazia assume a posição absoluta P M, e o sintagma nominal Tópico/Foco posiciona-se em P F.

Implante de tendão autógeno do músculo flexor superficial dos dedos no reparo de desvio do pênis de bovines

In this experiment, 18 hybrid steers, from the ages between 16 and 20 months, were submitted to the removal of the penis apical ligament to provoke a penis deviation. After passed 30 days, to be tested, the animals were put together to the females that were in heat. It was verified the formation of the ventral and right lateral penis deviation with incapacity of the copula. The bovines were submitted to the autologous implant from the one following superficial portion tendon from the superficial flexor fingers muscle, on the top of albuginea tunic, replacing the apical penis ligament removed. The material implanted was well tolerated, appearance a few inflammatory cells and formation of the fibrous conjunctive tissue, causing resistant adherence and capacity of sustaining the free penis extremity. The bovines, when were among the females that were in heat, did not show a penis deviation nor sexual behavior modification. The method revealed itself satisfactory and it can be indicated to correct ventral and right lateral penis deviation in bovines.