Analysis of population pharmacokinetic data using NONMEM and WinBUGS

The aim of this report is to describe the use of WinBUGS for two datasets that arise from typical population pharmacokinetic studies. The first dataset relates to gentamicin concentration-time data that arose as part of routine clinical care of 55 neonates. The second dataset incorporated data from 96 patients receiving enoxaparin. Both datasets were originally analyzed by using NONMEM. In the first instance, although NONMEM provided reasonable estimates of the fixed effects parameters it was unable to provide satisfactory estimates of the between-subject variance. In the second instance, the use of NONMEM resulted in the development of a successful model, albeit with limited available information on the between-subject variability of the pharmacokinetic parameters. WinBUGS was used to develop a model for both of these datasets. Model comparison for the enoxaparin dataset was performed by using the posterior distribution of the log-likelihood and a posterior predictive check. The use of WinBUGS supported the same structural models tried in NONMEM. For the gentamicin dataset a one-compartment model with intravenous infusion was developed, and the population parameters including the full between-subject variance-covariance matrix were available. Analysis of the enoxaparin dataset supported a two compartment model as superior to the one-compartment model, based on the posterior predictive check. Again, the full between-subject variance-covariance matrix parameters were available. Fully Bayesian approaches using MCMC methods, via WinBUGS, can offer added value for analysis of population pharmacokinetic data.

A Mixture model with random-effects components for clustering correlated gene-expression profiles

Motivation: The clustering of gene profiles across some experimental conditions of interest contributes significantly to the elucidation of unknown gene function, the validation of gene discoveries and the interpretation of biological processes. However, this clustering problem is not straightforward as the profiles of the genes are not all independently distributed and the expression levels may have been obtained from an experimental design involving replicated arrays. Ignoring the dependence between the gene profiles and the structure of the replicated data can result in important sources of variability in the experiments being overlooked in the analysis, with the consequent possibility of misleading inferences being made. We propose a random-effects model that provides a unified approach to the clustering of genes with correlated expression levels measured in a wide variety of experimental situations. Our model is an extension of the normal mixture model to account for the correlations between the gene profiles and to enable covariate information to be incorporated into the clustering process. Hence the model is applicable to longitudinal studies with or without replication, for example, time-course experiments by using time as a covariate, and to cross-sectional experiments by using categorical covariates to represent the different experimental classes. Results: We show that our random-effects model can be fitted by maximum likelihood via the EM algorithm for which the E(expectation) and M(maximization) steps can be implemented in closed form. Hence our model can be fitted deterministically without the need for time-consuming Monte Carlo approximations. The effectiveness of our model-based procedure for the clustering of correlated gene profiles is demonstrated on three real datasets, representing typical microarray experimental designs, covering time-course, repeated-measurement and cross-sectional data. In these examples, relevant clusters of the genes are obtained, which are supported by existing gene-function annotation. A synthetic dataset is considered too.

Ecological-economic modeling for biodiversity management: Potential, pitfalls, and prospects

Ecologists and economists both use models to help develop strategies for biodiversity management. The practical use of disciplinary models, however, can be limited because ecological models tend not to address the socioeconomic dimension of biodiversity management, whereas economic models tend to neglect the ecological dimension. Given these shortcomings of disciplinary models, there is a necessity to combine ecological and economic knowledge into ecological-economic models. It is insufficient if scientists work separately in their own disciplines and combine their knowledge only when it comes to formulating management recommendations. Such an approach does not capture feedback loops between the ecological and the socioeconomic systems. Furthermore, each discipline poses the management problem in its own way and comes up with its own most appropriate solution. These disciplinary solutions, however are likely to be so different that a combined solution considering aspects of both disciplines cannot be found. Preconditions for a successful model-based integration of ecology and economics include (1) an in-depth knowledge of the two disciplines, (2) the adequate identification and framing of the problem to be investigated, and (3) a common understanding between economists and ecologists of modeling and scale. To further advance ecological-economic modeling the development of common benchmarks, quality controls, and refereeing standards for ecological-economic models is desirable.

Leading coordinate analysis of reaction pathways in proton chain transfer: Application to a two-proton transfer model for the green fluorescent protein

