gems: An R Package for Simulating from Disease Progression Models

Mathematical models of disease progression predict disease outcomes and are useful epidemiological tools for planners and evaluators of health interventions. The R package gems is a tool that simulates disease progression in patients and predicts the effect of different interventions on patient outcome. Disease progression is represented by a series of events (e.g., diagnosis, treatment and death), displayed in a directed acyclic graph. The vertices correspond to disease states and the directed edges represent events. The package gems allows simulations based on a generalized multistate model that can be described by a directed acyclic graph with continuous transition-specific hazard functions. The user can specify an arbitrary hazard function and its parameters. The model includes parameter uncertainty, does not need to be a Markov model, and may take the history of previous events into account. Applications are not limited to the medical field and extend to other areas where multistate simulation is of interest. We provide a technical explanation of the multistate models used by gems, explain the functions of gems and their arguments, and show a sample application.

A Mathematical Model for Simplifying Representations of Objects in a Geographic Information System

The study of operations on representations of objects is well documented in the realm of spatial engineering. However, the mathematical structure and formal proof of these operational phenomena are not thoroughly explored. Other works have often focused on query-based models that seek to order classes and instances of objects in the form of semantic hierarchies or graphs. In some models, nodes of graphs represent objects and are connected by edges that represent different types of coarsening operators. This work, however, studies how the coarsening operator "simplification" can manipulate partitions of finite sets, independent from objects and their attributes. Partitions that are "simplified first have a collection of elements filtered (removed), and then the remaining partition is amalgamated (some sub-collections are unified). Simplification has many interesting mathematical properties. A finite composition of simplifications can also be accomplished with some single simplification. Also, if one partition is a simplification of the other, the simplified partition is defined to be less than the other partition according to the simp relation. This relation is shown to be a partial-order relation based on simplification. Collections of partitions can not only be proven to have a partial- order structure, but also have a lattice structure and are complete. In regard to a geographic information system (GIs), partitions related to subsets of attribute domains for objects are called views. Objects belong to different views based whether or not their attribute values lie in the underlying view domain. Given a particular view, objects with their attribute n-tuple codings contained in the view are part of the actualization set on views, and objects are labeled according to the particular subset of the view in which their coding lies. Though the scope of the work does not mainly focus on queries related directly to geographic objects, it provides verification for the existence of particular views in a system with this underlying structure. Given a finite attribute domain, one can say with mathematical certainty that different views of objects are partially ordered by simplification, and every collection of views has a greatest lower bound and least upper bound, which provides the validity for exploring queries in this regard.

Two-dimensional direct rainfall hydraulic models applied to direct calculation of hydrologic events in dams to prevent overtopping

One of the main concerns when conducting a dam test is the acute determination of the hydrograph for a specific flood event. The use of 2D direct rainfall hydraulic mathematical models on a finite elements mesh, combined with the efficiency of vector calculus that provides CUDA (Compute Unified Device Architecture) technology, enables nowadays the simulation of complex hydrological models without the need for terrain subbasin and transit splitting (as in HEC-HMS). Both the Spanish PNOA (National Plan of Aereal Orthophotography) Digital Terrain Model GRID with a 5 x 5 m accuracy and the CORINE GIS Land Cover (Coordination of INformation of the Environment) that allows assessment of the ground roughness, provide enough data to easily build these kind of models

Application of the dynamic models to the detection of imperfections in plate Structures

This paper is part of a set of publications related with the development of mathematical models aimed to simulate the dynamic input and output of experimental nondestructive tests in order to detect structural imperfections. The structures to be considered are composed by steel plates of thin thickness. The imperfections in these cases are cracks and they can penetrate either a significant part of the plate thickness or be micro cracks or superficial imperfections. The first class of cracks is related with structural safety and the second one is more connected to the structural protection to the environment, particularly if protective paintings can be deteriorated. Two mathematical groups of models have been developed. The first group tries to locate the position and extension of the imperfection of the first class of imperfections, i.e. cracks and it is the object of the present paper. Bending Kirchoff thin plate models belong to this first group and they are used to this respect. The another group of models is dealt with membrane structures under the superficial Rayleigh waves excitation. With this group of models the micro cracks detection is intended. In the application of the first group of models to the detection of cracks, it has been observed that the differences between the natural frequencies of the non cracked and the cracked structures are very small. However, geometry and crack position can be identified quite accurately if this comparison is carried out between first derivatives (mode rotations) of the natural modes are used instead. Finally, in relation with the analysis of the superficial crack existence the use of Rayleigh waves is very promising. The geometry and the penetration of the micro crack can be detected very accurately. The mathematical and numerical treatment of the generation of these Rayleigh waves present and a numerical application has been shown.

Models of immune memory: On the role of cross-reactive stimulation, competition, and homeostasis in maintaining immune memory

There has been much debate on the contribution of processes such as the persistence of antigens, cross-reactive stimulation, homeostasis, competition between different lineages of lymphocytes, and the rate of cell turnover on the duration of immune memory and the maintenance of the immune repertoire. We use simple mathematical models to investigate the contributions of these various processes to the longevity of immune memory (defined as the rate of decline of the population of antigen-specific memory cells). The models we develop incorporate a large repertoire of immune cells, each lineage having distinct antigenic specificities, and describe the dynamics of the individual lineages and total population of cells. Our results suggest that, if homeostatic control regulates the total population of memory cells, then, for a wide range of parameters, immune memory will be long-lived in the absence of persistent antigen (T1/2 > 1 year). We also show that the longevity of memory in this situation will be insensitive to the relative rates of cross-reactive stimulation, the rate of turnover of immune cells, and the functional form of the term for the maintenance of homeostasis.

