Perturbation expansions of multilocus fixation probabilities for frequency-dependent selection with applications to the Hill-Robertson effect and to the joint evolution of helping and punishment.

Natural populations are of finite size and organisms carry multilocus genotypes. There are, nevertheless, few results on multilocus models when both random genetic drift and natural selection affect the evolutionary dynamics. In this paper we describe a formalism to calculate systematic perturbation expansions of moments of allelic states around neutrality in populations of constant size. This allows us to evaluate multilocus fixation probabilities (long-term limits of the moments) under arbitrary strength of selection and gene action. We show that such fixation probabilities can be expressed in terms of selection coefficients weighted by mean first passages times of ancestral gene lineages within a single ancestor. These passage times extend the coalescence times that weight selection coefficients in one-locus perturbation formulas for fixation probabilities. We then apply these results to investigate the Hill-Robertson effect and the coevolution of helping and punishment. Finally, we discuss limitations and strengths of the perturbation approach. In particular, it provides accurate approximations for fixation probabilities for weak selection regimes only (Ns < or = 1), but it provides generally good prediction for the direction of selection under frequency-dependent selection.

Evaluation of two HbA1c point-of-care analyzers.

BACKGROUND Measurement of HbA1c is the most important parameter to assess glycemic control in diabetic patients. Different point-of-care devices for HbA1c are available. The aim of this study was to evaluate two point-of-care testing (POCT) analyzers (DCA Vantage from Siemens and Afinion from Axis-Shield). We studied the bias and precision as well as interference from carbamylated hemoglobin. METHODS Bias of the POCT analyzers was obtained by measuring 53 blood samples from diabetic patients with a wide range of HbA1c, 4%-14% (20-130 mmol/mol), and comparing the results with those obtained by the laboratory method: HPLC HA 8160 Menarini. Precision was performed by 20 successive determinations of two samples with low 4.2% (22 mmol/mol) and high 9.5% (80 mmol/mol) HbA1c values. The possible interference from carbamylated hemoglobin was studied using 25 samples from patients with chronic renal failure. RESULTS The means of the differences between measurements performed by each POCT analyzer and the laboratory method (95% confidence interval) were: 0.28% (p<0.005) (0.10-0.44) for DCA and 0.27% (p<0.001) (0.19-0.35) for Afinion. Correlation coefficients were: r=0.973 for DCA, and r=0.991 for Afinion. The mean bias observed by using samples from chronic renal failure patients were 0.2 (range -0.4, 0.4) for DCA and 0.2 (-0.2, 0.5) for Afinion. Imprecision results were: CV=3.1% (high HbA1c) and 2.97% (low HbA1c) for DCA, CV=1.95% (high HbA1c) and 2.66% (low HbA1c) for Afinion. CONCLUSIONS Both POCT analyzers for HbA1c show good correlation with the laboratory method and acceptable precision.

Space-time relatedness and Hamilton's rule for long-lasting behaviors in viscous populations.

Genes affect not only the behavior and fitness of their carriers but also that of other individuals. According to Hamilton's rule, whether a mutant gene will spread in the gene pool depends on the effects of its carrier on the fitness of all individuals in the population, each weighted by its relatedness to the carrier. However, social behaviors may affect not only recipients living in the generation of the actor but also individuals living in subsequent generations. In this note, I evaluate space-time relatedness coefficients for localized dispersal. These relatedness coefficients weight the selection pressures on long-lasting behaviors, which stem from a multigenerational gap between phenotypic expression by actors and the resulting environmental feedback on the fitness of recipients. Explicit values of space-time relatedness coefficients reveal that they can be surprisingly large for typical dispersal rates, even for hundreds of generations in the future.

Predictive value of Glasgow Coma Scale after brain trauma: change in trend over the past ten years.

