Influence of network topology and data collection on network inference.

We recently developed an approach for testing the accuracy of network inference algorithms by applying them to biologically realistic simulations with known network topology. Here, we seek to determine the degree to which the network topology and data sampling regime influence the ability of our Bayesian network inference algorithm, NETWORKINFERENCE, to recover gene regulatory networks. NETWORKINFERENCE performed well at recovering feedback loops and multiple targets of a regulator with small amounts of data, but required more data to recover multiple regulators of a gene. When collecting the same number of data samples at different intervals from the system, the best recovery was produced by sampling intervals long enough such that sampling covered propagation of regulation through the network but not so long such that intervals missed internal dynamics. These results further elucidate the possibilities and limitations of network inference based on biological data.

Evaluating functional network inference using simulations of complex biological systems.

MOTIVATION: Although many network inference algorithms have been presented in the bioinformatics literature, no suitable approach has been formulated for evaluating their effectiveness at recovering models of complex biological systems from limited data. To overcome this limitation, we propose an approach to evaluate network inference algorithms according to their ability to recover a complex functional network from biologically reasonable simulated data. RESULTS: We designed a simulator to generate data representing a complex biological system at multiple levels of organization: behaviour, neural anatomy, brain electrophysiology, and gene expression of songbirds. About 90% of the simulated variables are unregulated by other variables in the system and are included simply as distracters. We sampled the simulated data at intervals as one would sample from a biological system in practice, and then used the sampled data to evaluate the effectiveness of an algorithm we developed for functional network inference. We found that our algorithm is highly effective at recovering the functional network structure of the simulated system-including the irrelevance of unregulated variables-from sampled data alone. To assess the reproducibility of these results, we tested our inference algorithm on 50 separately simulated sets of data and it consistently recovered almost perfectly the complex functional network structure underlying the simulated data. To our knowledge, this is the first approach for evaluating the effectiveness of functional network inference algorithms at recovering models from limited data. Our simulation approach also enables researchers a priori to design experiments and data-collection protocols that are amenable to functional network inference.

Enhancing the Hodgkin-Huxley Equations: Simulations Based on the First Publication in the Biophysical Journal.

The experiments in the Cole and Moore article in the first issue of the Biophysical Journal provided the first independent experimental confirmation of the Hodgkin-Huxley (HH) equations. A log-log plot of the K current versus time showed that raising the HH variable n to the sixth power provided the best fit to the data. Subsequent simulations using n(6) and setting the resting potential at the in vivo value simplifies the HH equations by eliminating the leakage term. Our article also reported that the K current in response to a depolarizing step to ENa was delayed if the step was preceded by a hyperpolarization. While the interpretation of this phenomenon in the article was flawed, subsequent simulations show that the effect completely arises from the original HH equations.

Modeling the effect of lead-free soldering on flexible substrates

Flexible Circuit Boards (FPCs) are now being widely used in the electronic industries especially in the areas of electronic packages. Due to European lead-free legislation which has been implemented since July 2006, electronic packaging industries have to switch to use in the lead-free soldering technology. This change has posed a number of challenges in terms of development of lead-free solders and compatible substrates. An increase of at least 20-50 degrees in the reflow temperature is a concern and substantial research is required to investigate a sustainable design of flexible circuit boards as carrier substrates. This paper investigates a number of design variables such as copper conductor width, type of substrate materials, effect of insulating materials, etc. Computer modeling has been used to investigate thermo-mechanical behavior, and reliability, of flexible substrates after they have been subjected to a lead- free solder processing. Results will show particular designs that behave better for a particular rise in peak reflow temperature. Also presented will be the types of failures that can occur in these substrates and what particular materials are more reliable.

The impact of the passenger response time distribution on ship evacuation performance

The International Maritime Organisation (IMO) has adopted the use of computer simulation to assist in the assessment of the assembly time for passenger ships. A key parameter required for this analysis and specified as part of the IMO guidelines is the passenger response time distribution. It is demonstrated in this paper that the IMO specified response time distribution assumes an unrealistic mathematical form. This unrealistic mathematical form can lead to serious congestion issues being overlooked in the evacuation analysis and lead to incorrect conclusions concerning the suitability of vessel design. In light of these results, it is vital that IMO undertake research to generate passenger response time data suitable for use in evacuation analysis of passenger ships. Until this type of data becomes readily available, it is strongly recommended that rather than continuing to use the artificial and unrepresentative form of the response time distribution, IMO should adopt plausible and more realistic response time data derived from land based applications. © 2005: Royal Institution of Naval Architects.

Bandstructure approach to near edge structure

