Econometric analysis of microscopic simulation models

Microscopic simulation models are often evaluated based on visual inspection of the results. This paper presents formal econometric techniques to compare microscopic simulation (MS) models with real-life data. A related result is a methodology to compare different MS models with each other. For this purpose, possible parameters of interest, such as mean returns, or autocorrelation patterns, are classified and characterized. For each class of characteristics, the appropriate techniques are presented. We illustrate the methodology by comparing the MS model developed by He and Li [J. Econ. Dynam. Control, 2007, 31, 3396-3426, Quant. Finance, 2008, 8, 59-79] with actual data.

Integrating field sports, hare population management and conservation

Conflicts between field sports, animal welfare and species conservation are frequently contentious. In Ireland, the Irish Coursing Club (ICC) competitively tests the speed and agility of two greyhounds by using a live hare as a lure. Each coursing club is associated with a number of discrete localities, known as preserves, which are managed favourably for hares including predator control, prohibition of other forms of hunting such as shooting and poaching and the maintenance and enhancement of suitable hare habitat. We indirectly tested the efficacy of such management by comparing hare abundance within preserves to that in the wider countryside. In real terms, mean hare density was 18 times higher, and after controlling for variance in habitat remained 3 times higher, within ICC preserves than the wider countryside. Whilst we cannot rule out the role of habitat, our results suggest that hare numbers are maintained at high levels in ICC preserves either because clubs select areas of high hare density and subsequently have a negligible effect on numbers or that active population management positively increases hare abundance. The Irish hare Lepus timidus hibernicus Bell, 1837 is one of the highest priority species for conservation action in Ireland and without concessions for its role in conservation, any change in the legal status Of hare coursing under animal welfare grounds, may necessitate an increase in Government subsidies for conservation on private land together with a strengthened capacity for legislation enforcement.

Detailed analytical investigation of magnetic field line random walk in turbulent plasmas: I. Two-component slab/two-dimensional turbulence

The random displacement of magnetic field lines in the presence of magnetic turbulence in plasmas is investigated from first principles. A two-component (slab/two-dimensional composite) model for the turbulence spectrum is employes. An analytical investigation of the asymptotic behavior of the field-line mean square displacement (FL-MSD) is carried out. It is shown that the magnetic field lines behave superdifusively for every large values of the position variable z, since the FL-MSD sigma varies as sigma similar to z(4/3). An intermediate diffusive regime may also possible exist for finite values of z under conditions which are explicitly determined in terms of the intrinsic turbulent plasma parameters. The superdiffusie asymptotic result is confirmed numerically via an iterative algorithm. The relevance to previous resuslts is discussed.

Detailed analytical investigation of magnetic field line random walk in turbulent plasmas: II. Isotropic turbulence

The random walk of magnetic field lines in the presence of magnetic turbulence in plasmas is investigated from first principles. An isotropic model is employed for the magnetic turbulence spectrum. An analytical investigation of the asymptotic behavior of the field-line mean-square displacement is carried out. in terms of the position variable z. It is shown that varies as similar to z ln z for large distance z. This result corresponds to a superdiffusive behavior of field line wandering. This investigation complements previous work, which relied on a two-component model for the turbulence spectrum. Contrary to that model, quasilinear theory appears to provide an adequate description of the field line random walk for isotropic turbulence.

Strong-field approximation for Coulomb explosion of H-2(+) and D-2(+) by short intense laser pulses

We present a simple quantum mechanical model to describe Coulomb explosion of H-2(+) and D-2(+) by short, intense infrared laser pulses. The model is based on the length gauge version of the molecular strong-field approximation and is valid when the process of dissociation prior to ionization is negligible. The results are compared with recent experimental data for the proton kinetic energy spectrum [Th. Ergler , Phys. Rev. Lett. 95, 093001 (2005); D. S. Murphy , J. Phys. B 40, S359 (2007)]. Using a Franck-Condon distribution over initial vibrational states, the theory reproduces the overall shape of the spectrum with only a small overestimation of slow protons. The agreement between theory and experiment can be made perfect by using a non-Frank-Condon initial distribution characteristic for H-2(+) (D-2(+)) targets produced by strong-field ionization of H-2 (D-2). For comparison, we also present results obtained by two different tunneling models for this process.

Cross-validation of geo-electrical and hydrogeological models to evaluate seawater intrusion in coastal aquifers

The validation of variable-density flow models simulating seawater intrusion in coastal aquifers requires information about concentration distribution in groundwater. Electrical resistivity tomography (ERT) provides relevant data for this purpose. However, inverse modeling is not accurate because of the non-uniqueness of solutions. Such difficulties in evaluating seawater intrusion can be overcome by coupling geophysical data and groundwater modeling. First, the resistivity distribution obtained by inverse geo-electrical modeling is established. Second, a 3-D variable-density flow hydrogeological model is developed. Third, using Archie's Law, the electrical resistivity model deduced from salt concentration is compared to the formerly interpreted electrical model. Finally, aside from that usual comparison-validation, the theoretical geophysical response of concentrations simulated with the groundwater model can be compared to field-measured resistivity data. This constitutes a cross-validation of both the inverse geo-electrical model and the groundwater model.
[Comte, J.-C., and O. Banton (2007), Cross-validation of geo-electrical and hydrogeological models to evaluate seawater intrusion in coastal aquifers, Geophys. Res. Lett., 34, L10402, doi:10.1029/2007GL029981.]

