Stiffly accurate Runge-Kutta methods for stiff stochastic differential equations

In this paper we discuss implicit methods based on stiffly accurate Runge-Kutta methods and splitting techniques for solving Stratonovich stochastic differential equations (SDEs). Two splitting techniques: the balanced splitting technique and the deterministic splitting technique, are used in this paper. We construct a two-stage implicit Runge-Kutta method with strong order 1.0 which is corrected twice and no update is needed. The stability properties and numerical results show that this approach is suitable for solving stiff SDEs. (C) 2001 Elsevier Science B.V. All rights reserved.

Proximal arterial vasoconstriction precedes regression of the hyaloid vasculature

Purpose: To determine whether constriction of proximal arterial vessels precedes involution of the distal hyaloid vasculature in the mouse, under normal conditions, and whether this vasoconstriction is less pronounced when the distal hyaloid network persists, as it does in oxygen-induced retinopathy (OIR). Methods: Photomicrographs of the vasa hyaloidea propria were analysed from pre-term pups (1-2 days prior to birth), and on Days 1-11 post-birth. The OIR model involved exposing pups to similar to 90% O-2 from D1-5, followed by return to ambient air. At sampling times pups were anaesthetised and perfused with india ink. Retinal flatmounts were also incubated with FITC-lectin (BS-1, G. simplicifolia,); this labels all vessels, allowing identification of vessels not patent to the perfusate. Results: Mean diameter of proximal hyaloid vessels in preterm pups was 25.44 +/- 1.98 mum; +/-1 SEM). Within 3-12 hrs of birth, significant vasoconstriction was evident (diameter:12.45 +/- 0.88 mum), and normal hyaloid regression subsequently occurred. Similar vasoconstriction occurred in the O-2-treated group, but this was reversed upon return to room air, with significant dilation of proximal vessels by D7 (diameter: 31.75 +/- 11.99 mum) and distal hyaloid vessels subsequently became enlarged and tortuous. Conclusions: Under normal conditions, vasoconstriction of proximal hyaloid vessels occurs at birth, preceding attenuation of distal hyaloid vessels. Vasoconstriction also occurs in O-2-treated pups during treatment, but upon return to room air, the remaining hyaloid vessels dilate proximally, and the distal vessels become dilated and tortuous. These observations support the contention that regression of the hyaloid network is dependent, in the first instance, on proximal arterial vasoconstriction.

An integrative approach for studying the etiology of alcoholism and other addictions

Studies of alcoholism etiology often focus on genetic or psy-chosocial approaches, but not both. Greater understanding of the etiology of alcohol, tobacco and other addictions will come from integration of these research traditions. A research approach is outlined to test three models for the etiology of addictions — behavioral undercontrol, pharmacologic vulnerability, negative affect regulation — addressing key questions including (i) mediators of genetic effects, (ii) genotype-environment correlation effects, (iii) genotype x environment interaction effects, (iv) the developmental unfolding of genetic and environmental effects, (v) subtyping including identification of distinct trajectories of substance involvement, (vi) identification of individual genes that contribute to risk, and (vii) the consequences of excessive use. By using coordinated research designs, including prospective assessment of adolescent twins and their siblings and parents; of adult substance dependent and control twins and their MZ and DZ cotwins, the spouses of these pairs, and their adolescent offspring; and of regular families; by selecting for gene-mapping approaches sibships screened for extreme concordance or discordance on quantitative indices of substance use; and by using experimental (drug challenge) as well as survey approaches, a number of key questions concerning addiction etiology can be addressed. We discuss complementary strengths and weaknesses of different sampling strategies, as well as methods to implement such an integrated approach illustrated for the study of alcoholism etiology. A coordinated program of twin and family studies will allow a comprehensive dissection of the interplay of genetic and environmental risk-factors in the etiology of alcoholism and other addictions.

A retrospective analysis of mycophenolic acid and cyclosporin concentrations with acute rejection in renal transplant recipients

