Objectives. We introduce methods for the exploratory analysis of microarray data, especially focusing on cluster algorithms. Benefits and problems are discussed. Methods. We describe application and suitability of uns...
详细信息
Objectives. We introduce methods for the exploratory analysis of microarray data, especially focusing on cluster algorithms. Benefits and problems are discussed. Methods. We describe application and suitability of unsupervised learning methods for the classification of gene expression data. cluster algorithms are treated in more detail, including assessment of cluster quality. Results: When dealing with microarray data, most cluster algorithms must be applied with caution. As long as the structure of the true generating models of such data is not fully understood, the use of simple algorithms seems to be more appropriate than the application of complex block-box algorithms. New methods explicitly targeted to the analysis of microarray data are increasingly being developed in order to increase the amount of useful information extracted from the experiments. Conclusions: Unsupervised methods can be a helpful tool for the analysis of microarray data, but a critical choice of the algorithm and a careful interpretation of the results are required in order to avoid false conclusions.
The numerical simulation of various field theories at non-zero chemical potential suffers from severe complex action problems. In particular, QCD at non-zero quark density can presently not be simulated for that reaso...
详细信息
The numerical simulation of various field theories at non-zero chemical potential suffers from severe complex action problems. In particular, QCD at non-zero quark density can presently not be simulated for that reason. A similar complex action problem arises in the 2D O(3) model-a toy model for QCD. Here we construct the 2D 0(3) model at non-zero density via dimensional reduction of an antiferromagnetic quantum spin ladder in a magnetic field. The complex action problem of the 2D 0(3) model manifests itself as a sign problem of the ladder system. This sign problem is solved completely with a meron cluster algorithm. (C) 2002 Elsevier Science B.V. All rights reserved.
The numerical simulation of various field theories at non-zero chemical potential suffers from severe complex action problems. In particular, QCD at non-zero quark density can presently not be simulated for that reaso...
详细信息
The numerical simulation of various field theories at non-zero chemical potential suffers from severe complex action problems. In particular, QCD at non-zero quark density can presently not be simulated for that reason. A similar complex action problem arises in the 2D O(3) model-a toy model for QCD. Here we construct the 2D 0(3) model at non-zero density via dimensional reduction of an antiferromagnetic quantum spin ladder in a magnetic field. The complex action problem of the 2D 0(3) model manifests itself as a sign problem of the ladder system. This sign problem is solved completely with a meron cluster algorithm. (C) 2002 Elsevier Science B.V. All rights reserved.
A cluster methodology, motivated via density estimation, is proposed. It is based on the idea of estimating the population clusters, which, following Hartigan (1975), are defined as the connected parts of the "su...
详细信息
A cluster methodology, motivated via density estimation, is proposed. It is based on the idea of estimating the population clusters, which, following Hartigan (1975), are defined as the connected parts of the "substantial" support of the underlying density. The empirical clusters are defined by analogy in terms of the substantial support of a convolution (kemel-type) density estimator. The sample observations are grouped into data clusters, according to the empirical cluster they belong. An algorithm to implement the method, based on resampling ideas, is proposed. It allows either to automatically choose the number of clusters or to give this number as an input. Some theoretical and practical aspects are briefly discussed and a simulation study is given. The results show a good performance of our method, in terms of efficiency and robustness, when compared with two classical cluster algorithms: k-means and single linkage. Finally, a real-data example is discussed. (C) 2001 Elsevier Science B.V. All rights reserved.
We continue the study, initiated in Part I, of graphical representations and cluster algorithms for various models in (or related to) statistical mechanics. For certain models, e.g. the Blume-Emery-Griffths model and ...
详细信息
We continue the study, initiated in Part I, of graphical representations and cluster algorithms for various models in (or related to) statistical mechanics. For certain models, e.g. the Blume-Emery-Griffths model and various generalizations, we develop Fortuin Kasteleyn-type representations which lead immediately to Swendsen Wang-type algorithms. For other models, e.g. the random cluster model, that are defined by a graphical representation, we develop cluster algorithms without reference to an underlying spin system. In all cases, phase transitions are related to percolation (or incipient percolation) in the graphical representation which, via the IC algorithm, allows for the rapid simulation of these systems at the transition point. Pertinent examples include the (continuum) Widom-Rowlinson model, the restricted 1-step solid-on-solid model and the XY model. (C) 1998 Published by Elsevier Science B.V. All rights reserved.
After a brief general overview of Monte Carlo computer simulations in statistical physics, special emphasis is placed on applications to phase transitions and critical phenomena. Here, standard simulations employing l...
详细信息
After a brief general overview of Monte Carlo computer simulations in statistical physics, special emphasis is placed on applications to phase transitions and critical phenomena. Here, standard simulations employing local update algorithms are severely hampered by the problem of critical slowing down, that is by strong correlations between successively generated data It is shown that this problem can be greatly reduced by using nonlocal update techniques such as cluster and multigrid algorithms. The general ideas are illustrated for simple lattice spin models and Euclidean path integrals. (C) 1998 IMACS/ Elsevier Science B.V.
Damage spreading for Ising cluster dynamics :is investigated numerically by using random numbers in a way that conforms with the notion of submitting the two evolving replicas to the same thermal noise. Two damage spr...
详细信息
Damage spreading for Ising cluster dynamics :is investigated numerically by using random numbers in a way that conforms with the notion of submitting the two evolving replicas to the same thermal noise. Two damage spreading transitions are found;damage does not spread either at low or high temperatures. We determine some critical exponents at the high-temperature transition point, which seem consistent with directed percolation.
The Wolff single-cluster algorithm is the most efficient method known for Monte Carlo simulation of many spin models. Due to the irregular size, shape and position of the Wolff clusters, this method does not easily le...
详细信息
The Wolff single-cluster algorithm is the most efficient method known for Monte Carlo simulation of many spin models. Due to the irregular size, shape and position of the Wolff clusters, this method does not easily lend itself to efficient parallel implementation, so that simulations using this method have thus far been confined to workstations and vector machines. Here we present two parallel implementations of this algorithm, and show that one gives fairly good performance on a MIMD parallel computer.
We examine a number of models that generate random fractals. The models are studied using the tools of computational complexity theory from the perspective of parallel computation. Diffusion-limited aggregation and se...
详细信息
We examine a number of models that generate random fractals. The models are studied using the tools of computational complexity theory from the perspective of parallel computation. Diffusion-limited aggregation and several widely used algorithms for equilibrating the Ising model are shown to be highly sequential;it is unlikely they can be simulated efficiently in parallel. This is in contrast to Mandelbrot percolation, which can be simulated in constant parallel time. Our research helps shed light on the intrinsic complexity of these models relative to each other and to different growth processes that have been recently studied using complexity theory. In addition, the results may serve as a guide to simulation physics.
We present algorithmic procedures for generating systematically ideas and solutions to problems which are perceived as creative. Our method consists of identifying and characterizing the most creative ideas among a va...
详细信息
We present algorithmic procedures for generating systematically ideas and solutions to problems which are perceived as creative. Our method consists of identifying and characterizing the most creative ideas among a vast pool. We show that they fall within a few large classes (archetypes) which share the same conceptual structure (Macros). We prescribe well defined abstract algorithms which can act deterministically on arbitrary given objects. Each algorithm generates ideas with the same conceptual structure characteristic to one of the Macros. The resulting new ideas turn out to be perceived as highly creative. We support our claims by experiments in which senior advertising professionals graded advertisement ideas produced by our method according to their creativity. The marks (grade 4.6 +/- 0.2 on a 1-7 scale) obtained by laymen applying our algorithms (after being instructed for only two hours) were significantly better than the marks obtained by advertising professionals using standard methods (grade 3.6 +/- 0.2)). The method, which is currently taught in USA, Europe, and Israel and used by advertising agencies in Britain and Israel has received formal international recognition.
暂无评论