Distributions associated with general runs and patterns in hidden Markov models

This paper gives a method for computing distributions associated with patterns in the state sequence of a hidden Markov model, conditional on observing all or part of the observation sequence. Probabilities are computed for very general classes of patterns (competing patterns and generalized later patterns), and thus, the theory includes as special cases results for a large class of problems that have wide application. The unobserved state sequence is assumed to be Markovian with a general order of dependence. An auxiliary Markov chain is associated with the state sequence and is used to simplify the computations. Two examples are given to illustrate the use of the methodology. Whereas the first application is more to illustrate the basic steps in applying the theory, the second is a more detailed application to DNA sequences, and shows that the methods can be adapted to include restrictions related to biological knowledge.Comment: Published in at http://dx.doi.org/10.1214/07-AOAS125 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Dynamic filtering of static dipoles in magnetoencephalography

We consider the problem of estimating neural activity from measurements of the magnetic fields recorded by magnetoencephalography. We exploit the temporal structure of the problem and model the neural current as a collection of evolving current dipoles, which appear and disappear, but whose locations are constant throughout their lifetime. This fully reflects the physiological interpretation of the model. In order to conduct inference under this proposed model, it was necessary to develop an algorithm based around state-of-the-art sequential Monte Carlo methods employing carefully designed importance distributions. Previous work employed a bootstrap filter and an artificial dynamic structure where dipoles performed a random walk in space, yielding nonphysical artefacts in the reconstructions; such artefacts are not observed when using the proposed model. The algorithm is validated with simulated data, in which it provided an average localisation error which is approximately half that of the bootstrap filter. An application to complex real data derived from a somatosensory experiment is presented. Assessment of model fit via marginal likelihood showed a clear preference for the proposed model and the associated reconstructions show better localisation

A 6,000-year record of environmental change from the eastern Pacific margin of central Mexico

The transition from the mid- to late-Holocene in MesoAmerica saw increasing complexity in spatial patterns of change. Records from the western part of the region are sparse, with lacustrine sequences affected by long term anthropogenic disturbance or lacking chronological resolution. Here, we present a continuous palaeoecological and geochemical record from Laguna de Juanacatlán, a remote lake in the mountains of the western TMVB. Diatom assemblages, XRF scanning data and bulk organic geochemistry from a well-dated, 7.25-m laminated sequence were combined with summary pollen data from a 9-m partially laminated core to provide a continuous record of catchment and lake ecosystem changes over the last c. 6,000 years. Relatively humid conditions prevailed prior to c. 5.1 cal ka, which supported dense oak-pine forest cover around a deep, stratified lake. A trend towards drier conditions began c. 5.1 cal ka, intensifying after 4.0 cal ka, consistent with weakening of the North American Monsoon. Between 3.0 and 1.2 cal ka, lower lake levels and variable catchment run-off are consistent with increasing ENSO influence observed in the Late Holocene in the neotropics. From 1.2 to 0.9 cal ka, a marked change to catchment stability and more intense stratification reflected drier conditions and / or reduced rainfall variability and possibly warmer temperatures. After 0.9 cal ka, conditions were wetter, with an increase in catchment disturbance associated with the combined effects of climate and human activity. In recent decades, the lake ecosystem has changed markedly, possibly in response to recent climate change as well as local catchment dynamics

In silico approaches to study mass and energy flows in microbial consortia: a syntrophic case study

<p>Abstract</p> <p>Background</p> <p>Three methods were developed for the application of stoichiometry-based network analysis approaches including elementary mode analysis to the study of mass and energy flows in microbial communities. Each has distinct advantages and disadvantages suitable for analyzing systems with different degrees of complexity and <it>a priori </it>knowledge. These approaches were tested and compared using data from the thermophilic, phototrophic mat communities from Octopus and Mushroom Springs in Yellowstone National Park (USA). The models were based on three distinct microbial guilds: oxygenic phototrophs, filamentous anoxygenic phototrophs, and sulfate-reducing bacteria. Two phases, day and night, were modeled to account for differences in the sources of mass and energy and the routes available for their exchange.</p> <p>Results</p> <p>The <it>in silico </it>models were used to explore fundamental questions in ecology including the prediction of and explanation for measured relative abundances of primary producers in the mat, theoretical tradeoffs between overall productivity and the generation of toxic by-products, and the relative robustness of various guild interactions.</p> <p>Conclusion</p> <p>The three modeling approaches represent a flexible toolbox for creating cellular metabolic networks to study microbial communities on scales ranging from cells to ecosystems. A comparison of the three methods highlights considerations for selecting the one most appropriate for a given microbial system. For instance, communities represented only by metagenomic data can be modeled using the pooled method which analyzes a community's total metabolic potential without attempting to partition enzymes to different organisms. Systems with extensive <it>a priori </it>information on microbial guilds can be represented using the compartmentalized technique, employing distinct control volumes to separate guild-appropriate enzymes and metabolites. If the complexity of a compartmentalized network creates an unacceptable computational burden, the nested analysis approach permits greater scalability at the cost of more user intervention through multiple rounds of pathway analysis.</p

Discovery of Samarium, Europium, Gadolinium, and Terbium Isotopes

Currently, thirty-four samarium, thirty-four europium, thirty-one gadolinium, and thirty-one terbium isotopes have been observed and the discovery of these isotopes is discussed here. For each isotope a brief synopsis of the first refereed publication, including the production and identification method, is presented.Comment: To be published in At. Data Nucl. Data Table

Orbit size and estimated eye size in dinosaurs and other archosaurs and their implications for the evolution of visual capabilities

