A decision-making Fokker-Planck model in computational neuroscience

Minimal models for the explanation of decision-making in computational neuroscience are based on the analysis of the evolution for the average firing rates of two interacting neuron populations. While these models typically lead to multi-stable scenario for the basic derived dynamical systems, noise is an important feature of the model taking into account finite-size effects and robustness of the decisions. These stochastic dynamical systems can be analyzed by studying carefully their associated Fokker-Planck partial differential equation. In particular, we discuss the existence, positivity and uniqueness for the solution of the stationary equation, as well as for the time evolving problem. Moreover, we prove convergence of the solution to the the stationary state representing the probability distribution of finding the neuron families in each of the decision states characterized by their average firing rates. Finally, we propose a numerical scheme allowing for simulations performed on the Fokker-Planck equation which are in agreement with those obtained recently by a moment method applied to the stochastic differential system. Our approach leads to a more detailed analytical and numerical study of this decision-making model in computational neuroscience.

Evolutionary and convergence stability for continuous phenotypes in finite populations derived from two-allele models.

The evolution of a quantitative phenotype is often envisioned as a trait substitution sequence where mutant alleles repeatedly replace resident ones. In infinite populations, the invasion fitness of a mutant in this two-allele representation of the evolutionary process is used to characterize features about long-term phenotypic evolution, such as singular points, convergence stability (established from first-order effects of selection), branching points, and evolutionary stability (established from second-order effects of selection). Here, we try to characterize long-term phenotypic evolution in finite populations from this two-allele representation of the evolutionary process. We construct a stochastic model describing evolutionary dynamics at non-rare mutant allele frequency. We then derive stability conditions based on stationary average mutant frequencies in the presence of vanishing mutation rates. We find that the second-order stability condition obtained from second-order effects of selection is identical to convergence stability. Thus, in two-allele systems in finite populations, convergence stability is enough to characterize long-term evolution under the trait substitution sequence assumption. We perform individual-based simulations to confirm our analytic results.

Effects of Non-Susceptible Hosts on the Infection with Trypanosoma cruzi of the Vector Triatoma infestans: an Experimental Model

We tested experimentally the effects of the presence of non-susceptible hosts on the infection with Trypanosoma cruzi of the vector Triatoma infestans. The experiment consisted in two treatments: with chickens, including two chickens (non-susceptible hosts) and two infected guinea pigs (susceptible hosts), and without chickens, including only two infected guinea pigs. The hosts were held unrestrained in individual metal cages inside a closed tulle chamber. A total of 200 uninfected T. infestans third instar nymphs were liberated in each replica, collected on day 14, and examined for infection and blood meal sources on day 32-36. The additional presence of chickens relative to infected guinea pigs: (a) significantly modified the spatial distribution of bugs; (b) increased significantly the likelihoods of having a detectable blood meal on any host and molting to the next instar; (c) did not affect the bugs' probability of death by predation; and (d) decreased significantly the overall percentage of T. infestans infected with T. cruzi. The bugs collected from inside or close to the guinea pigs' cages showed a higher infection rate (71-88%) than those collected from the chickens' cages (22-32%). Mixed blood meals on chickens and guinea pigs were detected in 12-21% of bugs. Although the presence of chickens would decrease the overall percentage of infected bugs in short term experiments, the high rate of host change of T. infestans would make this difference fade out if longer exposure times had been provided.

Hybrid Multiscale Finite Volume method for two-phase flow in porous media

