Extreme value Birnbaum-Saunders regression models applied to environmental data

Extreme value models are widely used in different areas. The Birnbaum–Saunders distribution is receiving considerable attention due to its physical arguments and its good properties. We propose a methodology based on extreme value Birnbaum–Saunders regression models, which includes model formulation, estimation, inference and checking. We further conduct a simulation study for evaluating its performance. A statistical analysis with real-world extreme value environmental data using the methodology is provided as illustration.

Accounting for the dead in the longitudinal analysis of income-related health inequalities

This paper develops an accounting framework to consider the effect of deaths on the longitudinal analysis of income-related health inequalities. Ignoring deaths or using inverse probability weights (IPWs) to re-weight the sample for mortality-related attrition can produce misleading results, since to do so would be to disregard the most extreme of all health outcomes. Incorporating deaths into the longitudinal analysis of income-related health inequalities provides a more complete picture in terms of the evaluation of health changes in respect to socioeconomic status. We illustrate our work by investigating health mobility in Quality Adjusted Life Years (QALYs) as measured by the SF6D from 1999 till 2004 using the British Household Panel Survey (BHPS). We show that for Scottish males explicitly accounting for the dead, rather than using IPWs to account for mortality-related attrition, changes the direction of the relationship between relative health changes and initial income position, while for other population groups it increases the strength of this relationship by up to 14 times. When deaths are explicitly incorporated into the analysis it is found that over this five year period for both Scotland and England & Wales the relative health changes were significantly regressive such that the poor experienced a larger share of the health losses relative to their initial share of health and a large amount of this was related to mortality.

Bilateral symmetry of radial pulse in high-level tennis players: implications for the validity of central aortic pulse wave analysis.

BACKGROUND: Reconstruction of the central aortic pressure wave from the noninvasive recording of the radial pulse with applanation tonometry has become a standard tool in the field of hypertension. It is not presently known whether recording the radial pulse on the dominant or the nondominant side has any effect on such reconstruction. METHOD: We carried out radial applanation tonometry on both forearms in young, healthy, male volunteers, who were either sedentary (n = 11) or high-level tennis players (n = 10). The purpose of including tennis players was to investigate individuals with extreme asymmetry between the dominant and nondominant upper limb. RESULTS: In the sedentary individuals, forearm circumference and handgrip strength were slightly larger on the dominant (mean +/- SD respectively 27.9 +/- 1.5 cm and 53.8 +/- 10 kg) than on nondominant side (27.3 +/- 1.6 cm, P < 0.001 vs. dominant, and 52.1 +/- 11 kg, P = NS). In the tennis players, differences between sides were more conspicuous (forearm circumference: dominant 28.0 +/- 1.7 cm nondominant 26.4 +/- 1.5 cm, P < 0.001; handgrip strength 61.4 +/- 10.8 vs. 53.4 +/- 9.7 kg, P < 0.001). We found that in both sedentary individuals and tennis players, the radial pulse had identical shape on both sides and, consequently, the reconstructed central aortic pressure waveforms, as well as derived indices of central pulsatility, were not dependent on the side where applanation tonometry was carried out. CONCLUSION: Evidence from individuals with maximal asymmetry of dominant vs. nondominant upper limb indicates that laterality of measurement is not a methodological issue for central pulse wave analysis carried out with radial applanation tonometry.

Design and analysis of ChIP-Seq experiments for nuclear factor I DNA-binding proteins

