An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 1. Methodology for residual generation

For dynamic simulations to be credible, verification of the computer code must be an integral part of the modelling process. This two-part paper describes a novel approach to verification through program testing and debugging. In Part 1, a methodology is presented for detecting and isolating coding errors using back-to-back testing. Residuals are generated by comparing the output of two independent implementations, in response to identical inputs. The key feature of the methodology is that a specially modified observer is created using one of the implementations, so as to impose an error-dependent structure on these residuals. Each error can be associated with a fixed and known subspace, permitting errors to be isolated to specific equations in the code. It is shown that the geometric properties extend to multiple errors in either one of the two implementations. Copyright (C) 2003 John Wiley Sons, Ltd.

An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 2. Residuals analysis and a case study

In Part 1 of this paper a methodology for back-to-back testing of simulation software was described. Residuals with error-dependent geometric properties were generated. A set of potential coding errors was enumerated, along with a corresponding set of feature matrices, which describe the geometric properties imposed on the residuals by each of the errors. In this part of the paper, an algorithm is developed to isolate the coding errors present by analysing the residuals. A set of errors is isolated when the subspace spanned by their combined feature matrices corresponds to that of the residuals. Individual feature matrices are compared to the residuals and classified as 'definite', 'possible' or 'impossible'. The status of 'possible' errors is resolved using a dynamic subset testing algorithm. To demonstrate and validate the testing methodology presented in Part 1 and the isolation algorithm presented in Part 2, a case study is presented using a model for biological wastewater treatment. Both single and simultaneous errors that are deliberately introduced into the simulation code are correctly detected and isolated. Copyright (C) 2003 John Wiley Sons, Ltd.

Wurst: a protein threading server with a structural scoring function, sequence profiles and optimized substitution matrices