BACKGROUND: Age and the Glasgow Coma Scale (GCS) score on admission are considered important predictors of outcome after traumatic brain injury. We investigated the predictive value of the GCS in a large group of patients whose computerised multimodal bedside monitoring data had been collected over the previous 10 years. METHODS: Data from 358 subjects with head injury, collected between 1992 and 2001, were analysed retrospectively. Patients were grouped according to year of admission. Glasgow Outcome Scores (GOS) were determined at six months. Spearman's correlation coefficients between GCS and GOS scores were calculated for each year. RESULTS: On average 34 (SD: 7) patients were monitored every year. We found a significant correlation between the GCS and GOS for the first five years (overall 1992-1996: r = 0.41; p<0.00001; n = 183) and consistent lack of correlations from 1997 onwards (overall 1997-2001: r = 0.091; p = 0.226; n = 175). In contrast, correlations between age and GOS were in both time periods significant and similar (r = -0.24 v r = -0.24; p<0.002). CONCLUSIONS: The admission GCS lost its predictive value for outcome in this group of patients from 1997 onwards. The predictive value of the GCS should be carefully reconsidered when building prognostic models incorporating multimodality monitoring after head injury.

Training effects with the observer-rated cognitive errors and the coping action patterns scales

The use of observer-rated scales requires that raters be trained until they have become reliable in using the scales. However, few studies properly report how training in using a given rating scale is conducted or indeed how it should be conducted. This study examined progress in interrater reliability over 6 months of training with two observer-rated scales, the Cognitive Errors Rating Scale and the Coping Action Patterns Rating Scale. The evolution of the intraclass correlation coefficients was modeled using hierarchical linear modeling. Results showed an overall training effect as well as effects of the basic training phase and of the rater calibration phase, the latter being smaller than the former. The results are discussed in terms of implications for rater training in psychotherapy research.

Association between obesity indices and cardiovascular risk factors in late adolescence in the Seychelles.

ABSTRACT: BACKGROUND: The ability of different obesity indices to predict cardiovascular risk is still debated in youth and few data are available in sub Saharan Africa. We compared the associations between several indices of obesity and cardiovascular risk factors (CVRFs) in late adolescence in the Seychelles. METHODS: We measured body mass index (BMI), waist circumference, waist/hip ratio (WHiR), waist/height ratio (WHtR) and percent fat mass (by bioimpedance) and 6 CVRFs (blood pressure, LDL-cholesterol, HDL-cholesterol, triglycerides, fasting blood glucose and uric acid) in 423 youths aged 19-20 years from the general population. RESULTS: The prevalence of overweight/obesity and several CVRFs was high, with substantial sex differences. Except for glucose in males and LDL-cholesterol in females, all obesity indices were associated with CVRFs. BMI consistently predicted CVRFs at least as well as the other indices. Linear regression on BMI had standardized regression coefficients of 0.25-0.36 for most CVRFs (p<0.01) and ROC analysis had an AUC between 60%-75% for most CVRFs. BMI also predicted well various combinations of CVRFs: 36% of male and 16% of female lean subjects (BMI P90). CONCLUSION: There was an elevated prevalence of obesity and of several CVRFs in youths in Seychelles. BMI predicted single or combined CVRFs at least as well as other simple obesity indices.

Runaway Greenhouse Effect in a Semigray Radiative-Convective Model

The effects of the nongray absorption (i.e., atmospheric opacity varying with wavelength) on the possible upper bound of the outgoing longwave radiation (OLR) emitted by a planetary atmosphere have been examined. This analysis is based on the semigray approach, which appears to be a reasonable compromise between the complexity of nongray models and the simplicity of the gray assumption (i.e., atmospheric absorption independent of wavelength). Atmospheric gases in semigray atmospheres make use of constant absorption coefficients in finite-width spectral bands. Here, such a semigray absorption is introduced in a one-dimensional (1D) radiative– convective model with a stratosphere in radiative equilibrium and a troposphere fully saturated with water vapor, which is the semigray gas. A single atmospheric window in the infrared spectrum has been assumed. In contrast to the single absolute limit of OLR found in gray atmospheres, semigray ones may also show a relative limit. This means that both finite and infinite runaway effects may arise in some semigray cases. Of particular importance is the finding of an entirely new branch of stable steady states that does not appear in gray atmospheres. This new multiple equilibrium is a consequence of the nongray absorption only. It is suspected that this new set of stable solutions has not been previously revealed in analyses of radiative–convective models since it does not appear for an atmosphere with nongray parameters similar to those for the earth’s current state

Learning spatial filters for multispectral image segmentation

We present a novel filtering method for multispectral satellite image classification. The proposed method learns a set of spatial filters that maximize class separability of binary support vector machine (SVM) through a gradient descent approach. Regularization issues are discussed in detail and a Frobenius-norm regularization is proposed to efficiently exclude uninformative filters coefficients. Experiments carried out on multiclass one-against-all classification and target detection show the capabilities of the learned spatial filters.

Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

