A meshless method based on Point Interpolation Method (PIM) for the space fractional diffusion equation

This paper aims to develop a meshless approach based on the Point Interpolation Method (PIM) for numerical simulation of a space fractional diffusion equation. Two fully-discrete schemes for the one-dimensional space fractional diffusion equation are obtained by using the PIM and the strong-forms of the space diffusion equation. Numerical examples with different nodal distributions are studied to validate and investigate the accuracy and efficiency of the newly developed meshless approach.

Quantifying uncertainty in parameter estimates for stochastic models of collective cell spreading using approximate Bayesian computation

Wound healing and tumour growth involve collective cell spreading, which is driven by individual motility and proliferation events within a population of cells. Mathematical models are often used to interpret experimental data and to estimate the parameters so that predictions can be made. Existing methods for parameter estimation typically assume that these parameters are constants and often ignore any uncertainty in the estimated values. We use approximate Bayesian computation (ABC) to estimate the cell diffusivity, D, and the cell proliferation rate, λ, from a discrete model of collective cell spreading, and we quantify the uncertainty associated with these estimates using Bayesian inference. We use a detailed experimental data set describing the collective cell spreading of 3T3 fibroblast cells. The ABC analysis is conducted for different combinations of initial cell densities and experimental times in two separate scenarios: (i) where collective cell spreading is driven by cell motility alone, and (ii) where collective cell spreading is driven by combined cell motility and cell proliferation. We find that D can be estimated precisely, with a small coefficient of variation (CV) of 2–6%. Our results indicate that D appears to depend on the experimental time, which is a feature that has been previously overlooked. Assuming that the values of D are the same in both experimental scenarios, we use the information about D from the first experimental scenario to obtain reasonably precise estimates of λ, with a CV between 4 and 12%. Our estimates of D and λ are consistent with previously reported values; however, our method is based on a straightforward measurement of the position of the leading edge whereas previous approaches have involved expensive cell counting techniques. Additional insights gained using a fully Bayesian approach justify the computational cost, especially since it allows us to accommodate information from different experiments in a principled way.

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.

High order unconditionally stable difference schemes for the Riesz space-fractional telegraph equation

In this paper, a class of unconditionally stable difference schemes based on the Pad´e approximation is presented for the Riesz space-fractional telegraph equation. Firstly, we introduce a new variable to transform the original dfferential equation to an equivalent differential equation system. Then, we apply a second order fractional central difference scheme to discretise the Riesz space-fractional operator. Finally, we use (1, 1), (2, 2) and (3, 3) Pad´e approximations to give a fully discrete difference scheme for the resulting linear system of ordinary differential equations. Matrix analysis is used to show the unconditional stability of the proposed algorithms. Two examples with known exact solutions are chosen to assess the proposed difference schemes. Numerical results demonstrate that these schemes provide accurate and efficient methods for solving a space-fractional hyperbolic equation.

A fast semi-implicit difference method for a nonlinear two-sided space-fractional diffusion equation with variable diffusivity coefficients

In this paper, we derive a new nonlinear two-sided space-fractional diffusion equation with variable coefficients from the fractional Fick’s law. A semi-implicit difference method (SIDM) for this equation is proposed. The stability and convergence of the SIDM are discussed. For the implementation, we develop a fast accurate iterative method for the SIDM by decomposing the dense coefficient matrix into a combination of Toeplitz-like matrices. This fast iterative method significantly reduces the storage requirement of O(n2)O(n2) and computational cost of O(n3)O(n3) down to n and O(nlogn)O(nlogn), where n is the number of grid points. The method retains the same accuracy as the underlying SIDM solved with Gaussian elimination. Finally, some numerical results are shown to verify the accuracy and efficiency of the new method.

Machine learning approaches to analysing textual injury surveillance data: A systematic review

