A Quantile Variant of the EM Algorithm and Its Applications to Parameter Estimation with Interval Data

The expectation-maximization (EM) algorithm is a powerful computational technique for finding the maximum likelihood estimates for parametric models when the data are not fully observed. The EM is best suited for situations where the expectation in each E-step and the maximization in each M-step are straightforward. A difficulty with the implementation of the EM algorithm is that each E-step requires the integration of the log-likelihood function in closed form. The explicit integration can be avoided by using what is known as the Monte Carlo EM (MCEM) algorithm. The MCEM uses a random sample to estimate the integral at each E-step. However, the problem with the MCEM is that it often converges to the integral quite slowly and the convergence behavior can also be unstable, which causes a computational burden. In this paper, we propose what we refer to as the quantile variant of the EM (QEM) algorithm. We prove that the proposed QEM method has an accuracy of $O(1/K^2)$ while the MCEM method has an accuracy of $O_p(1/\sqrt{K})$. Thus, the proposed QEM method possesses faster and more stable convergence properties when compared with the MCEM algorithm. The improved performance is illustrated through the numerical studies. Several practical examples illustrating its use in interval-censored data problems are also provided

Anti-Allergic Cromones Inhibit Histamine and Eicosanoid Release from Activated Human and Murine Mast Cells by Releasing Annexin A1

PMCID: PMC3601088This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

A reference relative time-scale as an alternative to chronological age for cohorts with long follow-up

Background: Epidemiologists have debated the appropriate time-scale for cohort survival studies; chronological age or time-on-study being two such time-scales. Importantly, assessment of risk factors may depend on the choice of time-scale. Recently, chronological or attained age has gained support but a case can be made for a ‘reference relative time-scale’ as an alternative which circumvents difficulties that arise with this and other scales. The reference relative time of an individual participant is the integral of a reference population hazard function between time of entry and time of exit of the individual. The objective here is to describe the reference relative time-scale, illustrate its use, make comparison with attained age by simulation and explain its relationship to modern and traditional epidemiologic methods. Results: A comparison was made between two models; a stratified Cox model with age as the time-scale versus an un-stratified Cox model using the reference relative time-scale. The illustrative comparison used a UK cohort of cotton workers, with differing ages at entry to the study, with accrual over a time period and with long follow-up. Additionally, exponential and Weibull models were fitted since the reference relative time-scale analysis need not be restricted to the Cox model. A simulation study showed that analysis using the reference relative time-scale and analysis using chronological age had very similar power to detect a significant risk factor and both were equally unbiased. Further, the analysis using the reference relative time-scale supported fully-parametric survival modelling and allowed percentile predictions and mortality curves to be constructed. Conclusions: The reference relative time-scale was a viable alternative to chronological age, led to simplification of the modelling process and possessed the defined features of a good time-scale as defined in reliability theory. The reference relative time-scale has several interpretations and provides a unifying concept that links contemporary approaches in survival and reliability analysis to the traditional epidemiologic methods of Poisson regression and standardised mortality ratios. The community of practitioners has not previously made this connection

Emplacement of inflated Pāhoehoe flows in the Naude’s Nek Pass, Lesotho remnant, Karoo continental flood basalt province: use of flow-lobe tumuli in understanding flood basalt emplacement

Physical volcanological features are presented for a 710-m-thick section, of the Naude’s Nek Pass, within the lower part of the Lesotho remnant of the Karoo Large Igneous Province. The section consists of inflated pāhoehoe lava with thin, impersistent sedimentary interbeds towards the base. There are seven discreet packages of compound and hummocky pāhoehoe lobes containing flow-lobe tumuli, making up approximately 50% of the section. Approximately 45% of the sequence consists of 14 sheet lobes, between 10 and 52-m-thick. The majority of the sheet lobes are in two packages indicating prolonged periods of lava supply capable of producing thick sheet lobes. The other sheet lobes are as individual lobes or pairs, within compound flows, suggesting brief increases in lava supply rate. We suggest, contrary to current belief, that there is no evidence that compound flows are proximal to source and sheet lobes (simple flows) are distal to source and we propose that the presence of flow-lobe tumuli in compound flows could be an indicator that a flow is distal to source. We use detailed, previously published, studies of the Thakurvadi Formation (Deccan Traps) as an example. We show that the length of a lobe and therefore the sections that are ‘medial or distal to source’ are specific to each individual lobe and are dependent on the lava supply of each eruptive event, and as such flow lobe tumuli can be used as an indicator of relative distance from source

