A semi-alternating direction method for a 2-D fractional FitzHugh–Nagumo monodomain model on an approximate irregular domain

A FitzHugh-Nagumo monodomain model has been used to describe the propagation of the electrical potential in heterogeneous cardiac tissue. In this paper, we consider a two-dimensional fractional FitzHugh-Nagumo monodomain model on an irregular domain. The model consists of a coupled Riesz space fractional nonlinear reaction-diffusion model and an ordinary differential equation, describing the ionic fluxes as a function of the membrane potential. Secondly, we use a decoupling technique and focus on solving the Riesz space fractional nonlinear reaction-diffusion model. A novel spatially second-order accurate semi-implicit alternating direction method (SIADM) for this model on an approximate irregular domain is proposed. Thirdly, stability and convergence of the SIADM are proved. Finally, some numerical examples are given to support our theoretical analysis and these numerical techniques are employed to simulate a two-dimensional fractional Fitzhugh-Nagumo model on both an approximate circular and an approximate irregular domain.

A Crank--Nicolson ADI spectral method for a two-dimensional riesz space fractional nonlinear reaction-diffusion equation

In this paper, a new alternating direction implicit Galerkin--Legendre spectral method for the two-dimensional Riesz space fractional nonlinear reaction-diffusion equation is developed. The temporal component is discretized by the Crank--Nicolson method. The detailed implementation of the method is presented. The stability and convergence analysis is strictly proven, which shows that the derived method is stable and convergent of order $2$ in time. An optimal error estimate in space is also obtained by introducing a new orthogonal projector. The present method is extended to solve the fractional FitzHugh--Nagumo model. Numerical results are provided to verify the theoretical analysis.

A Novel High Order Space-Time Spectral Method for the Time Fractional Fokker--Planck Equation

The fractional Fokker-Planck equation is an important physical model for simulating anomalous diffusions with external forces. Because of the non-local property of the fractional derivative an interesting problem is to explore high accuracy numerical methods for fractional differential equations. In this paper, a space-time spectral method is presented for the numerical solution of the time fractional Fokker-Planck initial-boundary value problem. The proposed method employs the Jacobi polynomials for the temporal discretization and Fourier-like basis functions for the spatial discretization. Due to the diagonalizable trait of the Fourier-like basis functions, this leads to a reduced representation of the inner product in the Galerkin analysis. We prove that the time fractional Fokker-Planck equation attains the same approximation order as the time fractional diffusion equation developed in [23] by using the present method. That indicates an exponential decay may be achieved if the exact solution is sufficiently smooth. Finally, some numerical results are given to demonstrate the high order accuracy and efficiency of the new numerical scheme. The results show that the errors of the numerical solutions obtained by the space-time spectral method decay exponentially.

Experimental and numerical investigation of strain-rate dependent mechanical properties of single living cells

The objective of this project is to investigate the strain-rate dependent mechanical behaviour of single living cells using both experimental and numerical techniques. The results revealed that living cells behave as porohyperlastic materials and that both solid and fluid phases within the cells play important roles in their mechanical responses. The research reported in this thesis provides a better understanding of the mechanisms underlying the cellular responses to external mechanical loadings and of the process of mechanical signal transduction in living cells. It would help us to enhance knowledge of and insight into the role of mechanical forces in supporting tissue regeneration or degeneration.

A novel population balance model for the dilute acid hydrolysis of hemicellulose

Acid hydrolysis is a popular pretreatment for removing hemicellulose from lignocelluloses in order to produce a digestible substrate for enzymatic saccharification. In this work, a novel model for the dilute acid hydrolysis of hemicellulose within sugarcane bagasse is presented and calibrated against experimental oligomer profiles. The efficacy of mathematical models as hydrolysis yield predictors and as vehicles for investigating the mechanisms of acid hydrolysis is also examined. Experimental xylose, oligomer (degree of polymerisation 2 to 6) and furfural yield profiles were obtained for bagasse under dilute acid hydrolysis conditions at temperatures ranging from 110C to 170C. Population balance kinetics, diffusion and porosity evolution were incorporated into a mathematical model of the acid hydrolysis of sugarcane bagasse. This model was able to produce a good fit to experimental xylose yield data with only three unknown kinetic parameters ka, kb and kd. However, fitting this same model to an expanded data set of oligomeric and furfural yield profiles did not successfully reproduce the experimental results. It was found that a ``hard-to-hydrolyse'' parameter, $\alpha$, was required in the model to ensure reproducibility of the experimental oligomer profiles at 110C, 125C and 140C. The parameters obtained through the fitting exercises at lower temperatures were able to be used to predict the oligomer profiles at 155C and 170C with promising results. The interpretation of kinetic parameters obtained by fitting a model to only a single set of data may be ambiguous. Although these parameters may correctly reproduce the data, they may not be indicative of the actual rate parameters, unless some care has been taken to ensure that the model describes the true mechanisms of acid hydrolysis. It is possible to challenge the robustness of the model by expanding the experimental data set and hence limiting the parameter space for the fitting parameters. The novel combination of ``hard-to-hydrolyse'' and population balance dynamics in the model presented here appears to stand up to such rigorous fitting constraints.