Objective To synthesise recent research on the use of machine learning approaches to mining textual injury surveillance data. Design Systematic review. Data sources The electronic databases which were searched included PubMed, Cinahl, Medline, Google Scholar, and Proquest. The bibliography of all relevant articles was examined and associated articles were identified using a snowballing technique. Selection criteria For inclusion, articles were required to meet the following criteria: (a) used a health-related database, (b) focused on injury-related cases, AND used machine learning approaches to analyse textual data. Methods The papers identified through the search were screened resulting in 16 papers selected for review. Articles were reviewed to describe the databases and methodology used, the strength and limitations of different techniques, and quality assurance approaches used. Due to heterogeneity between studies meta-analysis was not performed. Results Occupational injuries were the focus of half of the machine learning studies and the most common methods described were Bayesian probability or Bayesian network based methods to either predict injury categories or extract common injury scenarios. Models were evaluated through either comparison with gold standard data or content expert evaluation or statistical measures of quality. Machine learning was found to provide high precision and accuracy when predicting a small number of categories, was valuable for visualisation of injury patterns and prediction of future outcomes. However, difficulties related to generalizability, source data quality, complexity of models and integration of content and technical knowledge were discussed. Conclusions The use of narrative text for injury surveillance has grown in popularity, complexity and quality over recent years. With advances in data mining techniques, increased capacity for analysis of large databases, and involvement of computer scientists in the injury prevention field, along with more comprehensive use and description of quality assurance methods in text mining approaches, it is likely that we will see a continued growth and advancement in knowledge of text mining in the injury field.

Towards strongly consistent online HMM parameter estimation using one-step Kerridge inaccuracy

In this paper, we propose a novel online hidden Markov model (HMM) parameter estimator based on the new information-theoretic concept of one-step Kerridge inaccuracy (OKI). Under several regulatory conditions, we establish a convergence result (and some limited strong consistency results) for our proposed online OKI-based parameter estimator. In simulation studies, we illustrate the global convergence behaviour of our proposed estimator and provide a counter-example illustrating the local convergence of other popular HMM parameter estimators.

Turbulence and Mixing in the Environment: Multi-Device Study in a Sub-tropical Estuary

In an estuary, mixing and dispersion result from a combination of large-scale advection and smallscale turbulence, which are complex to estimate. The predictions of scalar transport and mixing are often inferred and rarely accurate, due to inadequate understanding of the contributions of these difference scales to estuarine recirculation. A multi-device field study was conducted in a small sub-tropical estuary under neap tide conditions with near-zero fresh water discharge for about 48 hours. During the study, acoustic Doppler velocimeters (ADV) were sampled at high frequency (50 Hz), while an acoustic Doppler current profiler (ADCP) and global positioning system (GPS) tracked drifters were used to obtain some lower frequency spatial distribution of the flow parameters within the estuary. The velocity measurements were complemented with some continuous measurement of water depth, conductivity, temperature and some other physiochemical parameters. Thorough quality control was carried out by implementation of relevant error removal filters on the individual data set to intercept spurious data. A triple decomposition (TD) technique was introduced to access the contributions of tides, resonance and ‘true’ turbulence in the flow field. The time series of mean flow measurements for both the ADCP and drifter were consistent with those of the mean ADV data when sampled within a similar spatial domain. The tidal scale fluctuation of velocity and water level were used to examine the response of the estuary to tidal inertial current. The channel exhibited a mixed type wave with a typical phase-lag between 0.035π– 0.116π. A striking feature of the ADV velocity data was the slow fluctuations, which exhibited large amplitudes of up to 50% of the tidal amplitude, particularly in slack waters. Such slow fluctuations were simultaneously observed in a number of physiochemical properties of the channel. The ensuing turbulence field showed some degree of anisotropy. For all ADV units, the horizontal turbulence ratio ranged between 0.4 and 0.9, and decreased towards the bed, while the vertical turbulence ratio was on average unity at z = 0.32 m and approximately 0.5 for the upper ADV (z = 0.55 m). The result of the statistical analysis suggested that the ebb phase turbulence field was dominated by eddies that evolved from ejection type process, while that of the flood phase contained mixed eddies with significant amount related to sweep type process. Over 65% of the skewness values fell within the range expected of a finite Gaussian distribution and the bulk of the excess kurtosis values (over 70%) fell within the range of -0.5 and +2. The TD technique described herein allowed the characterisation of a broader temporal scale of fluctuations of the high frequency data sampled within the durations of a few tidal cycles. The study provides characterisation of the ranges of fluctuation required for an accurate modelling of shallow water dispersion and mixing in a sub-tropical estuary.

