Source author record

Tim van de Brug

Tim van de Brug appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.PR Machine Learning Methodology

Catalog footprint

What is connected

5works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2021arXiv

Fast marginal likelihood estimation of penalties for group-adaptive elastic net

Nowadays, clinical research routinely uses omics data, such as gene expression, for predicting clinical outcomes or selecting markers. Additionally, so-called co-data are often available, providing complementary information on the covariates, like p-values from previously published studies or groups of genes corresponding to pathways. Elastic net penalisation is widely used for prediction and covariate selection. Group-adaptive elastic net penalisation learns from co-data to improve the prediction and covariate selection, by penalising important groups of covariates less than other groups. Existing methods are, however, computationally expensive. Here we present a fast method for marginal likelihood estimation of group-adaptive elastic net penalties for generalised linear models. We first derive a low-dimensional representation of the Taylor approximation of the marginal likelihood and its first derivative for group-adaptive ridge penalties, to efficiently estimate these penalties. Then we show by using asymptotic normality of the linear predictors that the marginal likelihood for elastic net models may be approximated well by the marginal likelihood for ridge models. The ridge group penalties are then transformed to elastic net group penalties by using the variance function. The method allows for overlapping groups and unpenalised variables. We demonstrate the method in a model-based simulation study and an application to cancer genomics. The method substantially decreases computation time and outperforms or matches other methods by learning from co-data.

preprint2015arXiv

Random walk loop soups and conformal loop ensembles

The random walk loop soup is a Poissonian ensemble of lattice loops; it has been extensively studied because of its connections to the discrete Gaussian free field, but was originally introduced by Lawler and Trujillo Ferreras as a discrete version of the Brownian loop soup of Lawler and Werner, a conformally invariant Poissonian ensemble of planar loops with deep connections to conformal loop ensembles (CLEs) and the Schramm-Loewner evolution (SLE). Lawler and Trujillo Ferreras showed that, roughly speaking, in the continuum scaling limit, ``large'' lattice loops from the random walk loop soup converge to ``large'' loops from the Brownian loop soup. Their results, however, do not extend to clusters of loops, which are interesting because the connection between Brownian loop soup and CLE goes via cluster boundaries. In this paper, we study the scaling limit of clusters of ``large'' lattice loops, showing that they converge to Brownian loop soup clusters. In particular, our results imply that the collection of outer boundaries of outermost clusters composed of ``large'' lattice loops converges to CLE.

preprint2014arXiv

On central limit theorems in the random connection model

Consider a sequence of Poisson random connection models (X_n,lambda_n,g_n) on R^d, where lambda_n / n^d \to lambda > 0 and g_n(x) = g(nx) for some non-increasing, integrable connection function g. Let I_n(g) be the number of isolated vertices of (X_n,lambda_n,g_n) in some bounded Borel set K, where K has non-empty interior and boundary of Lebesgue measure zero. Roy and Sarkar [Phys. A 318 (2003), no. 1-2, 230-242] claim that (I_n(g) - E I_n(g)) / \sqrt Var I_n(g) converges in distribution to a standard normal random variable. However, their proof has errors. We correct their proof and extend the result to larger components when the connection function g has bounded support.

preprint2012arXiv

Fat fractal percolation and k-fractal percolation

We consider two variations on the Mandelbrot fractal percolation model. In the k-fractal percolation model, the d-dimensional unit cube is divided in N^d equal subcubes, k of which are retained while the others are discarded. The procedure is then iterated inside the retained cubes at all smaller scales. We show that the (properly rescaled) percolation critical value of this model converges to the critical value of ordinary site percolation on a particular d-dimensional lattice as N tends to infinity. This is analogous to the result of Falconer and Grimmett that the critical value for Mandelbrot fractal percolation converges to the critical value of site percolation on the same d-dimensional lattice. In the fat fractal percolation model, subcubes are retained with probability p_n at step n of the construction, where (p_n) is a non-decreasing sequence with \prod p_n > 0. The Lebesgue measure of the limit set is positive a.s. given non-extinction. We prove that either the set of connected components larger than one point has Lebesgue measure zero a.s. or its complement in the limit set has Lebesgue measure zero a.s.

preprint2011arXiv

Stochastic domination and weak convergence of conditioned Bernoulli random vectors

For n>=1 let X_n be a vector of n independent Bernoulli random variables. We assume that X_n consists of M "blocks" such that the Bernoulli random variables in block i have success probability p_i. Here M does not depend on n and the size of each block is essentially linear in n. Let X'_n be a random vector having the conditional distribution of X_n, conditioned on the total number of successes being at least k_n, where k_n is also essentially linear in n. Define Y'_n similarly, but with success probabilities q_i>=p_i. We prove that the law of X'_n converges weakly to a distribution that we can describe precisely. We then prove that sup Pr(X'_n <= Y'_n) converges to a constant, where the supremum is taken over all possible couplings of X'_n and Y'_n. This constant is expressed explicitly in terms of the parameters of the system.