News from compositions, the R package

The R-package “compositions”is a tool for advanced compositional analysis. Its basicfunctionality has seen some conceptual improvement, containing now some facilitiesto work with and represent ilr bases built from balances, and an elaborated subsys-tem for dealing with several kinds of irregular data: (rounded or structural) zeroes,incomplete observations and outliers. The general approach to these irregularities isbased on subcompositions: for an irregular datum, one can distinguish a “regular” sub-composition (where all parts are actually observed and the datum behaves typically)and a “problematic” subcomposition (with those unobserved, zero or rounded parts, orelse where the datum shows an erratic or atypical behaviour). Systematic classificationschemes are proposed for both outliers and missing values (including zeros) focusing onthe nature of irregularities in the datum subcomposition(s).To compute statistics with values missing at random and structural zeros, a projectionapproach is implemented: a given datum contributes to the estimation of the desiredparameters only on the subcompositon where it was observed. For data sets withvalues below the detection limit, two different approaches are provided: the well-knownimputation technique, and also the projection approach.To compute statistics in the presence of outliers, robust statistics are adapted to thecharacteristics of compositional data, based on the minimum covariance determinantapproach. The outlier classification is based on four different models of outlier occur-rence and Monte-Carlo-based tests for their characterization. Furthermore the packageprovides special plots helping to understand the nature of outliers in the dataset.Keywords: coda-dendrogram, lost values, MAR, missing data, MCD estimator,robustness, rounded zeros

Cerebral perfusion pressure and risk of brain hypoxia in severe head injury: a prospective observational study

INTRODUCTION Higher and lower cerebral perfusion pressure (CPP) thresholds have been proposed to improve brain tissue oxygen pressure (PtiO2) and outcome. We study the distribution of hypoxic PtiO2 samples at different CPP thresholds, using prospective multimodality monitoring in patients with severe traumatic brain injury. METHODS This is a prospective observational study of 22 severely head injured patients admitted to a neurosurgical critical care unit from whom multimodality data was collected during standard management directed at improving intracranial pressure, CPP and PtiO2. Local PtiO2 was continuously measured in uninjured areas and snapshot samples were collected hourly and analyzed in relation to simultaneous CPP. Other variables that influence tissue oxygen availability, mainly arterial oxygen saturation, end tidal carbon dioxide, body temperature and effective hemoglobin, were also monitored to keep them stable in order to avoid non-ischemic hypoxia. RESULTS Our main results indicate that half of PtiO2 samples were at risk of hypoxia (defined by a PtiO2 equal to or less than 15 mmHg) when CPP was below 60 mmHg, and that this percentage decreased to 25% and 10% when CPP was between 60 and 70 mmHg and above 70 mmHg, respectively (p < 0.01). CONCLUSION Our study indicates that the risk of brain tissue hypoxia in severely head injured patients could be really high when CPP is below the normally recommended threshold of 60 mmHg, is still elevated when CPP is slightly over it, but decreases at CPP values above it.

Measurement with an automated oscillometric wrist device with position sensor leads to lower values than measurements obtained with an automated oscillometric arm device from the same manufacturer in elderly persons

Introduction : Le monitoring de la tension artérielle à domicile est recommandé par plusieurs guidelines et a été montré être faisable chez la personne âgée. Les manomètres au poignet ont récemment été proposés pour la mesure de la tension artérielle à domicile, mais leur précision n'a pas été au préalable évaluée chez les patients âgés. Méthode : Quarante-huit participants (33 femmes et 15 hommes, moyenne d'âge 81.3±8.0 ans) ont leur tension artérielle mesurée avec un appareil au poignet avec capteur de position et un appareil au bras dans un ordre aléatoire et dans une position assise. Résultats : Les moyennes de mesures de tension artérielle étaient systématiquement plus basses avec l'appareil au poignet par rapport à celui du bras pour la pression systolique (120.1±2.2 vs. 130.5±2.2 mmHg, P < 0.001, moyenneidéviation standard) et pour la pression diastolique (66.011.3 vs. 69.7±1.3 mmHg, P < 0.001). De plus, une différence de lOmmHg ou plus grande entre l'appareil au bras et au poignet était observée dans 54.2 et 18,8% des mesures systoliques et diastoliques respectivement. Conclusion : Comparé à l'appareil au bras, l'appareil au poignet avec capteur de position sous-estimait systématiquement aussi bien la tension artérielle systolique que diastolique. L'ampleur de la différence est cliniquement significative et met en doute l'utilisation de l'appareil au poignet pour monitorer la tension artérielle chez la personne âgée. Cette étude indique le besoin de valider les appareils de mesures de la tension artérielle dans tous les groupes d'âge, y compris les personnes âgées.

Existential values transformation in all organ transplantation: A qualitative prospective study.

