Mutation analysis to verify feature matrices for isolating errors in simulation models

Software simulation models are computer programs that need to be verified and debugged like any other software. In previous work, a method for error isolation in simulation models has been proposed. The method relies on a set of feature matrices that can be used to determine which part of the model implementation is responsible for deviations in the output of the model. Currrently these feature matrices have to be generated by hand from the model implementation, which is a tedious and error-prone task. In this paper, a method based on mutation analysis, as well as prototype tool support for the verification of the manually generated feature matrices is presented. The application of the method and tool to a model for wastewater treatment shows that the feature matrices can be verified effectively using a minimal number of mutants.

Artificial pumping errors in the Kool-Parker scaling model of soil moisture hysteresis

Hysteresis models that eliminate the artificial pumping errors associated with the Kool-Parker (KP) soil moisture hysteresis model, such as the Parker-Lenhard (PL) method, can be computationally demanding in unsaturated transport models since they need to retain the wetting-drying history of the system. The pumping errors in these models need to be eliminated for correct simulation of cyclical systems (e.g. transport above a tidally forced watertable, infiltration and redistribution under periodic irrigation) if the soils exhibit significant hysteresis. A modification is made here to the PL method that allows it to be more readily applied to numerical models by eliminating the need to store a large number of soil moisture reversal points. The modified-PL method largely eliminates any artificial pumping error and so essentially retains the accuracy of the original PL approach. The modified-PL method is implemented in HYDRUS-1D (version 2.0), which is then used to simulate cyclic capillary fringe dynamics to show the influence of removing artificial pumping errors and to demonstrate the ease of implementation. Artificial pumping errors are shown to be significant for the soils and system characteristics used here in numerical experiments of transport above a fluctuating watertable. (c) 2005 Elsevier B.V. All rights reserved.

Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

A uniform white-noise model for fixed-point roundoff errors in digital systems

Fixed-point roundoff noise in digital implementation of linear systems arises due to overflow, quantization of coefficients and input signals, and arithmetical errors. In uniform white-noise models, the last two types of roundoff errors are regarded as uniformly distributed independent random vectors on cubes of suitable size. For input signal quantization errors, the heuristic model is justified by a quantization theorem, which cannot be directly applied to arithmetical errors due to the complicated input-dependence of errors. The complete uniform white-noise model is shown to be valid in the sense of weak convergence of probabilistic measures as the lattice step tends to zero if the matrices of realization of the system in the state space satisfy certain nonresonance conditions and the finite-dimensional distributions of the input signal are absolutely continuous.

Factors predictive of intravenous fluid administration errors in Australian surgical care wards

Background: Intravenous (IV) fluid administration is an integral component of clinical care. Errors in administration can cause detrimental patient outcomes and increase healthcare costs, although little is known about medication administration errors associated with continuous IV infusions. Objectives: ( 1) To ascertain the prevalence of medication administration errors for continuous IV infusions and identify the variables that caused them. ( 2) To quantify the probability of errors by fitting a logistic regression model to the data. Methods: A prospective study was conducted on three surgical wards at a teaching hospital in Australia. All study participants received continuous infusions of IV fluids. Parenteral nutrition and non-electrolyte containing intermittent drug infusions ( such as antibiotics) were excluded. Medication administration errors and contributing variables were documented using a direct observational approach. Results: Six hundred and eighty seven observations were made, with 124 (18.0%) having at least one medication administration error. The most common error observed was wrong administration rate. The median deviation from the prescribed rate was 247 ml/h (interquartile range 275 to + 33.8 ml/ h). Errors were more likely to occur if an IV infusion control device was not used and as the duration of the infusion increased. Conclusions: Administration errors involving continuous IV infusions occur frequently. They could be reduced by more common use of IV infusion control devices and regular checking of administration rates.

Potential errors in the application of mixture theory to multifrequency bioelectrical impedance analysis

Potential errors in the application of mixture theory to the analysis of multiple-frequency bioelectrical impedance data for the determination of body fluid volumes are assessed. Potential sources of error include: conductive length; tissue fluid resistivity; body density; weight and technical errors of measurement. Inclusion of inaccurate estimates of body density and weight introduce errors of typically < +/-3% but incorrect assumptions regarding conductive length or fluid resistivities may each incur errors of up to 20%.

Errors in lamina growth of primary olfactory axons in the rat and mouse olfactory bulb

In the adult olfactory nerve pathway of rodents, each primary olfactory axon forms a terminal arbor in a single glomerulus in the olfactory bulb. During development, axons are believed to project directly to and terminate precisely within a glomerulus without any exuberant growth or mistargeting. To gain insight into mechanisms underlying this process, the trajectories of primary olfactory axons during glomerular formation were studied in the neonatal period. Histochemical staining of mouse olfactory bulb sections with the lectin Dolichos biflorus-agglutinin revealed that many olfactory axons overshoot the glomerular layer and course into the deeper laminae of the bulb in the early postnatal period. Single primary olfactory axons were anterogradely labelled either with the lipophilic carbocyanine dye, 1,1'-dioctodecyl-3,3,3',3'-tetramethylindocarbocyanine perchlorate (DiI), or with horseradish peroxidase (HRP) by localized microinjections into the nerve fiber layer of the rat olfactory bulb. Five distinct trajectories of primary olfactory axons were observed in DLI-labelled preparations at postnatal day 1.5 (P1.5). Axons either coursed directly to and terminated specifically within a glomerulus, branched before terminating in a glomerulus, bypassed glomeruli and entered the underlying external plexiform layer, passed through the glomerular layer with side branches into glomeruli, or branched into more than one glomerulus. HRP-labelled axon arbors from eight postnatal ages were reconstructed by camera lucida and were used to determine arbor length, arbor area, and arbor branch number. Whereas primary olfactory axons display errors in laminar targeting in the mammalian olfactory bulb, axon arbors typically achieve their adult morphology without exuberant growth. Many olfactory axons appear not to recognize appropriate cues to terminate within the glomerular layer during the early postnatal period. However, primary olfactory axons exhibit precise targeting in the glomerular layer after P5.5, indicating temporal differences in either the presence of guidance cues or the ability of axons to respond to these cues. (C) 1999 Wiley-Liss, Inc.

