S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.

Implementing statistical learning methods through Bayesian networks. Part 1: a guide to Bayesian parameter estimation using forensic science data

As a thorough aggregation of probability and graph theory, Bayesian networks currently enjoy widespread interest as a means for studying factors that affect the coherent evaluation of scientific evidence in forensic science. Paper I of this series of papers intends to contribute to the discussion of Bayesian networks as a framework that is helpful for both illustrating and implementing statistical procedures that are commonly employed for the study of uncertainties (e.g. the estimation of unknown quantities). While the respective statistical procedures are widely described in literature, the primary aim of this paper is to offer an essentially non-technical introduction on how interested readers may use these analytical approaches - with the help of Bayesian networks - for processing their own forensic science data. Attention is mainly drawn to the structure and underlying rationale of a series of basic and context-independent network fragments that users may incorporate as building blocs while constructing larger inference models. As an example of how this may be done, the proposed concepts will be used in a second paper (Part II) for specifying graphical probability networks whose purpose is to assist forensic scientists in the evaluation of scientific evidence encountered in the context of forensic document examination (i.e. results of the analysis of black toners present on printed or copied documents).

Probabilistic evidential assessment of gunshot residue particle evidence (Part II) : Bayesian parameter estimation for experimental count data.

Part I of this series of articles focused on the construction of graphical probabilistic inference procedures, at various levels of detail, for assessing the evidential value of gunshot residue (GSR) particle evidence. The proposed models - in the form of Bayesian networks - address the issues of background presence of GSR particles, analytical performance (i.e., the efficiency of evidence searching and analysis procedures) and contamination. The use and practical implementation of Bayesian networks for case pre-assessment is also discussed. This paper, Part II, concentrates on Bayesian parameter estimation. This topic complements Part I in that it offers means for producing estimates useable for the numerical specification of the proposed probabilistic graphical models. Bayesian estimation procedures are given a primary focus of attention because they allow the scientist to combine (his/her) prior knowledge about the problem of interest with newly acquired experimental data. The present paper also considers further topics such as the sensitivity of the likelihood ratio due to uncertainty in parameters and the study of likelihood ratio values obtained for members of particular populations (e.g., individuals with or without exposure to GSR).

Parameter estimation for the distribution of single cell lag times.

In Quantitative Microbial Risk Assessment, it is vital to understand how lag times of individual cells are distributed over a bacterial population. Such identified distributions can be used to predict the time by which, in a growth-supporting environment, a few pathogenic cells can multiply to a poisoning concentration level. We model the lag time of a single cell, inoculated into a new environment, by the delay of the growth function characterizing the generated subpopulation. We introduce an easy-to-implement procedure, based on the method of moments, to estimate the parameters of the distribution of single cell lag times. The advantage of the method is especially apparent for cases where the initial number of cells is small and random, and the culture is detectable only in the exponential growth phase.

Frequency-dependent sexual selection with respect to offspring fitness returns is consistent with predictions from rock-paper-scissors dynamics in the European common lizard

Genetic polymorphism can be maintained over time by negative frequency-dependent (FD) selection induced by Rock-paper-scissors (RPS) social systems. RPS games produce cyclic dynamics, and have been suggested to exist in lizards, insects, isopods, plants, and bacteria. Sexual selection is predicted to accentuate the survival of the future progeny during negative FD survival selection. More specifically, females are predicted to select mates that produce progeny genotypes that exhibit highest survival during survival selection imposed by adult males. However, no empirical evidence demonstrates the existence of FD sexual selection with respect to fitness payoffs of genetic polymorphisms. Here we tested this prediction using the common lizard Zootoca vivipara, a species with three male color morphs (orange, white, yellow) that exhibit morph frequency cycles. In a first step we tested the congruence of the morph frequency change with the predicted change in three independent populations, differing in male color morph frequency and state of the FD morph cycle. Thereafter we ran standardized sexual selection assays in which we excluded alternative mechanisms that potentially induce negative FD selection, and we quantified inter-sexual behavior. The patterns of sexual selection and the observed behavior were in line with context-dependent female mate choice and male behavior played a minor role. Moreover, the strength of the sexual selection was within the magnitude of selection required to produce the observed 3-4-year and 6-8 year morph frequency cycles at low and high altitudes, respectively. In summary, the study provides the first experimental evidence that underpins the crucial assumption of the RPS games suggested to exist in lizards, insects, isopods, and plants; namely, that sexual selection produces negative-FD selection. This indicates that sexual selection, in our study exert by females, might be a crucial driver of the maintenance of genetic polymorphisms.