A cellular automata model to investigate immune cell–tumor cell interactions in growing tumors in two spatial dimensions

We develop a hybrid cellular automata model to describe the effect of the immune system and chemokines on a growing tumor. The hybrid cellular automata model consists of partial differential equations to model chemokine concentrations, and discrete cellular automata to model cell–cell interactions and changes. The computational implementation overlays these two components on the same spatial region. We present representative simulations of the model and show that increasing the number of immature dendritic cells (DCs) in the domain causes a decrease in the number of tumor cells. This result strongly supports the hypothesis that DCs can be used as a cancer treatment. Furthermore, we also use the hybrid cellular automata model to investigate the growth of a tumor in a number of computational “cancer patients.” Using these virtual patients, the model can explain that increasing the number of DCs in the domain causes longer “survival.” Not surprisingly, the model also reflects the fact that the parameter related to tumor division rate plays an important role in tumor metastasis.

A developmental cascade model of behavioral sleep problems and emotional and attentional self-regulation across early childhood

This paper documents the longitudinal and reciprocal relations among behavioral sleep problems, emotional and attentional self-regulation in a population sample of 4109 children participating in the Growing Up in Australia: The Longitudinal Study of Australian Children (LSAC) – Infant Cohort. Maternal reports of children’s sleep problems and self-regulation were collected at five time points from infancy to 8-9 years of age. Longitudinal structural equation modeling supported a developmental cascade model in which sleep problems have a persistent negative effect on emotional regulation, which in turn contributes to ongoing sleep problems and poorer attentional regulation in children over time. Findings suggest that sleep behaviors are a key target for interventions that aim to improve children’s self-regulatory capacities.

A comprehensive mechanistic kinetic model for dilute acid hydrolysis of switchgrass cellulose to glucose, 5-HMF and levulinic acid

Switchgrass was treated by 1% (w/w) H₂SO₄in batch tube reactors at temperatures ranging from 140–220°C for up to 60 minutes. In this study, release patterns of glucose, 5-hydroxymethylfurfural (5-HMF), and levulinic acid from switchgrass cellulose were investigated through a mechanistic kinetic model. The predictions were consistent with the measured products of interest when new parameters reflecting the effects of reaction limitations, such as cellulose crystallinity, acid soluble lignin–glucose complex (ASL–glucose) and humins that cannot be quantitatively analyzed, were included. The new mechanistic kinetic model incorporating these parameters simulated the experimental data with R² above 0.97. Results showed that glucose yield was most sensitive to variations in the parameter regarding the cellulose crystallinity at low temperatures (140–180°C), while the impact of crystallinity on the glucose yield became imperceptible at elevated temperatures (200–220 °C). Parameters related to the undesired products (e.g. ASL–glucose and humins) were the most sensitive factors compared with rate constants and other additional parameters in impacting the levulinic acid yield at elevated temperatures (200–220°C), while their impacts were negligible at 140–180°C. These new findings provide a more rational explanation for the kinetic changes in dilute acid pretreatment performance and suggest that the influences of cellulose crystallinity and undesired products including ASL–glucose and humins play key roles in determining the generation of glucose, 5-HMF and levulinic acid from biomass-derived cellulose.

Bivariate genome-wide association study of genetically correlated neuroimaging phenotypes from DTI and MRI through a seemingly unrelated regression model