The ‘leading coordinate’ approach to computing an approximate reaction pathway, with subsequent determination of the true minimum energy profile, is applied to a two-proton chain transfer model based on the chromophore and its surrounding moieties within the green fluorescent protein (GFP). Using an ab initio quantum chemical method, a number of different relaxed energy profiles are found for several plausible guesses at leading coordinates. The results obtained for different trial leading coordinates are rationalized through the calculation of a two-dimensional relaxed potential energy surface (PES) for the system. Analysis of the 2-D relaxed PES reveals that two of the trial pathways are entirely spurious, while two others contain useful information and can be used to furnish starting points for successful saddle-point searches. Implications for selection of trial leading coordinates in this class of proton chain transfer reactions are discussed, and a simple diagnostic function is proposed for revealing whether or not a relaxed pathway based on a trial leading coordinate is likely to furnish useful information.

Forest logging and institutional thresholds in developing south-east Asian economies: A conceptual model

Many developing south-east Asian governments are not capturing full rent from domestic forest logging operations. Such rent losses are commonly related to institutional failures, where informal institutions tend to dominate the control of forestry activity in spite of weakly enforced regulations. Our model is an attempt to add a new dimension to thinking about deforestation. We present a simple conceptual model, based on individual decisions rather than social or forest planning, which includes the human dynamics of participation in informal activity and the relatively slower ecological dynamics of changes in forest resources. We demonstrate how incumbent informal logging operations can be persistent, and that any spending aimed at replacing the informal institutions can only be successful if it pushes institutional settings past some threshold. (C) 2006 Elsevier B.V. All rights reserved.

Growth of Bose-Einstein condensates from thermal vapor

We report on a quantitative study of the growth process of 87Rb Bose-Einstein condensates. By continuous evaporative cooling we directly control the thermal cloud from which the condensate grows. We compare the experimental data with the results of a theoretical model based on quantum kinetic theory. We find quantitative agreement with theory for the situation of strong cooling, whereas in the weak cooling regime a distinctly different behavior is found in the experiment.

In vivo measurement of the effect of intra-abdominal pressure on the human spine

In humans, intra-abdominal pressure (IAP) is elevated during many everyday activities. This experiment aimed to investigate the extent to which increased IAP-without concurrent activity of the abdominal or back extensor muscles-produces an extensor torque. With subjects positioned in side lying on a swivel table with its axis at L3, moments about this vertebral level were measured when IAP was transiently increased by electrical stimulation of the diaphragm via the phrenic nerve. There was no electromyographic activity in abdominal and back extensor muscles. When IAP was increased artificially to similar to 15% of the maximum IAP amplitude that could be generated voluntarily with the trunk positioned in flexion, a trunk extensor moment (similar to6 Nm) was recorded. The size of the effect was proportional to the increase in pressure. The extensor moment was consistent with that predicted from a model based on measurements of abdominal cross-sectional area and IAP moment arm. When IAP was momentarily increased while the trunk was flexed passively at a constant velocity, the external torque required to maintain the velocity was increased. These results provide the first in vivo data of the amplitude of extensor moment that is produced by increased IAP. Although the net effect of this extensor torque in functional tasks would be dependent on the muscles used to increase the IAP and their associated flexion torque, the data do provide evidence that IAP contributes, at least in part, to spinal stability. (C) 2001 Elsevier Science Ltd. All rights reserved.

A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.

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.

Value congruence in leader-member exchange

In this study, hypotheses were tested that the quality of leader-member exchanges (LMX) depends on congruity of values between leader and member. Data on negotiating latitude and personal values were gathered from 160 members of 30 work groups in Australian organizations. Factor analysis revealed 5 value dimensions: Freedom, Achievement, Mateship, Obedience, and Coping. Analyses of variance supported the hypothesis that LMX quality is higher when leaders and members share achievement and obedience values. Subsequent exploratory analysis, however, indicated that a more complex model based on compatibility of leader authority and member affiliation values may provide a more complete representation.

Protecting Schrodinger cat states using feedback

A feedback model based on direct photodetection and micromaser-like atomic injection is proposed for the preservation of quantum coherence in a cavity. We show that in this way it is possible to slow down significantly the decoherence of Schrodinger cat states.

Features of colorectal cancers with high-level microsatellite instability occurring in familial and sporadic settings

