Search Results

Now showing 1 - 10 of 1319
  • Item
    Assessment of Stability in Partitional Clustering Using Resampling Techniques
    (Karlsruhe : KIT Scientific Publishing, 2016) Mucha, Hans-Joachim
    The assessment of stability in cluster analysis is strongly related to the main difficult problem of determining the number of clusters present in the data. The latter is subject of many investigations and papers considering different resampling techniques as practical tools. In this paper, we consider non-parametric resampling from the empirical distribution of a given dataset in order to investigate the stability of results of partitional clustering. In detail, we investigate here only the very popular K-means method. The estimation of the sampling distribution of the adjusted Rand index (ARI) and the averaged Jaccard index seems to be the most general way to do this. In addition, we compare bootstrapping with different subsampling schemes (i.e., with different cardinality of the drawn samples) with respect to their performance in finding the true number of clusters for both synthetic and real data.
  • Item
    Ultrashort optical pulse propagation in terms of analytic signal
    (New York, NY : Hindawi, 2011) Amiranashvili, Sh.; Demircan, A.
    We demonstrate that ultrashort optical pulses propagating in a nonlinear dispersive medium are naturally described through incorporation of analytic signal for the electric field. To this end a second-order nonlinear wave equation is first simplified using a unidirectional approximation. Then the analytic signal is introduced, and all nonresonant nonlinear terms are eliminated. The derived propagation equation accounts for arbitrary dispersion, resonant four-wave mixing processes, weak absorption, and arbitrary pulse duration. The model applies to the complex electric field and is independent of the slowly varying envelope approximation. Still the derived propagation equation posses universal structure of the generalized nonlinear Schrdinger equation (NSE). In particular, it can be solved numerically with only small changes of the standard split-step solver or more complicated spectral algorithms for NSE. We present exemplary numerical solutions describing supercontinuum generation with an ultrashort optical pulse.
  • Item
    A boundary control problem for the pure Cahn–Hilliard equation with dynamic boundary conditions
    (Berlin ; Boston, Mass. : de Gruyter, 2015) Colli, Pierluigi; Gilardi, Gianni; Sprekels, Jürgen
    A boundary control problem for the pure Cahn–Hilliard equations with possibly singular potentialsand dynamic boundary conditions is studied and rst-order necessary conditions for optimality are proved.
  • Item
    Simulation of microwave circuits and laser structures including PML by means of FIT
    (München : European Geopyhsical Union, 2004) Hebermehl, G.; Schefter, J.; Schlundt, R.; Tischler, Th.; Zscheile, H.; Heinrich, W.
    Field-oriented methods which describe the physical properties of microwave circuits and optical structures are an indispensable tool to avoid costly and time-consuming redesign cycles. Commonly the electromagnetic characteristics of the structures are described by the scattering matrix which is extracted from the orthogonal decomposition of the electric field. The electric field is the solution of an eigenvalue and a boundary value problem for Maxwell’s equations in the frequency domain. We discretize the equations with staggered orthogonal grids using the Finite Integration Technique (FIT). Maxwellian grid equations are formulated for staggered nonequidistant rectangular grids and for tetrahedral nets with corresponding dual Voronoi cells. The interesting modes of smallest attenuation are found solving a sequence of eigenvalue problems of modified matrices. To reduce the execution time for high-dimensional problems a coarse and a fine grid is used. The calculations are carried out, using two levels of parallelization. The discretized boundary value problem, a large-scale system of linear algebraic equations with different right-hand sides, is solved by a block Krylov subspace method with various preconditioning techniques. Special attention is paid to the Perfectly Matched Layer boundary condition (PML) which causes non physical modes and a significantly increased number of iterations in the iterative methods.
  • Item
    Distributed optimal control of a nonstandard nonlocal phase field system
    (Springfield, MO : AIMS Press, 2016) Colli, Pierluigi; Gilardi, Gianni; Sprekels, Jürgen
    We investigate a distributed optimal control problem for a nonlocal phase field model of viscous Cahn-Hilliard type. The model constitutes a nonlocal version of a model for two-species phase segregation on an atomic lattice under the presence of diffusion that has been studied in a series of papers by P. Podio-Guidugli and the present authors. The model consists of a highly nonlinear parabolic equation coupled to an ordinary differential equation. The latter equation contains both nonlocal and singular terms that render the analysis difficult. Standard arguments of optimal control theory do not apply directly, although the control constraints and the cost functional are of standard type. We show that the problem admits a solution, and we derive the first-order necessary conditions of optimality.
  • Item
    Large Deviations of Continuous Regular Conditional Probabilities
    (New York, NY [u.a.] : Springer Science + Business Media B.V., 2016) van Zuijlen, W.
    We study product regular conditional probabilities under measures of two coordinates with respect to the second coordinate that are weakly continuous on the support of the marginal of the second coordinate. Assuming that there exists a sequence of probability measures on the product space that satisfies a large deviation principle, we present necessary and sufficient conditions for the conditional probabilities under these measures to satisfy a large deviation principle. The arguments of these conditional probabilities are assumed to converge. A way to view regular conditional probabilities as a special case of product regular conditional probabilities is presented. This is used to derive conditions for large deviations of regular conditional probabilities. In addition, we derive a Sanov-type theorem for large deviations of the empirical distribution of the first coordinate conditioned on fixing the empirical distribution of the second coordinate.
  • Item
    Adaptive smoothing of digital images: The R package adimpro
    (Los Angeles, Calif. : UCLA, Dept. of Statistics, 2007) Polzehl, J.; Tabelow, K.
    Digital imaging has become omnipresent in the past years with a bulk of applications ranging from medical imaging to photography. When pushing the limits of resolution and sensitivity noise has ever been a major issue. However, commonly used non-adaptive filters can do noise reduction at the cost of a reduced effective spatial resolution only. Here we present a new package adimpro for R, which implements the propagationseparation approach by (Polzehl arid Spokoiriy 2006) for smoothing digital images. This method naturally adapts to different structures of different size in the image and thus avoids oversmoothing edges and fine structures. We extend the method for imaging data with spatial correlation. Furthermore we show how the estimation of the dependence between variance and mean value can be included. We illustrate the use of the package through some examples.
  • Item
    Statistical parametric maps for functional MRI experiments in R: The package fmri
    (Los Angeles : UCLA, 2011) Tabelow, K.; Polzehl, J.
    The purpose of the package fmri is the analysis of single subject functional magnetic resonance imaging (fMRI) data. It provides fMRI analysis from time series modeling by a linear model to signal detection and publication quality images. Specifically, it implements structural adaptive smoothing methods with signal detection for adaptive noise reduction which avoids blurring of activation areas. Within this paper we describe the complete pipeline for fMRI analysis using fmri. We describe data reading from various medical imaging formats and the linear modeling used to create the statistical parametric maps. We review the rationale behind the structural adaptive smoothing algorithms and explain their usage from the package fmri. We demonstrate the results of such analysis using two experimental datasets. Finally, we report on the usage of a graphical user interface for some of the package functions.
  • Item
    Hausdorff metric BV discontinuity of sweeping processes
    (Bristol : IOP Publ., 2016) Klein, Olaf; Recupero, Vincenzo
    Sweeping processes are a class of evolution differential inclusions arising in elastoplasticity and were introduced by J.J. Moreau in the early seventies. The solution operator of the sweeping processes represents a relevant example of rate independent operator. As a particular case we get the so called play operator, which is a typical example of a hysteresis operator. The continuity properties of these operators were studied in several works. In this note we address the continuity with respect to the strict metric in the space of functions of bounded variation with values in the metric space of closed convex subsets of a Hilbert space. We provide counterexamples showing that for all BV-formulations of the sweeping process the corresponding solution operator is not continuous when its domain is endowed with the strict topology of BV and its codomain is endowed with the L1-topology. This is at variance with the play operator which has a BV-extension that is continuous in this case.
  • Item
    Influence of cell shape, inhomogeneities and diffusion barriers in cell polarization models
    (Philadelphia, Pa. : IOP Publ., 2015) Giese, Wolfgang; Eigel, Martin; Westerheide, Sebastian; Engwer, Christian; Klipp, Edda
    In silico experiments bear the potential for further understanding of biological transport processes by allowing a systematic modification of any spatial property and providing immediate simulation results. Cell polarization and spatial reorganization of membrane proteins are fundamental for cell division, chemotaxis and morphogenesis. We chose the yeast Saccharomyces cerevisiae as an exemplary model system which entails the shuttling of small Rho GTPases such as Cdc42 and Rho, between an active membrane-bound form and an inactive cytosolic form. We used partial differential equations to describe the membrane-cytosol shuttling of proteins. In this study, a consistent extension of a class of 1D reaction-diffusion systems into higher space dimensions is suggested. The membrane is modeled as a thin layer to allow for lateral diffusion and the cytosol is modeled as an enclosed volume. Two well-known polarization mechanisms were considered. One shows the classical Turing-instability patterns, the other exhibits wave-pinning dynamics. For both models, we investigated how cell shape and diffusion barriers like septin structures or bud scars influence the formation of signaling molecule clusters and subsequent polarization. An extensive set of in silico experiments with different modeling hypotheses illustrated the dependence of cell polarization models on local membrane curvature, cell size and inhomogeneities on the membrane and in the cytosol. In particular, the results of our computer simulations suggested that for both mechanisms, local diffusion barriers on the membrane facilitate Rho GTPase aggregation, while diffusion barriers in the cytosol and cell protrusions limit spontaneous molecule aggregations of active Rho GTPase locally.