Large multisite efforts (e.g., the ENIGMA Consortium), have shown that neuroimaging traits including tract integrity (from DTI fractional anisotropy, FA) and subcortical volumes (from T1-weighted scans) are highly heritable and promising phenotypes for discovering genetic variants associated with brain structure. However, genetic correlations (rg) among measures from these different modalities for mapping the human genome to the brain remain unknown. Discovering these correlations can help map genetic and neuroanatomical pathways implicated in development and inherited risk for disease. We use structural equation models and a twin design to find rg between pairs of phenotypes extracted from DTI and MRI scans. When controlling for intracranial volume, the caudate as well as related measures from the limbic system - hippocampal volume - showed high rg with the cingulum FA. Using an unrelated sample and a Seemingly Unrelated Regression model for bivariate analysis of this connection, we show that a multivariate GWAS approach may be more promising for genetic discovery than a univariate approach applied to each trait separately.

Surgical fusion of early onset severe scoliosis increases survival in the child with rett syndrome: A population based study

Introduction. Rett Syndrome is a rare genetic neurodevelopmental disorder usually affecting females. Scoliosis is a common comorbidity and spinal fusion may be recommended if severe. Little is known about long term outcomes. We examined the impact of spinal fusion on survival and risk of severe lower respiratory tract infection (LRTI) in Rett Syndrome. Methods Data were ascertained from hospital medical records, the Australian Rett Syndrome Database, a longitudinal and population-based registry of Rett Syndrome cases established in 1993, and the Australian Institute of Health and Welfare National Death Index database. An extended Cox regression model was used to estimate the effect of spinal surgery on survival in females who developed severe scoliosis (Cobb angle > 45 degrees). Generalized estimating equation modelling was used to estimate the effect of spinal surgery on the odds of developing severe LRTI. Results Severe scoliosis was identified in 140 cases (60.3%) of whom slightly fewer than half (48.6%) developed scoliosis prior to eight years of age. Scoliosis surgery was performed in 98 (69.0%) of those at a median age of 13 years 3 months (IQR 11 years 5 months – 14 years 10 months). After adjusting for mutation type and age of scoliosis onset, the rate of death was lower in the surgery group (HR 0.30, 95% CI 0.12, 0.74, P = 0.009) compared to those without surgery. Rate of death was particularly reduced for those with early onset scoliosis (HR 0.17, 95% CI 0.06, 0.52, P = 0.002). Spinal fusion was not associated with reduction in the occurrence of a severe LRTI overall (OR 0.60, 95%CI 0.27, 1.33, P=0.206) but was associated with a large reduction in odds of severe LRTI among those with early onset scoliosis (OR 0.32, 95%CI 0.11, 0.93, P=0.036). Conclusion With appropriate cautions, spinal fusion confers an advantage to life expectancy in Rett syndrome.

Spatially explicit structural equation modeling

Structural equation modeling (SEM) is a powerful statistical approach for the testing of networks of direct and indirect theoretical causal relationships in complex data sets with intercorrelated dependent and independent variables. SEM is commonly applied in ecology, but the spatial information commonly found in ecological data remains difficult to model in a SEM framework. Here we propose a simple method for spatially explicit SEM (SE-SEM) based on the analysis of variance/covariance matrices calculated across a range of lag distances. This method provides readily interpretable plots of the change in path coefficients across scale and can be implemented using any standard SEM software package. We demonstrate the application of this method using three studies examining the relationships between environmental factors, plant community structure, nitrogen fixation, and plant competition. By design, these data sets had a spatial component, but were previously analyzed using standard SEM models. Using these data sets, we demonstrate the application of SE-SEM to regularly spaced, irregularly spaced, and ad hoc spatial sampling designs and discuss the increased inferential capability of this approach compared with standard SEM. We provide an R package, sesem, to easily implement spatial structural equation modeling.

Determinants of the citation rate of medical research publications from a developing country

