@article{3614,
  abstract     = {We analyze the changes in the mean and variance components of a quantitative trait caused by changes in allele frequencies, concentrating on the effects of genetic drift. We use a general representation of epistasis and dominance that allows an arbitrary relation between genotype and phenotype for any number of diallelic loci. We assume initial and final Hardy-Weinberg and linkage equilibrium in our analyses of drift-induced changes. Random drift generates transient linkage disequilibria that cause correlations between allele frequency fluctuations at different loci. However, we show that these have negligible effects, at least for interactions among small numbers of loci. Our analyses are based on diffusion approximations that summarize the effects of drift in terms of F, the inbreeding coefficient, interpreted as the expected proportional decrease in heterozygosity at each locus. For haploids, the variance of the trait mean after a population bottleneck is var(Δz̄) =inline imagewhere n is the number of loci contributing to the trait variance, VA(1)=VA is the additive genetic variance, and VA(k) is the kth-order additive epistatic variance. The expected additive genetic variance after the bottleneck, denoted (V*A), is closely related to var(Δz̄); (V*A) (1 –F)inline imageThus, epistasis inflates the expected additive variance above VA(1 –F), the expectation under additivity. For haploids (and diploids without dominance), the expected value of every variance component is inflated by the existence of higher order interactions (e.g., third-order epistasis inflates (V*AA)). This is not true in general with diploidy, because dominance alone can reduce (V*A) below VA(1 –F) (e.g., when dominant alleles are rare). Without dominance, diploidy produces simple expressions: var(Δz̄)=inline image=1 (2F) kVA(k) and (V*A) = (1 –F)inline imagek(2F)k-1VA(k) With dominance (and even without epistasis), var(Δz̄)and (V*A) no longer depend solely on the variance components in the base population. For small F, the expected additive variance simplifies to (V*A)(1 –F) VA+ 4FVAA+2FVD+2FCAD, where CAD is a sum of two terms describing covariances between additive effects and dominance and additive × dominance interactions. Whether population bottlenecks lead to expected increases in additive variance depends primarily on the ratio of nonadditive to additive genetic variance in the base population, but dominance precludes simple predictions based solely on variance components. We illustrate these results using a model in which genotypic values are drawn at random, allowing extreme and erratic epistatic interactions. Although our analyses clarify the conditions under which drift is expected to increase VA, we question the evolutionary importance of such increases.},
  author       = {Nicholas Barton and Turelli, Michael},
  journal      = {Evolution; International Journal of Organic Evolution},
  number       = {10},
  pages        = {2111 -- 2132},
  publisher    = {Wiley-Blackwell},
  title        = {{Effects of allele frequency changes on variance components under a general model of epistasis}},
  doi          = {10.1111/j.0014-3820.2004.tb01591.x},
  volume       = {58},
  year         = {2004},
}

@article{3615,
  abstract     = {We investigate three alternative selection-based scenarios proposed to maintain polygenic variation: pleiotropic balancing selection, G x E interactions (with spatial or temporal variation in allelic effects), and sex-dependent allelic effects. Each analysis assumes an additive polygenic trait with n diallelic loci under stabilizing selection. We allow loci to have different effects and consider equilibria at which the population mean departs from the stabilizing-selection optimum. Under weak selection, each model produces essentially identical, approximate allele-frequency dynamics. Variation is maintained under pleiotropic balancing selection only at loci for which the strength of balancing selection exceeds the effective strength of stabilizing selection. In addition, for all models, polymorphism requires that the population mean be close enough to the optimum that directional selection does not overwhelm balancing selection. This balance allows many simultaneously stable equilibria, and we explore their properties numerically. Both spatial and temporal G x E can maintain variation at loci for which the coefficient of variation (across environments) of the effect of a substitution exceeds a critical value greater than one. The critical value depends on the correlation between substitution effects at different loci. For large positive correlations (e.g., ρ2ij > 3/4), even extreme fluctuations in allelic effects cannot maintain variation. Surprisingly, this constraint on correlations implies that sex-dependent allelic effects cannot maintain polygenic variation. We present numerical results that support our analytical approximations and discuss our results in connection to relevant data and alternative variance-maintaining mechanisms.},
  author       = {Turelli, Michael and Nicholas Barton},
  journal      = {Genetics},
  number       = {2},
  pages        = {1053 -- 1079},
  publisher    = {Genetics Society of America},
  title        = {{Polygenic variation maintained by balancing selection: pleiotropy, sex-dependent allelic effects and GxE interactions}},
  doi          = {10.1534/genetics.166.2.1053},
  volume       = {166},
  year         = {2004},
}