We present a novel hybrid (or multiphysics) algorithm, which couples pore-scale and Darcy descriptions of two-phase flow in porous media. The flow at the pore-scale is described by the Navier?Stokes equations, and the Volume of Fluid (VOF) method is used to model the evolution of the fluid?fluid interface. An extension of the Multiscale Finite Volume (MsFV) method is employed to construct the Darcy-scale problem. First, a set of local interpolators for pressure and velocity is constructed by solving the Navier?Stokes equations; then, a coarse mass-conservation problem is constructed by averaging the pore-scale velocity over the cells of a coarse grid, which act as control volumes; finally, a conservative pore-scale velocity field is reconstructed and used to advect the fluid?fluid interface. The method relies on the localization assumptions used to compute the interpolators (which are quite straightforward extensions of the standard MsFV) and on the postulate that the coarse-scale fluxes are proportional to the coarse-pressure differences. By numerical simulations of two-phase problems, we demonstrate that these assumptions provide hybrid solutions that are in good agreement with reference pore-scale solutions and are able to model the transition from stable to unstable flow regimes. Our hybrid method can naturally take advantage of several adaptive strategies and allows considering pore-scale fluxes only in some regions, while Darcy fluxes are used in the rest of the domain. Moreover, since the method relies on the assumption that the relationship between coarse-scale fluxes and pressure differences is local, it can be used as a numerical tool to investigate the limits of validity of Darcy's law and to understand the link between pore-scale quantities and their corresponding Darcy-scale variables.

Experimental evaluation of an electromechanical artificial urinary sphincter in an animal model.

WHAT'S KNOWN ON THE SUBJECT? AND WHAT DOES THE STUDY ADD?: The AMS 800 urinary control system is the gold standard for the treatment of urinary incontinence due to sphincter insufficiency. Despite excellent functional outcome and latest technological improvements, the revision rate remains significant. To overcome the shortcomings of the current device, we developed a modern electromechanical artificial urinary sphincter. The results demonstrated that this new sphincter is effective and well tolerated up to 3 months. This preliminary study represents a first step in the clinical application of novel technologies and an alternative compression mechanism to the urethra. OBJECTIVES: To evaluate the effectiveness in continence achievement of a new electromechanical artificial urinary sphincter (emAUS) in an animal model. To assess urethral response and animal general response to short-term and mid-term activation of the emAUS. MATERIALS AND METHODS: The principle of the emAUS is electromechanical induction of alternating compression of successive segments of the urethra by a series of cuffs activated by artificial muscles. Between February 2009 and May 2010 the emAUS was implanted in 17 sheep divided into three groups. The first phase aimed to measure bladder leak point pressure during the activation of the device. The second and third phases aimed to assess tissue response to the presence of the device after 2-9 weeks and after 3 months respectively. Histopathological and immunohistochemistry evaluation of the urethra was performed. RESULTS: Bladder leak point pressure was measured at levels between 1091 ± 30.6 cmH2 O and 1244.1 ± 99 cmH2 O (mean ± standard deviation) depending on the number of cuffs used. At gross examination, the explanted urethra showed no sign of infection, atrophy or stricture. On microscopic examination no significant difference in structure was found between urethral structure surrounded by a cuff and control urethra. In the peripheral tissues, the implanted material elicited a chronic foreign body reaction. Apart from one case, specimens did not show significant presence of lymphocytes, polymorphonuclear leucocytes, necrosis or cell degeneration. Immunohistochemistry confirmed the absence of macrophages in the samples. CONCLUSIONS: This animal study shows that the emAUS can provide continence. This new electronic controlled sequential alternating compression mechanism can avoid damage to urethral vascularity, at least up to 3 months after implantation. After this positive proof of concept, long-term studies are needed before clinical application could be considered.

Random mating with a finite number of matings.

Random mating is the null model central to population genetics. One assumption behind random mating is that individuals mate an infinite number of times. This is obviously unrealistic. Here we show that when each female mates a finite number of times, the effective size of the population is substantially decreased.

How Low Can We Go in Contrast-Enhanced CT Imaging of the Chest?: A Dose-Finding Cadaver Study Using the Model-based Iterative Image Reconstruction Approach.

