PyRate: a new program to estimate speciation and extinction rates from incomplete fossil data

Despite the advancement of phylogenetic methods to estimate speciation and extinction rates, their power can be limited under variable rates, in particular for clades with high extinction rates and small number of extant species. Fossil data can provide a powerful alternative source of information to investigate diversification processes. Here, we present PyRate, a computer program to estimate speciation and extinction rates and their temporal dynamics from fossil occurrence data. The rates are inferred in a Bayesian framework and are comparable to those estimated from phylogenetic trees. We describe how PyRate can be used to explore different models of diversification. In addition to the diversification rates, it provides estimates of the parameters of the preservation process (fossilization and sampling) and the times of speciation and extinction of each species in the data set. Moreover, we develop a new birth-death model to correlate the variation of speciation/extinction rates with changes of a continuous trait. Finally, we demonstrate the use of Bayes factors for model selection and show how the posterior estimates of a PyRate analysis can be used to generate calibration densities for Bayesian molecular clock analysis. PyRate is an open-source command-line Python program available at http://sourceforge.net/projects/pyrate/.

Learning through Computer Model Improvisations

It has been convincingly argued that computer simulation modeling differs from traditional science. If we understand simulation modeling as a new way of doing science, the manner in which scientists learn about the world through models must also be considered differently. This article examines how researchers learn about environmental processes through computer simulation modeling. Suggesting a conceptual framework anchored in a performative philosophical approach, we examine two modeling projects undertaken by research teams in England, both aiming to inform flood risk management. One of the modeling teams operated in the research wing of a consultancy firm, the other were university scientists taking part in an interdisciplinary project experimenting with public engagement. We found that in the first context the use of standardized software was critical to the process of improvisation, the obstacles emerging in the process concerned data and were resolved through exploiting affordances for generating, organizing, and combining scientific information in new ways. In the second context, an environmental competency group, obstacles were related to the computer program and affordances emerged in the combination of experience-based knowledge with the scientists' skill enabling a reconfiguration of the mathematical structure of the model, allowing the group to learn about local flooding.

An objective anisotropic elastic plastic model and Algorithm applicable to bone mechanics

The implicit projection algorithm of isotropic plasticity is extended to an objective anisotropic elastic perfectly plastic model. The recursion formula developed to project the trial stress on the yield surface, is applicable to any non linear elastic law and any plastic yield function.A curvilinear transverse isotropic model based on a quadratic elastic potential and on Hill's quadratic yield criterion is then developed and implemented in a computer program for bone mechanics perspectives. The paper concludes with a numerical study of a schematic bone-prosthesis system to illustrate the potential of the model.

Physicians' estimates of the 10 year cardiovascular risk in hypertensive patients: an evaluation in primary care physicians in training.

To evaluate how young physicians in training perceive their patients' cardiovascular risk based on the medical charts and their clinical judgment. Cross sectional observational study. University outpatient clinic, Lausanne, Switzerland. Two hundred hypertensive patients and 50 non-hypertensive patients with at least one cardiovascular risk factor. Comparison of the absolute 10-year cardiovascular risk calculated by a computer program based on the Framingham score and adapted for physicians by the WHO/ISH with the perceived risk as assessed clinically by the physicians. Physicians underestimated the 10-year cardiovascular risk of their patients compared to that calculated with the Framingham score. Concordance between methods was 39% for hypertensive patients and 30% for non-hypertensive patients. Underestimation of cardiovascular risks for hypertensive patients was related to the fact they had a stabilized systolic blood pressure under 140 mm Hg (OR = 2.1 [1.1; 4.1]). These data show that young physicians in training often have an incorrect perception of the cardiovascular risk of their patients with a tendency to underestimate the risk. However, the calculated risk could also be slightly overestimated when applying the Framingham Heart Study model to a Swiss population. To implement a systematic evaluation of risk factors in primary care a greater emphasis should be placed on the teaching of cardiovascular risk evaluation and on the implementation of quality improvement programs.

The Early Toarcian anoxia, a synchronous event in the Western Tethys? An approach by quantitative biochronology (Unitary Associations), applied on calcareous nannofossils