Wurst is a protein threading program with an emphasis on high quality sequence to structure alignments (http://www.zbh.uni-hamburg.de/wurst). Submitted sequences are aligned to each of about 3000 templates with a conventional dynamic programming algorithm, but using a score function with sophisticated structure and sequence terms. The structure terms are a log-odds probability of sequence to structure fragment compatibility, obtained from a Bayesian classification procedure. A simplex optimization was used to optimize the sequence-based terms for the goal of alignment and model quality and to balance the sequence and structural contributions against each other. Both sequence and structural terms operate with sequence profiles.

Weibull model selection for reliability modelling

A large number of models have been derived from the two-parameter Weibull distribution and are referred to as Weibull models. They exhibit a wide range of shapes for the density and hazard functions, which makes them suitable for modelling complex failure data sets. The WPP and IWPP plot allows one to determine in a systematic manner if one or more of these models are suitable for modelling a given data set. This paper deals with this topic.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Structure determination of the three disulfide bond isomers of alpha-conotoxin GI: A model for the role of disulfide bonds in structural stability

The three possible disulfide bonded isomers of alpha-conotoxin GI have been selectively synthesised and their structures determined by H-1 NMR spectroscopy. alpha-Conotoxin GI derives from the venom of Conus geographus and is a useful neuropharmacological tool as it selectively binds to the nicotinic acetylcholine receptor (nAChR), a ligand-gated ion channel involved in nerve signal transmission. The peptide has the sequence ECCNPACGRHYSC-NH2, and the three disulfide bonded isomers are referred to as GI(2-7;3-13), GI(2-13;3-7) and GI(2-3;7-13). The NMR structure for the native isomer GI(2-7;3-13) is of excellent quality, with a backbone pairwise RMSD of 0.16 Angstrom for a family of 35 structures, and comprises primarily a distorted 3(10),, helix between residues 5 to 11. The two non-native isomers exhibit multiple conformers in solution, with the major populated forms being different in structure both from each other and from the native form. Structure-activity relationships for the native GI(2-7;3-13) as well as the role of the disulfide bonds on folding and stability of the three isomers are examined. It is concluded that the disulfide bonds in alpha-conotoxin GI play a crucial part in determining both the structure and stability of the peptide. A trend for increased conformational heterogeneity was observed in the order of GI(2-7;3-13) < GI(2-13;3-7) < GI(2-3;7-13). It was found that the peptide bond joining Cys2 to Cys3 in GI(2-3;7-13) is predominantly trans, rather than cis as theoretically predicted. These structural data are used to interpret the varying nAChR binding of the non-native forms. A model for the binding of native GI(2-7;3-13) to the mammalian nAChR is proposed, with an alpha-subunit binding face made up of Cys2, Asn4, Pro5, Ala6 and Cys7 and a selectivity face, comprised of Arg9 and His10. These two faces orient the molecule between the alpha and delta subunits of the receptor. The structure of the CCNPAC sequence of the native GI(2-7;3-13) is compared to the structure of the identical sequence from the toxic domain of heat-stable enterotoxins, which forms part of the receptor binding region of the enterotoxins, but which has a different disulfide connectivity. (C) 1998 Academic Press Limited.

Perceptual quality metrics applied to still image compression

We present a review of perceptual image quality metrics and their application to still image compression. The review describes how image quality metrics can be used to guide an image compression scheme and outlines the advantages, disadvantages and limitations of a number of quality metrics. We examine a broad range of metrics ranging from simple mathematical measures to those which incorporate full perceptual models. We highlight some variation in the models for luminance adaptation and the contrast sensitivity function and discuss what appears to be a lack of a general consensus regarding the models which best describe contrast masking and error summation. We identify how the various perceptual components have been incorporated in quality metrics, and identify a number of psychophysical testing techniques that can be used to validate the metrics. We conclude by illustrating some of the issues discussed throughout the paper with a simple demonstration. (C) 1998 Elsevier Science B.V. All rights reserved.

The clinical value compass: Achieving benchmarking and quality improvement in aged care

Quality measurement and benchmarking in aged cave presents several challenges. A model which addresses this by linking four dimensions of outcomes has been developed - the Clinical Value Compass (CVC). A CVC was developed for stroke rehabilitation and measured across four sites. The CVC teas well accepted by the treatment teams and proved practical to measure. The results revealed differences in practices and client groups that led to a closer analysis of process and subsequent changes in these processes. Remeasuring of the CVC is required to demonstrate improved outcomes arising from these process changes.

Image calibration to like-values in mapping shallow water quality from multitemporal data

The applicability of image calibration to like-values in mapping water quality parameters from multitemporal images is explored, Six sets of water samples were collected at satellite overpasses over Moreton Bay, Brisbane, Australia. Analysis of these samples reveals that waters in this shallow bay are mostly TSS-dominated, even though they are occasionally dominated by chlorophyll as well. Three of the images were calibrated to a reference image based on invariant targets. Predictive models constructed from the reference image were applied to estimating total suspended sediment (TSS) and Secchi depth from another image at a discrepancy of around 35 percent. Application of the predictive model for TSS concentration to another image acquired at a time of different water types resulted in a discrepancy of 152 percent. Therefore, image calibration to like-values could be used to reliably map certain water quality parameters from multitemporal TM images so long as the water type under study remains unchanged. This method is limited in that the mapped results could be rather inaccurate if the water type under study has changed considerably. Thus, the approach needs to be refined in shallow water from multitemporal satellite imagery.

When the going gets rough: effect of maternal size manipulation on larval quality

Variation in larval size has been shown to be an important factor for the post-metamorphic performance of marine invertebrates but, despite its importance, few sources of this variation have been identified. For a range of taxa, offspring size is positively correlated with maternal size but the reasons for this correlation remain unclear. We halved the size of colonies in the bryozoan Bugula neritina 1 wk prior to reproduction (but during embryogenesis) to determine if larval size is a fixed or plastic trait. We manipulated colonies in such a way that the ratio of feeding zooids to reproductive zooids was constant between treatment and control colonies. We found that manipulating colony size strongly affects larval size; halved colonies produced larvae that were similar to13% smaller than those produced by intact colonies. We entered these data into a simple model based on previous work to estimate the likely post-metamorphic consequences of this reduction in larval size. The model predicted that larvae that came from manipulated colonies would suffer similar to300% higher post-metamorphic mortality and similar to50% lower fecundity as adults. Colonies that are faced with a stress appear to be trading off current offspring fitness to maximize their own long-term fitness and this may explain previous observations of compensatory growth in colonial organisms. This study demonstrates that larval size is a surprisingly dynamic trait and strong links exist between the maternal phenotype and the fitness of the offspring. The performance of settling larvae may be determined not only by their larval experience but also by the experience of their mothers.

Developing a meta model for the Bunge-Wand-Weber ontological constructs

At the core of the analysis task in the development process is information systems requirements modelling, Modelling of requirements has been occurring for many years and the techniques used have progressed from flowcharting through data flow diagrams and entity-relationship diagrams to object-oriented schemas today. Unfortunately, researchers have been able to give little theoretical guidance only to practitioners on which techniques to use and when. In an attempt to address this situation, Wand and Weber have developed a series of models based on the ontological theory of Mario Bunge-the Bunge-Wand-Weber (BWW) models. Two particular criticisms of the models have persisted however-the understandability of the constructs in the BWW models and the difficulty in applying the models to a modelling technique. This paper addresses these issues by presenting a meta model of the BWW constructs using a meta language that is familiar to many IS professionals, more specific than plain English text, but easier to understand than the set-theoretic language of the original BWW models. Such a meta model also facilitates the application of the BWW theory to other modelling techniques that have similar meta models defined. Moreover, this approach supports the identification of patterns of constructs that might be common across meta models for modelling techniques. Such findings are useful in extending and refining the BWW theory. (C) 2002 Elsevier Science Ltd. All rights reserved.

A general model for host plant selection in phytophagous insects

We develop a general theoretical framework for exploring the host plant selection behaviour of herbivorous insects. This model can be used to address a number of questions, including the evolution of specialists, generalists, preference hierarchies, and learning. We use our model to: (i) demonstrate the consequences of the extent to which the reproductive success of a foraging female is limited by the rate at which they find host plants (host limitation) or the number of eggs they carry (egg limitation); (ii) emphasize the different consequences of variation in behaviour before and after landing on (locating) a host (termed pre- and post-alighting, respectively); (iii) show that, in contrast to previous predictions, learning can be favoured in post-alighting behaviour-in particular, individuals can be selected to concentrate oviposition on an abundant low-quality host, whilst ignoring a rare higher-quality host; (iv) emphasize the importance of interactions between mechanisms in favouring specialization or learning. (C) 2002 Elsevier Science Ltd.

Assessment of the ability to model proteins with leucine-rich repeats in light of the latest structural information

The three-dimensional structures of leucine-rich repeat (LRR) -containing proteins from five different families were previously predicted based on the crystal structure of the ribonuclease inhibitor. using an approach that combined homology-based modeling, structure-based sequence alignment of LRRs, and several rational assumptions. The structural models have been produced based on very limited sequence similarity, which, in general. cannot yield trustworthy predictions. Recently, the protein structures from three of these five families have been determined. In this report we estimate the quality of the modeling approach by comparing the models with the experimentally determined structures. The comparison suggests that the general architecture, curvature, interior/exterior orientations of side chains. and backbone conformation of the LRR structures can be predicted correctly. On the other hand. the analysis revealed that, in some cases. it is difficult to predict correctly the twist of the overall super-helical structure. Taking into consideration the conclusions from these comparisons, we identified a new family of bacterial LRR proteins and present its structural model. The reliability of the LRR protein modeling suggests that it would be informative to apply similar modeling approaches to other classes of solenoid proteins.