We review the current state of the art in EELS fingerprinting by computer simulation, focusing on the bandstructure approach to the problem. Currently calculations are made using a one electron theory, but we describe in principle the way to go beyond this to include final state effects. We include these effects within the one electron framework using the Slater transition state formula and assess the errors involved. Two examples are then given which illustrate the use of the one electron approximation within density functional theory. Our approach is to combine predicted atomic structure with predicted electronic structure to assist in fingerprinting of complex crystal structures.

Surface alloying and mixing at the Mn/Fe(001) interface: Real-time photoelectron spectroscopy and modified embedded atom simulations

Structural and magnetic properties of thin Mn films on the Fe(001) surface have been investigated by a combination of photoelectron spectroscopy and computer simulation in the temperature range 300 Kless than or equal toTless than or equal to750 K. Room-temperature as deposited Mn overlayers are found to be ferromagnetic up to 2.5-monolayer (ML) coverage, with a magnetic moment parallel to that of the iron substrate. The Mn atomic moment decreases with increasing coverage, and thicker samples (4-ML and 4.5-ML coverage) are antiferromagnetic. Photoemission measurements performed while the system temperature is rising at constant rate (dT/dtsimilar to0.5 K/s) detect the first signs of Mn-Fe interdiffusion at T=450 K, and reveal a broad temperature range (610 Kless than or equal toTless than or equal to680 K) in which the interface appears to be stable. Interdiffusion resumes at Tgreater than or equal to680 K. Molecular dynamics and Monte Carlo simulations allow us to attribute the stability plateau at 610 Kless than or equal toTless than or equal to680 K to the formation of a single-layer MnFe surface alloy with a 2x2 unit cell and a checkerboard distribution of Mn and Fe atoms. X-ray-absorption spectroscopy and analysis of the dichroic signal show that the alloy has a ferromagnetic spin structure, collinear with that of the substrate. The magnetic moments of Mn and Fe atoms in the alloy are estimated to be 0.8mu(B) and 1.1mu(B), respectively.

Turbo channel estimation and equalisation for a superposition-based cooperative system

This study investigates the superposition-based cooperative transmission system. In this system, a key point is for the relay node to detect data transmitted from the source node. This issued was less considered in the existing literature as the channel is usually assumed to be flat fading and a priori known. In practice, however, the channel is not only a priori unknown but subject to frequency selective fading. Channel estimation is thus necessary. Of particular interest is the channel estimation at the relay node which imposes extra requirement for the system resources. The authors propose a novel turbo least-square channel estimator by exploring the superposition structure of the transmission data. The proposed channel estimator not only requires no pilot symbols but also has significantly better performance than the classic approach. The soft-in-soft-out minimum mean square error (MMSE) equaliser is also re-derived to match the superimposed data structure. Finally computer simulation results are shown to verify the proposed algorithm.

Predicting revision risk for aseptic loosening of femoral components in total flip arthroplasty in individual patients - A finite element study