SUMMARY : Eukaryotic DNA interacts with the nuclear proteins using non-covalent ionic interactions. Proteins can recognize specific nucleotide sequences based on the sterical interactions with the DNA and these specific protein-DNA interactions are the basis for many nuclear processes, e.g. gene transcription, chromosomal replication, and recombination. New technology termed ChIP-Seq has been recently developed for the analysis of protein-DNA interactions on a whole genome scale and it is based on immunoprecipitation of chromatin and high-throughput DNA sequencing procedure. ChIP-Seq is a novel technique with a great potential to replace older techniques for mapping of protein-DNA interactions. In this thesis, we bring some new insights into the ChIP-Seq data analysis. First, we point out to some common and so far unknown artifacts of the method. Sequence tag distribution in the genome does not follow uniform distribution and we have found extreme hot-spots of tag accumulation over specific loci in the human and mouse genomes. These artifactual sequence tags accumulations will create false peaks in every ChIP-Seq dataset and we propose different filtering methods to reduce the number of false positives. Next, we propose random sampling as a powerful analytical tool in the ChIP-Seq data analysis that could be used to infer biological knowledge from the massive ChIP-Seq datasets. We created unbiased random sampling algorithm and we used this methodology to reveal some of the important biological properties of Nuclear Factor I DNA binding proteins. Finally, by analyzing the ChIP-Seq data in detail, we revealed that Nuclear Factor I transcription factors mainly act as activators of transcription, and that they are associated with specific chromatin modifications that are markers of open chromatin. We speculate that NFI factors only interact with the DNA wrapped around the nucleosome. We also found multiple loci that indicate possible chromatin barrier activity of NFI proteins, which could suggest the use of NFI binding sequences as chromatin insulators in biotechnology applications. RESUME : L'ADN des eucaryotes interagit avec les protéines nucléaires par des interactions noncovalentes ioniques. Les protéines peuvent reconnaître les séquences nucléotidiques spécifiques basées sur l'interaction stérique avec l'ADN, et des interactions spécifiques contrôlent de nombreux processus nucléaire, p.ex. transcription du gène, la réplication chromosomique, et la recombinaison. Une nouvelle technologie appelée ChIP-Seq a été récemment développée pour l'analyse des interactions protéine-ADN à l'échelle du génome entier et cette approche est basée sur l'immuno-précipitation de la chromatine et sur la procédure de séquençage de l'ADN à haut débit. La nouvelle approche ChIP-Seq a donc un fort potentiel pour remplacer les anciennes techniques de cartographie des interactions protéine-ADN. Dans cette thèse, nous apportons de nouvelles perspectives dans l'analyse des données ChIP-Seq. Tout d'abord, nous avons identifié des artefacts très communs associés à cette méthode qui étaient jusqu'à présent insoupçonnés. La distribution des séquences dans le génome ne suit pas une distribution uniforme et nous avons constaté des positions extrêmes d'accumulation de séquence à des régions spécifiques, des génomes humains et de la souris. Ces accumulations des séquences artéfactuelles créera de faux pics dans toutes les données ChIP-Seq, et nous proposons différentes méthodes de filtrage pour réduire le nombre de faux positifs. Ensuite, nous proposons un nouvel échantillonnage aléatoire comme un outil puissant d'analyse des données ChIP-Seq, ce qui pourraient augmenter l'acquisition de connaissances biologiques à partir des données ChIP-Seq. Nous avons créé un algorithme d'échantillonnage aléatoire et nous avons utilisé cette méthode pour révéler certaines des propriétés biologiques importantes de protéines liant à l'ADN nommés Facteur Nucléaire I (NFI). Enfin, en analysant en détail les données de ChIP-Seq pour la famille de facteurs de transcription nommés Facteur Nucléaire I, nous avons révélé que ces protéines agissent principalement comme des activateurs de transcription, et qu'elles sont associées à des modifications de la chromatine spécifiques qui sont des marqueurs de la chromatine ouverte. Nous pensons que lés facteurs NFI interagir uniquement avec l'ADN enroulé autour du nucléosome. Nous avons également constaté plusieurs régions génomiques qui indiquent une éventuelle activité de barrière chromatinienne des protéines NFI, ce qui pourrait suggérer l'utilisation de séquences de liaison NFI comme séquences isolatrices dans des applications de la biotechnologie.

Extreme genetic differentiation among the remnant populations of marble trout (Salmo marmoratus) in Slovenia.

