CMBfit: Rapid WMAP likelihood calculations with normal parameters

We present a method for ultra-fast confrontation of the WMAP cosmic microwave background observations with theoretical models, implemented as a publicly available software package called CMBfit, useful for anyone wishing to measure cosmological parameters by combining WMAP with other observations. The method takes advantage of the underlying physics by transforming into a set of parameters where the WMAP likelihood surface is accurately fit by the exponential of a quartic or sextic polynomial. Building on previous physics based approximations by Hu et.al., Kosowsky et.al. and Chu et.al., it combines their speed with precision cosmology grade accuracy. A Fortran code for computing the WMAP likelihood for a given set of parameters is provided, pre-calibrated against CMBfast, accurate to Delta lnL ~ 0.05 over the entire 2sigma region of the parameter space for 6 parameter ``vanilla'' Lambda CDM models. We also provide 7-parameter fits including spatial curvature, gravitational waves and a running spectral index.Comment: 14 pages, 8 figures, References added, accepted for publication in Phys.Rev.D., a Fortran code can be downloaded from http://space.mit.edu/home/tegmark/cmbfit

An Adaptive Markov Chain Monte Carlo Method for GARCH Model

We propose a method to construct a proposal density for the Metropolis-Hastings algorithm in Markov Chain Monte Carlo (MCMC) simulations of the GARCH model. The proposal density is constructed adaptively by using the data sampled by the MCMC metho d itself. It turns out that autocorrelations between the data generated with our adaptive proposal density are greatly reduced. Thus it is concluded that the adaptive construction method is very efficient and works well for the MCMC simulations of the GARCH model.Comment: 11 pages, 6 figure

Sensitivity analysis of reactive ecological dynamics

Author Posting. © Springer, 2008. This is the author's version of the work. It is posted here by permission of Springer for personal use, not for redistribution. The definitive version was published in Bulletin of Mathematical Biology 70 (2008): 1634-1659, doi:10.1007/s11538-008-9312-7.Ecological systems with asymptotically stable equilibria may exhibit significant transient dynamics following perturbations. In some cases, these transient dynamics include the possibility of excursions away from the equilibrium before the eventual return; systems that exhibit such amplification of perturbations are called reactive. Reactivity is a common property of ecological systems, and the amplification can be large and long-lasting. The transient response of a reactive ecosystem depends on the parameters of the underlying model. To investigate this dependence, we develop sensitivity analyses for indices of transient dynamics (reactivity, the amplification envelope, and the optimal perturbation) in both continuous- and discrete-time models written in matrix form. The sensitivity calculations require expressions, some of them new, for the derivatives of equilibria, eigenvalues, singular values, and singular vectors, obtained using matrix calculus. Sensitivity analysis provides a quantitative framework for investigating the mechanisms leading to transient growth. We apply the methodology to a predator-prey model and a size-structured food web model. The results suggest predator-driven and prey-driven mechanisms for transient amplification resulting from multispecies interactions.Financial support provided by NSF grant DEB-0343820, NOAA grant NA03-NMF4720491, the Ocean Life Institute of the Woods Hole Oceanographic Institution, and the Academic Programs Office of the MIT-WHOI Joint Program in Oceanography

Effects of dietary vegetable oil on atlantic salmon hepatocyte fatty acid desaturation and liver fatty acid compositions

Fatty acyl desaturase activities, involved in the conversion of the C18 EFA, 18:2n-6 and 18:3n-3, to the highly unsaturated fatty acids (HUFA) 20:4n-6, 20:5n-3 and 22:6n-3, are known to be under nutritional regulation. Specifically, the activity of the desaturation/elongation pathway is depressed when animals, including fish, are fed fish oils rich in n-3HUFA compared to animals fed vegetable oils rich in C18 EFA. The primary aims of the present study were a) to establish the relative importance of product inhibition (n-3HUFA) versus increased substrate concentration (C18 EFA) and, b) to determine whether 18:2n-6 and 18:3n-3 differ in their effects, on the hepatic fatty acyl desaturation/elongation pathway in Atlantic salmon (Salmo salar). Smolts were fed ten experimental diets containing blends of two vegetable oils, linseed (LO) and rapeseed oil (RO), and fish oil (FO) in a triangular mixture design for 50 weeks. Fish were sampled after 32 and 50 weeks, lipid and fatty acid composition of liver determined, fatty acyl desaturation/elongation activity estimated in hepatocytes using [1-14C]18:3n-3 as substrate, and the data subjected to regression analyses. Dietary 18:2n-6 was positively correlated, and n-3HUFA negatively correlated, with lipid content of liver. Dietary 20:5n-3 and 22:6n-3 were positively correlated with liver fatty acids with a slope greater than unity suggesting relative retention and deposition of these HUFA. In contrast, dietary 18:2n-6 and 18:3n-3 were positively correlated with liver fatty acids with a slope of less than unity suggesting metabolism via β-oxidation and/or desaturation/elongation. Consistent with this, fatty acyl desaturation/elongation in hepatocytes was significantly increased by feeding diets containing vegetable oils. Dietary 20:5n-3 and 22:6n-3 levels were negatively correlated with hepatocyte fatty acyl desaturation. At 32 weeks, 18:2n-6 but not 18:3n-3, was positively correlated with hepatocyte fatty acyl desaturation activity whereas the reverse was true at 50 weeks. The data indicate that both feedback inhibition through increased n-3HUFA and decreased C18 fatty acyl substrate concentration are probably important in determining hepatocyte fatty acyl desaturation activities, and that 18:2n-6 and 18:3n-3 may differ in their effects on this pathway