Steady-state modulation of voltage-gated K+ channels in rat arterial smooth muscle by cyclic AMP-dependent protein kinase and protein phosphatase 2B

Voltage-gated potassium channels (Kv) are important regulators of membrane potential in vascular smooth muscle cells, which is integral to controlling intracellular Ca2+ concentration and regulating vascular tone. Previous work indicates that Kv channels can be modulated by receptor-driven alterations of cyclic AMP-dependent protein kinase (PKA) activity. Here, we demonstrate that Kv channel activity is maintained by tonic activity of PKA. Whole-cell recording was used to assess the effect of manipulating PKA signalling on Kv and ATP-dependent K+ channels of rat mesenteric artery smooth muscle cells. Application of PKA inhibitors, KT5720 or H89, caused a significant inhibition of Kv currents. Tonic PKA-mediated activation of Kv appears maximal as application of isoprenaline (a β-adrenoceptor agonist) or dibutyryl-cAMP failed to enhance Kv currents. We also show that this modulation of Kv by PKA can be reversed by protein phosphatase 2B/calcineurin (PP2B). PKA-dependent inhibition of Kv by KT5720 can be abrogated by pre-treatment with the PP2B inhibitor cyclosporin A, or inclusion of a PP2B auto-inhibitory peptide in the pipette solution. Finally, we demonstrate that tonic PKA-mediated modulation of Kv requires intact caveolae. Pre-treatment of the cells with methyl-β-cyclodextrin to deplete cellular cholesterol, or adding caveolin-scaffolding domain peptide to the pipette solution to disrupt caveolae-dependent signalling each attenuated PKA-mediated modulation of the Kv current. These findings highlight a novel, caveolae-dependent, tonic modulatory role of PKA on Kv channels providing new insight into mechanisms and the potential for pharmacological manipulation of vascular tone

Unsupervised Bayesian linear unmixing of gene expression microarrays

Background: This paper introduces a new constrained model and the corresponding algorithm, called unsupervised Bayesian linear unmixing (uBLU), to identify biological signatures from high dimensional assays like gene expression microarrays. The basis for uBLU is a Bayesian model for the data samples which are represented as an additive mixture of random positive gene signatures, called factors, with random positive mixing coefficients, called factor scores, that specify the relative contribution of each signature to a specific sample. The particularity of the proposed method is that uBLU constrains the factor loadings to be non-negative and the factor scores to be probability distributions over the factors. Furthermore, it also provides estimates of the number of factors. A Gibbs sampling strategy is adopted here to generate random samples according to the posterior distribution of the factors, factor scores, and number of factors. These samples are then used to estimate all the unknown parameters. Results: Firstly, the proposed uBLU method is applied to several simulated datasets with known ground truth and compared with previous factor decomposition methods, such as principal component analysis (PCA), non negative matrix factorization (NMF), Bayesian factor regression modeling (BFRM), and the gradient-based algorithm for general matrix factorization (GB-GMF). Secondly, we illustrate the application of uBLU on a real time-evolving gene expression dataset from a recent viral challenge study in which individuals have been inoculated with influenza A/H3N2/Wisconsin. We show that the uBLU method significantly outperforms the other methods on the simulated and real data sets considered here. Conclusions: The results obtained on synthetic and real data illustrate the accuracy of the proposed uBLU method when compared to other factor decomposition methods from the literature (PCA, NMF, BFRM, and GB-GMF). The uBLU method identifies an inflammatory component closely associated with clinical symptom scores collected during the study. Using a constrained model allows recovery of all the inflammatory genes in a single factor

Localisation of RNAs into the germ plasm of vitellogenic xenopus oocytes

