Searching for convergence in phylogenetic Markov chain Monte Carlo

Markov chain Monte Carlo (MCMC) is a methodology that is gaining widespread use in the phylogenetics community and is central to phylogenetic software packages such as MrBayes. An important issue for users of MCMC methods is how to select appropriate values for adjustable parameters such as the length of the Markov chain or chains, the sampling density, the proposal mechanism, and, if Metropolis-coupled MCMC is being used, the number of heated chains and their temperatures. Although some parameter settings have been examined in detail in the literature, others are frequently chosen with more regard to computational time or personal experience with other data sets. Such choices may lead to inadequate sampling of tree space or an inefficient use of computational resources. We performed a detailed study of convergence and mixing for 70 randomly selected, putatively orthologous protein sets with different sizes and taxonomic compositions. Replicated runs from multiple random starting points permit a more rigorous assessment of convergence, and we developed two novel statistics, delta and epsilon, for this purpose. Although likelihood values invariably stabilized quickly, adequate sampling of the posterior distribution of tree topologies took considerably longer. Our results suggest that multimodality is common for data sets with 30 or more taxa and that this results in slow convergence and mixing. However, we also found that the pragmatic approach of combining data from several short, replicated runs into a metachain to estimate bipartition posterior probabilities provided good approximations, and that such estimates were no worse in approximating a reference posterior distribution than those obtained using a single long run of the same length as the metachain. Precision appears to be best when heated Markov chains have low temperatures, whereas chains with high temperatures appear to sample trees with high posterior probabilities only rarely. [Bayesian phylogenetic inference; heating parameter; Markov chain Monte Carlo; replicated chains.]

Determination of pore size distributions by regularization and finite element collocation

The problem of extracting pore size distributions from characterization data is solved here with particular reference to adsorption. The technique developed is based on a finite element collocation discretization of the adsorption integral, with fitting of the isotherm data by least squares using regularization. A rapid and simple technique for ensuring non-negativity of the solutions is also developed which modifies the original solution having some negativity. The technique yields stable and converged solutions, and is implemented in a package RIDFEC. The package is demonstrated to be robust, yielding results which are less sensitive to experimental error than conventional methods, with fitting errors matching the known data error. It is shown that the choice of relative or absolute error norm in the least-squares analysis is best based on the kind of error in the data. (C) 1998 Elsevier Science Ltd. All rights reserved.

X-ray computed tomography to quantify tree rooting spatial distributions

Poor root development due to constraining soil conditions could be an important factor influencing health of urban trees. Therefore, there is a need for efficient techniques to analyze the spatial distribution of tree roots. An analytical procedure for describing tree rooting patterns from X-ray computed tomography (CT) data is described and illustrated. Large irregularly shaped specimens of undisturbed sandy soil were sampled from Various positions around the base of trees using field impregnation with epoxy resin, to stabilize the cohesionless soil. Cores approximately 200 mm in diameter by 500 mm in height were extracted from these specimens. These large core samples were scanned with a medical X-ray CT device, and contiguous images of soil slices (2 mm thick) were thus produced. X-ray CT images are regarded as regularly-spaced sections through the soil although they are not actual 2D sections but matrices of voxels similar to 0.5 mm x 0.5 mm x 2 mm. The images were used to generate the equivalent of horizontal root contact maps from which three-dimensional objects, assumed to be roots, were reconstructed. The resulting connected objects were used to derive indices of the spatial organization of roots, namely: root length distribution, root length density, root growth angle distribution, root spatial distribution, and branching intensity. The successive steps of the method, from sampling to generation of indices of tree root organization, are illustrated through a case study examining rooting patterns of valuable urban trees. (C) 1999 Elsevier Science B.V. All rights reserved.

Habitat structure and egg distributions in the processionary caterpillar Ochrogaster lunifer: lessons for conservation and pest management