Out-of-field cell survival following exposure to intensity-modulated radiation fields

PURPOSE:
To determine the in-field and out-of-field cell survival of cells irradiated with either primary field or scattered radiation in the presence and absence of intercellular communication.
METHODS AND MATERIALS:
Cell survival was determined by clonogenic assay in human prostate cancer (DU145) and primary fibroblast (AGO1552) cells following exposure to different field configurations delivered using a 6-MV photon beam produced with a Varian linear accelerator.
RESULTS:
Nonuniform dose distributions were delivered using a multileaf collimator (MLC) in which half of the cell population was shielded. Clonogenic survival in the shielded region was significantly lower than that predicted from the linear quadratic model. In both cell lines, the out-of-field responses appeared to saturate at 40%-50% survival at a scattered dose of 0.70 Gy in DU-145 cells and 0.24 Gy in AGO1522 cells. There was an approximately eightfold difference in the initial slopes of the out-of-field response compared with the a-component of the uniform field response. In contrast, cells in the exposed part of the field showed increased survival. These observations were abrogated by direct physical inhibition of cellular communication and by the addition of the inducible nitric oxide synthase inhibitor aminoguanidine known to inhibit intercellular bystander effects. Additional studies showed the proportion of cells irradiated and dose delivered to the shielded and exposed regions of the field to impact on response.
CONCLUSIONS:
These data demonstrate out-of-field effects as important determinants of cell survival following exposure to modulated irradiation fields with cellular communication between differentially irradiated cell populations playing an important role. Validation of these observations in additional cell models may facilitate the refinement of existing radiobiological models and the observations considered important determinants of cell survival.

A method to reconstruct patient-specific proximal femur surface models from planar pre-operative radiographs

Three-dimensional reconstruction from volumetric medical images (e.g. CT, MRI) is a well-established technology used in patient-specific modelling. However, there are many cases where only 2D (planar) images may be available, e.g. if radiation dose must be limited or if retrospective data is being used from periods when 3D data was not available. This study aims to address such cases by proposing an automated method to create 3D surface models from planar radiographs. The method consists of (i) contour extraction from the radiograph using an Active Contour (Snake) algorithm, (ii) selection of a closest matching 3D model from a library of generic models, and (iii) warping the selected generic model to improve correlation with the extracted contour.

This method proved to be fully automated, rapid and robust on a given set of radiographs. Measured mean surface distance error values were low when comparing models reconstructed from matching pairs of CT scans and planar X-rays (2.57–3.74 mm) and within ranges of similar studies. Benefits of the method are that it requires a single radiographic image to perform the surface reconstruction task and it is fully automated. Mechanical simulations of loaded bone with different levels of reconstruction accuracy showed that an error in predicted strain fields grows proportionally to the error level in geometric precision. In conclusion, models generated by the proposed technique are deemed acceptable to perform realistic patient-specific simulations when 3D data sources are unavailable.

Validation of fixed speed induction generator models for inertial response using wind farm measurements

This paper outlines the use of phasor measurement unit (PMU) records to validate models of fixed speed induction generator (FSIG)-based wind farms during frequency transients. Wind turbine manufacturers usually create their own proprietary models which they can supply to power system utilities for stability studies, subject to confidentiality agreements. However, it is desirable to confirm the accuracy of supplied models with measurements from the particular installation, in order to assess their validity under real field conditions. This is prudent due to possible changes in control algorithms and design retrofits, not accurately reflected or omitted in the supplied model. One important aspect of such models, especially for smaller power systems with limited inertia, is their accuracy during system frequency transients. This paper, therefore, assesses the accuracy of FSIG models with regard to frequency stability, and hence validates a subset of the model dynamics. Such models can then be used with confidence to assess wider system stability implications. The measured and simulated response of a wind farm using doubly fed induction generator (DFIG) technology is also assessed.

Systemic psychotherapy as an intervention for post-traumatic stress responses: an introduction, theoretical rationale and overview of developments in an emerging field of interest

The purpose of this article is to critically examine the literature to provide a rationale for including systemic family therapy (SFT) in the psycho-social treatment of people suffering the impact of post-traumatic stress (PTS). Attention is drawn to the relatively underdeveloped academic literature on PTS and the family. The impact of PTS is conceptualized within a psycho-social framework and the current evidence base for psycho-social interventions for PTS responses is described, highlighting the opportunity and need to undergird this area of daily practice. The impact of PTS on the family at multiple levels is identified, emphasizing its recursive nature. The case for SFT is articulated and a range of models of family intervention for PTS briefly reviewed, concluding with an emphasis on Walsh's key processes in family resilience as a framework for practice.

Performance models of storage contention in cloud environments