Populations of the marble trout (Salmo marmoratus) have declined critically due to introgression by brown trout (Salmo trutta) strains. In order to define strategies for long-term conservation, we examined the genetic structure of the 8 known pure populations using 15 microsatellite loci. The analyses reveal extraordinarily strong genetic differentiation among populations separated by < 15 km, and extremely low levels of intrapopulation genetic variability. As natural recolonization seems highly unlikely, appropriate management and conservation strategies should comprise the reintroduction of pure populations from mixed stocks (translocation) to avoid further loss of genetic diversity.

Evolution of a cellular immune response in Drosophila: a phenotypic and genomic comparative analysis.

Understanding the genomic basis of evolutionary adaptation requires insight into the molecular basis underlying phenotypic variation. However, even changes in molecular pathways associated with extreme variation, gains and losses of specific phenotypes, remain largely uncharacterized. Here, we investigate the large interspecific differences in the ability to survive infection by parasitoids across 11 Drosophila species and identify genomic changes associated with gains and losses of parasitoid resistance. We show that a cellular immune defense, encapsulation, and the production of a specialized blood cell, lamellocytes, are restricted to a sublineage of Drosophila, but that encapsulation is absent in one species of this sublineage, Drosophila sechellia. Our comparative analyses of hemopoiesis pathway genes and of genes differentially expressed during the encapsulation response revealed that hemopoiesis-associated genes are highly conserved and present in all species independently of their resistance. In contrast, 11 genes that are differentially expressed during the response to parasitoids are novel genes, specific to the Drosophila sublineage capable of lamellocyte-mediated encapsulation. These novel genes, which are predominantly expressed in hemocytes, arose via duplications, whereby five of them also showed signatures of positive selection, as expected if they were recruited for new functions. Three of these novel genes further showed large-scale and presumably loss-of-function sequence changes in D. sechellia, consistent with the loss of resistance in this species. In combination, these convergent lines of evidence suggest that co-option of duplicated genes in existing pathways and subsequent neofunctionalization are likely to have contributed to the evolution of the lamellocyte-mediated encapsulation in Drosophila.

A compositional statistical analysis of capital stock

Most of economic literature has presented its analysis under the assumption of homogeneous capital stock.However, capital composition differs across countries. What has been the pattern of capital compositionassociated with World economies? We make an exploratory statistical analysis based on compositional datatransformed by Aitchinson logratio transformations and we use tools for visualizing and measuring statisticalestimators of association among the components. The goal is to detect distinctive patterns in the composition.As initial findings could be cited that:1. Sectorial components behaved in a correlated way, building industries on one side and , in a lessclear view, equipment industries on the other.2. Full sample estimation shows a negative correlation between durable goods component andother buildings component and between transportation and building industries components.3. Countries with zeros in some components are mainly low income countries at the bottom of theincome category and behaved in a extreme way distorting main results observed in the fullsample.4. After removing these extreme cases, conclusions seem not very sensitive to the presence ofanother isolated cases

Executive functions profile in extreme eating/weight conditions: from anorexia nervosa to obesity.

BACKGROUND Extreme weight conditions (EWC) groups along a continuum may share some biological risk factors and intermediate neurocognitive phenotypes. A core cognitive trait in EWC appears to be executive dysfunction, with a focus on decision making, response inhibition and cognitive flexibility. Differences between individuals in these areas are likely to contribute to the differences in vulnerability to EWC. The aim of the study was to investigate whether there is a common pattern of executive dysfunction in EWC while comparing anorexia nervosa patients (AN), obese subjects (OB) and healthy eating/weight controls (HC). METHODS Thirty five AN patients, fifty two OB and one hundred thirty seven HC were compared using the Wisconsin Card Sorting Test (WCST); Stroop Color and Word Test (SCWT); and Iowa Gambling Task (IGT). All participants were female, aged between 18 and 60 years. RESULTS There was a significant difference in IGT score (F(1.79); p<.001), with AN and OB groups showing the poorest performance compared to HC. On the WCST, AN and OB made significantly more errors than controls (F(25.73); p<.001), and had significantly fewer correct responses (F(2.71); p<.001). Post hoc analysis revealed that the two clinical groups were not significantly different from each other. Finally, OB showed a significant reduced performance in the inhibition response measured with the Stroop test (F(5.11); p<.001) compared with both AN and HC. CONCLUSIONS These findings suggest that EWC subjects (namely AN and OB) have similar dysfunctional executive profile that may play a role in the development and maintenance of such disorders.