USLE and WEPP Soil Erosion Models: Their Applicability and Use in Generating Cover Values to Stabilize Soil in the Los Alamos, New Mexico Area

Soil erosion is a naturally occurring process that involves the detachment, transport, and deposition of soil particles. Disturbances such as thinning and wildfire can reduce cover greatly and increase erosion rates. Forest managers may use erosion prediction tools, such as the Universal Soil Loss Equation (USLE) and Water Erosion Prediction Project (WEPP) to estimate erosion rates and develop techniques to manage erosion. However, it is important to understand the differences and the applications of each model. Erosion rates were generated by each model and the model most applicable to the study site, Los Alamos, New Mexico was determined. It was also used to find the amount of cover needed to stabilize soil. The USLE is a simpler model and less complicated than a computer model like WEPP, and thus easier to manipulate to estimate cover values. Predicted cover values were compared to field cover values. Cover is necessary to establish effective erosion control guidelines.

Integration of different models in the design of chemical processes: Application to the design of a power plant

With advances in the synthesis and design of chemical processes there is an increasing need for more complex mathematical models with which to screen the alternatives that constitute accurate and reliable process models. Despite the wide availability of sophisticated tools for simulation, optimization and synthesis of chemical processes, the user is frequently interested in using the ‘best available model’. However, in practice, these models are usually little more than a black box with a rigid input–output structure. In this paper we propose to tackle all these models using generalized disjunctive programming to capture the numerical characteristics of each model (in equation form, modular, noisy, etc.) and to deal with each of them according to their individual characteristics. The result is a hybrid modular–equation based approach that allows synthesizing complex processes using different models in a robust and reliable way. The capabilities of the proposed approach are discussed with a case study: the design of a utility system power plant that has been decomposed into its constitutive elements, each treated differently numerically. And finally, numerical results and conclusions are presented.

Saint Mathew Law and Bonini Paradox in Textual Theory of Complex Models

The mathematical models of the complex reality are texts belonging to a certain literature that is written in a semi-formal language, denominated L(MT) by the authors whose laws linguistic mathematics have been previously defined. This text possesses linguistic entropy that is the reflection of the physical entropy of the processes of real world that said text describes. Through the temperature of information defined by Mandelbrot, the authors begin a text-reality thermodynamic theory that drives to the existence of information attractors, or highly structured point, settling down a heterogeneity of the space text, the same one that of ontologic space, completing the well-known law of Saint Mathew, of the General Theory of Systems and formulated by Margalef saying: “To the one that has more he will be given, and to the one that doesn't have he will even be removed it little that it possesses.

Comparisons of models predicting ambient lake phosphorus concentrations /

Mode of access: Internet.

Comparison of numerical and physical hydraulic models, Masonboro Inlet, North Carolina /

"June 1977."

Mathematical instruction in France.

Appendices C and D "for later publication."

Models to illustrate the foundations of mathematics /

Mode of access: Internet.

Modelling of froth transportation in industrial flotation cells: Part I. Development of froth transportation models for attached particles

Modelling of froth transportation, as part of modelling of froth recovery, provides a scale-up procedure for flotation cell design. It can also assist in improving control of flotation operation. Mathematical models of froth velocity on the surface and froth residence time distribution in a cylindrical tank flotation cell are proposed, based on mass balance principle of the air entering the froth. The models take into account factors such as cell size, concentrate launder configuration, use of a froth crowder, cell operating conditions including froth height and air rate, and bubble bursting on the surface. (C) 2004 Elsevier Ltd. All rights reserved.

Using the canonical modelling approach to simplify the simulation of function in functional-structural plant models

Functional-structural plant models that include detailed mechanistic representation of underlying physiological processes can be expensive to construct and the resulting models can also be extremely complicated. On the other hand, purely empirical models are not able to simulate plant adaptability and response to different conditions. In this paper, we present an intermediate approach to modelling plant function that can simulate plant response without requiring detailed knowledge of underlying physiology. Plant function is modelled using a 'canonical' modelling approach, which uses compartment models with flux functions of a standard mathematical form, while plant structure is modelled using L-systems. Two modelling examples are used to demonstrate that canonical modelling can be used in conjunction with L-systems to create functional-structural plant models where function is represented either in an accurate and descriptive way, or in a more mechanistic and explanatory way. We conclude that canonical modelling provides a useful, flexible and relatively simple approach to modelling plant function at an intermediate level of abstraction.

A mathematical model of integrin-mediated haptotactic cell migration

Haptotactic cell migration, a directed response to gradients of cell—extracellular matrix adhesion, is an important process in a number of biological phenomena such as wound healing and tumour cell invasion. Previously, mathematical models of haptotaxis have been developed on the premise that cells migrate in response to gradients in the density of the extracellular matrix. In this paper, we develop a novel mathematical model of haptotaxis which includes the adhesion receptors known as integrins and a description of their functional activation, local recruitment and protrusion as part of lamellipodia. Through the inclusion of integrins, the modelled cell matter is able to respond to a true gradient of cell–matrix adhesion, represented by functionally active integrins. We also show that previous matrix-mediated models are in fact a subset of the novel integrin-mediated models, characterised by specific choices of diffusion and haptotaxis coefficients in their model equations. Numerical solutions suggest the existence of travelling waves of cell migration that are confirmed via a phase plane analysis of a simplified model.