Vision is one of the most important senses for animals, allowing them to interact with their environment and with further implications for evolutionary histories. However, relevant soft tissues, such as the eye and associated structures, are not preserved in fossil vertebrates, limiting our knowledge of their visual capabilities. Here, we quantified absolute and relative orbit size for 400 species of dinosaurs and other extinct archosaurs using linear measurements of the preserved skeletal elements as a proxy for visual capabilities. Our results demonstrate that the orbit makes up on average 20% of skull size with a strong and consistent correlation across all sampled groups. This trend is largely independent of temporal distribution, species richness, and phylogeny. In fact, relative orbit size is narrowly constrained and did not surpass 45% of skull size, suggesting physiological and functional controls. Estimated eye size was found to be absolutely larger in herbivores, whereas carnivores tended to have smaller eyes absolutely and compared with skull size. Relatively large eyes only occurred in small-bodied species and vice versa. However, eye size alone was not sufficient to discriminate between different activity patterns or to characterize visual capabilities in detail.</p

A proposal for a coordinated effort for the determination of brainwide neuroanatomical connectivity in model organisms at a mesoscopic scale

In this era of complete genomes, our knowledge of neuroanatomical circuitry remains surprisingly sparse. Such knowledge is however critical both for basic and clinical research into brain function. Here we advocate for a concerted effort to fill this gap, through systematic, experimental mapping of neural circuits at a mesoscopic scale of resolution suitable for comprehensive, brain-wide coverage, using injections of tracers or viral vectors. We detail the scientific and medical rationale and briefly review existing knowledge and experimental techniques. We define a set of desiderata, including brain-wide coverage; validated and extensible experimental techniques suitable for standardization and automation; centralized, open access data repository; compatibility with existing resources, and tractability with current informatics technology. We discuss a hypothetical but tractable plan for mouse, additional efforts for the macaque, and technique development for human. We estimate that the mouse connectivity project could be completed within five years with a comparatively modest budget.Comment: 41 page

Variants in KCNQ1 increase type II diabetes susceptibility in South Asians: A study of 3,310 subjects from India and the US

<p>Abstract</p> <p>Background</p> <p>Polymorphisms in intron 15 of potassium voltage-gated channel, KQT-like subfamily member 1 (<it>KCNQ1</it>) gene have been associated with type II diabetes (T2D) in Japanese genome-wide association studies (GWAS). More recently a meta-analysis of European GWAS has detected a new independent signal associated with T2D in intron 11 of the <it>KCNQ1 </it>gene. The purpose of this investigation is to examine the role of these variants with T2D in populations of Asian Indian descent from India and the US.</p> <p>Methods</p> <p>We examined the association between four variants in the <it>KCNQ1 </it>gene with T2D and related quantitative traits in a total of 3,310 Asian Indian participants from two different cohorts comprising 2,431 individuals of the Punjabi case-control cohort from the Sikh Diabetes Study and 879 migrant Asian Indians living in the US.</p> <p>Results</p> <p>Our data confirmed the association of a new signal at the <it>KCNQ1 </it>locus (rs231362) with T2D showing an allelic odds ratio (OR) of 1.24 95%CI [1.08-1.43], p = 0.002 in the Punjabi cohort. A moderate association with T2D was also seen for rs2237895 in the Punjabi (OR 1.14; p = 0.036) and combined cohorts (meta-analysis OR 1.14; p = 0.018). Three-site haplotype analysis of rs231362, rs2237892, rs2237895 exhibited considerably stronger evidence of association of the GCC haplotype with T2D showing OR of 1.24 95%CI [1.00-1.53], p = 0.001, permutation p = 8 × 10^-4in combined cohorts. The 'C' risk allele carriers of rs2237895 had significantly reduced measures of HOMA-B in the US cohort (p = 0.008) as well as in combined cohort in meta-analysis (p = 0.009).</p> <p>Conclusions</p> <p>Our investigation has confirmed that the variation within the <it>KCNQ1 </it>locus confers a significant risk to T2D among Asian Indians. Haplotype analysis further suggested that the T2D risk associated with <it>KCNQ1 </it>SNPs may be derived from 'G' allele of rs231362 and 'C' allele of rs2237895 and this appears to be mediated through β cell function.</p

Measurement of the Branching Fraction for B- --> D0 K*-

We present a measurement of the branching fraction for the decay B- --> D0 K*- using a sample of approximately 86 million BBbar pairs collected by the BaBar detector from e+e- collisions near the Y(4S) resonance. The D0 is detected through its decays to K- pi+, K- pi+ pi0 and K- pi+ pi- pi+, and the K*- through its decay to K0S pi-. We measure the branching fraction to be B.F.(B- --> D0 K*-)= (6.3 +/- 0.7(stat.) +/- 0.5(syst.)) x 10^{-4}.Comment: 7 pages, 1 postscript figure, submitted to Phys. Rev. D (Rapid Communications

Measurement of the quasi-elastic axial vector mass in neutrino-oxygen interactions

The weak nucleon axial-vector form factor for quasi-elastic interactions is determined using neutrino interaction data from the K2K Scintillating Fiber detector in the neutrino beam at KEK. More than 12,000 events are analyzed, of which half are charged-current quasi-elastic interactions nu-mu n to mu- p occurring primarily in oxygen nuclei. We use a relativistic Fermi gas model for oxygen and assume the form factor is approximately a dipole with one parameter, the axial vector mass M_A, and fit to the shape of the distribution of the square of the momentum transfer from the nucleon to the nucleus. Our best fit result for M_A = 1.20 \pm 0.12 GeV. Furthermore, this analysis includes updated vector form factors from recent electron scattering experiments and a discussion of the effects of the nucleon momentum on the shape of the fitted distributions.Comment: 14 pages, 10 figures, 6 table