Catalog footprint

What is connected

36works
24topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2022arXiv

Bayesian Semiparametric Hidden Markov Tensor Partition Models for Longitudinal Data with Local Variable Selection

We present a flexible Bayesian semiparametric mixed model for longitudinal data analysis in the presence of potentially high-dimensional categorical covariates. Building on a novel hidden Markov tensor decomposition technique, our proposed method allows the fixed effects components to vary between dependent random partitions of the covariate space at different time points. The mechanism not only allows different sets of covariates to be included in the model at different time points but also allows the selected predictors' influences to vary flexibly over time. Smooth time-varying additive random effects are used to capture subject specific heterogeneity. We establish posterior convergence guarantees for both function estimation and variable selection. We design a Markov chain Monte Carlo algorithm for posterior computation. We evaluate the method's empirical performances through synthetic experiments and demonstrate its practical utility through real world applications.

preprint2022arXiv

DT+GNN: A Fully Explainable Graph Neural Network using Decision Trees

We propose the fully explainable Decision Tree Graph Neural Network (DT+GNN) architecture. In contrast to existing black-box GNNs and post-hoc explanation methods, the reasoning of DT+GNN can be inspected at every step. To achieve this, we first construct a differentiable GNN layer, which uses a categorical state space for nodes and messages. This allows us to convert the trained MLPs in the GNN into decision trees. These trees are pruned using our newly proposed method to ensure they are small and easy to interpret. We can also use the decision trees to compute traditional explanations. We demonstrate on both real-world datasets and synthetic GNN explainability benchmarks that this architecture works as well as traditional GNNs. Furthermore, we leverage the explainability of DT+GNNs to find interesting insights into many of these datasets, with some surprising results. We also provide an interactive web tool to inspect DT+GNN's decision making.

preprint2022arXiv

Sound Automation of Magic Wands (extended version)

The magic wand $\mathbin{-\!\!*}$ (also called separating implication) is a separation logic connective commonly used to specify properties of partial data structures, for instance during iterative traversals. A footprint of a magic wand formula $A \mathbin{-\!\!*} B$ is a state that, combined with any state in which $A$ holds, yields a state in which $B$ holds. The key challenge of proving a magic wand (also called packaging a wand) is to find such a footprint. Existing package algorithms either have a high annotation overhead or, as we show in this paper, are unsound. We present a formal framework that precisely characterises a wide design space of possible package algorithms applicable to a large class of separation logics. We prove in Isabelle/HOL that our formal framework is sound and complete, and use it to develop a novel package algorithm that offers competitive automation and is sound. Moreover, we present a novel, restricted definition of wands and prove in Isabelle/HOL that it is possible to soundly combine fractions of such wands, which is not the case for arbitrary wands. We have implemented our techniques for the Viper language, and demonstrate that they are effective in practice.

preprint2020arXiv

Bayesian Nonparametric Bivariate Survival Regression for Current Status Data

We consider nonparametric inference for event time distributions based on current status data. We show that in this scenario conventional mixture priors, including the popular Dirichlet process mixture prior, lead to biologically uninterpretable results as they unnaturally skew the probability mass for the event times toward the extremes of the observed data. Simple assumptions on dependent censoring can fix the problem. We then extend the discussion to bivariate current status data with partial ordering of the two outcomes. In addition to dependent censoring, we also exploit some minimal known structure relating the two event times. We design a Markov chain Monte Carlo algorithm for posterior simulation. Applied to a recurrent infection study, the method provides novel insights into how symptoms-related hospital visits are affected by covariates.

preprint2020arXiv

Clustering and Prediction with Variable Dimension Covariates

In many applied fields incomplete covariate vectors are commonly encountered. It is well known that this can be problematic when making inference on model parameters, but its impact on prediction performance is less understood. We develop a method based on covariate dependent partition models that seamlessly handles missing covariates while completely avoiding any type of imputation. The method we develop allows in-sample predictions as well as out-of-sample prediction, even if the missing pattern in the new subjects' incomplete covariate vector was not seen in the training data. Any data type, including categorical or continuous covariates are permitted. In simulation studies the proposed method compares favorably. We illustrate the method in two application examples.

preprint2020arXiv