Tonotopic gradients in human primary auditory cortex: concurring evidence from high-resolution 7 t and 3 t FMRI.

The tonotopic representations within the primary auditory cortex (PAC) have been successfully mapped with ultra-high field fMRI. Here, we compared the reliability of this tonotopic mapping paradigm at 7 T with 1.5 mm spatial resolution with maps acquired at 3 T with the same stimulation paradigm, but with spatial resolutions of 1.8 and 2.4 mm. For all subjects, the mirror-symmetric gradients within PAC were highly similar at 7 T and 3 T and across renderings at different spatial resolutions; albeit with lower percent signal changes at 3 T. In contrast, the frequency maps outside PAC tended to suffer from a reduced BOLD contrast-to-noise ratio at 3 T for a 1.8 mm voxel size, while robust at 2.4 mm and at 1.5 mm at 7 T. Overall, our results showed the robustness of the phase-encoding paradigm used here to map tonotopic representations across scanners.

Fat signal suppression for coronary MRA at 3T using a water-selective adiabatic T2 -preparation technique.

PURPOSE: To improve fat saturation in coronary MRA at 3T by using a spectrally selective adiabatic T2 -Prep (WSA-T2 -Prep). METHODS: A conventional adiabatic T2 -Prep (CA-T2 -Prep) was modified, such that the excitation and restoration pulses were of differing bandwidths. On-resonance spins are T2 -Prepared, whereas off-resonance spins, such as fat, are spoiled. This approach was combined with a CHEmically Selective Saturation (CHESS) pulse to achieve even greater fat suppression. Numerical simulations were followed by phantom validation and in vivo coronary MRA. RESULTS: Numerical simulations demonstrated that augmenting a CHESS pulse with a WSA-T2 -Prep improved robustness to B1 inhomogeneities and that this combined fat suppression was effective over a broader spectral range than that of a CHESS pulse in a conventional T2 -Prepared sequence. Phantom studies also demonstrated that the WSA-T2 -Prep+CHESS combination produced greater fat suppression across a range of B1 values than did a CA-T2 -Prep+CHESS combination. Lastly, in vivo measurements demonstrated that the contrast-to-noise ratio between blood and myocardium was not adversely affected by using a WSA-T2 -Prep, despite the improved abdominal and epicardial fat suppression. Additionally, vessel sharpness improved. CONCLUSION: The proposed WSA-T2 -Prep method was shown to improve fat suppression and vessel sharpness as compared to a CA-T2 -Prep technique, and to also increase fat suppression when combined with a CHESS pulse.

Off-resonance angiography: a new method to depict vessels--phantom and rabbit studies.

PURPOSE: To evaluate the utility of inversion recovery with on-resonant water suppression (IRON) in combination with injection of the long-circulating monocrystalline iron oxide nanoparticle (MION)-47 for contrast material-enhanced magnetic resonance (MR) angiography. MATERIALS AND METhods: Experiments were approved by the institutional animal care committee. Eleven rabbits were imaged at baseline before injection of a contrast agent and then serially 5-30 minutes, 2 hours, 1 day, and 3 days after a single intravenous bolus injection of 80 micromol of MION-47 per kilogram of body weight (n = 6) or 250 micromol/kg MION-47 (n = 5). Conventional T1-weighted MR angiography and IRON MR angiography were performed on a clinical 3.0-T imager. Signal-to-noise and contrast-to-noise ratios were measured in the aorta of rabbits in vivo. Venous blood was obtained from the rabbits before and after MION-47 injection for use in phantom studies. RESULTS: In vitro blood that contained MION-47 appeared signal attenuated on T1-weighted angiograms, while characteristic signal-enhanced dipolar fields were observed on IRON angiograms. In vivo, the vessel lumen was signal attenuated on T1-weighted MR angiograms after MION-47 injection, while IRON supported high intravascular contrast by simultaneously providing positive signal within the vessels and suppressing background tissue (mean contrast-to-noise ratio, 61.9 +/- 12.4 [standard deviation] after injection vs 1.1 +/- 0.4 at baseline, P < .001). Contrast-to-noise ratio was higher on IRON MR angiograms than on conventional T1-weighted MR angiograms (9.0 +/- 2.5, P < .001 vs IRON MR angiography) and persisted up to 24 hours after MION-47 injection (76.2 +/- 15.9, P < .001 vs baseline). CONCLUSION: IRON MR angiography in conjunction with superparamagnetic nanoparticle administration provides high intravascular contrast over a long time and without the need for image subtraction.