In vivo and in vitro models demonstrate a role for caveolin-1 in the pathogenesis of ischaemic acute renal failure

Caveolae and their proteins, the caveolins, transport macromolecules; compartmentalize signalling molecules; and are involved in various repair processes. There is little information regarding their role in the pathogenesis of significant renal syndromes such as acute renal failure (ARF). In this study, an in vivo rat model of 30 min bilateral renal ischaemia followed by reperfusion times from 4 h to 1 week was used to map the temporal and spatial association between caveolin-1 and tubular epithelial damage (desquamation, apoptosis, necrosis). An in vitro model of ischaemic ARF was also studied, where cultured renal tubular epithelial cells or arterial endothelial cells were subjected to injury initiators modelled on ischaemia-reperfusion (hypoxia, serum deprivation, free radical damage or hypoxia-hyperoxia). Expression of caveolin proteins was investigated using immunohistochemistry, immunoelectron microscopy, and immunoblots of whole cell, membrane or cytosol protein extracts. In vivo, healthy kidney had abundant caveolin-1 in vascular endothelial cells and also some expression in membrane surfaces of distal tubular epithelium. In the kidneys of ARF animals, punctate cytoplasmic localization of caveolin-1 was identified, with high intensity expression in injured proximal tubules that were losing basement membrane adhesion or were apoptotic, 24 h to 4 days after ischaemia-reperfusion. Western immunoblots indicated a marked increase in caveolin-1 expression in the cortex where some proximal tubular injury was located. In vitro, the main treatment-induced change in both cell types was translocation of caveolin-1 from the original plasma membrane site into membrane-associated sites in the cytoplasm. Overall, expression levels did not alter for whole cell extracts and the protein remained membrane-bound, as indicated by cell fractionation analyses. Caveolin-1 was also found to localize intensely within apoptotic cells. The results are indicative of a role for caveolin-1 in ARF-induced renal injury. Whether it functions for cell repair or death remains to be elucidated. Copyright (C) 2003 John Wiley Sons, Ltd.

Pseudo-randomness of round-off errors in discretized linear maps on the plane

We analyze the sequences of round-off errors of the orbits of a discretized planar rotation, from a probabilistic angle. It was shown [Bosio & Vivaldi, 2000] that for a dense set of parameters, the discretized map can be embedded into an expanding p-adic dynamical system, which serves as a source of deterministic randomness. For each parameter value, these systems can generate infinitely many distinct pseudo-random sequences over a finite alphabet, whose average period is conjectured to grow exponentially with the bit-length of the initial condition (the seed). We study some properties of these symbolic sequences, deriving a central limit theorem for the deviations between round-off and exact orbits, and obtain bounds concerning repetitions of words. We also explore some asymptotic problems computationally, verifying, among other things, that the occurrence of words of a given length is consistent with that of an abstract Bernoulli sequence.

Evaluation of individual protein errors in silver-stained two-dimensional gels

The relationship between spot volume and variation for all protein spots observed on large format 2D gels when utilising silver stain technology and a model system based on mammalian NSO cell extracts is reported. By running multiple gels we have shown that the reproducibility of data generated in this way is dependent on individual protein spot volumes, which in turn are directly correlated with the coefficient of variation. The coefficients of variation across all observed protein spots were highest for low abundant proteins which are the primary contributors to process error, and lowest for more abundant proteins. Using the relationship between spot volume and coefficient of variation we show it is necessary to calculate variation for individual protein spot volumes. The inherent limitations of silver staining therefore mean that errors in individual protein spot volumes must be considered when assessing significant changes in protein spot volume and not global error. (C) 2003 Elsevier Science (USA). All rights reserved.

On parameter estimation in population models

We describe methods for estimating the parameters of Markovian population processes in continuous time, thus increasing their utility in modelling real biological systems. A general approach, applicable to any finite-state continuous-time Markovian model, is presented, and this is specialised to a computationally more efficient method applicable to a class of models called density-dependent Markov population processes. We illustrate the versatility of both approaches by estimating the parameters of the stochastic SIS logistic model from simulated data. This model is also fitted to data from a population of Bay checkerspot butterfly (Euphydryas editha bayensis), allowing us to assess the viability of this population. (c) 2006 Elsevier Inc. All rights reserved.

Elasticity, yielding and episodicity in simple models of mantle convection

We explore the implications of refinements in the mechanical description of planetary constituents on the convection modes predicted by finite-element simulations. The refinements consist in the inclusion of incremental elasticity, plasticity (yielding) and multiple simultaneous creep mechanisms in addition to the usual visco-plastic models employed in the context of unified plate-mantle models. The main emphasis of this paper rests on the constitutive and computational formulation of the model. We apply a consistent incremental formulation of the non-linear governing equations avoiding the computationally expensive iterations that are otherwise necessary to handle the onset of plastic yield. In connection with episodic convection simulations, we point out the strong dependency of the results on the choice of the initial temperature distribution. Our results also indicate that the inclusion of elasticity in the constitutive relationships lowers the mechanical energy associated with subduction events.