High-level microsatellite instability (AISI-H) is demonstrated in 10 to 15% of sporadic colorectal cancers and in most cancers presenting In the inherited condition hereditary nonpolyposis colorectal cancer (HNPCC). Distinction between these categories of MSI-H cancer is of clinical importance and the aim of this study was to assess clinical, pathological, and molecular features that might he discriminatory. One hundred and twelve MSI-H colorectal cancers from families fulfilling the Bethesda criteria were compared with 57 sporadic MSI-H colorectal cancers. HNPCC cancers presented at a lower age (P < 0.001) with no sporadic MSI-H cancer being diagnosed before the age of 57 years. MSI was less extensive in HNPCC cancers with 72% microsatellite markers showing band shifts compared with 87% in sporadic tumors (P < 0.001). Absent immunostaining for hMSH2 was only found in HNPCC tumors. Methylation of bMLH1 was observed in 87% of sporadic cancers but also in 55% of HNPCC tumors that showed loss of expression of hMLH1 (P = 0.02). HNPCC cancers were more frequently characterized by aberrant beta -catenin immunostaining as evidenced by nuclear positivity (P < 0.001). Aberrant p53 immunostaining was infrequent in both groups. There were no differences with respect to 5q loss of heterozygosity or codon 12 K-ras mutation, which were infrequent in both groups. Sporadic MSI-H cancers were more frequently heterogeneous (P < 0.001), poorly differentiated (P = 0.02), mucinous (P = 0.02), and proximally located (P = 0.04) than RNPCC tumors. In sporadic MSI-H cancers, contiguous adenomas were likely to be serrated whereas traditional adenomas were dominant in HNPCC. Lymphocytic infiltration was more pronounced in HNPCC but the results did not reach statistical significance. Overall, HNPCC cancers were more like common colorectal cancer in terms of morphology and expression of beta -catenin whereas sporadic MSI-H cancers displayed features consistent with a different morphogenesis. No individual feature was discriminatory for all RN-PCC cancers. However, a model based on four features was able to classify 94.5% of tumors as sporadic or HNPCC. The finding of multiple differences between sporadic and familial MSI-H colorectal cancer with respect to both genotype and phenotype is consistent with tumorigenesis through parallel evolutionary pathways and emphasizes the importance of studying the two groups separately.

Role of circumventricular organs in pro-inflammatory cytokine-induced activation of the hypothalamic-pituitary-adrenal axis

1. The past 15 years has seen the emergence of a new field of neuroscience research based primarily on how the immune system and the central nervous system can interact. A notable example of this interaction occurs when peripheral inflammation, infection or tissue injury activates the hypothalamic- pituitary-adrenal axis (HPA). 2. During such assaults, immune cells release the pro- inflammatory cytokines interleukin (IL)-1, IL-6 and tumour necrosis factor-alpha into the general circulation. 3. These cytokines are believed to act as mediators for HPA axis activation. However, physical limitations of cytokines impede their movement across the blood-brain barrier and, consequently, it has been unclear as to precisely how and where IL-1beta signals cross into the brain to trigger HPA axis activation. 4. Evidence from recent anatomical and functional studies suggests two neuronal networks may be involved in triggering HPA axis activity in response to circulating cytokines. These are catecholamine cells of the medulla oblongata and the circumventricular organs (CVO). 5. The present paper examines the role of CVO in generating HPA axis responses to pro-inflammatory cytokines and culminates with a proposed model based on cytokine signalling primarily involving the area postrema and catecholamine cells in the ventrolateral and dorsal medulla.

MapScript: A map algebra programming language incorporating neighborhood analysis

Map algebra is a data model and simple functional notation to study the distribution and patterns of spatial phenomena. It uses a uniform representation of space as discrete grids, which are organized into layers. This paper discusses extensions to map algebra to handle neighborhood operations with a new data type called a template. Templates provide general windowing operations on grids to enable spatial models for cellular automata, mathematical morphology, and local spatial statistics. A programming language for map algebra that incorporates templates and special processing constructs is described. The programming language is called MapScript. Example program scripts are presented to perform diverse and interesting neighborhood analysis for descriptive, model-based and processed-based analysis.

Probabilistic modelling of CO2 corrosion laboratory data using neural networks

The present paper addresses two major concerns that were identified when developing neural network based prediction models and which can limit their wider applicability in the industry. The first problem is that it appears neural network models are not readily available to a corrosion engineer. Therefore the first part of this paper describes a neural network model of CO2 corrosion which was created using a standard commercial software package and simple modelling strategies. It was found that such a model was able to capture practically all of the trends noticed in the experimental data with acceptable accuracy. This exercise has proven that a corrosion engineer could readily develop a neural network model such as the one described below for any problem at hand, given that sufficient experimental data exist. This applies even in the cases when the understanding of the underlying processes is poor. The second problem arises from cases when all the required inputs for a model are not known or can be estimated with a limited degree of accuracy. It seems advantageous to have models that can take as input a range rather than a single value. One such model, based on the so-called Monte Carlo approach, is presented. A number of comparisons are shown which have illustrated how a corrosion engineer might use this approach to rapidly test the sensitivity of a model to the uncertainities associated with the input parameters. (C) 2001 Elsevier Science Ltd. All rights reserved.