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

Occupational balance: What tips the scales for new students?

The open question, ‘What prevents you from reaching occupational balance?’, was posed within a questionnaire aimed at exploring the meanings of occupation, health and wellbeing with a cohort of first-year occupational therapy students during their initial few weeks at university. Their written responses to the question about occupational balance were analysed and are discussed in this paper. Not surprisingly, occupational balance appeared to be achieved by only a few and more by chance than design. People, time and money factors were identified as the main impediments to achieving occupational balance, with psychological and emotional pressures being at the forefront. Interestingly, despite these barriers, the overall educational benefit of considering the occupational balance question in this way raised the students’ awareness of its relationship to health and wellbeing. This increased awareness might have longer-term health benefits, both personally and professionally, which would be worthy of further research

Simonsenia aveniformis sp nov (Bacillariophyceae), molecular phylogeny and systematics of the genus, and a new type of canal raphe system

The genus Simonsenia is reviewed and S. aveniformis described as new for science by light and electron microscopy. The new species originated from estuarine environments in southern Iberia (Atlantic coast) and was isolated into culture. In LM, Simonsenia resembles Nitzschia, with bridges (fibulae) beneath the raphe, which is marginal. It is only electron microscope (EM) examination that reveals the true structure of the raphe system, which consists of a raphe canal raised on a keel (wing), supported by rib like braces (fenestral bars) and tube-like portulae; between the portulae the keel is perforated by open windows (fenestrae). Based on the presence of portulae and a fenestrated keel, Simonsenia has been proposed to be intermediate between Bacillariaceae and Surirellaceae. However, an rbcL phylogeny revealed that Simonsenia belongs firmly in the Bacillariaceae, with which it shares a similar chloroplast arrangement, rather than in the Surirellaceae. Lack of homology between the surirelloid and simonsenioid keels is reflected in subtle differences in the morphology and ontogeny of the portulae and fenestrae. The diversity of Simonsenia has probably been underestimated, particularly in the marine environment.Polish National Science Centre in Cracow within the Maestro program [N 2012/04/A/ST10/00544]; Sciences and Technologies Foundation-FCT (Portugal) [SFRH/BD/62405/2009]info:eu-repo/semantics/publishedVersio

The demand for sports and exercise: Results from an illustrative survey

Funding from the Department of Health policy research programme was used in this study.There is a paucity of empirical evidence on the extent to which price and perceived benefits affect the level of participation in sports and exercise. Using an illustrative sample of 60 adults at Brunel University, West London, we investigate the determinants of demand for sports and exercise. The data were collected through face-to-face interviews that covered indicators of sports and exercise behaviour; money/time price and perceived benefits of participation; and socio- economic/demographic details. Count, linear and probit regression models were fitted as appropriate. Seventy eight per cent of the sample participated in sports and exercise and spent an average of £27 per month and an average of 20 min travelling per occasion of sports and exercise. The demand for sport and exercise was negatively associated with time (travel or access time) and ‘variable’ price and positively correlated with ‘fixed’ price. Demand was price inelastic, except in the case of meeting the UK government’s recommended level of participation, which is time price elastic (elasticity = −2.2). The implications of data from a larger nationally representative sample as well as the role of economic incentives in influencing uptake of sports and exercise are discussed.This article is available through the Brunel Open Access Publishing Fund

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

Charting the landscape of N=4 flux compactifications

We analyse the vacuum structure of isotropic Z_2 x Z_2 flux compactifications, allowing for a single set of sources. Combining algebraic geometry with supergravity techniques, we are able to classify all vacua for both type IIA and IIB backgrounds with arbitrary gauge and geometric fluxes. Surprisingly, geometric IIA compactifications lead to a unique theory with four different vacua. In this case we also perform the general analysis allowing for sources compatible with minimal supersymmetry. Moreover, some relevant examples of type IIB non-geometric compactifications are studied. The computation of the full N=4 mass spectrum reveals the presence of a number of non-supersymmetric and nevertheless stable AdS_4 vacua. In addition we find a novel dS_4 solution based on a non-semisimple gauging.Comment: Minor corrections and references added. Version published in JHE

Challenges of Profile Likelihood Evaluation in Multi-Dimensional SUSY Scans

Statistical inference of the fundamental parameters of supersymmetric theories is a challenging and active endeavor. Several sophisticated algorithms have been employed to this end. While Markov-Chain Monte Carlo (MCMC) and nested sampling techniques are geared towards Bayesian inference, they have also been used to estimate frequentist confidence intervals based on the profile likelihood ratio. We investigate the performance and appropriate configuration of MultiNest, a nested sampling based algorithm, when used for profile likelihood-based analyses both on toy models and on the parameter space of the Constrained MSSM. We find that while the standard configuration is appropriate for an accurate reconstruction of the Bayesian posterior, the profile likelihood is poorly approximated. We identify a more appropriate MultiNest configuration for profile likelihood analyses, which gives an excellent exploration of the profile likelihood (albeit at a larger computational cost), including the identification of the global maximum likelihood value. We conclude that with the appropriate configuration MultiNest is a suitable tool for profile likelihood studies, indicating previous claims to the contrary are not well founded.Comment: 21 pages, 9 figures, 1 table; minor changes following referee report. Matches version accepted by JHE

Semiparametric Multivariate Accelerated Failure Time Model with Generalized Estimating Equations

The semiparametric accelerated failure time model is not as widely used as the Cox relative risk model mainly due to computational difficulties. Recent developments in least squares estimation and induced smoothing estimating equations provide promising tools to make the accelerate failure time models more attractive in practice. For semiparametric multivariate accelerated failure time models, we propose a generalized estimating equation approach to account for the multivariate dependence through working correlation structures. The marginal error distributions can be either identical as in sequential event settings or different as in parallel event settings. Some regression coefficients can be shared across margins as needed. The initial estimator is a rank-based estimator with Gehan's weight, but obtained from an induced smoothing approach with computation ease. The resulting estimator is consistent and asymptotically normal, with a variance estimated through a multiplier resampling method. In a simulation study, our estimator was up to three times as efficient as the initial estimator, especially with stronger multivariate dependence and heavier censoring percentage. Two real examples demonstrate the utility of the proposed method