Organ transplantation offers a treatment of choice for patients suffering from end stage illnesses. The aim of this IRB approved prospective qualitative study was to analyze patients’ psychological concerns from their inclusion on the waiting list for first organ transplantation (TX) (T1; N=71; kidney, K=30; liver, Li=11; lung, Lu=15; heart, H=15) and six months after TX (T2; N=49; K=15; Li=10; Lu=14; H=10). Semi-structured interviews were conducted at home or in a place selected by patients. Qualitative pattern analysis (QUAPA) of the verbatim transcriptions was applied. T1 (K) Patients maintained an apparent normality (87%), building emotional protection (23%), and developing a fatalist attitude towards life (43%). (Li) Physical limits were set to spare energy until TX (73%). Illness led to reevaluation of life values (66%). (Lu) Physical and psychological self-protection was prioritized when health declined (67%). Modified life values, fatalism (33%) and spirituality (27 %) were mentioned. (H) Patients husbanded physical (80%) and psychological (67%) resources and self-protection. Modified life values and fatalist attitude towards life were reported (40%). T2 (K) New perspective on life was described, with increase of empathy towards others (20%). (Li) Positive identity and life values modifications (60%), greater openness towards others, closeness to significant ones (30%) and a more self-centered attitude (30%) prioritizing the essential (20%) were reported. Lack of respect of life values generated anger (40%). (Lu) Setting existential priorities and increase in spirituality (64%), along with the development of new life values, greater openness to others (57%) and closeness to significant ones (21%) were underlined. Lack of respect of human values induced negative feelings (36%). Self-centered attitudes, setting limits to other people were mentioned (29%). (H) Change in life values with setting life priorities was reported (70%) with increase in spirituality, and the lack of respect of life values generated anger (50%). Self-centered attitudes were reported (60%). TX not only comes with positive physical benefits, but also with positive existential values and psychological transformation, and the development of a more altruistic attitude and humanistic values.

Incongruence between test statistics and P values in medical papers

Given an observed test statistic and its degrees of freedom, one may compute the observed P value with most statistical packages. It is unknown to what extent test statistics and P values are congruent in published medical papers. Methods:We checked the congruence of statistical results reported in all the papers of volumes 409–412 of Nature (2001) and a random sample of 63 results from volumes 322–323 of BMJ (2001). We also tested whether the frequencies of the last digit of a sample of 610 test statistics deviated from a uniform distribution (i.e., equally probable digits).Results: 11.6% (21 of 181) and 11.1% (7 of 63) of the statistical results published in Nature and BMJ respectively during 2001 were incongruent, probably mostly due to rounding, transcription, or type-setting errors. At least one such error appeared in 38% and 25% of the papers of Nature and BMJ, respectively. In 12% of the cases, the significance level might change one or more orders of magnitude. The frequencies of the last digit of statistics deviated from the uniform distribution and suggested digit preference in rounding and reporting.Conclusions: this incongruence of test statistics and P values is another example that statistical practice is generally poor, even in the most renowned scientific journals, and that quality of papers should be more controlled and valued

Hexyl-aminolevulinate mediated photodynamic therapy: how to spare normal urothelium. An in vitro approach