Stability of the enhanced area law of the entanglement entropy

We consider a multi-dimensional continuum Schrödinger operator which is given by a perturbation of the negative Laplacian by a compactly supported potential. We establish both an upper and a lower bound on the bipartite entanglement entropy of the ground state of the corresponding quasi-free Fermi gas. The bounds prove that the scaling behaviour of the entanglement entropy remains a logarithmically enhanced area law as in the unperturbed case of the free Fermi gas. The central idea for the upper bound is to use a limiting absorption principle for such kinds of Schrödinger operators.

preprint2019arXiv

Perturbations of continuum random Schrödinger operators with applications to Anderson orthogonality and the spectral shift function

We study effects of a bounded and compactly supported perturbation on multi-dimensional continuum random Schrödinger operators in the region of complete localisation. Our main emphasis is on Anderson orthogonality for random Schrödinger operators. Among others, we prove that Anderson orthogonality does occur for Fermi energies in the region of complete localisation with a non-zero probability. This partially confirms recent non-rigorous findings [V. Khemani et al., Nature Phys. 11, 560-565 (2015)]. The spectral shift function plays an important role in our analysis of Anderson orthogonality. We identify it with the index of the corresponding pair of spectral projections and explore the consequences thereof. All our results rely on the main technical estimate of this paper which guarantees separate exponential decay of the disorder-averaged Schatten $p$-norm of $χ_{a}(f(H) - f(H^τ)) χ_{b}$ in $a$ and $b$. Here, $H^τ$ is a perturbation of the random Schrödinger operator $H$, $χ_{a}$ is the multiplication operator corresponding to the indicator function of a unit cube centred about $a\in\mathbb{R}^{d}$, and $f$ is in a suitable class of functions of bounded variation with distributional derivative supported in the region of complete localisation for $H$.

preprint2016arXiv

Automatic Verification of Iterated Separating Conjunctions using Symbolic Execution

In permission logics such as separation logic, the iterated separating conjunction is a quantifier denoting access permission to an unbounded set of heap locations. In contrast to recursive predicates, iterated separating conjunctions do not prescribe a structure on the locations they range over, and so do not restrict how to traverse and modify these locations. This flexibility is important for the verification of random-access data structures such as arrays and data structures that can be traversed in multiple ways such as graphs. Despite its usefulness, no automatic program verifier natively supports iterated separating conjunctions; they are especially difficult to incorporate into symbolic execution engines, the prevalent technique for building verifiers for these logics. In this paper, we present the first symbolic execution technique to support general iterated separating conjunctions. We propose a novel representation of symbolic heaps and flexible support for logical specifications that quantify over heap locations. Our technique exhibits predictable and fast performance despite employing quantifiers at the SMT level, by carefully controlling quantifier instantiations. It is compatible with other features of permission logics such as fractional permissions, recursive predicates, and abstraction functions. Our technique is implemented as an extension of the Viper verification infrastructure.

preprint2016arXiv

The IDE as a Scriptable Information System (extended version)

Software engineering is extremely information-intensive. Every day developers work with source code, version repositories, issue trackers, documentation, web-based and other information resources. However, three key aspects of information work lack good support: (i) combining information from different sources; (ii) flexibly presenting collected information to enable easier comprehension; and (iii) automatically acting on collected information, for example to perform a refactoring. Poor support for these activities makes many common development tasks time-consuming and error-prone. We propose an approach that directly addresses these three issues by integrating a flexible query mechanism into the development environment. Our approach enables diverse ways to process and visualize information and can be extended via scripts. We demonstrate how an implementation of the approach can be used to rapidly write queries that meet a wide range of information needs.

preprint2015arXiv

A Bayesian feature allocation model for tumor heterogeneity