During the Early Toarcian, major paleoenvironnemental and paleoceanographical changes occurred, leading to an oceanic anoxic event (OAE) and to a perturbation of the carbon isotope cycle. Although the standard biochronology of the Lower Jurassic is essentially based upon ammonites, in recent years biostratigraphy based on calcareous nannofossils and dinoflagellate cysts is increasingly used to date Jurassic rocks. However, the precise dating and correlation of the Early Toarcian OAE, and of the associated delta C-13 anomaly in different settings of the western Tethys, are still partly problematic, and it is still unclear whether these events are synchronous or not. In order to allow more accurate correlations of the organic rich levels recorded in the Lower Toarcian OAE, this account proposes a new biozonation based on a quantitative biochronology approach, the Unitary Associations (UA), applied to calcareous nannofossils. This study represents the first attempt to apply the UA method to Jurassic nannofossils. The study incorporates eighteen sections distributed across western Tethys and ranging from the Pliensbachian to Aalenian, comprising 1220 samples and 72 calcareous nannofossil taxa. The BioGraph [Savary, J., Guex, J., 1999. Discrete biochronological scales and unitary associations: description of the Biograph Computer program. Memoires de Geologie de Lausanne 34, 282 pp] and UA-Graph (Copyright Hammer O., Guex and Savary, 2002) softwares provide a discrete biochronological framework based upon multi-taxa concurrent range zones in the different sections. The optimized dataset generates nine UAs using the co-occurrences of 56 taxa. These UAs are grouped into six Unitary Association Zones (UA-Z), which constitute a robust biostratigraphic synthesis of all the observed or deduced biostratigraphic relationships between the analysed taxa. The UA zonation proposed here is compared to ``classic'' calcareous nannofossil biozonations, which are commonly used for the southern and the northern sides of Tethys. The biostratigraphic resolution of the UA-Zones varies from one nannofossil subzone or part of it to several subzones, and can be related to the pattern of calcareous nannoplankton originations and extinctions during the studied time interval. The Late Pliensbachian - Early Toarcian interval (corresponding to the UA-Z II) represents a major step in the Jurassic nannoplankton radiation. The recognized UA-Zones are also compared to the carbon isotopic negative excursion and TOC maximum in five sections of central Italy, Germany and England, with the aim of providing a more reliable correlation tool for the Early Toarcian OAE, and of the associated isotopic anomaly, between the southern and northern part of western Tethys. The results of this work show that the TOC maximum and delta C-13 negative excursion correspond to the upper part of the UA-Z II (i.e., UA 3) in the sections analysed. This suggests that the Early Toarcian OAE was a synchronous event within the western Tethys. (c) 2006 Elsevier B.V. All rights reserved.

Technical note: Preliminary estimation of rockfall runout zones

Rockfall propagation areas can be determined using a simple geometric rule known as shadow angle or energy line method based on a simple Coulomb frictional model implemented in the CONEFALL computer program. Runout zones are estimated from a digital terrain model (DTM) and a grid file containing the cells representing rockfall potential source areas. The cells of the DTM that are lowest in altitude and located within a cone centered on a rockfall source cell belong to the potential propagation area associated with that grid cell. In addition, the CONEFALL method allows estimation of mean and maximum velocities and energies of blocks in the rockfall propagation areas. Previous studies indicate that the slope angle cone ranges from 27° to 37° depending on the assumptions made, i.e. slope morphology, probability of reaching a point, maximum run-out, field observations. Different solutions based on previous work and an example of an actual rockfall event are presented here.

Concomitant chemoradiotherapy versus acceleration of radiotherapy with or without concomitant chemotherapy in locally advanced head and neck carcinoma (GORTEC 99-02): an open-label phase 3 randomised trial.