Résumé Objectifs : La thérapie photodynamique a pour but la destruction sélective du tissu néoplasique par interaction de lumière, d'oxygène et d'une substance photosensibilisatrice (la Protoporphyrine IX dans notre étude). Malgré une accumulation sélective du photosensibilisateur dans le tissu tumoral, la thérapie photodynamique du carcinome urothélial de la vessie peut endommager les cellules normales de l'épithélium urinaire. La prévention de ces lésions est importante pour la régénération de la muqueuse. Notre étude sur un modèle in vitro d'urothélium porcin étudie l'influence de la concentration du photosensibilisateur, des paramètres d'irradiation et de la production d'intermédiaires réactifs de l'oxygène (ROS) sur les effets photodynamique. Le but était de déterminer les conditions seuil pour épargner l'urothélium sain. Méthode: Dans une chambre de culture transparente à deux compartiments, des muqueuses vésicales de porc maintenues en vie ont été incubées avec une solution d'hexyl-aminolévulinate (HAL), le précurseur de la Protoporphyrine IX. Ces muqueuses ont ensuite été irradiées avec des doses lumineuses croissantes en lumière bleue et en lumière blanche, et les altérations cellulaires ont été évaluées par microscopie électronique à balayage et par un colorant fluorescent, le Sytox green. Nous avons également évalué la production d'intermédiaires réactifs de l'oxygène parla mesure de la fluorescence intracellulaire de Rhodamine 123 (R123), produit de l'oxydation de la Dihydrorhodamine 123 (DHR123) non fluorescente. Ces valeurs ont été corrélées avec celles du photo blanchiment de la PAIX. Résultats : Le taux de mortalité cellulaire était dépendant de la concentration de PAIX. Après 3 heures d'incubation, la valeur seuil de dose lumineuse pour la lumière bleu était de 0.15 et 0.75 J/cm2 (irradiance 30 et 75 mW/cm2, respectivement) et pour la lumière blanche de 0.55 J/cm2 (irradiante 30 mW/cm2). Le taux de photo blanchiment était inversement proportionnel à l'irradiante. Le système de détection des intermédiaires réactifs de l'oxygène DHR123/R123 a démontré une bonne corrélation avec les valeurs seuil pour toutes les conditions d'irradiation utilisées. Conclusions : Nous avons déterminé les doses lumineuses permettant d'épargner 50% des cellules urothéliales saines. L'utilisation d'une faible irradiante associée à des systèmes permettant de mesurer la production d'intermédiaires réactifs de l'oxygène dans les tissus irradiés pourrait améliorer la dosimétrie in vivo et l'efficacité de la thérapie photodynamique. Abstract Background and Objectives: Photodynamic therapy of superficial bladder cancer may cause damages to the normal surrounding bladder wall. Prevention of these is important for bladder healing. We studied the influence of photosensitizes concentration, irradiation parameters and production of reactive oxygen species (ROS) on the photodynamically induced damage in the porcine urothelium in vitro. The aim was to determine the threshold conditions for the cell survival. Methods: Living porcine bladder mucosae were incubated with solution of hexylester of 5-aminolevulinic acid (HAL). The mucosae were irradiated with increasing doses and cell alterations were evaluated by scanning electron microscopy and by Sytox green fluorescence. The urothelial survival score was correlated with Protoporphyrin IX (PpIX) photobleaching and intracellular fluorescence of Rhodamine 123 reflecting the ROS production. Results: The mortality ratio was dependent on PpIX concentration. After 3 hours of incubation, the threshold radiant exposures for blue light were 0.15 and 0.75 J/cm2 (irradiance 30 and 75 mW/cm2, respectively) and for white light 0.55 J/cm2 (irradiance 30 mW/cm2). Photobleaching rate increased with decreasing irradiance. Interestingly, the DHR123/R123 reporter system correlated well with the threshold exposures under all conditions used. Conclusions: we have determined radiant exposures sparing half of normal urothelial cells. We propose that the use of low irradiance combined with systems reporting the ROS production in the irradiated tissue could improve the in vivo dosimetry and optimize the PDT.

Sensitivity of genome-wide-association signals to phenotyping strategy: the PROP-TAS2R38 taste association as a benchmark.

Natural genetic variation can have a pronounced influence on human taste perception, which in turn may influence food preference and dietary choice. Genome-wide association studies represent a powerful tool to understand this influence. To help optimize the design of future genome-wide-association studies on human taste perception we have used the well-known TAS2R38-PROP association as a tool to determine the relative power and efficiency of different phenotyping and data-analysis strategies. The results show that the choice of both data collection and data processing schemes can have a very substantial impact on the power to detect genotypic variation that affects chemosensory perception. Based on these results we provide practical guidelines for the design of future GWAS studies on chemosensory phenotypes. Moreover, in addition to the TAS2R38 gene past studies have implicated a number of other genetic loci to affect taste sensitivity to PROP and the related bitter compound PTC. None of these other locations showed genome-wide significant associations in our study. To facilitate further, target-gene driven, studies on PROP taste perception we provide the genome-wide list of p-values for all SNPs genotyped in the current study.

From Leadership to Citizenship Behavior in Public Organizations: When Values Matter

After decades of management reforms in the public sector, questions on the impact of leader-ship behavior in public organizations have been attracting increasing attention. This paper investigates the relationship between transformational leadership behavior and organizational citizenship behavior as one major extra-role outcome of transformational leadership. Refer-ring to a growing body of research that shows the importance of public service values and employee identification in public administration research, we include public service motiva-tion and organizational goal clarification as mediating variables in our analysis. Structural equation modeling is applied as the method of analysis for a sample of 569 public managers at the local level of Switzerland. The findings of our study support the assumed indirect relation-ship between leadership and employee behavior and emphasize the relevance of public ser-vice values when analyzing leadership behavior in public sector organizations.

Tuberculin skin test and interferon-gamma release assay values are associated with antimicrobial peptides expression in polymorphonuclear cells during latent tuberculous infection

It has been reported that patients with progressive tuberculosis (TB) express abundant amounts of the antimicrobial peptides (AMPs) cathelicidin (LL-37) and human neutrophil peptide-1 (HNP-1) in circulating cells, whereas latent TB infected donors showed no differences when compared with purified protein derivative (PPD) and QuantiFERON®-TB Gold (QFT)-healthy individuals. The aim of this study was to determine whether LL-37 and HNP-1 production correlates with higher tuberculin skin test (TST) and QFT values in TB household contacts. Twenty-six TB household contact individuals between 26-58 years old TST and QFT positive with at last two years of latent TB infection were recruited. AMPs production by polymorphonuclear cells was determined by flow cytometry and correlation between TST and QFT values was analysed. Our results showed that there is a positive correlation between levels of HNP-1 and LL-37 production with reactivity to TST and/or QFT levels. This preliminary study suggests the potential use of the expression levels of these peptides as biomarkers for progression in latent infected individuals.