@misc{3616,
  author       = {Nicholas Barton},
  booktitle    = {Current Biology},
  number       = {15},
  pages        = {R603 -- R604},
  publisher    = {Cell Press},
  title        = {{Speciation: Why, how, where and when?}},
  doi          = {10.1016/j.cub.2004.07.037},
  volume       = {14},
  year         = {2004},
}

@article{3617,
  abstract     = {The coalescent process can describe the effects of selection at linked loci only if selection is so strong that genotype frequencies evolve deterministically. Here, we develop methods proposed by Kaplan, Darden, and Hudson to find the effects of weak selection. We show that the overall effect is given by an extension to Price's equation: the change in properties such as moments of coalescence times is equal to the covariance between those properties and the fitness of the sample of genes. The distribution of coalescence times differs substantially between allelic classes, even in the absence of selection. However, the average coalescence time between randomly chosen genes is insensitive to the current allele frequency and is affected significantly by purifying selection only if deleterious mutations are common and selection is strong (i.e., the product of population size and selection coefficient, Ns > 3). Balancing selection increases mean coalescence times, but the effect becomes large only when mutation rates between allelic classes are low and when selection is extremely strong. Our analysis supports previous simulations that show that selection has surprisingly little effect on genealogies. Moreover, small fluctuations in allele frequency due to random drift can greatly reduce any such effects. This will make it difficult to detect the action of selection from neutral variation alone.},
  author       = {Nicholas Barton and Etheridge, Alison M},
  journal      = {Genetics},
  number       = {2},
  pages        = {1115 -- 1131},
  publisher    = {Genetics Society of America},
  title        = {{The effect of selection on genealogies}},
  doi          = {10.1534/genetics.166.2.1115},
  volume       = {166},
  year         = {2004},
}

@inproceedings{3688,
  abstract     = {Capturing images of documents using handheld digital cameras has a variety of applications in academia, research, knowledge management, retail, and office settings. The ultimate goal of such systems is to achieve image quality comparable to that currently achieved with flatbed scanners even for curved, warped, or curled pages. This can be achieved by high-accuracy 3D modeling of the page surface, followed by a "flattening" of the surface. A number of previous systems have either assumed only perspective distortions, or used techniques like structured lighting, shading, or side-imaging for obtaining 3D shape. This paper describes a system for handheld camera-based document capture using general purpose stereo vision methods followed by a new document dewarping technique. Examples of shape modeling and dewarping of book images is shown.},
  author       = {Ulges, Adrian and Christoph Lampert and Breuel,Thomas M},
  pages        = {198 -- 200},
  publisher    = {ACM},
  title        = {{Document capture using stereo vision}},
  doi          = {10.1145/1030397.1030434},
  year         = {2004},
}

@article{3805,
  abstract     = {The operation of neuronal networks crucially depends on a fast time course of signaling in inhibitory interneurons. Synapses that excite interneurons generate fast currents, owing to the expression of glutamate receptors of specific subunit composition. Interneurons generate brief action potentials in response to transient synaptic activation and discharge repetitively at very high frequencies during sustained stimulation. The ability to generate short-duration action potentials at high frequencies depends on the expression of specific voltage-gated K+ channels. Factors facilitating fast action potential initiation following synaptic excitation include depolarized interneuron resting potential, subthreshold conductances and active dendrites. Finally, GABA release at interneuron output synapses is rapid and highly synchronized, leading to a faster inhibition in postsynaptic interneurons than in principal cells. Thus, the expression of distinct transmitter receptors and voltage-gated ion channels ensures that interneurons operate with high speed and temporal precision.},
  author       = {Peter Jonas and Bischofberger, Josef and Fricker, Desdemona and Miles, Richard},
  journal      = {Trends in Neurosciences},
  number       = {1},
  pages        = {30 -- 40},
  publisher    = {Elsevier},
  title        = {{Interneuron Diversity series: Fast in, fast out--temporal and spatial signal processing in hippocampal interneurons}},
  doi          = {doi:10.1016/j.tins.2003.10.010},
  volume       = {27},
  year         = {2004},
}

@article{3807,
  abstract     = {The time course of Mg(2+) block and unblock of NMDA receptors (NMDARs) determines the extent they are activated by depolarization. Here, we directly measure the rate of NMDAR channel opening in response to depolarizations at different times after brief (1 ms) and sustained (4.6 s) applications of glutamate to nucleated patches from neocortical pyramidal neurons. The kinetics of Mg(2+) unblock were found to be non-instantaneous and complex, consisting of a prominent fast component (time constant approximately 100 micros) and slower components (time constants 4 and approximately 300 ms), the relative amplitudes of which depended on the timing of the depolarizing pulse. Fitting a kinetic model to these data indicated that Mg(2+) not only blocks the NMDAR channel, but reduces both the open probability and affinity for glutamate, while enhancing desensitization. These effects slow the rate of NMDAR channel opening in response to depolarization in a time-dependent manner such that the slower components of Mg(2+) unblock are enhanced during depolarizations at later times after glutamate application. One physiological consequence of this is that brief depolarizations occurring earlier in time after glutamate application are better able to open NMDAR channels. This finding has important implications for spike-timing-dependent synaptic plasticity (STDP), where the precise (millisecond) timing of action potentials relative to synaptic inputs determines the magnitude and sign of changes in synaptic strength. Indeed, we find that STDP timing curves of NMDAR channel activation elicited by realistic dendritic action potential waveforms are narrower than expected assuming instantaneous Mg(2+) unblock, indicating that slow Mg(2+) unblock of NMDAR channels makes the STDP timing window more precise.},
  author       = {Kampa, Bjorn M and Clements, John and Peter Jonas and Stuart, Greg J},
  journal      = {Journal of Physiology},
  number       = {Pt 2},
  pages        = {337 -- 45},
  publisher    = {Wiley-Blackwell},
  title        = {{Kinetics of Mg(2+) unblock of NMDA receptors: implications for spike-timing dependent synaptic plasticity}},
  doi          = {10.1113/jphysiol.2003.058842 },
  volume       = {556},
  year         = {2004},
}

@article{3809,
  abstract     = {Neural stem cells in various regions of the vertebrate brain continuously generate neurons throughout life. In the mammalian hippocampus, a region important for spatial and episodic memory, thousands of new granule cells are produced per day, with the exact number depending on environmental conditions and physical exercise. The survival of these neurons is improved by learning and conversely learning may be promoted by neurogenesis. Although it has been suggested that newly generated neurons may have specific properties to facilitate learning, the cellular and synaptic mechanisms of plasticity in these neurons are largely unknown. Here we show that young granule cells in the adult hippocampus differ substantially from mature granule cells in both active and passive membrane properties. In young neurons, T-type Ca2+ channels can generate isolated Ca2+ spikes and boost fast Na+ action potentials, contributing to the induction of synaptic plasticity. Associative long-term potentiation can be induced more easily in young neurons than in mature neurons under identical conditions. Thus, newly generated neurons express unique mechanisms to facilitate synaptic plasticity, which may be important for the formation of new memories.},
  author       = {Schmidt-Hieber, Christoph and Peter Jonas and Bischofberger, Josef},
  journal      = {Nature},
  number       = {6988},
  pages        = {184 -- 7},
  publisher    = {Nature Publishing Group},
  title        = {{Enhanced synaptic plasticity in newly generated granule cells of the adult hippocampus}},
  doi          = {10.1038/nature02553},
  volume       = {429},
  year         = {2004},
}

@article{3810,
  abstract     = {Voltage-gated potassium (Kv) channels control action potential repolarization, interspike membrane potential, and action potential frequency in excitable cells. It is thought that the combinatorial association between distinct alpha and beta subunits determines whether Kv channels function as non-inactivating delayed rectifiers or as rapidly inactivating A-type channels. We show that membrane lipids can convert A-type channels into delayed rectifiers and vice versa. Phosphoinositides remove N-type inactivation from A-type channels by immobilizing the inactivation domains. Conversely, arachidonic acid and its amide anandamide endow delayed rectifiers with rapid voltage-dependent inactivation. The bidirectional control of Kv channel gating by lipids may provide a mechanism for the dynamic regulation of electrical signaling in the nervous system.},
  author       = {Oliver, Dominik and Lien, Cheng-Chang and Soom, Malle and Baukrowitz, Thomas and Peter Jonas and Fakler, Bernd},
  journal      = {Science},
  number       = {5668},
  pages        = {265 -- 70},
  publisher    = {American Association for the Advancement of Science},
  title        = {{Functional conversion between A-type and delayed rectifier K+ channels by membrane lipids}},
  doi          = {10.1126/science.1094113},
  volume       = {304},
  year         = {2004},
}

@inproceedings{3894,
  abstract     = {We study infinite stochastic games played by n-players on a finite graph with goals given by sets of infinite traces. The games are stochastic (each player simultaneously and independently chooses an action at each round, and the next state is determined by a probability distribution depending on the current state and the chosen actions), infinite (the game continues for an infinite number of rounds), nonzero sum (the players' goals are not necessarily conflicting), and undiscounted. We show that if each player has a reachability objective, that is, if the goal for each player i is to visit some subset R-i of the states, then there exists an epsilon-Nash equilibrium in memoryless strategies, for every epsilon > 0. However, exact Nash equilibria need not exist. We study the complexity of finding such Nash equilibria, and show that the payoff of some epsilon-Nash equilibrium in memoryless strategies can be epsilon-approximated in NP. We study the important subclass of n-player turn-based probabilistic games, where at each state at most one player has a nontrivial choice of moves. For turn-based probabilistic games, we show the existence of epsilon-Nash equilibria in pure strategies for games where the objective of player i is a Borel set B-i of infinite traces. However, exact Nash equilibria may not exist. For the special case of omega-regular objectives, we show exact Nash equilibria exist, and can be computed in NP when the omega-regular objectives are expressed as parity objectives.},
  author       = {Krishnendu Chatterjee and Majumdar, Ritankar S and Jurdziński, Marcin},
  pages        = {26 -- 40},
  publisher    = {Springer},
  title        = {{On Nash equilibria in stochastic games}},
  doi          = {10.1007/978-3-540-30124-0_6},
  volume       = {3210},
  year         = {2004},
}

@inproceedings{3895,
  abstract     = {In 2-player non-zero-sum games, Nash equilibria capture the options for rational behavior if each player attempts to maximize her payoff. In contrast to classical game theory, we consider lexicographic objectives: first, each player tries to maximize her own payoff, and then, the player tries to minimize the opponent's payoff. Such objectives arise naturally in the verification of systems with multiple components. There, instead of proving that each component satisfies its specification no matter how the other components behave, it often suffices to prove that each component satisfies its specification provided that the other components satisfy their specifications. We say that a Nash equilibrium is secure if it is an equilibrium with respect to the lexicographic objectives of both players. We prove that in graph games with Borel objectives, which include the games that arise in verification, there may be several Nash equilibria, but there is always a unique maximal payoff profile of secure equilibria. We show how this equilibrium can be computed in the case of omega-regular objectives, and we characterize the memory requirements of strategies that achieve the equilibrium.},
  author       = {Krishnendu Chatterjee and Thomas Henzinger and Jurdziński, Marcin},
  pages        = {160 -- 169},
  publisher    = {IEEE},
  title        = {{Games with secure equilibria}},
  doi          = {10.1109/LICS.2004.1319610},
  year         = {2004},
}

@article{3918,
  abstract     = {Wingless (ergatoid) males of the tramp ant Cardiocondyla minutior attack and kill their young ergatoid rivals and thus attempt to monopolize mating with female sexuals reared in the colony. Because of the different strength of local mate competition in colonies with one or several reproductive queens, we expected the production of new ergatoid males to vary with queen number. Sex ratios were mostly female-biased, but in contrast to the sympatric species C. obscurior (Cremer and Heinze, 2002) neither the percentage of ergatoid males nor of female sexuals among the first 20 sexuals produced varied considerably with queen number. As in C. obscurior, experimental colony fragmentation led to the production of winged males, whereas in unfragmented control colonies only ergatoid males eclosed.},
  author       = {Heinze, Jürgen and Böttcher, A. and Cremer, Sylvia},
  journal      = {Insectes Sociaux},
  number       = {3},
  pages        = {275 -- 278},
  publisher    = {Springer},
  title        = {{Production of winged and wingless males in the ant, Cardiocondyla minutior}},
  doi          = {10.1007/s00040-004-0740-6},
  volume       = {51},
  year         = {2004},
}

@article{3929,
  abstract     = {The Nef protein of human and simian immunodeficiency virus (HIV/SIV) is believed to interfere with T cell activation signals by forming a signaling complex at the plasma membrane. Composition and function of the complex are not fully understood. Here we report that Nef recruits the Polycomb Group (PcG) protein Eed, so far known as a nuclear factor and repressor of transcription, to the membrane of cells. The Nef-induced translocation of Eed led to a potent stimulation of Tat-dependent HIV transcription, implying that Eed removal from the nucleus is required for optimal Tat function. Similar to Nef action, activation of integrin receptors recruited Eed to the plasma membrane, also leading to enhanced Tat/Nef-mediated transcription. Our results suggest a link between membrane-associated activation processes and transcriptional derepression and demonstrate how HIV exploits this mechanism.},
  author       = {Witte, Vanessa and Laffert, Bernd and Rosorius, Olaf and Lischka, Peter and Blume, Katja and Galler, Gunther and Stilper, Andrea and Willbold, Dieter and D'Aloja, Paola and Michael Sixt and Kolanus, Johanna and Ott, Melanie and Kolanus, Waldemar and Schuler, Gerold and Baur, Andreas S},
  journal      = {Molecular Cell},
  number       = {2},
  pages        = {179 -- 190},
  publisher    = {Cell Press},
  title        = {{HIV-1 Nef mimics an integrin receptor signal that recruits the polycomb group protein Eed to the plasma membrane}},
  doi          = {10.1016/S1097-2765(04)00004-8},
  volume       = {13},
  year         = {2004},
}

@article{3931,
  abstract     = {Hyaluronan is an unsulfated glycosaminoglycan (GAG) that is ubiquitously expressed in the extracellular matrix (ECM) of all vertebrates, where hyaluronan rich matrices constitute a particular permissive environment for the development of complex biological structures and also for tumor progression. Because of its conserved structure and ubiquitous expression, antibodies for its histochemical detection cannot be produced. We have engineered a fusion protein, neurocan-GFP, and expressed it as a secreted molecule in mammalian cells. Neurocan-GFP fusion protein specifically binds to hyaluronan and directly visualizes hyaluronan on tissue sections, revealing a very detailed picture of hyaluronan distribution. The fluorescent fusion protein can be used in combination with antibodies and nuclear markers for double or triple staining. In addition, it is suitable to visualize hyaluronan on living cells by time-lapse video microscopy. The successful production and application of the neurocan-GFP fusion protein opens up new perspectives for using GFP fusion proteins as detection tools in histological and cytological studies complementing conventional antibody and biotin/avidin techniques.},
  author       = {Zhang, Hui and Baader, Stephan L and Michael Sixt and Kappler, Joachim and Rauch, Uwe},
  journal      = {Journal of Histochemistry and Cytochemistry},
  number       = {7},
  pages        = {915 -- 922},
  publisher    = {Histochemical Society},
  title        = {{Neurocan-GFP fusion protein: a new approach to detect hyaluronan on tissue sections and living cells}},
  doi          = {10.1369/jhc.3A6221.2004},
  volume       = {52},
  year         = {2004},
}

@article{3984,
  abstract     = {We combine topological and geometric methods to construct a multiresolution representation for a function over a two-dimensional domain. In a preprocessing stage, we create the Morse-Smale complex of the function and progressively simplify its topology by cancelling pairs of critical points. Based on a simple notion of dependency among these cancellations, we construct a hierarchical data structure supporting traversal and reconstruction operations similarly to traditional geometry-based representations. We use this data structure to extract topologically valid approximations that satisfy error bounds provided at runtime.},
  author       = {Bremer, Peer-Timo and Herbert Edelsbrunner and Hamann, Bernd and Pascucci, Valerio},
  journal      = {IEEE Transactions on Visualization and Computer Graphics},
  number       = {4},
  pages        = {385 -- 396},
  publisher    = {IEEE},
  title        = {{A topological hierarchy for functions on triangulated surfaces}},
  doi          = {10.1109/TVCG.2004.3},
  volume       = {10},
  year         = {2004},
}

@article{3985,
  abstract     = {Given a Morse function f over a 2-manifold with or without boundary, the Reeb graph is obtained by contracting the connected components of the level sets to points. We prove tight upper and lower bounds on the number of loops in the Reeb graph that depend on the genus, the number of boundary components, and whether or not the 2-manifold is orientable. We also give an algorithm that constructs the Reeb graph in time O(n log n), where n is the number of edges in the triangulation used to represent the 2-manifold and the Morse function.},
  author       = {Cole-McLaughlin, Kree and Herbert Edelsbrunner and Harer, John and Natarajan, Vijay and Pascucci, Valerio},
  journal      = {Discrete & Computational Geometry},
  number       = {2},
  pages        = {231 -- 244},
  publisher    = {Springer},
  title        = {{Loops in Reeb graphs of 2-manifolds}},
  doi          = {10.1007/s00454-004-1122-6},
  volume       = {32},
  year         = {2004},
}

@article{3986,
  abstract     = {The motion of a biomolecule greatly depends on the engulfing solution, which is mostly water. Instead of representing individual water molecules, it is desirable to develop implicit solvent models that nevertheless accurately represent the contribution of the solvent interaction to the motion. In such models, hydrophobicity is expressed as a weighted sum of atomic surface areas. The derivatives of these weighted areas contribute to the force that drives the motion. In this paper we give formulas for the weighted and unweighted area derivatives of a molecule modeled as a space-filling diagram made up of balls in motion. Other than the radii and the centers of the balls, the formulas are given in terms of the sizes of circular arcs of the boundary and edges of the power diagram. We also give inclusion-exclusion formulas for these sizes.},
  author       = {Bryant, Robert and Herbert Edelsbrunner and Koehl, Patrice and Levitt, Michael},
  journal      = {Discrete & Computational Geometry},
  number       = {3},
  pages        = {293 -- 308},
  publisher    = {Springer},
  title        = {{The area derivative of a space-filling diagram}},
  doi          = {10.1007/s00454-004-1099-1},
  volume       = {32},
  year         = {2004},
}

@article{3987,
  abstract     = {We consider scientific data sets that describe density functions over three-dimensional geometric domains. Such data sets are often large and coarsened representations are needed for visualization and analysis. Assuming a tetrahedral mesh representation, we construct such representations with a simplification algorithm that combines three goals: the approximation of the function, the preservation of the mesh topology, and the improvement of the mesh quality. The third goal is achieved with a novel extension of the well-known quadric error metric. We perform a number of computational experiments to understand the effect of mesh quality improvement on the density map approximation. In addition, we study the effect of geometric simplification on the topological features of the function by monitoring its critical points.},
  author       = {Natarajan, Vijay and Herbert Edelsbrunner},
  journal      = {IEEE Transactions on Visualization and Computer Graphics},
  number       = {5},
  pages        = {587 -- 597},
  publisher    = {IEEE},
  title        = {{Simplification of three-dimensional density maps}},
  doi          = {10.1109/TVCG.2004.32},
  volume       = {10},
  year         = {2004},
}

@inproceedings{3988,
  abstract     = {We give an algorithm that locally improves the fit between two proteins modeled as space-filling diagrams. The algorithm defines the fit in purely geometric terms and improves by applying a rigid motion to one of the two proteins. Our implementation of the algorithm takes between three and ten seconds and converges with high likelihood to the correct docked configuration, provided it starts at a position away from the correct one by at most 18 degrees of rotation and at most 3.0Angstrom of translation. The speed and convergence radius make this an attractive algorithm to use in combination with a coarse sampling of the six-dimensional space of rigid motions.},
  author       = {Choi, Vicky and Agarwal, Pankaj K and Herbert Edelsbrunner and Rudolph, Johannes},
  pages        = {218 -- 229},
  publisher    = {Springer},
  title        = {{Local search heuristic for rigid protein docking}},
  doi          = {10.1007/978-3-540-30219-3_19},
  volume       = {3240},
  year         = {2004},
}

@inproceedings{3989,
  abstract     = {We introduce local and global comparison measures for a collection of k less than or equal to d real-valued smooth functions on a common d-dimensional Riemannian manifold. For k = d = 2 we relate the measures to the set of critical points of one function restricted to the level sets of the other. The definition of the measures extends to piecewise linear functions for which they ace easy to compute. The computation of the measures forms the centerpiece of a software tool which we use to study scientific datasets.},
  author       = {Herbert Edelsbrunner and Harer, John and Natarajan, Vijay and Pascucci, Valerio},
  pages        = {275 -- 280},
  publisher    = {IEEE},
  title        = {{Local and global comparison of continuous functions}},
  doi          = {10.1109/VISUAL.2004.68},
  year         = {2004},
}