BACKGROUND: Concomitant chemoradiotherapy and accelerated radiotherapy independently improve outcomes for patients with locally advanced head and neck squamous-cell carcinoma (HNSCC). We aimed to assess the efficacy and safety of a combination of these approaches. METHODS: In our open-label phase 3 randomised trial, we enrolled patients with locally advanced, stage III and IV (non-metastatic) HNSCC and an Eastern Cooperative Oncology Group performance status of 0-2. We randomly allocated patients centrally with a computer program (with centre, T stage, N stage, and localisation as minimisation factors) in a 1:1:1 ratio to receive conventional chemoradiotherapy (70 Gy in 7 weeks plus three cycles of 4 days' concomitant carboplatin-fluorouracil), accelerated radiotherapy-chemotherapy (70 Gy in 6 weeks plus two cycles of 5 days' concomitant carboplatin-fluorouracil), or very accelerated radiotherapy alone (64·8 Gy [1·8 Gy twice daily] in 3·5 weeks). The primary endpoint, progression-free survival (PFS), was assessed in all enrolled patients. This trial is completed. The trial is registered with ClinicalTrials.gov, number NCT00828386. FINDINGS: Between Feb 29, 2000, and May 9, 2007, we randomly allocated 279 patients to receive conventional chemoradiotherapy, 280 to accelerated radiotherapy-chemotherapy, and 281 to very accelerated radiotherapy. Median follow-up was 5·2 years (IQR 4·9-6·2); rates of chemotherapy and radiotherapy compliance were good in all groups. Accelerated radiotherapy-chemotherapy offered no PFS benefit compared with conventional chemoradiotherapy (HR 1·02, 95% CI 0·84-1·23; p=0·88) or very accelerated radiotherapy (0·83, 0·69-1·01; p=0·060); conventional chemoradiotherapy improved PFS compared with very accelerated radiotherapy (0·82, 0·67-0·99; p=0·041). 3-year PFS was 37·6% (95% CI 32·1-43·4) after conventional chemoradiotherapy, 34·1% (28·7-39·8) after accelerated radiotherapy-chemotherapy, and 32·2% (27·0-37·9) after very accelerated radiotherapy. More patients in the very accelerated radiotherapy group had RTOG grade 3-4 acute mucosal toxicity (226 [84%] of 268 patients) compared with accelerated radiotherapy-chemotherapy (205 [76%] of 271 patients) or conventional chemoradiotherapy (180 [69%] of 262; p=0·0001). 158 (60%) of 265 patients in the conventional chemoradiotherapy group, 176 (64%) of 276 patients in the accelerated radiotherapy-chemotherapy group, and 190 (70%) of 272 patients in the very accelerated radiotherapy group were intubated with feeding tubes during treatment (p=0·045). INTERPRETATION: Chemotherapy has a substantial treatment effect given concomitantly with radiotherapy and acceleration of radiotherapy cannot compensate for the absence of chemotherapy. We noted the most favourable outcomes for conventional chemoradiotherapy, suggesting that acceleration of radiotherapy is probably not beneficial in concomitant chemoradiotherapy schedules. FUNDING: French Ministry of Health.

Limited genetic diversity and high differentiation among the remnant adder (Vipera berus) populations in the Swiss and French Jura Mountains

Although the adder (Vipera berus) has a large distribution area, this species is particularly threatened in Western Europe due to high habitat fragmentation and human persecution. We developed 13 new microsatellite markers in order to evaluate population structure and genetic diversity in the Swiss and French Jura Mountains, where the species is limited to only a few scattered populations. We found that V. berus exhibits a considerable genetic differentiation among populations (global F-ST = 0.269), even if these are not geographically isolated. Moreover, the genetic diversity within populations in the Jura Mountains and in the less perturbed Swiss Alps is significantly lower than in other French populations, possibly due to post-glacial recolonisation processes. Finally, in order to minimize losses of genetic diversities within isolated populations, suggestions for the conservation of this species in fragmented habitats are proposed.

Apprentissage automatique de la morphologie : le cas des structures racine-schème