We develop a feature allocation model for inference on genetic tumor variation using next-generation sequencing data. Specifically, we record single nucleotide variants (SNVs) based on short reads mapped to human reference genome and characterize tumor heterogeneity by latent haplotypes defined as a scaffold of SNVs on the same homologous genome. For multiple samples from a single tumor, assuming that each sample is composed of some sample-specific proportions of these haplotypes, we then fit the observed variant allele fractions of SNVs for each sample and estimate the proportions of haplotypes. Varying proportions of haplotypes across samples is evidence of tumor heterogeneity since it implies varying composition of cell subpopulations. Taking a Bayesian perspective, we proceed with a prior probability model for all relevant unknown quantities, including, in particular, a prior probability model on the binary indicators that characterize the latent haplotypes. Such prior models are known as feature allocation models. Specifically, we define a simplified version of the Indian buffet process, one of the most traditional feature allocation models. The proposed model allows overlapping clustering of SNVs in defining latent haplotypes, which reflects the evolutionary process of subclonal expansion in tumor samples.

preprint2015arXiv

Ergodicity and dynamical localization for Delone-Anderson operators

We study the ergodic properties of Delone-Anderson operators, using the framework of randomly coloured Delone sets and Delone dynamical systems. In particular, we show the existence of the integrated density of states and, under some assumptions on the geometric complexity of the underlying Delone sets, we obtain information on the almost-sure spectrum of the family of random operators. We then exploit these results to study the Lifshitz-tail behaviour of the integrated density of states of a Delone-Anderson operator at the bottom of the spectrum. Furthermore, we use Lifshitz-tail estimates as an input for the multi-scale analysis to prove dynamical localization.

preprint2014arXiv

Ac-conductivity and electromagnetic energy absorption for the Anderson model in linear response theory

We continue our study of the ac-conductivity in linear response theory for the Anderson model using the conductivity measure. We establish further properties of the conductivity measure, including nontriviality at nonzero temperature, the high temperature limit, and asymptotics with respect to the disorder. We also calculate the electromagnetic energy absorption in linear response theory in terms of the conductivity measure.

preprint2014arXiv

Anderson's orthogonality catastrophe

We give an upper bound on the modulus of the ground-state overlap of two non-interacting fermionic quantum systems with $N$ particles in a large but finite volume $L^d$ of $d$-dimensional Euclidean space. The underlying one-particle Hamiltonians of the two systems are standard Schrödinger operators that differ by a non-negative compactly supported scalar potential. In the thermodynamic limit, the bound exhibits an asymptotic power-law decay in the system size $L$, showing that the ground-state overlap vanishes for macroscopic systems. The decay exponent can be interpreted in terms of the total scattering cross section averaged over all incident directions. The result confirms and generalises P. W. Anderson's informal computation [Phys. Rev. Lett. 18, 1049--1051 (1967)].

preprint2014arXiv

Characterization of an in-vacuum PILATUS 1M detector

A dedicated in-vacuum X-ray detector based on the hybrid pixel PILATUS 1M detector has been installed at the four-crystal monochromator beamline of PTB at the electron storage ring BESSY II in Berlin. Due to its windowless operation, the detector can be used in the entire photon energy range of the beamline from 10 keV down to 1.75 keV for small-angle X-ray scattering (SAXS) experiments and anomalous SAXS (ASAXS) at absorption edges of light elements. The radiometric and geometric properties of the detector like quantum efficiency, pixel pitch and module alignment have been determined with low uncertainties. The first grazing incidence SAXS (GISAXS) results demonstrate the superior resolution in momentum transfer achievable at low photon energies.

preprint2014arXiv

Subgroup-Based Adaptive (SUBA) Designs for Multi-Arm Biomarker Trials

Targeted therapies based on biomarker profiling are becoming a mainstream direction of cancer research and treatment. Depending on the expression of specific prognostic biomarkers, targeted therapies assign different cancer drugs to subgroups of patients even if they are diagnosed with the same type of cancer by traditional means, such as tumor location. For example, Herceptin is only indicated for the subgroup of patients with HER2+ breast cancer, but not other types of breast cancer. However, subgroups like HER2+ breast cancer with effective targeted therapies are rare and most cancer drugs are still being applied to large patient populations that include many patients who might not respond or benefit. Also, the response to targeted agents in human is usually unpredictable. To address these issues, we propose SUBA, subgroup-based adaptive designs that simultaneously search for prognostic subgroups and allocate patients adaptively to the best subgroup-specific treatments throughout the course of the trial. The main features of SUBA include the continuous reclassification of patient subgroups based on a random partition model and the adaptive allocation of patients to the best treatment arm based on posterior predictive probabilities. We compare the SUBA design with three alternative designs including equal randomization, outcome-adaptive randomization and a design based on a probit regression. In simulation studies we find that SUBA compares favorably against the alternatives.

