A Krylov-based finite state projection algorithm for solving the chemical master equation arising in the discrete modelling of biological systems

Biochemical reactions underlying genetic regulation are often modelled as a continuous-time, discrete-state, Markov process, and the evolution of the associated probability density is described by the so-called chemical master equation (CME). However the CME is typically difficult to solve, since the state-space involved can be very large or even countably infinite. Recently a finite state projection method (FSP) that truncates the state-space was suggested and shown to be effective in an example of a model of the Pap-pili epigenetic switch. However in this example, both the model and the final time at which the solution was computed, were relatively small. Presented here is a Krylov FSP algorithm based on a combination of state-space truncation and inexact matrix-vector product routines. This allows larger-scale models to be studied and solutions for larger final times to be computed in a realistic execution time. Additionally the new method computes the solution at intermediate times at virtually no extra cost, since it is derived from Krylov-type methods for computing matrix exponentials. For the purpose of comparison the new algorithm is applied to the model of the Pap-pili epigenetic switch, where the original FSP was first demonstrated. Also the method is applied to a more sophisticated model of regulated transcription. Numerical results indicate that the new approach is significantly faster and extendable to larger biological models.

Accelerated leap methods for simulating discrete stochastic chemical kinetics

Biologists are increasingly conscious of the critical role that noise plays in cellular functions such as genetic regulation, often in connection with fluctuations in small numbers of key regulatory molecules. This has inspired the development of models that capture this fundamentally discrete and stochastic nature of cellular biology - most notably the Gillespie stochastic simulation algorithm (SSA). The SSA simulates a temporally homogeneous, discrete-state, continuous-time Markov process, and of course the corresponding probabilities and numbers of each molecular species must all remain positive. While accurately serving this purpose, the SSA can be computationally inefficient due to very small time stepping so faster approximations such as the Poisson and Binomial τ-leap methods have been suggested. This work places these leap methods in the context of numerical methods for the solution of stochastic differential equations (SDEs) driven by Poisson noise. This allows analogues of Euler-Maruyuma, Milstein and even higher order methods to be developed through the Itô-Taylor expansions as well as similar derivative-free Runge-Kutta approaches. Numerical results demonstrate that these novel methods compare favourably with existing techniques for simulating biochemical reactions by more accurately capturing crucial properties such as the mean and variance than existing methods.

Sequential Monte Carlo for Bayesian sequentially designed experiments for discrete data

In this paper we present a sequential Monte Carlo algorithm for Bayesian sequential experimental design applied to generalised non-linear models for discrete data. The approach is computationally convenient in that the information of newly observed data can be incorporated through a simple re-weighting step. We also consider a flexible parametric model for the stimulus-response relationship together with a newly developed hybrid design utility that can produce more robust estimates of the target stimulus in the presence of substantial model and parameter uncertainty. The algorithm is applied to hypothetical clinical trial or bioassay scenarios. In the discussion, potential generalisations of the algorithm are suggested to possibly extend its applicability to a wide variety of scenarios

Solder joint defects classification using the Log-Gabor Filter, the Discrete Wavelet Transform and the Discrete Cosine Transform

Inspection of solder joints has been a critical process in the electronic manufacturing industry to reduce manufacturing cost, improve yield, and ensure project quality and reliability. This paper proposes the use of the Log-Gabor filter bank, Discrete Wavelet Transform and Discrete Cosine Transform for feature extraction of solder joint images on Printed Circuit Boards (PCBs). A distance based on the Mahalanobis Cosine metric is also presented for classification of five different types of solder joints. From the experimental results, this methodology achieved high accuracy and a well generalised performance. This can be an effective method to reduce cost and improve quality in the production of PCBs in the manufacturing industry.

Profiling and characterisation of key weather variables and their implication on building design