RATIONALE AND OBJECTIVES: Dose reduction may compromise patients because of a decrease of image quality. Therefore, the amount of dose savings in new dose-reduction techniques needs to be thoroughly assessed. To avoid repeated studies in one patient, chest computed tomography (CT) scans with different dose levels were performed in corpses comparing model-based iterative reconstruction (MBIR) as a tool to enhance image quality with current standard full-dose imaging. MATERIALS AND METHODS: Twenty-five human cadavers were scanned (CT HD750) after contrast medium injection at different, decreasing dose levels D0-D5 and respectively reconstructed with MBIR. The data at full-dose level, D0, have been additionally reconstructed with standard adaptive statistical iterative reconstruction (ASIR), which represented the full-dose baseline reference (FDBR). Two radiologists independently compared image quality (IQ) in 3-mm multiplanar reformations for soft-tissue evaluation of D0-D5 to FDBR (-2, diagnostically inferior; -1, inferior; 0, equal; +1, superior; and +2, diagnostically superior). For statistical analysis, the intraclass correlation coefficient (ICC) and the Wilcoxon test were used. RESULTS: Mean CT dose index values (mGy) were as follows: D0/FDBR = 10.1 ± 1.7, D1 = 6.2 ± 2.8, D2 = 5.7 ± 2.7, D3 = 3.5 ± 1.9, D4 = 1.8 ± 1.0, and D5 = 0.9 ± 0.5. Mean IQ ratings were as follows: D0 = +1.8 ± 0.2, D1 = +1.5 ± 0.3, D2 = +1.1 ± 0.3, D3 = +0.7 ± 0.5, D4 = +0.1 ± 0.5, and D5 = -1.2 ± 0.5. All values demonstrated a significant difference to baseline (P < .05), except mean IQ for D4 (P = .61). ICC was 0.91. CONCLUSIONS: Compared to ASIR, MBIR allowed for a significant dose reduction of 82% without impairment of IQ. This resulted in a calculated mean effective dose below 1 mSv.

Neonatal and adult ICU ventilators to provide ventilation in neonates, infants, and children: a bench model study.

BACKGROUND: Using a bench test model, we investigated the hypothesis that neonatal and/or adult ventilators equipped with neonatal/pediatric modes currently do not reliably administer pressure support (PS) in neonatal or pediatric patient groups in either the absence or presence of air leaks. METHODS: PS was evaluated in 4 neonatal and 6 adult ventilators using a bench model to evaluate triggering, pressurization, and cycling in both the absence and presence of leaks. Delivered tidal volumes were also assessed. Three patients were simulated: a preterm infant (resistance 100 cm H2O/L/s, compliance 2 mL/cm H2O, inspiratory time of the patient [TI] 400 ms, inspiratory effort 1 and 2 cm H2O), a full-term infant (resistance 50 cm H2O/L/s, compliance 5 mL/cm H2O, TI 500 ms, inspiratory effort 2 and 4 cm H2O), and a child (resistance 30 cm H2O/L/s, compliance 10 mL/cm H2O, TI 600 ms, inspiratory effort 5 and 10 cm H2O). Two PS levels were tested (10 and 15 cm H2O) with and without leaks and with and without the leak compensation algorithm activated. RESULTS: Without leaks, only 2 neonatal ventilators and one adult ventilator had trigger delays under a given predefined acceptable limit (1/8 TI). Pressurization showed high variability between ventilators. Most ventilators showed TI in excess high enough to seriously impair patient-ventilator synchronization (> 50% of the TI of the subject). In some ventilators, leaks led to autotriggering and impairment of ventilation performance, but the influence of leaks was generally lower in neonatal ventilators. When a noninvasive ventilation algorithm was available, this was partially corrected. In general, tidal volume was calculated too low by the ventilators in the presence of leaks; the noninvasive ventilation algorithm was able to correct this difference in only 2 adult ventilators. CONCLUSIONS: No ventilator performed equally well under all tested conditions for all explored parameters. However, neonatal ventilators tended to perform better in the presence of leaks. These findings emphasize the need to improve algorithms for assisted ventilation modes to better deal with situations of high airway resistance, low pulmonary compliance, and the presence of leaks.