1. The spatial and temporal distribution of eggs laid by herbivorous insects is a crucial component of herbivore population stability, as it influences overall mortality within the population. Thus an ecologist studying populations of an endangered butterfly can do little to increase its numbers through habitat management without knowledge of its egg-laying patterns across individual host-plants under different habitat management regimes. At the other end of the spectrum, a knowledge of egg-laying behaviour can do much to control pest outbreaks by disrupting egg distributions that lead to rapid population growth. 2. The distribution of egg batches of the processionary caterpillar Ochrogaster lunifer on acacia trees was monitored in 21 habitats during 2 years in coastal Australia. The presence of egg batches on acacias was affected by host-tree 'quality' (tree size and foliar chemistry that led to increased caterpillar survival) and host-tree 'apparency' (the amount of vegetation surrounding host-trees). 3. In open homogeneous habitats, more egg batches were laid on high-quality trees, increasing potential population growth. In diverse mixed-species habitats, more egg batches were laid on low-quality highly apparent trees, reducing population growth and so reducing the potential for unstable population dynamics. The aggregation of batches on small apparent trees in diverse habitats led to outbreaks on these trees year after year, even when population levels were low, while site-wide outbreaks were rare. 4. These results predict that diverse habitats with mixed plant species should increase insect aggregation and increase population stability. In contrast, in open disturbed habitats or in regular plantations, where egg batches are more evenly distributed across high-quality hosts, populations should be more unstable, with site-wide outbreaks and extinctions being more common. 5. Mixed planting should be used on habitat regeneration sites to increase the population stability of immigrating or reintroduced insect species. Mixed planting also increases the diversity of resources, leading to higher herbivore species richness. With regard to the conservation of single species, different practices of habitat management will need to be employed depending on whether a project is concerned with methods of rapidly increasing the abundance of an endangered insect or concerned with the maintenance of a stable, established insect population that is perhaps endemic to an area. Suggestions for habitat management in these different cases are discussed. 6. Finally, intercropping can be highly effective in reducing pest outbreaks, although the economic gains of reduced pest attack may be outweighed by reduced crop yields in mixed-crop systems.

Further results on the relationship between mu-invariant measures and quasi-stationary distributions for absorbing continuous-time Markov chains

This note considers continuous-time Markov chains whose state space consists of an irreducible class, C, and an absorbing state which is accessible from C. The purpose is to provide results on mu-invariant and mu-subinvariant measures where absorption occurs with probability less than one. In particular, the well-known premise that the mu-invariant measure, m, for the transition rates be finite is replaced by the more natural premise that m be finite with respect to the absorption probabilities. The relationship between mu-invariant measures and quasi-stationary distributions is discussed. (C) 2000 Elsevier Science Ltd. All rights reserved.

New methods for determining quasi-stationary distributions for Markov chains

We shall be concerned with the problem of determining quasi-stationary distributions for Markovian models directly from their transition rates Q. We shall present simple conditions for a mu-invariant measure m for Q to be mu-invariant for the transition function, so that if m is finite, it can be normalized to produce a quasi-stationary distribution. (C) 2000 Elsevier Science Ltd. All rights reserved.

Characterization of pore size distributions of mesoporous materials from adsorption isotherms

In this article, a new hybrid model for estimating the pore size distribution of micro- and mesoporous materials is developed, and tested with the adsorption data of nitrogen, oxygen, and argon on ordered mesoporous materials reported in the literature. For the micropore region, the model uses the Dubinin-Rudushkevich (DR) isotherm with the Chen-Yang modification. A recent isotherm model of the authors for nonporous materials, which uses a continuum-mechanical model for the multilayer region and the Unilan model for the submonolayer region, has been extended for adsorption in mesopores. The experimental data is inverted using regularization to obtain the pore size distribution. The present model was found to be successful in predicting the pore size distribution of pure as well as binary physical mixtures of MCM-41 synthesized with different templates, with results in agreement with those from the XRD method and nonlocal density functional theory. It was found that various other recent methods, as well as the classical Broekhoff and de Beer method, underpredict the pore diameter of MCM-41. The present model has been successfully applied to MCM-48, SBA's, CMK, KIT, HMS, FSM, MTS, mesoporous fly ash, and a large number of other regular mesoporous materials.

Determination of activation energy distributions for chemisorption of oxygen on carbon: an improved approach

The present paper proposes an approach to obtaining the activation energy distribution for chemisorption of oxygen onto carbon surfaces, while simultaneously allowing for the activation energy dependence of the pre-exponential factor of the rate constant. Prior studies in this area have considered this factor to be uniform, thereby biasing estimated distributions. The results show that the derived activation energy distribution is not sensitive to the chemisorption mechanism because of the step function like property of the coverage. The activation energy distribution is essentially uniform for some carbons, and has two or possibly more discrete stages, suggestive of at least two types of sites, each with its own uniform distribution. The pre-exponential factors of the reactions are determined directly from the experimental data, and are found not to be constant as assumed in earlier work, but correlated with the activation energy. The latter results empirically follow an exponential function, supporting some earlier statistical and experimental work. The activation energy distribution obtained in the present paper permits improved correlation of chemisorption data in comparison to earlier studies. (C) 2000 Elsevier Science Ltd. All rights reserved.