RÉSUMÉ Cette thèse porte sur le développement de méthodes algorithmiques pour découvrir automatiquement la structure morphologique des mots d'un corpus. On considère en particulier le cas des langues s'approchant du type introflexionnel, comme l'arabe ou l'hébreu. La tradition linguistique décrit la morphologie de ces langues en termes d'unités discontinues : les racines consonantiques et les schèmes vocaliques. Ce genre de structure constitue un défi pour les systèmes actuels d'apprentissage automatique, qui opèrent généralement avec des unités continues. La stratégie adoptée ici consiste à traiter le problème comme une séquence de deux sous-problèmes. Le premier est d'ordre phonologique : il s'agit de diviser les symboles (phonèmes, lettres) du corpus en deux groupes correspondant autant que possible aux consonnes et voyelles phonétiques. Le second est de nature morphologique et repose sur les résultats du premier : il s'agit d'établir l'inventaire des racines et schèmes du corpus et de déterminer leurs règles de combinaison. On examine la portée et les limites d'une approche basée sur deux hypothèses : (i) la distinction entre consonnes et voyelles peut être inférée sur la base de leur tendance à alterner dans la chaîne parlée; (ii) les racines et les schèmes peuvent être identifiés respectivement aux séquences de consonnes et voyelles découvertes précédemment. L'algorithme proposé utilise une méthode purement distributionnelle pour partitionner les symboles du corpus. Puis il applique des principes analogiques pour identifier un ensemble de candidats sérieux au titre de racine ou de schème, et pour élargir progressivement cet ensemble. Cette extension est soumise à une procédure d'évaluation basée sur le principe de la longueur de description minimale, dans- l'esprit de LINGUISTICA (Goldsmith, 2001). L'algorithme est implémenté sous la forme d'un programme informatique nommé ARABICA, et évalué sur un corpus de noms arabes, du point de vue de sa capacité à décrire le système du pluriel. Cette étude montre que des structures linguistiques complexes peuvent être découvertes en ne faisant qu'un minimum d'hypothèses a priori sur les phénomènes considérés. Elle illustre la synergie possible entre des mécanismes d'apprentissage portant sur des niveaux de description linguistique distincts, et cherche à déterminer quand et pourquoi cette coopération échoue. Elle conclut que la tension entre l'universalité de la distinction consonnes-voyelles et la spécificité de la structuration racine-schème est cruciale pour expliquer les forces et les faiblesses d'une telle approche. ABSTRACT This dissertation is concerned with the development of algorithmic methods for the unsupervised learning of natural language morphology, using a symbolically transcribed wordlist. It focuses on the case of languages approaching the introflectional type, such as Arabic or Hebrew. The morphology of such languages is traditionally described in terms of discontinuous units: consonantal roots and vocalic patterns. Inferring this kind of structure is a challenging task for current unsupervised learning systems, which generally operate with continuous units. In this study, the problem of learning root-and-pattern morphology is divided into a phonological and a morphological subproblem. The phonological component of the analysis seeks to partition the symbols of a corpus (phonemes, letters) into two subsets that correspond well with the phonetic definition of consonants and vowels; building around this result, the morphological component attempts to establish the list of roots and patterns in the corpus, and to infer the rules that govern their combinations. We assess the extent to which this can be done on the basis of two hypotheses: (i) the distinction between consonants and vowels can be learned by observing their tendency to alternate in speech; (ii) roots and patterns can be identified as sequences of the previously discovered consonants and vowels respectively. The proposed algorithm uses a purely distributional method for partitioning symbols. Then it applies analogical principles to identify a preliminary set of reliable roots and patterns, and gradually enlarge it. This extension process is guided by an evaluation procedure based on the minimum description length principle, in line with the approach to morphological learning embodied in LINGUISTICA (Goldsmith, 2001). The algorithm is implemented as a computer program named ARABICA; it is evaluated with regard to its ability to account for the system of plural formation in a corpus of Arabic nouns. This thesis shows that complex linguistic structures can be discovered without recourse to a rich set of a priori hypotheses about the phenomena under consideration. It illustrates the possible synergy between learning mechanisms operating at distinct levels of linguistic description, and attempts to determine where and why such a cooperation fails. It concludes that the tension between the universality of the consonant-vowel distinction and the specificity of root-and-pattern structure is crucial for understanding the advantages and weaknesses of this approach.

Urinary saturation and nephrocalcinosis in preterm infants: effect of parenteral nutrition.