Editorial: Research reporting: Clarity and accuracy improve quality

Research publication is one of the final steps in the research process, which begins with development of a research idea. Moving through the process of bringing together collaborators, design of the study protocol, securing of grant or study funding, and obtaining ethic(s) approval to conduct the research, and implementation of the research, analysis and drawing of conclusions based on the data leads to publication of the study results. Although a final step in the research process entails dissemination of the results, many studies go unreported or are improperly reported. Indeed, reviewers have suggested that many randomized controlled trials, observational studies, and qualitative studies lack crucial methodological features or details that lend credibility to study results (Simera et al., 2010).

Reliability control of GNSS carrier-phase integer ambiguity resolution

This research investigates how to obtain accurate and reliable positioning results with global navigation satellite systems (GNSS). The work provides a theoretical framework for reliability control in GNSS carrier phase ambiguity resolution, which is the key technique for precise GNSS positioning in centimetre levels. The proposed approach includes identification and exclusion procedures of unreliable solutions and hypothesis tests, allowing the reliability of solutions to be controlled in the aspects of mathematical models, integer estimation and ambiguity acceptance tests. Extensive experimental results with both simulation and observed data sets effectively demonstrate the reliability performance characteristics based on the proposed theoretical framework and procedures.

Association mapping

Association mapping seeks to identify marker alleles present at significantly different frequencies in cases carrying a particular disease or trait compared with controls. Genome-wide association studies are increasingly replacing candidate gene-based association studies for complex diseases, where a number of loci are likely to contribute to disease risk and the effect size of each particular risk allele is typically modest or low. Good study design is essential to the success of an association study, and factors such as the heritability of the disease under investigation, the choice of controls, statistical power, multiple testing and whether the association can be replicated need to be considered before beginning. Likewise, thorough quality control of the genotype data needs to be undertaken prior to running any association analyses. Finally, it should be kept in mind that a significant genetic association is not proof positive that a particular genetic locus causes a disease, but rather an important first step in discovering the genetic variants underlying a complex disease.

Genome-wide association studies of quantitative traits with related individuals: Little (power) lost but much to be gained

For complex disease genetics research in human populations, remarkable progress has been made in recent times with the publication of a number of genome-wide association scans (GWAS) and subsequent statistical replications. These studies have identified new genes and pathways implicated in disease, many of which were not known before. Given these early successes, more GWAS are being conducted and planned, both for disease and quantitative phenotypes. Many researchers and clinicians have DNA samples available on collections of families, including both cases and controls. Twin registries around the world have facilitated the collection of large numbers of families, with DNA and multiple quantitative phenotypes collected on twin pairs and their relatives. In the design of a new GWAS with a fixed budget for the number of chips, the question arises whether to include or exclude related individuals. It is commonly believed to be preferable to use unrelated individuals in the first stage of a GWAS because relatives are 'over-matched' for genotypes. In this study, we quantify that for GWAS of a quantitative phenotype, relative to a sample of unrelated individuals surprisingly little power is lost when using relatives. The advantages of using relatives are manifold, including the ability to perform more quality control, the choice to perform within-family tests of association that are robust to population stratification, and the ability to perform joint linkage and association analysis. Therefore, the advantages of using relatives in GWAS for quantitative traits may well outweigh the small disadvantage in terms of statistical power.

Electronic states of trans-polyacetylene, poly(p-phenylene vinylene) and sp-hybridised carbon species in amorphous hydrogenated carbon probed by resonant Raman scattering