Evaluation of a linear spectral mixture model and vegetation indices (NDVI and EVI) in a study of schistosomiasis mansoni and Biomphalaria glabrata distribution in the state of Minas Gerais, Brazil

This paper analyses the associations between Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) on the prevalence of schistosomiasis and the presence of Biomphalaria glabrata in the state of Minas Gerais (MG), Brazil. Additionally, vegetation, soil and shade fraction images were created using a Linear Spectral Mixture Model (LSMM) from the blue, red and infrared channels of the Moderate Resolution Imaging Spectroradiometer spaceborne sensor and the relationship between these images and the prevalence of schistosomiasis and the presence of B. glabrata was analysed. First, we found a high correlation between the vegetation fraction image and EVI and second, a high correlation between soil fraction image and NDVI. The results also indicate that there was a positive correlation between prevalence and the vegetation fraction image (July 2002), a negative correlation between prevalence and the soil fraction image (July 2002) and a positive correlation between B. glabrata and the shade fraction image (July 2002). This paper demonstrates that the LSMM variables can be used as a substitute for the standard vegetation indices (EVI and NDVI) to determine and delimit risk areas for B. glabrata and schistosomiasis in MG, which can be used to improve the allocation of resources for disease control.

Explaining Fiscal Balances with a Simultaneous Equation Model of Revenue and Expenditure : A Case Study of Swiss Cantons Using Panel Data

Empirical literature on the analysis of the efficiency of measures for reducing persistent government deficits has mainly focused on the direct explanation of deficit. By contrast, this paper aims at modeling government revenue and expenditure within a simultaneous framework and deriving the fiscal balance (surplus or deficit) equation as the difference between the two variables. This setting enables one to not only judge how relevant the explanatory variables are in explaining the fiscal balance but also understand their impact on revenue and/or expenditure. Our empirical results, obtained by using a panel data set on Swiss Cantons for the period 1980-2002, confirm the relevance of the approach followed here, by providing unambiguous evidence of a simultaneous relationship between revenue and expenditure. They also reveal strong dynamic components in revenue, expenditure, and fiscal balance. Among the significant determinants of public fiscal balance we not only find the usual business cycle elements, but also and more importantly institutional factors such as the number of administrative units, and the ease with which people can resort to political (direct democracy) instruments, such as public initiatives and referendum.

Inoculum effect on the efficacies of amoxicillin-clavulanate, piperacillin-tazobactam, and imipenem against extended-spectrum β-lactamase (ESBL)-producing and non-ESBL-producing Escherichia coli in an experimental murine sepsis model.

Escherichia coli is commonly involved in infections with a heavy bacterial burden. Piperacillin-tazobactam and carbapenems are among the recommended empirical treatments for health care-associated complicated intra-abdominal infections. In contrast to amoxicillin-clavulanate, both have reduced in vitro activity in the presence of high concentrations of extended-spectrum β-lactamase (ESBL)-producing and non-ESBL-producing E. coli bacteria. Our goal was to compare the efficacy of these antimicrobials against different concentrations of two clinical E. coli strains, one an ESBL-producer and the other a non-ESBL-producer, in a murine sepsis model. An experimental sepsis model {~5.5 log10 CFU/g [low inoculum concentration (LI)] or ~7.5 log(10) CFU/g [high inoculum concentration (HI)]} using E. coli strains ATCC 25922 (non-ESBL producer) and Ec1062 (CTX-M-14 producer), which are susceptible to the three antimicrobials, was used. Amoxicillin-clavulanate (50/12.5 mg/kg given intramuscularly [i.m.]), piperacillin-tazobactam (25/3.125 mg/kg given intraperitoneally [i.p.]), and imipenem (30 mg/kg i.m.) were used. Piperacillin-tazobactam and imipenem reduced spleen ATCC 25922 strain concentrations (-2.53 and -2.14 log10 CFU/g [P < 0.05, respectively]) in the HI versus LI groups, while amoxicillin-clavulanate maintained its efficacy (-1.01 log10 CFU/g [no statistically significant difference]). Regarding the Ec1062 strain, the antimicrobials showed lower efficacy in the HI than in the LI groups: -0.73, -1.89, and -1.62 log10 CFU/g (P < 0.05, for piperacillin-tazobactam, imipenem, and amoxicillin-clavulanate, respectively, although imipenem and amoxicillin-clavulanate were more efficacious than piperacillin-tazobactam). An adapted imipenem treatment (based on the time for which the serum drug concentration remained above the MIC obtained with a HI of the ATCC 25922 strain) improved its efficacy to -1.67 log10 CFU/g (P < 0.05). These results suggest that amoxicillin-clavulanate could be an alternative to imipenem treatment of infections caused by ESBL- and non-ESBL-producing E. coli strains in patients with therapeutic failure with piperacillin-tazobactam.

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

