491 research outputs found

    Accurate Profiling of Microbial Communities from Massively Parallel Sequencing using Convex Optimization

    Full text link
    We describe the Microbial Community Reconstruction ({\bf MCR}) Problem, which is fundamental for microbiome analysis. In this problem, the goal is to reconstruct the identity and frequency of species comprising a microbial community, using short sequence reads from Massively Parallel Sequencing (MPS) data obtained for specified genomic regions. We formulate the problem mathematically as a convex optimization problem and provide sufficient conditions for identifiability, namely the ability to reconstruct species identity and frequency correctly when the data size (number of reads) grows to infinity. We discuss different metrics for assessing the quality of the reconstructed solution, including a novel phylogenetically-aware metric based on the Mahalanobis distance, and give upper-bounds on the reconstruction error for a finite number of reads under different metrics. We propose a scalable divide-and-conquer algorithm for the problem using convex optimization, which enables us to handle large problems (with 106\sim10^6 species). We show using numerical simulations that for realistic scenarios, where the microbial communities are sparse, our algorithm gives solutions with high accuracy, both in terms of obtaining accurate frequency, and in terms of species phylogenetic resolution.Comment: To appear in SPIRE 1

    Systematic model behavior of adsorption on flat surfaces

    Full text link
    A low density film on a flat surface is described by an expansion involving the first four virial coefficients. The first coefficient (alone) yields the Henry's law regime, while the next three correct for the effects of interactions. The results permit exploration of the idea of universal adsorption behavior, which is compared with experimental data for a number of systems

    Lombardi Drawings of Graphs

    Full text link
    We introduce the notion of Lombardi graph drawings, named after the American abstract artist Mark Lombardi. In these drawings, edges are represented as circular arcs rather than as line segments or polylines, and the vertices have perfect angular resolution: the edges are equally spaced around each vertex. We describe algorithms for finding Lombardi drawings of regular graphs, graphs of bounded degeneracy, and certain families of planar graphs.Comment: Expanded version of paper appearing in the 18th International Symposium on Graph Drawing (GD 2010). 13 pages, 7 figure

    Effect of pre-cardiac and adult stages of Dirofilaria immitis in pulmonary disease of cats: CBC, bronchial lavage cytology, serology, radiographs, CT images, bronchial reactivity, and histopathology

    Get PDF
    AbstractA controlled, blind study was conducted to define the initial inflammatory response and lung damage associated with the death of precardiac stages of Dirofilaria immitis in cats as compared to adult heartworm infections and normal cats. Three groups of six cats each were used: UU: uninfected untreated controls; PreS I: infected with 100 D. immitis L3 by subcutaneous injection and treated topically with selamectin 32 and 2 days pre-infection and once monthly for 8 months); IU: infected with 100 D. immitis L3 and left untreated. Peripheral blood, serum, bronchial lavage, and thoracic radiographic images were collected from all cats on Days 0, 70, 110, 168, and 240. CT images were acquired on Days 0, 110, and 240. Cats were euthanized, and necropsies were conducted on Day 240 to determine the presence of heartworms. Bronchial rings were collected for in vitro reactivity. Lung, heart, brain, kidney, and liver tissues were collected for histopathology. Results were compared for changes within each group. Pearson and Spearman correlations were performed for association between histologic, radiographic, serologic, hematologic and bronchoalveolar lavage (BAL) results. Infected cats treated with selamectin did not develop radiographically evident changes throughout the study, were heartworm antibody negative, and were free of adult heartworms and worm fragments at necropsy. Histologic lung scores and CT analysis were not significantly different between PreS I cats and UU controls. Subtle alveolar myofibrosis was noted in isolated areas of several PreS I cats and an eosinophilic BAL cytology was noted on Days 75 and 120. Bronchial ring reactivity was blunted in IU cats but was normal in PreS I and UU cats. The IU cats became antibody positive, and five cats developed adult heartworms. All cats with heartworms were antigen positive at one time point; but one cat was antibody positive, antigen negative, with viable adult females at necropsy. The CT revealed early involvement of all pulmonary arteries and a random pattern of parenchymal disease with severe lesions immediately adjacent to normal areas. Analysis of CT 3D reconstruction and Hounsfield units demonstrated lung disease consistent with restrictive pulmonary fibrosis with an interstitial infiltrate, absence of air trapping, and decrease in total lung volume in Group IU as compared to Groups UU and PreS I. The clinical implications of this study are that cats pretreated with selamectin 1 month before D. immitis L3 infection did not become serologically positive and did not develop pulmonary arterial hypertrophy and myofibrosis

    The acceleration of the universe and the physics behind it

    Full text link
    Using a general classification of dark enegy models in four classes, we discuss the complementarity of cosmological observations to tackle down the physics beyond the acceleration of our universe. We discuss the tests distinguishing the four classes and then focus on the dynamics of the perturbations in the Newtonian regime. We also exhibit explicitely models that have identical predictions for a subset of observations.Comment: 18 pages, 18 figure

    Cosmological distance indicators

    Full text link
    We review three distance measurement techniques beyond the local universe: (1) gravitational lens time delays, (2) baryon acoustic oscillation (BAO), and (3) HI intensity mapping. We describe the principles and theory behind each method, the ingredients needed for measuring such distances, the current observational results, and future prospects. Time delays from strongly lensed quasars currently provide constraints on H0H_0 with < 4% uncertainty, and with 1% within reach from ongoing surveys and efforts. Recent exciting discoveries of strongly lensed supernovae hold great promise for time-delay cosmography. BAO features have been detected in redshift surveys up to z <~ 0.8 with galaxies and z ~ 2 with Ly-α\alpha forest, providing precise distance measurements and H0H_0 with < 2% uncertainty in flat Λ\LambdaCDM. Future BAO surveys will probe the distance scale with percent-level precision. HI intensity mapping has great potential to map BAO distances at z ~ 0.8 and beyond with precisions of a few percent. The next years ahead will be exciting as various cosmological probes reach 1% uncertainty in determining H0H_0, to assess the current tension in H0H_0 measurements that could indicate new physics.Comment: Review article accepted for publication in Space Science Reviews (Springer), 45 pages, 10 figures. Chapter of a special collection resulting from the May 2016 ISSI-BJ workshop on Astronomical Distance Determination in the Space Ag

    Charge Transport in one Dimension:Dissipative and Non-Dissipative Space-Charge Limited Currents

    Full text link
    We consider charge transport in nanopores where the dielectric constant inside the nanopore is much greater than in the surrounding material, so that the flux of the electric fields due to the charges is almost entirely confined to the nanopore. That means that we may model the electric fields due to charge densities in the nanopore in terms of average properties across the nanopore as solutions of one dimensional Poisson equations. We develop basic equations for an M component system using equations of continuity to relate concentrations to currents, and flux equations relating currents to concentration gradients and conductivities. We then derive simplified scaled versions of the equations. We develop exact solutions for the one component case in a variety of boundary conditions using a Hopf-Cole transformation, Fourier series, and periodic solutions of the Burgers equation. These are compared with a simpler model in which the scaled diffusivity is zero so that all charge motion is driven by the electric field. In this non-dissipative case, recourse to an admissibility condition is utilised to obtain the physically relevant weak solution of a Riemann problem concerning the electric field. It is shown that the admissibility condition is Poynting's theorem.Comment: v2 contains the minor updates and corrections that have been incorporated into the published articl

    Velocity-space sensitivity of the time-of-flight neutron spectrometer at JET

    Get PDF
    The velocity-space sensitivities of fast-ion diagnostics are often described by so-called weight functions. Recently, we formulated weight functions showing the velocity-space sensitivity of the often dominant beam-target part of neutron energy spectra. These weight functions for neutron emission spectrometry (NES) are independent of the particular NES diagnostic. Here we apply these NES weight functions to the time-of-flight spectrometer TOFOR at JET. By taking the instrumental response function of TOFOR into account, we calculate time-of-flight NES weight functions that enable us to directly determine the velocity-space sensitivity of a given part of a measured time-of-flight spectrum from TOFOR

    The solution of Burgers' equation for sinusoidal excitation at the upstream boundary

    Full text link
    This paper generates an exact solution to Burgers' nonlinear diffusion equation on a convective stream with sinusoidal excitation applied at the upstream boundary, x =0. The downstream boundary, effectively at x =∞, is assumed to always be far enough ahead of the convective front at x=V s t that no disturbance is felt there. The Hopf-Cole transformation is applied in achieving the analytical solution, but only after integrating the equation and its conditions in x to avoid a nonlinearity in the transformed upstream boundary condition.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/42700/1/10665_2006_Article_BF02383570.pd
    corecore