Highly unsaturated fatty acid synthesis in marine fish: Cloning, functional characterization, and nutritional regulation of fatty acyl delta6 desaturase of Atlantic cod (Gadus morhua L.)

Fish contain high levels of the n-3 highly unsaturated fatty acids (HUFA), eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids that are crucial to the health of higher vertebrates. Biosynthesis of HUFA requires enzyme-mediated desaturation of fatty acids. Here we report cloning and functional characterisation of a ∆6 fatty acyl desaturase of Atlantic cod (Gadus morhua), and describe its tissue expression and nutritional regulation. PCR primers were designed based on the sequences of conserved motifs in available fish desaturases and used to isolate a cDNA fragment from liver of cod. The full-length cDNA was obtained by Rapid Amplification of cDNA Ends (RACE). The cDNA for the putative fatty acyl desaturase was shown to comprise 1980bp which included a 5’-UTR of 261bp and a 3’-UTR of 375bp. Sequencing revealed that the cDNA included an ORF of 1344 bp that specified a protein of 447 amino acids. The protein sequence included three histidine boxes, two transmembrane regions, and an N-terminal cytochrome b5 domain containing the haem-binding motif HPGG, all of which are characteristic of microsomal fatty acid desaturases. The cDNA displayed Δ6 desaturase activity in a heterologous yeast expression system. Quantitative real time PCR assay of gene expression in cod showed that the ∆6 desaturase gene, was highly expressed in brain, relatively highly expressed in liver, kidney, intestine, red muscle and gill, and expressed at much lower levels in white muscle, spleen and heart. In contrast, the abundance of a cod fatty acyl elongase transcript was high in brain and gill, with intermediate levels in kidney, spleen, intestine and heart, and relatively low expression in liver. The expression of the Δ6 desaturase gene and the PUFA elongase gene may be under a degree of nutritional regulation, with levels being marginally increased in livers and intestine of fish fed a vegetable oil blend by comparison with levels in fish fed fish oil. However, this was not reflected in increased Δ6 desaturase activity in hepatocytes or enterocytes, which showed very little highly unsaturated fatty acid biosynthesis activity irrespective of diet. The study described has demonstrated that Atlantic cod express a fatty acid desaturase gene with functional Δ6 activity in a yeast expression system. This is consistent with an established hypothesis that the poor ability of marine fish to synthesise HUFA is not due to lack of a Δ6 desaturase, but rather to deficiencies in other parts of the biosynthetic pathway. However, further studies are required to determine why the Δ6 desaturase appears to be barely functional in cod under the conditions tested

Cosmological parameters from SDSS and WMAP

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.Comment: Minor revisions to match accepted PRD version. SDSS data and ppt figures available at http://www.hep.upenn.edu/~max/sdsspars.htm

Online Interactive Teaching Modules Enhance Quantitative Proficiency of Introductory Biology Students

There is widespread agreement within the scientific and education communities that undergraduate biology curricula fall short in providing students with the quantitative and interdisciplinary problem-solving skills they need to obtain a deep understanding of biological phenomena and be prepared fully to contribute to future scientific inquiry. MathBench Biology Modules were designed to address these needs through a series of interactive, Web-based modules that can be used to supplement existing course content across the biological sciences curriculum. The effect of the modules was assessed in an introductory biology course at the University of Maryland. Over the course of the semester, students showed significant increases in quantitative skills that were independent of previous math course work. Students also showed increased comfort with solving quantitative problems, whether or not they ultimately arrived at the correct answer. A survey of spring 2009 graduates indicated that those who had experienced MathBench in their course work had a greater appreciation for the role of mathematics in modern biology than those who had not used MathBench. MathBench modules allow students from diverse educational backgrounds to hone their quantitative skills, preparing them for more complex mathematical approaches in upper-division courses