Local climate is a critical element in the design of buildings. In this paper, ten years of historical weather data in Australia's all eight capital cities are analyzed to characterize the variation profiles of climatic variables. The method of descriptive statistics is employed. Either the pattern of cumulative distribution and/or the profile of percentage distribution are used to graphically illustrate the similarity and difference between different study locations. It is found that although the weather variables vary with different locations, except for the extreme parts, there is often a good, nearly linear relation between weather variable and its cumulative percentage for the majority of middle part. The implication of these extreme parts and the slopes of the middle parts on building design is also discussed.

On the significance of omitted variables in intersection crash modeling

Advances in safety research—trying to improve the collective understanding of motor vehicle crash causes and contributing factors—rest upon the pursuit of numerous lines of research inquiry. The research community has focused considerable attention on analytical methods development (negative binomial models, simultaneous equations, etc.), on better experimental designs (before-after studies, comparison sites, etc.), on improving exposure measures, and on model specification improvements (additive terms, non-linear relations, etc.). One might logically seek to know which lines of inquiry might provide the most significant improvements in understanding crash causation and/or prediction. It is the contention of this paper that the exclusion of important variables (causal or surrogate measures of causal variables) cause omitted variable bias in model estimation and is an important and neglected line of inquiry in safety research. In particular, spatially related variables are often difficult to collect and omitted from crash models—but offer significant opportunities to better understand contributing factors and/or causes of crashes. This study examines the role of important variables (other than Average Annual Daily Traffic (AADT)) that are generally omitted from intersection crash prediction models. In addition to the geometric and traffic regulatory information of intersection, the proposed model includes many spatial factors such as local influences of weather, sun glare, proximity to drinking establishments, and proximity to schools—representing a mix of potential environmental and human factors that are theoretically important, but rarely used. Results suggest that these variables in addition to AADT have significant explanatory power, and their exclusion leads to omitted variable bias. Provided is evidence that variable exclusion overstates the effect of minor road AADT by as much as 40% and major road AADT by 14%.

Analytical solutions for the two- and three-dimensional time-fractional telegraph equations by the method of separating variables

In this paper, a method of separating variables is effectively implemented for solving a time-fractional telegraph equation (TFTE) in two and three dimensions. We discuss and derive the analytical solution of the TFTE in two and three dimensions with nonhomogeneous Dirichlet boundary condition. This method can be extended to other kinds of the boundary conditions.

Practical stability of approximating discrete-time filters with respect to model mismatch

This paper establishes practical stability results for an important range of approximate discrete-time filtering problems involving mismatch between the true system and the approximating filter model. Practical stability is established in the sense of an asymptotic bound on the amount of bias introduced by the model approximation. Our analysis applies to a wide range of estimation problems and justifies the common practice of approximating intractable infinite dimensional nonlinear filters by simpler computationally tractable filters.

A history of heart interventions moderates the relationship between psychological variables and the presence of chest pain in older women with self-reported coronary heart disease

Objectives: This study examines the hypothesis that a past history of heart interventions will moderate the relationship between psychosocial factors (stressful life events, social support, perceived stress, having a current partner, having a past diagnosis of depression or anxiety over the past 3 years, time pressure, education level, and the mental health index) and the presence of chest pain in a sample of older women. Design: Longitudinal survey over a 3-year period. Methods: The sample was taken from a prospective cohort study of 10,432 women initially aged between 70 and 75 years, who were surveyed in 1996 and then again in 1999. Two groups of women were identified: those reporting to have heart disease but no past history of heart interventions (i.e., coronary artery bypass graft/angioplasty) and those reporting to have heart disease with a past history of heart interventions. Results: Binary logistic regression analysis was used to show that for the women with self-reported coronary heart disease but without a past history of heart intervention, feelings of time pressure as well as the number of stressful life events experienced in the 12 months prior to 1996 were independent risk factors for the presence of chest pain, even after accounting for a range of traditional risk factors. In comparison, for the women with self-reported coronary heart disease who did report a past history of heart interventions, a diagnosis of depression in the previous 3 years was the significant independent risk factor for chest pain even after accounting for traditional risk factors. Conclusion: The results indicate that it is important to consider a history of heart interventions as a moderator of the associations between psychosocial variables and the frequency of chest pain in older women. Statement of Contribution: What is already known on this subject? Psychological factors have been shown to be independent predictors of a range of health outcomes in individuals with coronary heart disease, including the presence of chest pain. Most research has been conducted with men or with small samples of women; however, the evidence does suggest that these relationships exist in women as well as in men. What does this study add? Most studies have looked at overall relationships between psychological variables and health outcomes. The few studies that have looked at moderators have mainly examined gender as a moderator. To our knowledge, this is the first published study to examine a history of heart interventions as a moderator of the relationship between psychological variables and the presence of chest pain.