We propose simple models to predict the performance degradation of disk requests due to storage device contention in consolidated virtualized environments. Model parameters can be deduced from measurements obtained inside Virtual Machines (VMs) from a system where a single VM accesses a remote storage server. The parameterized model can then be used to predict the effect of storage contention when multiple VMs are consolidated on the same server. We first propose a trace-driven approach that evaluates a queueing network with fair share scheduling using simulation. The model parameters consider Virtual Machine Monitor level disk access optimizations and rely on a calibration technique. We further present a measurement-based approach that allows a distinct characterization of read/write performance attributes. In particular, we define simple linear prediction models for I/O request mean response times, throughputs and read/write mixes, as well as a simulation model for predicting response time distributions. We found our models to be effective in predicting such quantities across a range of synthetic and emulated application workloads. 

A Monte Carlo model of DNA double-strand break clustering and rejoining kinetics for the analysis of pulsed-field gel electrophoresis data

In studies of radiation-induced DNA fragmentation and repair, analytical models may provide rapid and easy-to-use methods to test simple hypotheses regarding the breakage and rejoining mechanisms involved. The random breakage model, according to which lesions are distributed uniformly and independently of each other along the DNA, has been the model most used to describe spatial distribution of radiation-induced DNA damage. Recently several mechanistic approaches have been proposed that model clustered damage to DNA. In general, such approaches focus on the study of initial radiation-induced DNA damage and repair, without considering the effects of additional (unwanted and unavoidable) fragmentation that may take place during the experimental procedures. While most approaches, including measurement of total DNA mass below a specified value, allow for the occurrence of background experimental damage by means of simple subtractive procedures, a more detailed analysis of DNA fragmentation necessitates a more accurate treatment. We have developed a new, relatively simple model of DNA breakage and the resulting rejoining kinetics of broken fragments. Initial radiation-induced DNA damage is simulated using a clustered breakage approach, with three free parameters: the number of independently located clusters, each containing several DNA double-strand breaks (DSBs), the average number of DSBs within a cluster (multiplicity of the cluster), and the maximum allowed radius within which DSBs belonging to the same cluster are distributed. Random breakage is simulated as a special case of the DSB clustering procedure. When the model is applied to the analysis of DNA fragmentation as measured with pulsed-field gel electrophoresis (PFGE), the hypothesis that DSBs in proximity rejoin at a different rate from that of sparse isolated breaks can be tested, since the kinetics of rejoining of fragments of varying size may be followed by means of computer simulations. The problem of how to account for background damage from experimental handling is also carefully considered. We have shown that the conventional procedure of subtracting the background damage from the experimental data may lead to erroneous conclusions during the analysis of both initial fragmentation and DSB rejoining. Despite its relative simplicity, the method presented allows both the quantitative and qualitative description of radiation-induced DNA fragmentation and subsequent rejoining of double-stranded DNA fragments. (C) 2004 by Radiation Research Society.

A New Gradient Descent Approach for Local Learning of Fuzzy Neural Models

The majority of reported learning methods for Takagi-Sugeno-Kang fuzzy neural models to date mainly focus on the improvement of their accuracy. However, one of the key design requirements in building an interpretable fuzzy model is that each obtained rule consequent must match well with the system local behaviour when all the rules are aggregated to produce the overall system output. This is one of the distinctive characteristics from black-box models such as neural networks. Therefore, how to find a desirable set of fuzzy partitions and, hence, to identify the corresponding consequent models which can be directly explained in terms of system behaviour presents a critical step in fuzzy neural modelling. In this paper, a new learning approach considering both nonlinear parameters in the rule premises and linear parameters in the rule consequents is proposed. Unlike the conventional two-stage optimization procedure widely practised in the field where the two sets of parameters are optimized separately, the consequent parameters are transformed into a dependent set on the premise parameters, thereby enabling the introduction of a new integrated gradient descent learning approach. A new Jacobian matrix is thus proposed and efficiently computed to achieve a more accurate approximation of the cost function by using the second-order Levenberg-Marquardt optimization method. Several other interpretability issues about the fuzzy neural model are also discussed and integrated into this new learning approach. Numerical examples are presented to illustrate the resultant structure of the fuzzy neural models and the effectiveness of the proposed new algorithm, and compared with the results from some well-known methods.

Selectivity fields: comparative molecular field analysis (CoMFA) of the glycine/NMDA and AMPA receptors

An approach for evaluation of binding selectivity was suggested and exemplified using glycine/NMDA and AMPA receptors. For analyzing the pairwise selectivity, we propose to use the difference between biological activities (expressed as -log Ki) of ligands with respect to different receptor subtypes as a dependent variable for building comparative molecular field analysis (CoMFA) models. The resulting fields (which will be referred to as the "selectivity fields") indicate the ways of increasing selectivity of binding, inhibition, etc. As an example, CoMFA of a set of pyrazolo[1,5-c]quinazolines and triazolo[1,5-c]quinazolines was used for considering the binding selectivity with respect to glycine/NMDA and AMPA receptors. In addition, the mapping of these fields onto the molecular models of the corresponding receptors makes it possible to reveal the reasons for experimentally observed selectivity as well as to suggest additional ways of increasing selectivity.