Background The number of citations received by an article is considered as an objective marker judging the importance and the quality of the research work. The present study aims to study the determinants of citations for research articles published by Sri Lankan authors. Methods Papers were selectively retrieved from the SciVerse Scopus® (Elsevier Properties S.A, USA) database for 10 years from 1st January 1997 to 31st December 2006, of which 50% were selected for inclusion by simple random sampling. The primary outcome measure was citation rate (defined as the number of citations during the 2 subsequent years after publication). Citation data was collected using the SciVerse Scopus® Citation Analyzer and self citations were excluded. A linear regression analysis was performed with ‘number of citations’ as the continuous dependent variable and other independent variables. Result The number of publications has steadily increased during the period of study. Over three quarter of papers were published in international journals. More than half of publications were research studies (55.3%), and most of the research studies were descriptive cross-sectional studies (27.1%). The mean number of citations within 2 years of publication was 1.7 and 52.1% of papers were not cited within the first two years of publication. The mean number of citations for collaborative studies (2.74) was significantly higher than that of non-collaborative studies (0.66). The mean number of citations did not significantly change depending on whether the publication had a positive result (2.08) or not (2.92) and was also not influenced by the presence (2.30) or absence (1.99) of the main study conclusion in the title of the article. In the linear regression model, the journal rank, number of authors, conducting the study abroad, being a research study or systematic review/meta-analysis and having regional and/or international collaboration all significantly increased the number of citations. Conclusion The journal rank, number of authors, conducting the study abroad, being a research study or systematic review/meta-analysis and having regional and/or international collaboration all significantly increased the number of citations. However, the presence of a positive result in the study did not influence the citation rate.

Physical mechanisms underlying the strain-rate-dependent mechanical behavior of kangaroo shoulder cartilage

Due to anatomical and biomechanical similarities to human shoulder, kangaroo was chosen as a model to study shoulder cartilage. Comprehensive enzymatic degradation and indentation tests were applied on kangaroo shoulder cartilage to study mechanisms underlying its strain-rate-dependent mechanical behavior. We report that superficial collagen plays a more significant role than proteoglycans in facilitating strain-rate-dependent behavior of kangaroo shoulder cartilage. By comparing the mechanical properties of degraded and normal cartilages it was noted that proteoglycan and collagen degradation significantly compromised strain-rate-dependent mechanical behavior of the cartilage. Superficial collagen contributed equally to the tissue behavior at all strain-rates. This is different to studies reported on knee cartilage and confirms the importance of superficial collagen on shoulder cartilage mechanical behavior. A porohyperelastic numerical model also indicated that collagen disruption would lead to faster damage of the shoulder cartilage than when proteoglycans are depleted.

Prediction of the fate of oxytetracycline and oxolinic acid in a fish pond using simulation model: A preliminary study

The fate of two popular antibiotics, oxytetracycline and oxolinic acid, in a fish pond were simulated using a computational model. The VDC model, which is designed based on a model for predicting pesticide fate and transport in paddy fields, was modified to take into account the differences between the pond and the paddies as well as those between the fish and the rice plant behaviors. The pond conditions were set following the typical practice in South East Asia aquaculture. The two antibiotics were administered to the animal in the pond through medicated feed during a period of 5 days as in actual practice. Concentrations of oxytetracycline in pond water were higher than those of oxolinic acid at the beginning of the simulation. Dissipation rate of oxytetracycline is also higher as it is more readily available for degradation in the water. For the long term, oxolinic acid was present at higher concentration than oxytetracycline in pond water as well as pond sediment. The simulated results were expected to be conservative and can be useful for the lower tier assessment of exposure risk of veterinary medicine in aquaculture industry but more data are needed for the complete validation of the model.

Predictors of condom use behaviour among male street labourers in urban Vietnam using a modified Information-Motivation-Behavioral Skills (IMB) model

HIV risk in vulnerable groups such as itinerant male street labourers is often examined via a focus on individual determinants. This study provides a test of a modified Information-Motivation-Behavioral Skills (IMB) model to predict condom use behaviour among male street workers in urban Vietnam. In a cross-sectional survey using a social mapping technique, 450 male street labourers from 13 districts of Hanoi, Vietnam were recruited and interviewed. Collected data were first examined for completeness; structural equation modelling was then employed to test the model fit. Condoms were used inconsistently by many of these men, and usage varied in relation to a number of factors. A modified IMB model had a better fit than the original IMB model in predicting condom use behaviour. This modified model accounted for 49% of the variance, versus 10% by the original version. In the modified model, the influence of psychosocial factors was moderately high, whilst the influence of HIV prevention information, motivation and perceived behavioural skills was moderately low, explaining in part the limited level of condom use behaviour. This study provides insights into social factors that should be taken into account in public health planning to promote safer sexual behaviour among Asian male street labourers.