preprint2013arXiv

Defining Predictive Probability Functions for Species Sampling Models

We review the class of species sampling models (SSM). In particular, we investigate the relation between the exchangeable partition probability function (EPPF) and the predictive probability function (PPF). It is straightforward to define a PPF from an EPPF, but the converse is not necessarily true. In this paper we introduce the notion of putative PPFs and show novel conditions for a putative PPF to define an EPPF. We show that all possible PPFs in a certain class have to define (unnormalized) probabilities for cluster membership that are linear in cluster size. We give a new necessary and sufficient condition for arbitrary putative PPFs to define an EPPF. Finally, we show posterior inference for a large class of SSMs with a PPF that is not linear in cluster size and discuss a numerical method to derive its PPF.

preprint2012arXiv

Bayesian Nonparametric Shrinkage Applied to Cepheid Star Oscillations

Bayesian nonparametric regression with dependent wavelets has dual shrinkage properties: there is shrinkage through a dependent prior put on functional differences, and shrinkage through the setting of most of the wavelet coefficients to zero through Bayesian variable selection methods. The methodology can deal with unequally spaced data and is efficient because of the existence of fast moves in model space for the MCMC computation. The methodology is illustrated on the problem of modeling the oscillations of Cepheid variable stars; these are a class of pulsating variable stars with the useful property that their periods of variability are strongly correlated with their absolute luminosity. Once this relationship has been calibrated, knowledge of the period gives knowledge of the luminosity. This makes these stars useful as "standard candles" for estimating distances in the universe.

preprint2012arXiv

Localization for random block operators

We continue the investigations of Kirsch, Metzger and the second-named author [J. Stat. Phys. 143, 1035--1054 (2011)] on spectral properties of a certain type of random block operators. In particular, we establish an alternative version of a Wegner estimate and an improved result on Lifschitz tails at the internal band edges. Using these ingredients and the bootstrap multi-scale analysis, we also prove dynamical localization in a neighbourhood of the internal band edges.

preprint2012arXiv

Modeling dependent gene expression

In this paper we propose a Bayesian approach for inference about dependence of high throughput gene expression. Our goals are to use prior knowledge about pathways to anchor inference about dependence among genes; to account for this dependence while making inferences about differences in mean expression across phenotypes; and to explore differences in the dependence itself across phenotypes. Useful features of the proposed approach are a model-based parsimonious representation of expression as an ordinal outcome, a novel and flexible representation of prior information on the nature of dependencies, and the use of a coherent probability model over both the structure and strength of the dependencies of interest. We evaluate our approach through simulations and in the analysis of data on expression of genes in the Complement and Coagulation Cascade pathway in ovarian cancer.

preprint2011arXiv

Galois groups of multivariate Tutte polynomials

The multivariate Tutte polynomial $\hat Z_M$ of a matroid $M$ is a generalization of the standard two-variable version, obtained by assigning a separate variable $v_e$ to each element $e$ of the ground set $E$. It encodes the full structure of $M$. Let $\bv = \{v_e\}_{e\in E}$, let $K$ be an arbitrary field, and suppose $M$ is connected. We show that $\hat Z_M$ is irreducible over $K(\bv)$, and give three self-contained proofs that the Galois group of $\hat Z_M$ over $K(\bv)$ is the symmetric group of degree $n$, where $n$ is the rank of $M$. An immediate consequence of this result is that the Galois group of the multivariate Tutte polynomial of any matroid is a direct product of symmetric groups. Finally, we conjecture a similar result for the standard Tutte polynomial of a connected matroid.

preprint2011arXiv

Synthesis and Characterization of Single Crystal Samples of Spin-$1/2$ Kagome Lattice Antiferromagnets in the Zn-Paratacamite Family Zn$_{x}$Cu$_{4-x}$(OH)$_{6}$Cl$_{2}$