Measurement of νˉμ\bar{\nu}_{\mu} and νμ\nu_{\mu} charged current inclusive cross sections and their ratio with the T2K off-axis near detector

We report a measurement of cross section $\sigma(\nu_{\mu}+{\rm nucleus}\rightarrow\mu^{-}+X)$ and the first measurements of the cross section $\sigma(\bar{\nu}_{\mu}+{\rm nucleus}\rightarrow\mu^{+}+X)$ and their ratio $R(\frac{\sigma(\bar \nu)}{\sigma(\nu)})$ at (anti-)neutrino energies below 1.5 GeV. We determine the single momentum bin cross section measurements, averaged over the T2K $\bar{\nu}/\nu$-flux, for the detector target material (mainly Carbon, Oxygen, Hydrogen and Copper) with phase space restricted laboratory frame kinematics of $\theta_{\mu}$500 MeV/c. The results are $\sigma(\bar{\nu})=\left( 0.900\pm0.029{\rm (stat.)}\pm0.088{\rm (syst.)}\right)\times10^{-39}$ and $\sigma(\nu)=\left( 2.41\ \pm0.022{\rm{(stat.)}}\pm0.231{\rm (syst.)}\ \right)\times10^{-39}$in units of cm$^{2}$/nucleon and$R\left(\frac{\sigma(\bar{\nu})}{\sigma(\nu)}\right)= 0.373\pm0.012{\rm (stat.)}\pm0.015{\rm (syst.)}$.Comment: 18 pages, 8 figure

Disorder-assisted error correction in Majorana chains

It was recently realized that quenched disorder may enhance the reliability of topological qubits by reducing the mobility of anyons at zero temperature. Here we compute storage times with and without disorder for quantum chains with unpaired Majorana fermions - the simplest toy model of a quantum memory. Disorder takes the form of a random site-dependent chemical potential. The corresponding one-particle problem is a one-dimensional Anderson model with disorder in the hopping amplitudes. We focus on the zero-temperature storage of a qubit encoded in the ground state of the Majorana chain. Storage and retrieval are modeled by a unitary evolution under the memory Hamiltonian with an unknown weak perturbation followed by an error-correction step. Assuming dynamical localization of the one-particle problem, we show that the storage time grows exponentially with the system size. We give supporting evidence for the required localization property by estimating Lyapunov exponents of the one-particle eigenfunctions. We also simulate the storage process for chains with a few hundred sites. Our numerical results indicate that in the absence of disorder, the storage time grows only as a logarithm of the system size. We provide numerical evidence for the beneficial effect of disorder on storage times and show that suitably chosen pseudorandom potentials can outperform random ones.Comment: 50 pages, 7 figure

Harmonising and linking biomedical and clinical data across disparate data archives to enable integrative cross-biobank research

A wealth of biospecimen samples are stored in modern globally distributed biobanks. Biomedical researchers worldwide need to be able to combine the available resources to improve the power of large-scale studies. A prerequisite for this effort is to be able to search and access phenotypic, clinical and other information about samples that are currently stored at biobanks in an integrated manner. However, privacy issues together with heterogeneous information systems and the lack of agreed-upon vocabularies have made specimen searching across multiple biobanks extremely challenging. We describe three case studies where we have linked samples and sample descriptions in order to facilitate global searching of available samples for research. The use cases include the ENGAGE (European Network for Genetic and Genomic Epidemiology) consortium comprising at least 39 cohorts, the SUMMIT (surrogate markers for micro- and macro-vascular hard endpoints for innovative diabetes tools) consortium and a pilot for data integration between a Swedish clinical health registry and a biobank. We used the Sample avAILability (SAIL) method for data linking: first, created harmonised variables and then annotated and made searchable information on the number of specimens available in individual biobanks for various phenotypic categories. By operating on this categorised availability data we sidestep many obstacles related to privacy that arise when handling real values and show that harmonised and annotated records about data availability across disparate biomedical archives provide a key methodological advance in pre-analysis exchange of information between biobanks, that is, during the project planning phase