Objectives: Although monitoring of cyclosporin (CsA) is standard clinical practice postrenal transplantation. mycophenolic acid (MPA) concentrations are not routinely measured. There is evidence that a relationship exists between MPA area under the concentration-time curve (AUC) and rejection. In this study, a retrospective analysis was undertaken of 27 adult renal transplant recipients. Methods: Patients received CsA and MPA therapy and had a four-point MPA AUC investigation. The relationship between MPA AUC performed in the first week after transplantation, as well as median trough cyclosporin concentrations, and clinical outcomes in the first month posttransplant were evaluated. Results: A total of 12 patients experienced biopsy proven rejection (44.4%) and 4 patients had gastrointestinal adverse events (14.8%). A statistically significant relationship was observed between the incidence of biopsy proven rejection and both MPA AUC (p = 0.02) and median trough CsA concentration (p = 0.008). No relationship between trough MPA concentration and rejection was observed (p = 0.21). Only 3 of 11 (27%) patients with an MPA AUC > 30 mg.h/L and a median trough CsA > 175 mug/L experienced acute rejection, compared with a 56% incidence of rejection for the remaining 16 patients. Patients who experienced adverse gastrointestinal events had significantly lower MPA AUC (p = 0.04), but median trough CsA concentrations were not significantly different (p = 0.24). Further, 3 of these 4 patients had rejection episodes. Conclusions: in addition to standard CsA monitoring, we propose further investigation of the use of a 4-point sampling strategy to predict MPA AUC in the first week posttransplant, which may facilitate optimization of mycophenolate mofetil dose at a rime when patients are most vulnerable to acute rejection. (C) 2001 The Canadian Society of Clinical Chemists. All rights reserved.

Corneal endothelial response to induced myopia in the chicken

Purpose: To investigate the role of corneal endothelial surface enlargement in the chicken myopia model in inducing corneal endothelial changes. Methods: Lid suture was performed on one eye of 1-day-old cockerels. Five chickens were killed at 1 week, and four chickens killed at each of 3 weeks, 6 weeks, and 10 weeks postnatal. The endothelial morphology was obtained by flat mounting the endothelial surface and the subsequent digitisation. Comparisons were undertaken between the control unsutured eye and the lid-sutured eye endothelium, and between the central endothelial areas compared to the peripheral endothelial areas in both the myopic and the normal corneas. Calculation of the contribution to the endothelial change by hypertrophy and mitosis were calculated using Bahn's formula. Results: Total endothelial surface area increased significantly over time in the myopic model compared to control eyes but the mean cell area of endothelial cells remained the same for both the enlarged myopic endothelial surface area and in the normal controls. Sampling from the central and the peripheral corneal endothelial surface also disclosed no difference. The mean cell area did increase steadily with age but was the same for both normal and myopic corneas. Conclusions: It would appear that there are equal contributions from hypertrophy and mitosis in the myopic group and the normal corneal group with a slightly increasing trend towards mitotic activity in the myopic corneal endothelial layer.

Time evolution in the unimolecular master equation at low temperatures: full spectral solution with scalable iterative methods and high precision

A new method is presented to determine an accurate eigendecomposition of difficult low temperature unimolecular master equation problems. Based on a generalisation of the Nesbet method, the new method is capable of achieving complete spectral resolution of the master equation matrix with relative accuracy in the eigenvectors. The method is applied to a test case of the decomposition of ethane at 300 K from a microcanonical initial population with energy transfer modelled by both Ergodic Collision Theory and the exponential-down model. The fact that quadruple precision (16-byte) arithmetic is required irrespective of the eigensolution method used is demonstrated. (C) 2001 Elsevier Science B.V. All rights reserved.

Design of clinical pharmacology trials

1. There are a variety of methods that could be used to increase the efficiency of the design of experiments. However, it is only recently that such methods have been considered in the design of clinical pharmacology trials. 2. Two such methods, termed data-dependent (e.g. simulation) and data-independent (e.g. analytical evaluation of the information in a particular design), are becoming increasingly used as efficient methods for designing clinical trials. These two design methods have tended to be viewed as competitive, although a complementary role in design is proposed here. 3. The impetus for the use of these two methods has been the need for a more fully integrated approach to the drug development process that specifically allows for sequential development (i.e. where the results of early phase studies influence later-phase studies). 4. The present article briefly presents the background and theory that underpins both the data-dependent and -independent methods with the use of illustrative examples from the literature. In addition, the potential advantages and disadvantages of each method are discussed.

Application of Petrov-Galerkin methods to transient boundary value problems in chemical engineering: adsorption with steep gradients in bidisperse solids