DNA multigene sequencing of topotypic specimens of the fascioliasis vector Lymnaea diaphana and phylogenetic analysis of the genus Pectinidens (Gastropoda)

Freshwater lymnaeid snails are crucial in defining transmission and epidemiology of fascioliasis. In South America, human endemic areas are related to high altitudes in Andean regions. The species Lymnaea diaphana has, however, been involved in low altitude areas of Chile, Argentina and Peru where human infection also occurs. Complete nuclear ribosomal DNA 18S, internal transcribed spacer (ITS)-2 and ITS-1 and fragments of mitochondrial DNA 16S and cytochrome c oxidase (cox)1 genes of L. diaphana specimens from its type locality offered 1,848, 495, 520, 424 and 672 bp long sequences. Comparisons with New and Old World Galba/Fossaria, Palaearctic stagnicolines, Nearctic stagnicolines, Old World Radix and Pseudosuccinea allowed to conclude that (i) L. diaphana shows sequences very different from all other lymnaeids, (ii) each marker allows its differentiation, except cox1 amino acid sequence, and (iii) L. diaphana is not a fossarine lymnaeid, but rather an archaic relict form derived from the oldest North American stagnicoline ancestors. Phylogeny and large genetic distances support the genus Pectinidens as the first stagnicoline representative in the southern hemisphere, including colonization of extreme world regions, as most southern Patagonia, long time ago. The phylogenetic link of L. diaphana with the stagnicoline group may give light to the aforementioned peculiar low altitude epidemiological scenario of fascioliasis.

A Factor analysis of hidrochemical composition of Llobregat river basin

Hydrogeological research usually includes some statistical studies devised to elucidate mean background state, characterise relationships among different hydrochemical parameters, and show the influence of human activities. These goals are achieved either by means of a statistical approach or by mixing modelsbetween end-members. Compositional data analysis has proved to be effective with the first approach, but there is no commonly accepted solution to the end-member problem in a compositional framework.We present here a possible solution based on factor analysis of compositions illustrated with a case study.We find two factors on the compositional bi-plot fitting two non-centered orthogonal axes to the most representative variables. Each one of these axes defines a subcomposition, grouping those variables thatlay nearest to it. With each subcomposition a log-contrast is computed and rewritten as an equilibrium equation. These two factors can be interpreted as the isometric log-ratio coordinates (ilr) of three hiddencomponents, that can be plotted in a ternary diagram. These hidden components might be interpreted as end-members.We have analysed 14 molarities in 31 sampling stations all along the Llobregat River and its tributaries, with a monthly measure during two years. We have obtained a bi-plot with a 57% of explained totalvariance, from which we have extracted two factors: factor G, reflecting geological background enhanced by potash mining; and factor A, essentially controlled by urban and/or farming wastewater. Graphicalrepresentation of these two factors allows us to identify three extreme samples, corresponding to pristine waters, potash mining influence and urban sewage influence. To confirm this, we have available analysisof diffused and widespread point sources identified in the area: springs, potash mining lixiviates, sewage, and fertilisers. Each one of these sources shows a clear link with one of the extreme samples, exceptfertilisers due to the heterogeneity of their composition.This approach is a useful tool to distinguish end-members, and characterise them, an issue generally difficult to solve. It is worth note that the end-member composition cannot be fully estimated but only characterised through log-ratio relationships among components. Moreover, the influence of each endmember in a given sample must be evaluated in relative terms of the other samples. These limitations areintrinsic to the relative nature of compositional data

Preservation of an extreme transient geotherm in the Raft River detachment shear zone