Advances in surgical procedure, prosthesis design, and biomaterials performance have considerably increased the longevity of total joint replacements. Preoperative planning is another step in joint replacement that may have the potential to improve clinical outcome for the individual patient, but has remained relatively consistent for a longtime. One means of advancing this aspect of joint replacement surgery may be to include predictive computer simulation into the planning process. In this article, the potential of patient-specific finite element analysis in preoperative assessment is investigated. Seventeen patient-specific finite element models of cemented Charnley reconstructions were created, of which six were early (

Two-dimensional imaging of low temperature laser produced plasmas

Laser induced fluorescence images of a low temperature laser-produced plasma expanding into vacuum are presented and compared to a computer simulation. The complex nature of a plume expanding into background gas is highlighted, along with a potential means of simplifying the study of such systems.

A probabilistic modelling scheme for analysis of long-term failure of cemented femoral joint replacements

Reliable prediction of long-term medical device performance using computer simulation requires consideration of variability in surgical procedure, as well as patient-specific factors. However, even deterministic simulation of long-term failure processes for such devices is time and resource consuming so that including variability can lead to excessive time to achieve useful predictions. This study investigates the use of an accelerated probabilistic framework for predicting the likely performance envelope of a device and applies it to femoral prosthesis loosening in cemented hip arthroplasty.
A creep and fatigue damage failure model for bone cement, in conjunction with an interfacial fatigue model for the implant–cement interface, was used to simulate loosening of a prosthesis within a cement mantle. A deterministic set of trial simulations was used to account for variability of a set of surgical and patient factors, and a response surface method was used to perform and accelerate a Monte Carlo simulation to achieve an estimate of the likely range of prosthesis loosening. The proposed framework was used to conceptually investigate the influence of prosthesis selection and surgical placement on prosthesis migration.
Results demonstrate that the response surface method is capable of dramatically reducing the time to achieve convergence in mean and variance of predicted response variables. A critical requirement for realistic predictions is the size and quality of the initial training dataset used to generate the response surface and further work is required to determine the recommendations for a minimum number of initial trials. Results of this conceptual application predicted that loosening was sensitive to the implant size and femoral width. Furthermore, different rankings of implant performance were predicted when only individual simulations (e.g. an average condition) were used to rank implants, compared with when stochastic simulations were used. In conclusion, the proposed framework provides a viable approach to predicting realistic ranges of loosening behaviour for orthopaedic implants in reduced timeframes compared with conventional Monte Carlo simulations.

Using temporal analysis of products and flux response technology to determine diffusion coefficients in catalytic monoliths

The importance of accurately measuring gas diffusivity in porous materials has led to a number of methods being developed. In this study the Temporal Analysis of Products (TAP) reactor and Flux Response Technology (FRT) have been used to examine the diffusivity in the washcoat supported on cordierite monoliths. Herein, the molecular diffusion of propane within four monoliths with differently prepared alumina/CeZrOx washcoats was investigated as a function of temperature. Moment-based analysis of the observed TAP responses led to the calculation of the apparent intermediate gas constant, Kp, that characterises adsorption into the mesoporous network and apparent time delay, tapp, that characterises residence time in the mesoporous network. Additionally, FRT has been successfully adapted as an extensive in situ perturbation technique in measuring intraphase diffusion coefficients in the washcoats of the same four monolith samples. The diffusion coefficients obtained by moment-based analysis of TAP responses are larger than the coefficients determined by zero length column (ZLC) analysis of flux response profiles with measured values of the same monolith samples between 20 and 100 °C ranging from 2–5×10^-9 m² s^-1 to 4–8×10^-10 m² s^-1, respectively. The TAP and FRT data, therefore, provide a range of the lower and upper limits of diffusivity, respectively. The reported activation energies and diffusivities clearly correlate with the difference in the washcoat structure of different monolith samples.

Mechanical characterization of the varian Exact-arm and R-arm support systems for eight aS500 electronic portal imaging devices

The aim of this study is to compare the positioning accuracy at different gantry angles of two electronic portal imaging devices (EPIDs) support arm systems by using EPID difference images as a measure for displacement. This work presents a comparison of the mechanical performance of eight Varian aS500 (Varian Medical Systems, Palo Alto, CA) EPIDs, mounted using either the Varian Exact-arm or R-arm.

Pressure effect on vibrational frequency and dephasing of 1-alkyl-3-methylimidazolium hexafluorophosphate ionic liquids

Raman spectra in the range of the totally symmetric stretching mode of the [PF6]− anion, νs(PF6), have been measured for 1-alkyl-3-methylimidazolium ionic liquids [CnC1im][PF6], for n = 4, 6, and 8, as a function of pressure at room temperature. The ionic liquids [C6C1im][PF6] and [C8C1im][PF6] remain in an amorphous phase up to 3.5 GPa, in contrast to [C4C1im][PF6], whichcrystallizes above ∼0.5 GPa. Equations of state based either on a group contribution model or Carnahan-Starling-van der Waals model have been used to estimate the densities of the ionic liquids at high pressures. The shifts of the vibrational frequency of νs(PF6) with density observed in [C6C1im][PF6] and in [C8C1im][PF6] have been calculated by a hard-sphere model of a pseudo-diatomic solute under short-range repulsive interactions with the neighboring particles. The stochastic model of Kubo for vibrational dephasing has been used to obtain the amplitude of vibrational frequency fluctuation, ⟨Δω 2⟩, and the relaxation time of frequency fluctuation, τ c , as a function of density by Raman band shape analysis of the νs(PF6) mode of [C6C1im][PF6] and [C8C1im][PF6].

Estimation of fetal and effective dose for CT examinations

Doses from CT examinations are difficult to estimate. However, they are requested more frequently due to the increase in CT examinations. In particular, fetal dose estimations are frequently required for patients who have discovered, subsequent to the examination, that they were pregnant when the examination was conducted. A computer model has been developed to facilitate such dose calculations. This model combines empirical beam data with anatomical information. The model has been verified using thermoluminescent dosemeter (TLD) readings of internal and surface dose from both phantoms and patients, including intrauterine doses for patients undergoing afterloading gynaecological intracavitary treatment. Although only limited experimental data were available, the results indicate that the model accurately predicts uterine doses within acceptable errors. This approach has been validated for fetal dose estimation. The model was also used in a comparison with the nationally available CT dose data from the National Radiological Protection Board (NRPB). The two models were found to be in agreement for fetal dose estimations.