The Zn-paratacamite family, Zn$_{x}$Cu$_{4-x}$(OH)$_{6}$Cl$_{2}$ for $x \, \geq$ 0.33, is an ideal system for studying spin-1/2 frustrated magnetism in the form of antiferromagnetic Cu$^{2+}$ kagome planes. Here we report a new synthesis method by which high quality millimeter-sized single crystals of Zn-paratacamite have been produced. These crystals have been characterized by metal analysis, x-ray diffraction, neutron diffraction, and thermodynamic measurements. The $x$ = 1 member of the series displays a magnetic susceptibility that is slightly anisotropic at high temperatures with $χ_{c} \, > \, χ_{ab}$. Neutron and synchrotron x-ray diffraction experiments confirm the quality of these $x$ = 1 single crystals and indicate no obvious structural transition down to temperatures of T=2 K.

preprint2010arXiv

Random Block Operators

We study fundamental spectral properties of random block operators that are common in the physical modelling of mesoscopic disordered systems such as dirty superconductors. Our results include ergodic properties, the location of the spectrum, existence and regularity of the integrated density of states, as well as Lifshits tails. Special attention is paid to the peculiarities arising from the block structure such as the occurrence of a robust gap in the middle of the spectrum. Without randomness in the off-diagonal blocks the density of states typically exhibits an inverse square-root singularity at the edges of the gap. In the presence of randomness we establish a Wegner estimate that is valid at all energies. It implies that the singularities are smeared out by randomness, and the density of states is bounded. We also show Lifshits tails at these band edges. Technically, one has to cope with a non-monotone dependence on the random couplings.

preprint2010arXiv

The spectral shift function for compactly supported perturbations of Schrödinger operators on large bounded domains

We study the asymptotic behavior as L \to \infty of the finite-volume spectral shift function for a positive, compactly-supported perturbation of a Schrödinger operator in d-dimensional Euclidean space, restricted to a cube of side length L with Dirichlet boundary conditions. The size of the support of the perturbation is fixed and independent of L. We prove that the Cesàro mean of finite-volume spectral shift functions remains pointwise bounded along certain sequences L_n \to \infty for Lebesgue-almost every energy. In deriving this result, we give a short proof of the vague convergence of the finite-volume spectral shift functions to the infinite-volume spectral shift function as L \to\infty . Our findings complement earlier results of W. Kirsch [Proc. Amer. Math. Soc. 101, 509 - 512 (1987), Int. Eqns. Op. Th. 12, 383 - 391 (1989)] who gave examples of positive, compactly-supported perturbations of finite-volume Dirichlet Laplacians for which the pointwise limit of the spectral shift function does not exist for any given positive energy. Our methods also provide a new proof of the Birman--Solomyak formula for the spectral shift function that may be used to express the measure given by the infinite-volume spectral shift function directly in terms of the potential.

preprint2010arXiv

Uniform convergence of spectral shift functions

The spectral shift function ξ_{L}(E) for a Schrödinger operator restricted to a finite cube of length L in multi-dimensional Euclidean space, with Dirichlet boundary conditions, counts the number of eigenvalues less than or equal to E \in \RR created by a perturbation potential V. We study the behavior of this function ξ_{L}(E) as L to infinity for the case of a compactly-supported and bounded potential V. After reviewing results of Kirsch [Proc. Amer. Math. Soc. 101, 509-512 (1987)], and our recent pointwise convergence result for the Cesàro mean [Proc. Amer. Math. Soc. 138, 2141-2150 (2010)], we present a new result on the convergence of the energy-averaged spectral shift function that is uniform with respect to the location of the potential V within the finite box.

preprint2007arXiv

A lower bound for the density of states of the lattice Anderson model

We consider the Anderson model on the multi-dimensional cubic lattice and prove a positive lower bound on the density of states under certain conditions. For example, if the random variables are independently and identically distributed and the probability measure has a bounded Lebesgue density with compact support, and if this density is essentially bounded away from zero on its support, then we prove that the density of states is strictly positive for Lebesgue-almost every energy in the deterministic spectrum.

preprint2007arXiv

Spectral asymptotics of the Laplacian on supercritical bond-percolation graphs

We investigate Laplacians on supercritical bond-percolation graphs with different boundary conditions at cluster borders. The integrated density of states of the Dirichlet Laplacian is found to exhibit a Lifshits tail at the lower spectral edge, while that of the Neumann Laplacian shows a van Hove asymptotics, which results from the percolating cluster. At the upper spectral edge, the behaviour is reversed.

