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.

Dynamical analysis in growth models: Blumberg’s equation

We present a new dynamical approach to the Blumberg's equation, a family of unimodal maps. These maps are proportional to Beta(p, q) probability densities functions. Using the symmetry of the Beta(p, q) distribution and symbolic dynamics techniques, a new concept of mirror symmetry is defined for this family of maps. The kneading theory is used to analyze the effect of such symmetry in the presented models. The main result proves that two mirror symmetric unimodal maps have the same topological entropy. Different population dynamics regimes are identified, when the intrinsic growth rate is modified: extinctions, stabilities, bifurcations, chaos and Allee effect. To illustrate our results, we present a numerical analysis, where are demonstrated: monotonicity of the topological entropy with the variation of the intrinsic growth rate, existence of isentropic sets in the parameters space and mirror symmetry.

From weak Allee effect to no Allee effect in Richards’ growth models

Population dynamics have been attracting interest since many years. Among the considered models, the Richards’ equations remain one of the most popular to describe biological growth processes. On the other hand, Allee effect is currently a major focus of ecological research, which occurs when positive density dependence dominates at low densities. In this chapter, we propose the dynamical study of classes of functions based on Richards’ models describing the existence or not of Allee effect. We investigate bifurcation structures in generalized Richards’ functions and we look for the conditions in the (β, r) parameter plane for the existence of a weak Allee effect region. We show that the existence of this region is related with the existence of a dovetail structure. When the Allee limit varies, the weak Allee effect region disappears when the dovetail structure also disappears. Consequently, we deduce the transition from the weak Allee effect to no Allee effect to this family of functions. To support our analysis, we present fold and flip bifurcation curves and numerical simulations of several bifurcation diagrams.

Landslide susceptibility assessment in Karanganyar regency - Indonesia - Comparison of knowledge-based and Data-driven Models

Dissertation submitted in partial fulfillment of the requirements for the Degree of Master of Science in Geospatial Technologies.

Entrepreneurial Orientation, Intangible Assets and Firm Growth: the impact of ‘Spirit and Material’ on the growth of Chinese private firms

This paper has three contributions. First, it shows how field work within small firms in PR Chinese has provided new evidence which enables us to measure and calibrate Entrepreneurial Orientation (EO), as ‘spirit’, and Intangible Assets (IA), as ‘material’, for use in models of small firm growth. Second, it uses inter-item correlation analysis and both exploratory and confirmatory factor analysis to provide new measures of EO and IA, in index and in vector form, for use in econometric models of firm growth. Third, it estimates two new econometric models of small firm employment growth in PR China, under the null hypothesis of Gibrat’s Law, using our two new index-based and vector-based measures of EO and IA. Estimation is by OLS with adjustment for heteroscedasticity, and for sample selectivity. Broadly, it finds that EO attributes have had little significant impact on small firm growth, and indeed innovativeness and pro-activity paradoxically may even dampen growth. However, IA attributes have had a positive and significant impact on growth, with networking, and technological knowledge being of prime importance, and intellectual property and human capital being of lesser but still significant importance. In the light of these results, Gibrat’s Law is generalized, and Jovanovic’s learning theory is extended, to emphasise the importance of IA to growth. These findings cast new empirical light on the oft-quoted national slogan in PR China of “spirit and material”. So far as small firms are concerned, this paper suggests that their contribution to PR China’s remarkable economic growth is not so much attributable to the ‘spirit’ of enterprise (as suggested by propaganda) as, more prosaically, to the pursuit of the ‘material’.

Two different 3D biogeochemical models for the NW mediterranean sea: validation with Meris data and intercomparison

Report for the scientific sojourn carried out at the University of New South Wales from February to June the 2007. Two different biogeochemical models are coupled to a three dimensional configuration of the Princeton Ocean Model (POM) for the Northwestern Mediterranean Sea (Ahumada and Cruzado, 2007). The first biogeochemical model (BLANES) is the three-dimensional version of the model described by Bahamon and Cruzado (2003) and computes the nitrogen fluxes through six compartments using semi-empirical descriptions of biological processes. The second biogeochemical model (BIOMEC) is the biomechanical NPZD model described in Baird et al. (2004), which uses a combination of physiological and physical descriptions to quantify the rates of planktonic interactions. Physical descriptions include, for example, the diffusion of nutrients to phytoplankton cells and the encounter rate of predators and prey. The link between physical and biogeochemical processes in both models is expressed by the advection-diffusion of the non-conservative tracers. The similarities in the mathematical formulation of the biogeochemical processes in the two models are exploited to determine the parameter set for the biomechanical model that best fits the parameter set used in the first model. Three years of integration have been carried out for each model to reach the so called perpetual year run for biogeochemical conditions. Outputs from both models are averaged monthly and then compared to remote sensing images obtained from sensor MERIS for chlorophyll.