Development of model observers applied to 3D breast tomosynthesis microcalcifications and masses

The development of model observers for mimicking human detection strategies has followed from symmetric signals in simple noise to increasingly complex backgrounds. In this study we implement different model observers for the complex task of detecting a signal in a 3D image stack. The backgrounds come from real breast tomosynthesis acquisitions and the signals were simulated and reconstructed within the volume. Two different tasks relevant to the early detection of breast cancer were considered: detecting an 8 mm mass and detecting a cluster of microcalcifications. The model observers were calculated using a channelized Hotelling observer (CHO) with dense difference-of-Gaussian channels, and a modified (Partial prewhitening [PPW]) observer which was adapted to realistic signals which are not circularly symmetric. The sustained temporal sensitivity function was used to filter the images before applying the spatial templates. For a frame rate of five frames per second, the only CHO that we calculated performed worse than the humans in a 4-AFC experiment. The other observers were variations of PPW and outperformed human observers in every single case. This initial frame rate was a rather low speed and the temporal filtering did not affect the results compared to a data set with no human temporal effects taken into account. We subsequently investigated two higher speeds at 5, 15 and 30 frames per second. We observed that for large masses, the two types of model observers investigated outperformed the human observers and would be suitable with the appropriate addition of internal noise. However, for microcalcifications both only the PPW observer consistently outperformed the humans. The study demonstrated the possibility of using a model observer which takes into account the temporal effects of scrolling through an image stack while being able to effectively detect a range of mass sizes and distributions.

Methodology for studying the numerical speed of sound in finite differences schemes

The space and time discretization inherent to all FDTD schemesintroduce non-physical dispersion errors, i.e. deviations ofthe speed of sound from the theoretical value predicted bythe governing Euler differential equations. A generalmethodologyfor computing this dispersion error via straightforwardnumerical simulations of the FDTD schemes is presented.The method is shown to provide remarkable accuraciesof the order of 1/1000 in a wide variety of twodimensionalfinite difference schemes.

The use of the quality model of Parasuraman, Zeithaml and Berry in health services

This is an article about the theoretical model for assessing quality in health services proposed by Parasuraman, Zheitaml and Berry, in order to measure the degree of satisfaction of users. This model is based on the analysis of expectations and perceptions of users of health services, by means of five dimensions: tangibility, reliability, responsiveness, assurance and empathy. From the difference between what is expected by the user and the service offered, gaps or shortcomings are derived that may be the main obstacle for users to perceive the provision of such services with quality. It was observed that the use of the psychometric scale called Service Quality (SERVQUAL) in some studies about satisfaction, obtained very favorable results in the institutions in which it was employed. The analysis revealed the need to improve the existing models of evaluation, as well as the importance of measuring user satisfaction in health institutions.