preprint2007arXiv

Variational bounds for the shear viscosity of gelling melts

We study shear stress relaxation for a gelling melt of randomly crosslinked, interacting monomers. We derive a lower bound for the static shear viscosity $η$, which implies that it diverges algebraically with a critical exponent $k\ge 2ν-β$. Here, $ν$ and $β$ are the critical exponents of percolation theory for the correlation length and the gel fraction. In particular, the divergence is stronger than in the Rouse model, proving the relevance of excluded-volume interactions for the dynamic critical behaviour at the gel transition. Precisely at the critical point, our exact results imply a Mark-Houwink relation for the shear viscosity of isolated clusters of fixed size.

preprint2006arXiv

Lifshitz tails for spectra of Erdős--Rényi random graphs

We consider the discrete Laplace operator $Δ^{(N)}$ on Erdős--Rényi random graphs with $N$ vertices and edge probability $p/N$. We are interested in the limiting spectral properties of $Δ^{(N)}$ as $N\to\infty$ in the subcritical regime $0<p<1$ where no giant cluster emerges. We prove that in this limit the expectation value of the integrated density of states of $Δ^{(N)}$ exhibits a Lifshitz-tail behavior at the lower spectral edge E=0.

preprint2006arXiv

On Mott's formula for the ac-conductivity in the Anderson model

We study the ac-conductivity in linear response theory in the general framework of ergodic magnetic Schrödinger operators. For the Anderson model, if the Fermi energy lies in the localization regime, we prove that the ac-conductivity is bounded by $C ν^2 (\log \frac 1 ν)^{d+2}$ at small frequencies $ν$. This is to be compared to Mott's formula, which predicts the leading term to be $C ν^2 (\log \frac 1 ν)^{d+1}$.

preprint2006arXiv

Persistence of Anderson localization in Schrödinger operators with decaying random potentials

We show persistence of both Anderson and dynamical localization in Schrödinger operators with non-positive (attractive) random decaying potential. We consider an Anderson-type Schrödinger operator with a non-positive ergodic random potential, and multiply the random potential by a decaying envelope function. If the envelope function decays slower than $|x|^{-2}$ at infinity, we prove that the operator has infinitely many eigenvalues below zero. For envelopes decaying as $|x|^{-α}$ at infinity, we determine the number of bound states below a given energy $E<0$, asymptotically as $α\downarrow 0$. To show that bound states located at the bottom of the spectrum are related to the phenomenon of Anderson localization in the corresponding ergodic model, we prove: (a) these states are exponentially localized with a localization length that is uniform in the decay exponent $α$; (b)~ dynamical localization holds uniformly in $α$.

preprint1999arXiv

Dynamical properties of random Schrödinger operators

We study dynamical properties of random Schrödinger operators $H^{(ω)}$ defined on the Hilbert space $\ell^2(\bbZ^d)$ or $L^2(\bbR^d)$. Building on results from existing multi-scale analyses, we give sufficient conditions on $H^{(ω)}$ to obtain the vanishing of the diffusion exponent $$ σ_{\rm diff}^+ := \limsup_{T\rightarrow\infty } \frac{\log \bbE \left(\la\la\vert X \vert^2\ra\ra_{T,f_I(H^{(ω)})ψ}\right) }{\log T}=0. $$ Here $\bbE$ is the expectation over randomness, $f_{I}$ is any smooth characteristic function of a bounded energy-interval $I$ and $ψ$ is a state vector in the domain of $H^{(ω)}$ with compact spatial support. The quantity $\la\la |X|^2 \ra\ra_{T,φ}$ denotes the Cesaro mean up to time $T$ of the second moment of position $\la |X|^2\ra_{t,φ}$ at times $0\le t\le T$ of an initial state vector $φ$. If the Hilbert space is $\ell^2(\bbZ^d)$, the method of proof can be strengthened to yield dynamical localization. Under weaker assumptions, we also prove a theorem on the absence of diffusion. The results are applied to a randomly perturbed periodic Schrödinger operator on $L^2(\bbR^d)$, to a simple Anderson-type model on the lattice and to a model with a correlated random potential in continuous space.