Select-divide-and-conquer method for large-scale configuration interaction

A select-divide-and-conquer variational method to approximate configuration interaction (CI) is presented. Given an orthonormal set made up of occupied orbitals (Hartree-Fock or similar) and suitable correlation orbitals (natural or localized orbitals), a large N-electron target space S is split into subspaces S0,S1,S2,...,SR. S0, of dimension d0, contains all configurations K with attributes (energy contributions, etc.) above thresholds T0={T0egy, T0etc.}; the CI coefficients in S0 remain always free to vary. S1 accommodates KS with attributes above T1≤T0. An eigenproblem of dimension d0+d1 for S0+S 1 is solved first, after which the last d1 rows and columns are contracted into a single row and column, thus freezing the last d1 CI coefficients hereinafter. The process is repeated with successive Sj(j≥2) chosen so that corresponding CI matrices fit random access memory (RAM). Davidson's eigensolver is used R times. The final energy eigenvalue (lowest or excited one) is always above the corresponding exact eigenvalue in S. Threshold values {Tj;j=0, 1, 2,...,R} regulate accuracy; for large-dimensional S, high accuracy requires S 0+S1 to be solved outside RAM. From there on, however, usually a few Davidson iterations in RAM are needed for each step, so that Hamiltonian matrix-element evaluation becomes rate determining. One μhartree accuracy is achieved for an eigenproblem of order 24 × 106, involving 1.2 × 1012 nonzero matrix elements, and 8.4×109 Slater determinants

A fractionally integrated approach to monetary policy and inflation dynamics

This paper relaxes the standard I(0) and I(1) assumptions typically stated in the monetary VAR literature by considering a richer framework that encompasses the previous two processes as well as other fractionally integrated possibilities. First, a timevarying multivariate spectrum is estimated for post WWII US data. Then, a structural fractionally integrated VAR (VARFIMA) is fitted to each of the resulting time dependent spectra. In this way, both the coefficients of the VAR and the innovation variances are allowed to evolve freely. The model is employed to analyze inflation persistence and to evaluate the stance of US monetary policy. Our findings indicate a strong decline in the innovation variances during the great disinflation, consistent with the view that the good performance of the economy during the 80’s and 90’s is in part a tale of good luck. However, we also find evidence of a decline in inflation persistence together with a stronger monetary response to inflation during the same period. This last result suggests that the Fed may still play a role in accounting for the observed differences in the US inflation history. Finally, we conclude that previous evidence against drifting coefficients could be an artifact of parameter restriction towards the stationary region. Keywords: monetary policy, inflation persistence, fractional integration, timevarying coefficients, VARFIMA. JEL Classification: E52, C32

The evolution of XY recombination: sexually antagonistic selection versus deleterious mutation load.

Recombination arrest between X and Y chromosomes, driven by sexually antagonistic genes, is expected to induce their progressive differentiation. However, in contrast to birds and mammals (which display the predicted pattern), most cold-blooded vertebrates have homomorphic sex chromosomes. Two main hypotheses have been proposed to account for this, namely high turnover rates of sex-determining systems and occasional XY recombination. Using individual-based simulations, we formalize the evolution of XY recombination (here mediated by sex reversal; the "fountain-of-youth" model) under the contrasting forces of sexually antagonistic selection and deleterious mutations. The shift between the domains of elimination and accumulation occurs at much lower selection coefficients for the Y than for the X. In the absence of dosage compensation, mildly deleterious mutations accumulating on the Y depress male fitness, thereby providing incentives for XY recombination. Under our settings, this occurs via "demasculinization" of the Y, allowing recombination in XY (sex-reversed) females. As we also show, this generates a conflict with the X, which coevolves to oppose sex reversal. The resulting rare events of XY sex reversal are enough to purge the Y from its load of deleterious mutations. Our results support the "fountain of youth" as a plausible mechanism to account for the maintenance of sex-chromosome homomorphy.

Epidemiology of inflammatory bowel disease: Is there a shift towards onset at a younger age?

