1,277 research outputs found
Clustering by compression
We present a new method for clustering based on compression. The method
doesn't use subject-specific features or background knowledge, and works as
follows: First, we determine a universal similarity distance, the normalized
compression distance or NCD, computed from the lengths of compressed data files
(singly and in pairwise concatenation). Second, we apply a hierarchical
clustering method. The NCD is universal in that it is not restricted to a
specific application area, and works across application area boundaries. A
theoretical precursor, the normalized information distance, co-developed by one
of the authors, is provably optimal but uses the non-computable notion of
Kolmogorov complexity. We propose precise notions of similarity metric, normal
compressor, and show that the NCD based on a normal compressor is a similarity
metric that approximates universality. To extract a hierarchy of clusters from
the distance matrix, we determine a dendrogram (binary tree) by a new quartet
method and a fast heuristic to implement it. The method is implemented and
available as public software, and is robust under choice of different
compressors. To substantiate our claims of universality and robustness, we
report evidence of successful application in areas as diverse as genomics,
virology, languages, literature, music, handwritten digits, astronomy, and
combinations of objects from completely different domains, using statistical,
dictionary, and block sorting compressors. In genomics we presented new
evidence for major questions in Mammalian evolution, based on
whole-mitochondrial genomic analysis: the Eutherian orders and the Marsupionta
hypothesis against the Theria hypothesis.Comment: LaTeX, 27 pages, 20 figure
Algorithmic Clustering of Music
We present a fully automatic method for music classification, based only on
compression of strings that represent the music pieces. The method uses no
background knowledge about music whatsoever: it is completely general and can,
without change, be used in different areas like linguistic classification and
genomics. It is based on an ideal theory of the information content in
individual objects (Kolmogorov complexity), information distance, and a
universal similarity metric. Experiments show that the method distinguishes
reasonably well between various musical genres and can even cluster pieces by
composer.Comment: 17 pages, 11 figure
An Approximation Algorithm for Ancestral Maximum-Likelihood and Phylogeography Inference Problems under Time Reversible Markov Evolutionary Models
The ancestral maximum-likelihood and phylogeography problems are two
fundamental problems involving evolutionary studies. The ancestral
maximum-likelihood problem involves identifying a rooted tree alongside
internal node sequences that maximizes the probability of observing a given set
of sequences as leaves. The phylogeography problem extends the ancestral
maximum-likelihood problem to incorporate geolocation of leaf and internal
nodes. While a constant factor approximation algorithm has been established for
the ancestral maximum-likelihood problem concerning two-state sequences, no
such algorithm has been devised for any generalized instances of the problem.
In this paper, we focus on a generalization of the two-state model, the time
reversible Markov evolutionary models for sequences and geolocations. Under
this evolutionary model, we present a -approximation algorithm,
where is the number of input samples, addressing both the ancestral
maximum-likelihood and phylogeography problems. This is the first approximation
algorithm for the phylogeography problem. Furthermore, we show how to apply the
algorithm on popular evolutionary models like generalized time-reversible (GTR)
model and its specialization Jukes and Cantor 69 (JC69)
Computational Molecular Biology
Computational Biology is a fairly new subject that arose in response to the computational problems posed by the analysis and the processing of biomolecular sequence and structure data. The field was initiated in the late 60's and early 70's largely by pioneers working in the life sciences. Physicists and mathematicians entered the field in the 70's and 80's, while Computer Science became involved with the new biological problems in the late 1980's. Computational problems have gained further importance in molecular biology through the various genome projects which produce enormous amounts of data. For this bibliography we focus on those areas of computational molecular biology that involve discrete algorithms or discrete optimization. We thus neglect several other areas of computational molecular biology, like most of the literature on the protein folding problem, as well as databases for molecular and genetic data, and genetic mapping algorithms. Due to the availability of review papers and a bibliography this bibliography
Bayesian Statistical Methods for Genetic Association Studies with Case-Control and Cohort Design
Large-scale genetic association studies are carried out with the hope of discovering single
nucleotide polymorphisms involved in the etiology of complex diseases. We propose a
coalescent-based model for association mapping which potentially increases the power to
detect disease-susceptibility variants in genetic association studies with case-control and cohort
design. The approach uses Bayesian partition modelling to cluster haplotypes with
similar disease risks by exploiting evolutionary information. We focus on candidate gene
regions and we split the chromosomal region of interest into sub-regions or windows of high
linkage disequilibrium (LD) therein assuming a perfect phylogeny. The haplotype space is
then partitioned into disjoint clusters within which the phenotype-haplotype association is
assumed to be the same. The novelty of our approach consists in the fact that the distance
used for clustering haplotypes has an evolutionary interpretation, as haplotypes are clustered
according to the time to their most recent common mutation. Our approach is fully
Bayesian and we develop Markov Chain Monte Carlo algorithms to sample efficiently over
the space of possible partitions. We have also developed a Bayesian survival regression model
for high-dimension and small sample size settings. We provide a Bayesian variable selection
procedure and shrinkage tool by imposing shrinkage priors on the regression coefficients. We
have developed a computationally efficient optimization algorithm to explore the posterior
surface and find the maximum a posteriori estimates of the regression coefficients. We compare
the performance of the proposed methods in simulation studies and using real datasets
to both single-marker analyses and recently proposed multi-marker methods and show that
our methods perform similarly in localizing the causal allele while yielding lower false positive
rates. Moreover, our methods offer computational advantages over other multi-marker
approaches
Phylogenetic Trees and Their Analysis
Determining the best possible evolutionary history, the lowest-cost phylogenetic tree, to fit a given set of taxa and character sequences using maximum parsimony is an active area of research due to its underlying importance in understanding biological processes. As several steps in this process are NP-Hard when using popular, biologically-motivated optimality criteria, significant amounts of resources are dedicated to both both heuristics and to making exact methods more computationally tractable. We examine both phylogenetic data and the structure of the search space in order to suggest methods to reduce the number of possible trees that must be examined to find an exact solution for any given set of taxa and associated character data. Our work on four related problems combines theoretical insight with empirical study to improve searching of the tree space. First, we show that there is a Hamiltonian path through tree space for the most common tree metrics, answering Bryant\u27s Challenge for the minimal such path. We next examine the topology of the search space under various metrics, showing that some metrics have local maxima and minima even with perfect data, while some others do not. We further characterize conditions for which sequences simulated under the Jukes-Cantor model of evolution yield well-behaved search spaces. Next, we reduce the search space needed for an exact solution by splitting the set of characters into mutually-incompatible subsets of compatible characters, building trees based on the perfect phylogenies implied by these sets, and then searching in the neighborhoods of these trees. We validate this work empirically. Finally, we compare two approaches to the generalized tree alignment problem, or GTAP: Sequence alignment followed by tree search vs. Direct Optimization, on both biological and simulated data
Automated Gene Classification using Nonnegative Matrix Factorization on Biomedical Literature
Understanding functional gene relationships is a challenging problem for biological applications. High-throughput technologies such as DNA microarrays have inundated biologists with a wealth of information, however, processing that information remains problematic. To help with this problem, researchers have begun applying text mining techniques to the biological literature. This work extends previous work based on Latent Semantic Indexing (LSI) by examining Nonnegative Matrix Factorization (NMF). Whereas LSI incorporates the singular value decomposition (SVD) to approximate data in a dense, mixed-sign space, NMF produces a parts-based factorization that is directly interpretable. This space can, in theory, be used to augment existing ontologies and annotations by identifying themes within the literature. Of course, performing NMF does not come without a priceânamely, the large number of parameters. This work attempts to analyze the effects of some of the NMF parameters on both convergence and labeling accuracy. Since there is a dearth of automated label evaluation techniques as well as âgold standardâ hierarchies, a method to produce âcorrectâ trees is proposed as well as a technique to label trees and to evaluate those labels
- âŠ