The Role of Institutions in Cross-Section Income and Panel Data Growth Models: A Deeper Investigation on the Weakness and Proliferation of Instruments

This paper investigates the role of institutions in determining per capita income levels and growth. It contributes to the empirical literature by using different variables as proxies for institutions and by developing a deeper analysis of the issues arising from the use of weak and too many instruments in per capita income and growth regressions. The cross-section estimation suggests that institutions seem to matter, regardless if they are the only explanatory variable or are combined with geographical and integration variables, although most models suffer from the issue of weak instruments. The results from the growth models provides some interesting results: there is mixed evidence on the role of institutions and such evidence is more likely to be associated with law and order and investment profile; government spending is an important policy variable; collapsing the number of instruments results in fewer significant coefficients for institutions.

Stability of Growth Models with Generalised Lag Structures

This paper considers the lag structures of dynamic models in economics, arguing that the standard approach is too simple to capture the complexity of actual lag structures arising, for example, from production and investment decisions. It is argued that recent (1990s) developments in the the theory of functional differential equations provide a means to analyse models with generalised lag structures. The stability and asymptotic stability of two growth models with generalised lag structures are analysed. The paper concludes with some speculative discussion of time-varying parameters.

Improvement in growth after two years of growth hormone therapy in very young children born small for gestational age and without spontaneous catch-up growth: results of a multicenter, controlled, randomized, open clinical trial.

CONTEXT GH treatment is effective in children born small for gestational age (SGA); however, its effectiveness and safety in very young SGA children is unknown. OBJECTIVE The aim was to analyze the outcome of very young SGA children treated with GH and followed for 2 yr. The results after 24 months of treatment, compared with a control group without treatment during 12 months followed by 12 months of treatment, are shown. DESIGN We performed a multicenter, controlled, randomized, open trial. SETTINGS The pediatric endocrinology departments of 14 public hospitals in Spain participated in the study. PATIENTS Seventy-six children, aged 2-5 yr born SGA and without catch-up growth, were studied. INTERVENTION Children received GH at 0.06 mg/kg.d for 2 yr (group I) or were followed for 12 months with no treatment and then treated for 12 months (group II). MAIN OUTCOME MEASURES Age, general health status, pubertal stage, bone age, height, weight, biochemical and hormonal analyses, and adverse side effects were determined at biannual check-ups. RESULTS The mean height sd score gain for chronological age in children treated for 24 months (group I) was 2.10, whereas in those treated only during the last 12 months (group II) was 1.43. In both groups, children under 4 yr of age had the greatest gain in growth velocity. No significant acceleration of bone age or side effects related to treatment was seen. CONCLUSION Very young SGA children without spontaneous catch-up growth could benefit from GH treatment because growth was accelerated and no negative side effects were observed.

The chemokine receptor CXCR4 strongly promotes neuroblastoma primary tumour and metastatic growth, but not invasion.

Neuroblastoma (NB) is a heterogeneous, and particularly malignant childhood neoplasm in its higher stages, with a propensity to form metastasis in selected organs, in particular liver and bone marrow, and for which there is still no efficient treatment available beyond surgery. Recent evidence indicates that the CXCR4/CXCL12 chemokine/receptor axis may be involved in promoting NB invasion and metastasis. In this study, we explored the potential role of CXCR4 in the malignant behaviour of NB, using a combination of in vitro functional analyses and in vivo growth and metastasis assessment in an orthotopic NB mouse model. We show here that CXCR4 overexpression in non-metastatic CXCR4-negative NB cells IGR-NB8 and in moderately metastatic, CXCR4 expressing NB cells IGR-N91, strongly increased tumour growth of primary tumours and liver metastases, without altering the frequency or the pattern of metastasis. Moreover shRNA-mediated knock-down experiments confirmed our observations by showing that silencing CXCR4 in NB cells impairs in vitro and almost abrogates in vivo growth. High levels of CXCL12 were detected in the mouse adrenal gland (the primary tumour site), and in the liver suggesting a paracrine effect of host-derived CXCL12 on NB growth. In conclusion, this study reveals a yet unreported NB-specific predominant growth and survival-promoting role of CXCR4, which warrants a critical reconsideration of the role of CXCR4 in the malignant behaviour of NB and other cancers.

One, two, (three), infinity: Newspaper and lab beauty-contest experiments

"Beauty-contest" is a game in which participants have to choose, typically, a number in [0,100], the winner being the person whose number is closest to a proportion of the average of all chosen numbers. We describe and analyze Beauty-contest experiments run in newspapers in UK, Spain, and Germany and find stable patterns of behavior across them, despite the uncontrollability of these experiments. These results are then compared with lab experiments involving undergraduates and game theorists as subjects, in what must be one of the largest empirical corroborations of interactive behavior ever tried. We claim that all observed behavior, across a wide variety of treatments and subject pools, can be interpretedas iterative reasoning. Level-1 reasoning, Level-2 reasoning and Level-3 reasoning are commonly observed in all the samples, while the equilibrium choice (Level-Maximum reasoning) is only prominently chosen by newspaper readers and theorists. The results show the empirical power of experiments run with large subject-pools, and open the door for more experimental work performed on the rich platform offered by newspapers and magazines.