A critical evaluation and comparison of two formative measures of system quality using criterion variables

Evaluating the validity of formative variables has presented ongoing challenges for researchers. In this paper we use global criterion measures to compare and critically evaluate two alternative formative measures of System Quality. One model is based on the ISO-9126 software quality standard, and the other is based on a leading information systems research model. We find that despite both models having a strong provenance, many of the items appear to be non-significant in our study. We examine the implications of this by evaluating the quality of the criterion variables we used, and the performance of PLS when evaluating formative models with a large number of items. We find that our respondents had difficulty distinguishing between global criterion variables measuring different aspects of overall System Quality. Also, because formative indicators “compete with one another” in PLS, it may be difficult to develop a set of measures which are all significant for a complex formative construct with a broad scope and a large number of items. Overall, we suggest that there is cautious evidence that both sets of measures are valid and largely equivalent, although questions still remain about the measures, the use of criterion variables, and the use of PLS for this type of model evaluation.

On the probability density function of the product of Rayleigh distributed random variables

In this paper we investigate the distribution of the product of Rayleigh distributed random variables. Considering the Mellin-Barnes inversion formula and using the saddle point approach we obtain an upper bound for the product distribution. The accuracy of this tail-approximation increases as the number of random variables in the product increase.

Synthesis of discrete allyl-functional polyethylene glycol for efficient hydrophilization of dihydr (polydimethylsiloxane)

Amphiphilic poly(ethylene glycol)-block-pol (dimethylsiloxane)-block-poly(ethylene glycol)(PEG-block-PDMS block-PEG) triblock copolymers have been successfully prepared via hydrosilylation using discrete and polydisperse PEG of various chain lengths. Facile synthesis of discrete PEG (dPEG) is achieved via systematic tosylation and etherification of lower glycols. Amphiphilicity of the dPEG block-PDMS-block-dPEG triblock copolymer is illustrated by dynamic light scattering (DLS) and measurement of the critical micelle concentration (CMC).

Effects of bovine colostrum supplementation on immune variables of highly-trained cyclists

The aim of this study was to investigate the influence of low-dose bovine colostrum protein concentrate (CPC) supplementation on selected immune variables in cyclists. Twenty-nine highly trained male road cyclists completed an initial 40-km time trial (TT(40)) and were then randomly assigned to either a supplement (n = 14, 10 g bovine CPC/day) or placebo group (n = 15, 10 g whey protein concentrate/day). After 5 wk of supplementation, the cyclists completed a second TT(40). They then completed 5 consecutive days of high-intensity training (HIT) that included a TT(40), followed by a final TT(40) in the following week. Venous blood and saliva samples were collected immediately before and after each TT(40), and upper respiratory illness symptoms were recorded over the experimental period. Compared with the placebo group, bovine CPC supplementation significantly increased preexercise serum soluble TNF receptor 1 during the HIT period (bovine CPC = 882 +/- 233 pg/ml, placebo = 468 +/- 139 pg/ml; P = 0.039). Supplementation also suppressed the postexercise decrease in cytotoxic/suppressor T cells during the HIT period (bovine CPC = -1.0 +/- 2.7%, placebo = -9.2 +/- 2.8%; P = 0.017) and during the following week (bovine CPC = 1.4 +/- 2.9%, placebo = -8.2 +/- 2.8%; P = 0.004). Bovine CPC supplementation prevented a postexercise decrease in serum IgG(2) concentration at the end of the HIT period (bovine CPC = 4.8 +/- 6.8%, P = 0.88; placebo = -9.7 +/- 6.9%, P = 0.013). There was a trend toward reduced incidence of upper respiratory illness symptoms in the bovine CPC group (P = 0.055). In summary, low-dose bovine CPC supplementation modulates immune parameters during normal training and after an acute period of intense exercise, which may have contributed to the trend toward reduced upper respiratory illness in the bovine CPC group.