Petrov-Galerkin methods are known to be versatile techniques for the solution of a wide variety of convection-dispersion transport problems, including those involving steep gradients. but have hitherto received little attention by chemical engineers. We illustrate the technique by means of the well-known problem of simultaneous diffusion and adsorption in a spherical sorbent pellet comprised of spherical, non-overlapping microparticles of uniform size and investigate the uptake dynamics. Solutions to adsorption problems exhibit steep gradients when macropore diffusion controls or micropore diffusion controls, and the application of classical numerical methods to such problems can present difficulties. In this paper, a semi-discrete Petrov-Galerkin finite element method for numerically solving adsorption problems with steep gradients in bidisperse solids is presented. The numerical solution was found to match the analytical solution when the adsorption isotherm is linear and the diffusivities are constant. Computed results for the Langmuir isotherm and non-constant diffusivity in microparticle are numerically evaluated for comparison with results of a fitted-mesh collocation method, which was proposed by Liu and Bhatia (Comput. Chem. Engng. 23 (1999) 933-943). The new method is simple, highly efficient, and well-suited to a variety of adsorption and desorption problems involving steep gradients. (C) 2001 Elsevier Science Ltd. All rights reserved.

On surrogate methods for detecting lateral gene transfer

Surrogate methods for detecting lateral gene transfer are those that do not require inference of phylogenetic trees. Herein I apply four such methods to identify open reading frames (ORFs) in the genome of Escherichia coli K12 that may have arisen by lateral gene transfer. Only two of these methods detect the same ORFs more frequently than expected by chance, whereas several intersections contain many fewer ORFs than expected. Each of the four methods detects a different non-random set of ORFs. The methods may detect lateral ORFs of different relative ages; testing this hypothesis will require rigorous inference of trees. (C) 2001 Federation of European Microbiological Societies. Published by Elsevier Science BN. All rights reserved.

Stickiness in foods: A review of mechanisms and test methods

Problems associated with the stickiness of food in processing and storage practices along with its causative factors are outlined. Fundamental mechanisms that explain why and how food products become sticky are discussed. Methods currently in use for characterizing and overcoming stickiness problems in food processing and storage operations are described. The use of glass transition temperature-based model, which provides a rational basis for understanding and characterizing the stickiness of many food products, is highlighted.

'Neighbourhood' size, dispersal and density estimates in the prickly forest skink (Gnypetoscincus queenslandiae) using individual genetic and demographic methods

Dispersal, or the amount of dispersion between an individual's birthplace and that of its offspring, is of great importance in population biology, behavioural ecology and conservation, however, obtaining direct estimates from field data on natural populations can be problematic. The prickly forest skink, Gnypetoscincus queenslandiae, is a rainforest endemic skink from the wet tropics of Australia. Because of its log-dwelling habits and lack of definite nesting sites, a demographic estimate of dispersal distance is difficult to obtain. Neighbourhood size, defined as 4 piD sigma (2) (where D is the population density and sigma (2) the mean axial squared parent-offspring dispersal rate), dispersal and density were estimated directly and indirectly for this species using mark-recapture and microsatellite data, respectively, on lizards captured at a local geographical scale of 3 ha. Mark-recapture data gave a dispersal rate of 843 m(2)/generation (assuming a generation time of 6.5 years), a time-scaled density of 13 635 individuals * generation/km(2) and, hence, a neighbourhood size of 144 individuals. A genetic method based on the multilocus (10 loci) microsatellite genotypes of individuals and their geographical location indicated that there is a significant isolation by distance pattern, and gave a neighbourhood size of 69 individuals, with a 95% confidence interval between 48 and 184. This translates into a dispersal rate of 404 m(2)/generation when using the mark-recapture density estimation, or an estimate of time-scaled population density of 6520 individuals * generation/km(2) when using the mark-recapture dispersal rate estimate. The relationship between the two categories of neighbourhood size, dispersal and density estimates and reasons for any disparities are discussed.

Surface diffusion of hydrocarbons in activated carbon: Comparison between constant molar flow, differential permeation and differential adsorption bed methods

This paper presents the comparison of surface diffusivities of hydrocarbons in activated carbon. The surface diffusivities are obtained from the analysis of kinetic data collected using three different kinetics methods- the constant molar flow, the differential adsorption bed and the differential permeation methods. In general the values of surface diffusivity obtained by these methods agree with each other, and it is found that the surface diffusivity increases very fast with loading. Such a fast increase can not be accounted for by a thermodynamic Darken factor, and the surface heterogeneity only partially accounts for the fast rise of surface diffusivity versus loading. Surface diffusivities of methane, ethane, propane, n-butane, n-hexane, benzene and ethanol on activated carbon are reported in this paper.