Local and global error models to improve uncertainty quantification

In groundwater applications, Monte Carlo methods are employed to model the uncertainty on geological parameters. However, their brute-force application becomes computationally prohibitive for highly detailed geological descriptions, complex physical processes, and a large number of realizations. The Distance Kernel Method (DKM) overcomes this issue by clustering the realizations in a multidimensional space based on the flow responses obtained by means of an approximate (computationally cheaper) model; then, the uncertainty is estimated from the exact responses that are computed only for one representative realization per cluster (the medoid). Usually, DKM is employed to decrease the size of the sample of realizations that are considered to estimate the uncertainty. We propose to use the information from the approximate responses for uncertainty quantification. The subset of exact solutions provided by DKM is then employed to construct an error model and correct the potential bias of the approximate model. Two error models are devised that both employ the difference between approximate and exact medoid solutions, but differ in the way medoid errors are interpolated to correct the whole set of realizations. The Local Error Model rests upon the clustering defined by DKM and can be seen as a natural way to account for intra-cluster variability; the Global Error Model employs a linear interpolation of all medoid errors regardless of the cluster to which the single realization belongs. These error models are evaluated for an idealized pollution problem in which the uncertainty of the breakthrough curve needs to be estimated. For this numerical test case, we demonstrate that the error models improve the uncertainty quantification provided by the DKM algorithm and are effective in correcting the bias of the estimate computed solely from the MsFV results. The framework presented here is not specific to the methods considered and can be applied to other combinations of approximate models and techniques to select a subset of realizations

Monte Carlo study of the relation between vacancy diffusion and domain growth in two-dimensional binary alloys

Domain growth in a two-dimensional binary alloy is studied by means of Monte Carlo simulation of an ABV model. The dynamics consists of exchanges of particles with a small concentration of vacancies. The influence of changing the vacancy concentration and finite-size effects has been analyzed. Features of the vacancy diffusion during domain growth are also studied. The anomalous character of the diffusion due to its correlation with local order is responsible for the obtained fast-growth behavior.

Individual-tree basal area growth models for scots pine, pubescent birch and Norway spruce on drained peatlands in Finland

Summary

On the problematic of reserving and stochastic loss models

Abstract Traditionally, the common reserving methods used by the non-life actuaries are based on the assumption that future claims are going to behave in the same way as they did in the past. There are two main sources of variability in the processus of development of the claims: the variability of the speed with which the claims are settled and the variability between the severity of the claims from different accident years. High changes in these processes will generate distortions in the estimation of the claims reserves. The main objective of this thesis is to provide an indicator which firstly identifies and quantifies these two influences and secondly to determine which model is adequate for a specific situation. Two stochastic models were analysed and the predictive distributions of the future claims were obtained. The main advantage of the stochastic models is that they provide measures of variability of the reserves estimates. The first model (PDM) combines one conjugate family Dirichlet - Multinomial with the Poisson distribution. The second model (NBDM) improves the first one by combining two conjugate families Poisson -Gamma (for distribution of the ultimate amounts) and Dirichlet Multinomial (for distribution of the incremental claims payments). It was found that the second model allows to find the speed variability in the reporting process and development of the claims severity as function of two above mentioned distributions' parameters. These are the shape parameter of the Gamma distribution and the Dirichlet parameter. Depending on the relation between them we can decide on the adequacy of the claims reserve estimation method. The parameters have been estimated by the Methods of Moments and Maximum Likelihood. The results were tested using chosen simulation data and then using real data originating from the three lines of business: Property/Casualty, General Liability, and Accident Insurance. These data include different developments and specificities. The outcome of the thesis shows that when the Dirichlet parameter is greater than the shape parameter of the Gamma, resulting in a model with positive correlation between the past and future claims payments, suggests the Chain-Ladder method as appropriate for the claims reserve estimation. In terms of claims reserves, if the cumulated payments are high the positive correlation will imply high expectations for the future payments resulting in high claims reserves estimates. The negative correlation appears when the Dirichlet parameter is lower than the shape parameter of the Gamma, meaning low expected future payments for the same high observed cumulated payments. This corresponds to the situation when claims are reported rapidly and fewer claims remain expected subsequently. The extreme case appears in the situation when all claims are reported at the same time leading to expectations for the future payments of zero or equal to the aggregated amount of the ultimate paid claims. For this latter case, the Chain-Ladder is not recommended.