OBJECTIVES: Increasing numbers of paediatric and adolescent patients with Crohn disease (CD) and ulcerative colitis (UC) are reported. To determine whether this observation is a consequence of a shift towards onset at a younger age, we analysed retrospective data from patients enrolled in the Swiss IBD Cohort Study (SIBDCS). PATIENTS AND METHODS: The SIBDCS is a disease-based cohort in Switzerland, which collects retrospective and prospective data on a large sample of patients with inflammatory bowel disease (IBD). Patients, diagnosed from 1980, were stratified according to diagnosis of CD and UC. Age at disease onset (age at first symptoms and age at diagnosis) was analysed in relation to calendar year of disease onset. Data were extracted from physician and patient questionnaires. Linear regressions of age at disease onset by calendar year of disease onset adjusted by sex, country of birth, and education were performed. RESULTS: Adjusted regression coefficients for CD and UC were significantly positive, that is, age at disease onset has increased with time. Male sex was associated with an increase in age at disease onset, and birth in Switzerland with a decrease. These associations were statistically significant. CONCLUSIONS: The results from the SIBDCS do not support the hypothesis that disease onset of both CD and UC occur today at a younger age. On the contrary, our results show that there is a significant trend for age at disease onset occurring at an older age today as compared with recent decades. We conclude that the observation of increasing numbers of paediatric and adolescent patients with IBD is not caused by a trend towards disease onset at a younger age, but that this may rather be a consequence of the overall increasing incidence of these conditions.

Evidence for validity and reliability of a french version of the FAAM.

BACKGROUND: The Foot and Ankle Ability Measure (FAAM) is a self reported questionnaire for patients with foot and ankle disorders available in English, German, and Persian. This study plans to translate the FAAM from English to French (FAAM-F) and assess the validity and reliability of this new version.METHODS: The FAAM-F Activities of Daily Living (ADL) and sports subscales were completed by 105 French-speaking patients (average age 50.5 years) presenting various chronic foot and ankle disorders. Convergent and divergent validity was assessed by Pearson's correlation coefficients between the FAAM-F subscales and the SF-36 scales: Physical Functioning (PF), Physical Component Summary (PCS), Mental Health (MH) and Mental Component Summary (MCS). Internal consistency was calculated by Cronbach's Alpha (CA). To assess test re-test reliability, 22 patients filled out the questionnaire a second time to estimate minimal detectable changes (MDC) and intraclass correlation coefficients (ICC).RESULTS: Correlations for FAAM-F ADL subscale were 0.85 with PF, 0.81 with PCS, 0.26 with MH, 0.37 with MCS. Correlations for FAAM-F Sports subscale were 0.72 with PF, 0.72 with PCS, 0.21 with MH, 0.29 with MCS. CA estimates were 0.97 for both subscales. Respectively for the ADL and Sports subscales, ICC were 0.97 and 0.94, errors for a single measure were 8 and 10 points at 95% confidence and the MDC values at 95% confidence were 7 and 18 points.CONCLUSION: The FAAM-F is valid and reliable for the self-assessment of physical function in French-speaking patients with a wide range of chronic foot and ankle disorders.

Simulation methods with extended stability for stiff biochemical kinetics

Background: With increasing computer power, simulating the dynamics of complex systems in chemistry and biology is becoming increasingly routine. The modelling of individual reactions in (bio)chemical systems involves a large number of random events that can be simulated by the stochastic simulation algorithm (SSA). The key quantity is the step size, or waiting time, τ, whose value inversely depends on the size of the propensities of the different channel reactions and which needs to be re-evaluated after every firing event. Such a discrete event simulation may be extremely expensive, in particular for stiff systems where τ can be very short due to the fast kinetics of some of the channel reactions. Several alternative methods have been put forward to increase the integration step size. The so-called τ-leap approach takes a larger step size by allowing all the reactions to fire, from a Poisson or Binomial distribution, within that step. Although the expected value for the different species in the reactive system is maintained with respect to more precise methods, the variance at steady state can suffer from large errors as τ grows. Results: In this paper we extend Poisson τ-leap methods to a general class of Runge-Kutta (RK) τ-leap methods. We show that with the proper selection of the coefficients, the variance of the extended τ-leap can be well-behaved, leading to significantly larger step sizes.Conclusions: The benefit of adapting the extended method to the use of RK frameworks is clear in terms of speed of calculation, as the number of evaluations of the Poisson distribution is still one set per time step, as in the original τ-leap method. The approach paves the way to explore new multiscale methods to simulate (bio)chemical systems.