Frequency and determinants of urgent requests to home care agencies for community-dwelling elderly.

We analyzed a one-year case series and performed a longitudinal (4 month) cohort analysis of urgent requests made to home care agencies by and for their > or = 65 years old clients in order to estimate the frequency of unscheduled services delivered by home care agencies and to identify risk factors. All 40 home care agencies located in a Swiss region were included in the study and we registered 3,816 urgent requests (75/1,000 > or = 65 years residents per year). Among home care users, the presence of a urinary catheter, incontinence and the need for assistance in bathing were predictors of unscheduled services. Resources should be planned in order to help home care teams to handle unexpected, disruptive clusters of urgent requests that may compromise their scheduled activities.

Reproducibility of free-breathing cardiovascular magnetic resonance coronary angiography.

OBJECTIVE: Contemporary free-breathing non contrast enhanced cardiovascular magnetic resonance angiography (CMRA) was qualitatively and quantitatively evaluated to ascertain the reproducibility of the method for coronary artery luminal dimension measurements. SUBJECTS AND METHODS: Twenty-two healthy volunteers (mean age 32 +/- 7 years, 12 males) without coronary artery disease were imaged at 2 centers (1 each in Europe and North America) using navigator-gated and corrected SSFP CMRA on a commercial whole body 1.5T System. Repeat images of right (RCA, n = 21), left anterior descending (LAD, n = 14) and left circumflex (LCX, n = 14) coronary arteries were obtained in separate sessions using identical scan protocol and imaging parameters. True visible vessel length, signal-to-noise (SNR), contrast-to-noise ratios (CNR) and the average luminal diameter over the first 4 cm of the vessel were measured. Intra-observer, inter-observer and inter-scan reproducibility of coronary artery luminal diameter were determined using Pearson's correlation, Bland-Altman analysis and intraclass correlation coefficients (ICC). RESULTS: CNR, SNR and the mean length of the RCA, LAD and LCX imaged for original and repeat scans were not significantly different (all p > 0.30). There was a high degree of intra-observer, inter-observer and inter-scan agreements for RCA, LAD and LCX luminal diameter respectively on Bland-Altman and ICC analysis (ICC's for RCA: 0.98. 0.98 and 0.86; LAD: 0.89, 0.89 and 0.63; LCX: 0.95, 0.94 and 0.79). CONCLUSION: In a 2-center study, we demonstrate that free-breathing 3D SSFP CMRA can visualize long continuous segments of coronary vessels with highly reproducible measurements of luminal diameter.

Does Semantic Context Benefit Speech Understanding through "Top-Down" Processes? Evidence from Time-resolved Sparse fMRI.

When speech is degraded, word report is higher for semantically coherent sentences (e.g., her new skirt was made of denim) than for anomalous sentences (e.g., her good slope was done in carrot). Such increased intelligibility is often described as resulting from "top-down" processes, reflecting an assumption that higher-level (semantic) neural processes support lower-level (perceptual) mechanisms. We used time-resolved sparse fMRI to test for top-down neural mechanisms, measuring activity while participants heard coherent and anomalous sentences presented in speech envelope/spectrum noise at varying signal-to-noise ratios (SNR). The timing of BOLD responses to more intelligible speech provides evidence of hierarchical organization, with earlier responses in peri-auditory regions of the posterior superior temporal gyrus than in more distant temporal and frontal regions. Despite Sentence content × SNR interactions in the superior temporal gyrus, prefrontal regions respond after auditory/perceptual regions. Although we cannot rule out top-down effects, this pattern is more compatible with a purely feedforward or bottom-up account, in which the results of lower-level perceptual processing are passed to inferior frontal regions. Behavioral and neural evidence that sentence content influences perception of degraded speech does not necessarily imply "top-down" neural processes.