Residence time distributions of gas flowing through rotating drum bioreactors

Residence time distribution studies of gas through a rotating drum bioreactor for solid-state fermentation were performed using carbon monoxide as a tracer gas. The exit concentration as a function of time differed considerably from profiles expected for plug flow, plug flow with axial dispersion, and continuous stirred tank reactor (CSTR) models. The data were then fitted by least-squares analysis to mathematical models describing a central plug flow region surrounded by either one dead region (a three-parameter model) or two dead regions (a five-parameter model). Model parameters were the dispersion coefficient in the central plug flow region, the volumes of the dead regions, and the exchange rates between the different regions. The superficial velocity of the gas through the reactor has a large effect on parameter values. Increased superficial velocity tends to decrease dead region volumes, interregion transfer rates, and axial dispersion. The significant deviation from CSTR, plug flow, and plug flow with axial dispersion of the residence time distribution of gas within small-scale reactors can lead to underestimation of the calculation of mass and heat transfer coefficients and hence has implications for reactor design and scaleup. (C) 2001 John Wiley & Sons, Inc.

Limiting conditional distributions for birth-death processes

In a recent paper [16], one of us identified all of the quasi-stationary distributions for a non-explosive, evanescent birth-death process for which absorption is certain, and established conditions for the existence of the corresponding limiting conditional distributions. Our purpose is to extend these results in a number of directions. We shall consider separately two cases depending on whether or not the process is evanescent. In the former case we shall relax the condition that absorption is certain. Furthermore, we shall allow for the possibility that the minimal process might be explosive, so that the transition rates alone will not necessarily determine the birth-death process uniquely. Although we shall be concerned mainly with the minimal process, our most general results hold for any birth-death process whose transition probabilities satisfy both the backward and the forward Kolmogorov differential equations.

Proposal for the interaction of non-conventional partial agonists and catecholamines with the 'putative beta(4)-adrenoceptor' in mammalian heart

1. Evidence for a 'putative beta(4)-adrenoceptor' originated over 20 years ago when cardiostimulant effects were observed to nonconventional partial agonists, These agonists were originally described as beta(1)- and beta(2)-adrenoceptor antagonists; however, they cause cardiostimulant effects at much higher concentrations than those required to block beta(1)- and beta(2)-adrenoceptors. Cardiostimulant effects of non-conventional partial agonists have been observed in mouse, rat, guinea-pig, cat, ferret and human heart tissues, 2. The receptor is expressed in several heart regions, including the sinoatrial node, atrium and ventricle, 3. The receptor is resistant to blockade by most antagonists that possess high affinity for beta(1)- and beta(2)- adrenoceptors, but is blocked with moderate affinity by (-)-bupranolol and CGP 20712A. 4. The receptor is pharmacologically distinct from the beta(3)-adrenoceptor. Micromolar concentrations of beta(3)-adrenoceptor agonists have no agonist or blocking activity, The receptor is also resistant to blockade by a beta(3)-adrenoceptor-selective antagonist. 5. The receptor mediates increases in cAMP levels and cAMP-dependent protein kinase (PK) A activity in cardiac tissues. Phosphodiesterase inhibition potentiates the positive chronotropic and inotropic effects of non-conventional partial agonists. 6. The receptor mediates hastening of atrial and ventricular relaxation, which is consistent with involvement of a cAMP-dependent pathway. 7. The non-conventional partial agonist (-)-[H-3]-CGP 12177A labels the cardiac putative beta(4)-adrenoceptor, Non-conventional partial agonists compete for binding with affinities that are closely similar to their agonist potencies, Catecholamines compete for binding in a stereoselective manner with a rank order of affinity of (-)-R0363 > (-)-isoprenaline > (-)-noradrenaline greater than or equal to (-)-adrenaline much greater than (-)-isoprenaline, suggesting that catecholamines can interact with the receptor. 8. The putative beta(4)-adrenoceptor appears to be coupled to the G(s)-adenylyl cyclase system, which could serve as a guide to its future cloning, Activation of the receptor may plausibly improve diastolic function but could also mediate arrhythmias.