Discrete loading of the foot during hurdle drills: A pilot investigation

Introduction: Lower limb function in hurdling is patently asymmetrical. The lead limb undertakes the preparatory and landing steps while the trail limb contends with the hurdle and recovery steps. Discrete loading profiles of these steps will reflect the asymmetrical function and may provide useful insight into injury mechanisms. A pilot study was undertaken to determine the loading profiles of the hurdle, landing and recovery steps of elite male hurdlers. Equivalent data for steps between hurdles, where the running action is more symmetrical, were used for the purpose of comparison, simultaneously minimising the confounding effect of speed. Methodology: In-shoe pressures were recorded (FScan, 200 Hz) for four elite male hurdlers while they completed a routine hurdle drill at a self-selected fast but sub-race speed. The drill comprised of three consecutive hurdles. Data for the hurdle, landing and recovery steps of the first and second hurdles, along with data for the running steps between hurdles 1 and 2, and 2 and 3, were used for the purpose of analysis. Peak pressures within 1cm2 masks were determined for the hallux, first, central and fifth metatarsals (T1, M1, M2–4 and M5 respectively). Peak pressure (kPa) and loading duration (ms) for the hurdle, landing and recovery steps are reported as a percentage of the respective limb-matched values for between-hurdle steps. Results/discussion: For between-hurdle steps, T1, M1 and M2–4 peak pressures were 312/357, 356/306 and 362/368 kPa, lead/trail limbs respectively. For the hurdle, landing and recovery steps, pressures at T1 and M1 increased. For T1 the increases were in the order of 17%, 36% and 8% (hurdle, landing and recovery steps, respectively) while the corresponding increases at M1 were 7%, 54% and 20%. Pressures at M2–4 were similar for all steps, while M5 loaded erratically. For the between-hurdle steps, the loading durations at T1, M1 and M2–4, were 160/162, 170/142 and 190/191 ms, respectively. For the landing step, loading duration decreased for T1, M1and M2–4 (−8%, −19% and −18%, respectively). In the hurdle step, loading duration decreased for the metatarsals but not for T1. Conclusions: The hurdling action leads to regional pressure increases that act for shorter durations in comparison to the between-hurdle running steps. These changes are most notable at the first metatarsal, a common site of foot injury.

Lattice-free models of cell invasion : discrete simulations and travelling waves

Invasion waves of cells play an important role in development, disease and repair. Standard discrete models of such processes typically involve simulating cell motility, cell proliferation and cell-to-cell crowding effects in a lattice-based framework. The continuum-limit description is often given by a reaction–diffusion equation that is related to the Fisher–Kolmogorov equation. One of the limitations of a standard lattice-based approach is that real cells move and proliferate in continuous space and are not restricted to a predefined lattice structure. We present a lattice-free model of cell motility and proliferation, with cell-to-cell crowding effects, and we use the model to replicate invasion wave-type behaviour. The continuum-limit description of the discrete model is a reaction–diffusion equation with a proliferation term that is different from lattice-based models. Comparing lattice based and lattice-free simulations indicates that both models lead to invasion fronts that are similar at the leading edge, where the cell density is low. Conversely, the two models make different predictions in the high density region of the domain, well behind the leading edge. We analyse the continuum-limit description of the lattice based and lattice-free models to show that both give rise to invasion wave type solutions that move with the same speed but have very different shapes. We explore the significance of these differences by calibrating the parameters in the standard Fisher–Kolmogorov equation using data from the lattice-free model. We conclude that estimating parameters using this kind of standard procedure can produce misleading results.