Hierarchic species-area relationships and the management of forest habitat islands in intensive farmland

Habitat loss and fragmentation due to land use changes are major threats to biodiversity in forest ecosystems, and they are expected to have important impacts on many taxa and at various spatial scales. Species richness and area relationships (SARs) have been used to assess species diversity patterns and drivers, and thereby in the establishment of conservation and management strategies. Here we propose a hierarchical approach to achieve deeper insights on SARs in small forest islets in intensive farmland and to address the impacts of decreasing naturalness on such relationships. In the intensive dairy landscapes of Northwest Portugal, where small forest stands (dominated by pines, eucalypts or both) represent semi-natural habitat islands, 50 small forest stands were selected and surveyed for vascular plant diversity. A hierarchical analytical framework was devised to determine species richness and inter- and intra-patch SARs for the whole set of forest patches (general patterns) and for each type of forest (specific patterns). Differences in SARs for distinct groups were also tested by considering subsets of species (native, alien, woody, and herbaceous). Overall, values for species richness were confirmed to be different between forest patches exhibiting different levels of naturalness. Whereas higher values of plant diversity were found in pine stands, higher values for alien species were observed in eucalypt stands. Total area of forest (inter-patch SAR) was found not to have a significant impact on species richness for any of the targeted groups of species. However, significant intra-patch SARs were obtained for all groups of species and forest types. A hierarchical approach was successfully applied to scrutinise SARs along a gradient of forest naturalness in intensively managed landscapes. Dominant canopy tree and management intensity were found to reflect differently on distinct species groups as well as to compensate for increasing stand area, buffering SARs among patches, but not within patches. Thus, the maintenance of small semi-natural patches dominated by pines, under extensive practices of forest management, will promote native plant diversity while at the same time contributing to limit the expansion of problematic alien invasive species.

Modelling world investment markets using threshold conditional correlation models

In this paper we propose a parsimonious regime-switching approach to model the correlations between assets, the threshold conditional correlation (TCC) model. This method allows the dynamics of the correlations to change from one state (or regime) to another as a function of observable transition variables. Our model is similar in spirit to Silvennoinen and Teräsvirta (2009) and Pelletier (2006) but with the appealing feature that it does not suffer from the course of dimensionality. In particular, estimation of the parameters of the TCC involves a simple grid search procedure. In addition, it is easy to guarantee a positive definite correlation matrix because the TCC estimator is given by the sample correlation matrix, which is positive definite by construction. The methodology is illustrated by evaluating the behaviour of international equities, govenrment bonds and major exchange rates, first separately and then jointly. We also test and allow for different parts in the correlation matrix to be governed by different transition variables. For this, we estimate a multi-threshold TCC specification. Further, we evaluate the economic performance of the TCC model against a constant conditional correlation (CCC) estimator using a Diebold-Mariano type test. We conclude that threshold correlation modelling gives rise to a significant reduction in portfolio´s variance.

Markov chain montecarlo method applied to rounding zeros of compositional data: first approach

As stated in Aitchison (1986), a proper study of relative variation in a compositional data set should be based on logratios, and dealing with logratios excludes dealing with zeros. Nevertheless, it is clear that zero observations might be present in real data sets, either because the corresponding part is completelyabsent –essential zeros– or because it is below detection limit –rounded zeros. Because the second kind of zeros is usually understood as “a trace too small to measure”, it seems reasonable to replace them by a suitable small value, and this has been the traditional approach. As stated, e.g. by Tauber (1999) and byMartín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000), the principal problem in compositional data analysis is related to rounded zeros. One should be careful to use a replacement strategy that does not seriously distort the general structure of the data. In particular, the covariance structure of the involvedparts –and thus the metric properties– should be preserved, as otherwise further analysis on subpopulations could be misleading. Following this point of view, a non-parametric imputation method isintroduced in Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000). This method is analyzed in depth by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2003) where it is shown that thetheoretical drawbacks of the additive zero replacement method proposed in Aitchison (1986) can be overcome using a new multiplicative approach on the non-zero parts of a composition. The new approachhas reasonable properties from a compositional point of view. In particular, it is “natural” in the sense thatit recovers the “true” composition if replacement values are identical to the missing values, and it is coherent with the basic operations on the simplex. This coherence implies that the covariance structure of subcompositions with no zeros is preserved. As a generalization of the multiplicative replacement, in thesame paper a substitution method for missing values on compositional data sets is introduced