Urinary lithogenic and inhibitory factors were studied in 27 preterm infants; 16 had total parenteral nutrition (TPN) and 11 had breastmilk with an additional glucose-sodium chloride infusion. Urines were collected for 24 hours on day 2 (period A), day 3 (B), and once between days 4 and 10 (C). Urinary calcium oxalate saturation was calculated by the computer program EQUIL 2. Renal ultrasonography was performed every second week until discharge. The calcium/creatinine ratio increased in infants on TPN (A 0.91; C 1.68 mol/mol) and was significantly higher at period C than that in infants on breastmilk/infusion (A 0.52; C 0.36). The oxalate/creatinine ratio was persistently higher with TPN (203 mmol/mol) than with breastmilk/infusion (98; 137). The citrate/creatinine remained constant with TPN (0.44 mol/mol), whereas it increased significantly with breastmilk/infusion (0.26; 0.49). Calcium/citrate rose considerably with TPN, but decreased with breastmilk/infusion to a significantly lower level than with TPN. The urinary calcium oxalate saturation increased with TPN (2.4; 4.5) and decreased with breastmilk/infusion (2.1; 1.5) to a significantly lower value than with TPN. Nephrocalcinosis developed in two infants on TPN. Mean daily calcium intake was similar in both groups, whereas protein, sodium, and phosphorus intake were significantly higher on TPN. It is concluded that the increase in urinary calcium oxalate saturation observed with TPN is due to the combined effect of an increased urinary calcium excretion and higher urinary oxalate/creatinine and calcium/citrate ratios. The changes observed are likely to be caused by TPN itself, which differs in several respects from breastmilk feeding.

Illite "crystallinity'' revisited

The Kubler Index (KI) is defined as the full width at half-maximum height (FWHM) of the 10-Angstrom X-ray diffraction peak of illite-smectite interstratified (I-S) clay minerals, The only parameters controlling the Kubler index are assumed to be the mean number of layers (N) in the coherent scattering domains (CSD). the variance of the distribution of the number of layers of the CSB, the mean percentage of smectite layers in I-S (%S), and the probability of layer stocking (Reichweite). The Kubler-Index measurements on ah-dried (KIAD) and ethylene-glycolated (KIEG) samples were compared to N and %S using the NEWMOD computer program to simulate X-ray diffraction patterns. Charts of KIAD versus KIEG corrected for instrumental broadening were made and isolines were mapped for constant N and CLS. Isolines allow a direct and rapid determination of N and RS from KI measurements. The method allows quantification of the metamorphic anchizone limits by considering mean thickness of fundamental particles in MacEwan crystallites. The transition from diagenesis to the anchizone and from the anchizone to the epizone of low-grade metamorphism corresponds to thicknesses of 20- and 70-layer fundamental particles, respectively.

A randomized controlled trial of the effectiveness of a Computer-Assisted Cognitive Remediation (CACR) program in adolescents with psychosis or at high risk of psychosis: Short term and long term outcomes

The present study investigates the short- and long-term outcomes of a computer-assisted cognitive remediation (CACR) program in adolescents with psychosis or at high risk. 32 adolescents participated in a blinded 8-week randomized controlled trial of CACR treatment compared to computer games (CG). Clinical and neuropsychological evaluations were undertaken at baseline, at the end of the program and at 6-month. At the end of the program (n = 28), results indicated that visuospatial abilities (Repeatable Battery for the Assessment of Neuropsychological Status, RBANS; P = .005) improved signifi cantly more in the CACR group compared to the CG group. Furthermore, other cognitive functions (RBANS), psychotic symptoms (Positive and Negative Symptom Scale) and psychosocial functioning (Social and Occupational Functioning Assessment Scale) improved signifi cantly, but at similar rates, in the two groups. At long term (n = 22), cognitive abilities did not demonstrated any amelioration in the control group while, in the CACR group, signifi cant long-term improvements in inhibition (Stroop; P = .040) and reasoning (Block Design Test; P = .005) were observed. In addition, symptom severity (Clinical Global Improvement) decreased signifi cantly in the control group (P = .046) and marginally in the CACR group (P = .088). To sum up, CACR can be successfully administered in this population. CACR proved to be effective over and above CG for the most intensively trained cognitive ability. Finally, on the long-term, enhanced reasoning and inhibition abilities, which are necessary to execute higher-order goals or to adapt behavior to the ever-changing environment, were observed in adolescents benefi ting from a CACR.