Extensional detachment systems separate hot footwalls from cool hanging walls, but the degree to which this thermal gradient is the product of ductile or brittle deformation or a preserved original transient geotherm is unclear. Oxygen isotope thermometry using recrystallized quartz-muscovite pairs indicates a smooth thermal gradient (140 degrees C/100 m) across the gently dipping, quartzite-dominated detachment zone that bounds the Raft River core complex in northwest Utah (United States). Hydrogen isotope values of muscovite (delta D-Ms similar to-100 parts per thousand) and fluid inclusions in quartz (delta D-Fluid similar to-85 parts per thousand) indicate the presence of meteoric fluids during detachment dynamics. Recrystallized grain-shape fabrics and quartz c-axis fabric patterns reveal a large component of coaxial strain (pure shear), consistent with thinning of the detachment section. Therefore, the high thermal gradient preserved in the Raft River detachment reflects the transient geotherm that developed owing to shearing, thinning, and the potentially prominent role of convective flow of surface fluids.

A systematic review and meta-analysis of perineural dexamethasone for peripheral nerve blocks.

We systematically reviewed the safety and efficacy of perineural dexamethasone as an adjunct for peripheral nerve blockade in 29 controlled trials of 1695 participants. We grouped trials by the duration of local anaesthetic action (short- or medium- vs long-term). Dexamethasone increased the mean (95% CI) duration of analgesia by 233 (172-295) min when injected with short- or medium-term action local anaesthetics and by 488 (419-557) min when injected with long-term action local anaesthetics, p < 0.00001 for both. However, these results should be interpreted with caution due to the extreme heterogeneity of results, with I2 exceeding 90% for both analyses. Meta-regression did not show an interaction between dose of perineural dexamethasone (4-10 mg) and duration of analgesia (r2 = 0.02, p = 0.54). There were no differences between 4 and 8 mg dexamethasone on subgroup analysis.

Large-system analysis of a CDMA dynamic channel under a Markovian input process

We study the minimum mean square error (MMSE) and the multiuser efficiency η of large dynamic multiple access communication systems in which optimal multiuser detection is performed at the receiver as the number and the identities of active users is allowed to change at each transmission time. The system dynamics are ruled by a Markov model describing the evolution of the channel occupancy and a large-system analysis is performed when the number of observations grow large. Starting on the equivalent scalar channel and the fixed-point equation tying multiuser efficiency and MMSE, we extend it to the case of a dynamic channel, and derive lower and upper bounds for the MMSE (and, thus, for η as well) holding true in the limit of large signal–to–noise ratios and increasingly large observation time T.

Demographic effects of extreme winter weather in the barn owl.

Extreme weather events can lead to immediate catastrophic mortality. Due to their rare occurrence, however, the long-term impacts of such events for ecological processes are unclear. We examined the effect of extreme winters on barn owl (Tyto alba) survival and reproduction in Switzerland over a 68-year period (approximately 20 generations). This long-term data set allowed us to compare events that occurred only once in several decades to more frequent events. Winter harshness explained 17 and 49% of the variance in juvenile and adult survival, respectively, and the two harshest winters were associated with major population crashes caused by simultaneous low juvenile and adult survival. These two winters increased the correlation between juvenile and adult survival from 0.63 to 0.69. Overall, survival decreased non-linearly with increasing winter harshness in adults, and linearly in juveniles. In contrast, brood size was not related to the harshness of the preceding winter. Our results thus reveal complex interactions between climate and demography. The relationship between weather and survival observed during regular years is likely to underestimate the importance of climate variation for population dynamics.

Spanish monetary policy: A structural VAR analysis

Some past studies analyzed Spanish monetary policy with the standard VAR. Their problem is that this method obliges researchers to impose a certain extreme form of the short run policy rule on their models. Hence, it does not allow researchers to study the possibility of structural changes in this rule, either. This paper overcomes these problems by using the structural VAR. I find that the rule has always been that of partial accommodation. Prior to 1984, it was quite close to money targeting. After 1984, it became closer to the interest rate targeting, with more emphasis on the exchange rate.