Inclusions of sp-hybridised, trans-polyacetylene [trans-(CH)x] and poly(p-phenylene vinylene) (PPV) chains are revealed using resonant Raman scattering (RRS) investigation of amorphous hydrogenated carbon (a-C:H) films in the near IR – UV range. The RRS spectra of trans-(CH)x core Ag modes and the PPV CC-H phenylene mode are found to transform and disperse as the laser excitation energy ћωL is increased from near IR through visible to UV, whereas sp-bonded inclusions only become evident in UV. This is attributed to ћωL probing of trans-(CH)x chain inhomogeneity and the distribution of chains with varying conjugation length; for PPV to the resonant probing of phelynene ring disorder; and for sp segments, to ћωL probing of a local band gap of end-terminated polyynes. The IR spectra analysis confirmed the presence of sp, trans-(CH)x and PPV inclusions. The obtained RRS results for a-C:H denote differentiation between the core Ag trans-(CH)x modes and the PPV phenylene mode. Furthermore, it was found that at various laser excitation energies the changes in Raman spectra features for trans-(CH)x segments included in an amorphous carbon matrix are the same as in bulk trans-polyacetylene. The latter finding can be used to facilitate identification of trans-(CH)x in the spectra of complex carbonaceous materials.

Creating a therapeutic environment : a non-randomised controlled trial of a quiet time intervention for patients in acute care

Background: Noise is a significant barrier to sleep for acute care hospital patients, and sleep has been shown to be therapeutic for health, healing and recovery. Scheduled quiet time interventions to promote inpatient rest and sleep have been successfully trialled in critical care but not in acute care settings. Objectives: The study aim was to evaluate as cheduled quiet time intervention in an acute care setting. The study measured the effect of a scheduled quiet time on noise levels, inpatients’ rest and sleep behaviour, and wellbeing. The study also examined the impact of the intervention on patients’, visitors’ and health professionals’ satisfaction, and organisational functioning. Design: The study was a multi-centred non-randomised parallel group trial. Settings: The research was conducted in the acute orthopaedic wards of two major urban public hospitals in Brisbane, Australia. Participants: All patientsadmitted to the two wards in the5-month period of the study were invited to participate, withafinalsample of 299 participants recruited. This sample produced an effect size of 0.89 for an increase in the number of patients asleep during the quiet time. Methods: Demographic data were collected to enable comparison between groups. Data for noise level, sleep status, sleepiness and well being were collected using previously validated instruments: a Castle Model 824 digital sound level indicator; a three point sleep status scale; the Epworth Sleepiness Scale; and the SF12 V2 questionnaire. The staff, patient and visitor surveys on the experimental ward were adapted from published instruments. Results: Significant differences were found between the two groups in mean decibel level and numbers of patients awake and asleep. The difference in mean measured noise levels between the two environments corresponded to a ‘perceived’ difference of 2 to 1. There were significant correlations between average decibel level and number of patients awake and asleep in the experimental group, and between average decibel level and number of patients awake in the control group. Overall, patients, visitors and health professionals were satisfied with the quiet time intervention. Conclusions: The findings show that a quiet time intervention on an acute care hospital ward can affect noise level and patient sleep/wake patterns during the intervention period. The overall strongly positive response from surveys suggests that scheduled quiet time would be a positively perceived intervention with therapeutic benefit.

Sound radiation characteristics of a box-type structure

The finite element and boundary element methods are employed in this study to investigate the sound radiation characteristics of a box-type structure. It has been shown [T.R. Lin, J. Pan, Vibration characteristics of a box-type structure, Journal of Vibration and Acoustics, Transactions of ASME 131 (2009) 031004-1–031004-9] that modes of natural vibration of a box-type structure can be classified into six groups according to the symmetry properties of the three panel pairs forming the box. In this paper, we demonstrate that such properties also reveal information about sound radiation effectiveness of each group of modes. The changes of radiation efficiencies and directivity patterns with the wavenumber ratio (the ratio between the acoustic and the plate bending wavenumbers) are examined for typical modes from each group. Similar characteristics of modal radiation efficiencies between a box structure and a corresponding simply supported panel are observed. The change of sound radiation patterns as a function of the wavenumber ratio is also illustrated. It is found that the sound radiation directivity of each box mode can be correlated to that of elementary sound sources (monopole, dipole, etc.) at frequencies well below the critical frequency of the plates of the box. The sound radiation pattern on the box surface also closely related to the vibration amplitude distribution of the box structure at frequencies above the critical frequency. In the medium frequency range, the radiated sound field is dominated by the edge vibration pattern of the box. The radiation efficiency of all box modes reaches a peak at frequencies above the critical frequency, and gradually approaches unity at higher frequencies.