We have studied the localisation of mRNAs in full-grown Xenopus laevis oocytes by injecting fluorescent RNAs, followed by confocal microscopy of the oocyte cortex. Concentrating on RNA encoding the Xenopus Nanos homologue, nanos1 (formerly Xcat2), we find that it consistently localised into aggregated germ plasm ribonucleoprotein (RNP) particles, independently of cytoskeletal integrity. This implies that a diffusion/entrapment-mediated mechanism is active, as previously reported for previtellogenic oocytes. Sometimes this was accompanied by localisation into scattered particles of the “late”, Vg1/VegT pathway; occasionally only late pathway localisation was seen. The Xpat RNA behaved in an identical fashion and for neither RNA was the localisation changed by any culture conditions tested. The identity of the labelled RNP aggregates as definitive germ plasm was confirmed by their inclusion of abundant mitochondria and co-localisation with the germ plasm protein Hermes. Further, the nanos1/Hermes RNP particles are interspersed with those containing the germ plasm protein Xpat. These aggregates may be followed into the germ plasm of unfertilized eggs, but with a notable reduction in its quantity, both in terms of injected molecules and endogenous structures. Our results conflict with previous reports that there is no RNA localisation in large oocytes, and that during mid-oogenesis even germ plasm RNAs localise exclusively by the late pathway. We find that in mid oogenesis nanos1 RNA also localises to germ plasm but also by the late pathway. Late pathway RNAs, Vg1 and VegT, also may localise into germ plasm. Our results support the view that mechanistically the two modes of localisation are extremely similar, and that in an injection experiment RNAs might utilise either pathway, the distinction in fates being very subtle and subject to variation. We discuss these results in relation to their biological significance and the results of others

Critical animal and media studies: Expanding the understanding of oppression in communication research

Critical and communication studies have traditionally neglected the oppression conducted by humans towards other animals. However, our (mis)treatment of other animals is the result of public consent supported by a morally speciesist-anthropocentric system of values. Speciesism or anthroparchy, as much as any other mainstream ideologies, feeds the media and at the same time is perpetuated by them. The goal of this article is to remedy this neglect by introducing the subdiscipline of Critical Animal and Media Studies. Critical Animal and Media Studies takes inspiration both from critical animal studies – which is so far the most consolidated critical field of research in the social sciences addressing our exploitation of other animals – and from the normative-moral stance rooted in the cornerstones of traditional critical media studies. The authors argue that the Critical Animal and Media Studies approach is an unavoidable step forward for critical media and communication studies to engage with the expanded circle of concerns of contemporary ethical thinking

Ears of the Armadillo: Global Health Research and Neglected Diseases in Texas

Neglected tropical diseases (NTDs) have\ud been recently identified as significant public\ud health problems in Texas and elsewhere in\ud the American South. A one-day forum on the\ud landscape of research and development and\ud the hidden burden of NTDs in Texas\ud explored the next steps to coordinate advocacy,\ud public health, and research into a\ud cogent health policy framework for the\ud American NTDs. It also highlighted how\ud U.S.-funded global health research can serve\ud to combat these health disparities in the\ud United States, in addition to benefiting\ud communities abroad

The influence of soil communities on the temperature sensitivity of soil respiration

Soil respiration represents a major carbon flux between terrestrial ecosystems and the atmosphere, and is expected to accelerate under climate warming. Despite its importance in climate change forecasts, however, our understanding of the effects of temperature on soil respiration (RS) is incomplete. Using a metabolic ecology approach we link soil biota metabolism, community composition and heterotrophic activity, to predict RS rates across five biomes. We find that accounting for the ecological mechanisms underpinning decomposition processes predicts climatological RS variations observed in an independent dataset (n = 312). The importance of community composition is evident because without it RS is substantially underestimated. With increasing temperature, we predict a latitudinal increase in RS temperature sensitivity, with Q10 values ranging between 2.33 ±0.01 in tropical forests to 2.72 ±0.03 in tundra. This global trend has been widely observed, but has not previously been linked to soil communities