Repositioning precision of coronary arteries measured on X-ray angiography and its implications for coronary MR angiography.

BACKGROUND: To test the hypothesis that intervals with superior beat-to-beat coronary artery repositioning precision exist in the cardiac cycle, to design a coronary MR angiography (MRA) methodology in response, and to ascertain its performance. METHODS: Coronary repositioning precision in consecutive heartbeats was measured on x-ray coronary angiograms of 17 patients and periods with the highest repositioning precision were identified. In response, the temporal order of coronary MRA pulse sequence elements required modification and the T2 -prep now follows (T2 -post) rather than precedes the imaging part of the sequence. The performance of T2 -post was quantitatively compared (signal-to-noise [SNR], contrast-to-noise [CNR], vessel sharpness) to that of T2 -prep in vivo. RESULTS: Coronary repositioning precision is <1 mm at peak systole and in mid diastole. When comparing systolic T2 -post to diastolic T2 -prep, CNR and vessel sharpness remained unchanged (both P = NS) but SNR for muscle and blood increased by 104% and 36% (both P < 0.05), respectively. CONCLUSION: Windows with improved coronary repositioning precision exist in the cardiac cycle: one in peak systole and one in mid diastole. Peak-systolic imaging necessitates a re-design of conventional coronary MRA pulse sequences and leads to image quality very similar to that of conventional mid-diastolic data acquisition but improved SNR. J. Magn. Reson. Imaging 2015;41:1251-1258. © 2014 Wiley Periodicals, Inc.

Using a NPWE model observer to assess suitable image quality for a digital mammography quality assurance programme.

A method of objectively determining imaging performance for a mammography quality assurance programme for digital systems was developed. The method is based on the assessment of the visibility of a spherical microcalcification of 0.2 mm using a quasi-ideal observer model. It requires the assessment of the spatial resolution (modulation transfer function) and the noise power spectra of the systems. The contrast is measured using a 0.2-mm thick Al sheet and Polymethylmethacrylate (PMMA) blocks. The minimal image quality was defined as that giving a target contrast-to-noise ratio (CNR) of 5.4. Several evaluations of this objective method for evaluating image quality in mammography quality assurance programmes have been considered on computed radiography (CR) and digital radiography (DR) mammography systems. The measurement gives a threshold CNR necessary to reach the minimum standard image quality required with regards to the visibility of a 0.2-mm microcalcification. This method may replace the CDMAM image evaluation and simplify the threshold contrast visibility test used in mammography quality.

Navigator-gated free-breathing three-dimensional balanced fast field echo (TrueFISP) coronary magnetic resonance angiography.

RATIONALE AND OBJECTIVES: Recent developments of MR imaging equipment enabled high-quality steady state-free-precession (Balanced FFE, True-FISP) MR-imaging with a substantial 'T2 like' contrast, resulting in a high signal intensity of the blood-pool without the application of exogenous contrast agents. It is hypothesized that Balanced-FFE may be valuable for contrast enhancement in 3D free-breathing coronary MRA. MATERIALS AND METHODS: Navigator-gated free-breathing cardiac triggered coronary MRA was performed in 10 healthy adult subjects and three patients with radiograph defined coronary artery disease using a segmented k-space 3D Balanced FFE imaging sequence. RESULTS: High contrast-to-noise ratio between the blood-pool and the myocardium (29 +/- 8) and long segment visualization of both coronary arteries could be obtained in about 5 minutes during free breathing using the present navigator-gated Balanced-FFE coronary MRA approach. First patient results demonstrated successful display of coronary artery stenoses. CONCLUSION: Balanced FFE offers a potential alternative for endogenous contrast enhancement in navigator-gated free-breathing 3D coronary MRA. The obtained results together with the relatively short scanning time warrant further studies in larger patient collectives.