Source author record

Lasse Leskelä

Lasse Leskelä appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
14topics
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

15 published item(s)

preprint2022arXiv

Community recovery in non-binary and temporal stochastic block models

This article studies the estimation of latent community memberships from pairwise interactions in a network of $N$ nodes, where the observed interactions can be of arbitrary type, including binary, categorical, and vector-valued, and not excluding even more general objects such as time series or spatial point patterns. As a generative model for such data, we introduce a stochastic block model with a general measurable interaction space $\mathcal S$, for which we derive information-theoretic bounds for the minimum achievable error rate. These bounds yield sharp criteria for the existence of consistent and strongly consistent estimators in terms of data sparsity, statistical similarity between intra- and inter-block interaction distributions, and the shape and size of the interaction space. The general framework makes it possible to study temporal and multiplex networks with $\mathcal S = \{0,1\}^T$, in settings where both $N \to \infty$ and $T \to \infty$, and the temporal interaction patterns are correlated over time. For temporal Markov interactions, we derive sharp consistency thresholds. We also present fast online estimation algorithms which fully utilise the non-binary nature of the observed data. Numerical experiments on synthetic and real data show that these algorithms rapidly produce accurate estimates even for very sparse data arrays.

preprint2022arXiv

Consistent Bayesian community recovery in multilayer networks

Revealing underlying relations between nodes in a network is one of the most important tasks in network analysis. Using tools and techniques from a variety of disciplines, many community recovery methods have been developed for different scenarios. Despite the recent interest on community recovery in multilayer networks, theoretical results on the accuracy of the estimates are few and far between. Given a multilayer, e.g. temporal, network and a multilayer stochastic block model, we derive bounds for sufficient separation between intra- and inter-block connectivity parameters to achieve posterior exact and almost exact community recovery. These conditions are comparable to a well known threshold for community detection by a single-layer stochastic block model. A simulation study shows that the derived bounds translate to classification accuracy that improves as the number of observed layers increases.

preprint2021arXiv

Adaptive and optimized COVID-19 vaccination strategies across geographical regions and age groups

We evaluate the efficiency of various heuristic strategies for allocating vaccines against COVID-19 and compare them to strategies found using optimal control theory. Our approach is based on a mathematical model which tracks the spread of disease among different age groups and across different geographical regions, and we introduce a method to combine age-specific contact data to geographical movement data. As a case study, we model the epidemic in the population of mainland Finland utilizing mobility data from a major telecom operator. Our approach allows to determine which geographical regions and age groups should be targeted first in order to minimize the number of deaths. In the scenarios that we test, we find that distributing vaccines demographically and in an age-descending order is not optimal for minimizing deaths and the burden of disease. Instead, more lives could potentially be saved by using strategies which emphasize high-incidence regions and distribute vaccines in parallel to multiple age groups. The level of emphasis that high-incidence regions should be given depends on the overall transmission rate in the population. This observation highlights the importance of updating the vaccination strategy when the effective reproduction number changes due to the general contact patterns changing and new virus variants entering.

preprint2020arXiv

Towards analyzing large graphs with quantum annealing and quantum gate computers

The use of quantum computing in graph community detection and regularity checking related to Szemeredi's Regularity Lemma (SRL) are demonstrated with D-Wave Systems' quantum annealer and simulations. We demonstrate the capability of quantum computing in solving hard problems relevant to big data. A new community detection algorithm based on SRL is also introduced and tested. In worst case scenario of regularity check we use Grover's algorithm and quantum phase estimation algorithm, in order to speed-up computations using a quantum gate computers.

preprint2016arXiv

Diclique clustering in a directed random graph

We discuss a notion of clustering for directed graphs, which describes how likely two followers of a node are to follow a common target. The associated network motifs, called dicliques or bi-fans, have been found to be key structural components in various real-world networks. We introduce a two-mode statistical network model consisting of actors and auxiliary attributes, where an actor i decides to follow an actor j whenever i demands an attribute supplied by j. We show that the digraph admits nontrivial clustering properties of the aforementioned type, as well as power-law indegree and outdegree distributions.

preprint2015arXiv

Geometric juggling with q-analogues

We derive a combinatorial equilibrium for bounded juggling patterns with a random, $q$-geometric throw distribution. The dynamics are analyzed via rook placements on staircase Ferrers boards, which leads to a steady-state distribution containing $q$-rook polynomial coefficients and $q$-Stirling numbers of the second kind. We show that the equilibrium probabilities of the bounded model can be uniformly approximated with the equilibrium probabilities of a corresponding unbounded model. This observation leads to new limit formulae for $q$-analogues. Keywords: juggling pattern; $q$-Stirling number of the second kind; Ferrers board; Markov process; combinatorial equilibrium

preprint2014arXiv

Flow coupling and stochastic ordering of throughputs in linear networks

Robust estimates for the performance of complicated queueing networks can be obtained by showing that the number of jobs in the network is stochastically comparable to a simpler, analytically tractable reference network. Classical coupling results on stochastic ordering of network populations require strong monotonicity assumptions which are often violated in practice. However, in most real-world applications we care more about what goes through a network than what sits inside it. This paper describes a new approach for ordering flows instead of populations by augmenting network states with their associated flow counting processes and deriving Markov couplings of the augmented state-flow processes.

preprint2012arXiv

Hard-core thinnings of germ-grain models with power-law grain sizes

Random sets with long-range dependence can be generated using a Boolean model with power-law grain sizes. We study thinnings of such Boolean models which have the hard-core property that no grains overlap in the resulting germ-grain model. A fundamental question is whether long-range dependence is preserved under such thinnings. To answer this question we study four natural thinnings of a Poisson germ-grain model where the grains are spheres with a regularly varying size distribution. We show that a thinning which favors large grains preserves the slow correlation decay of the original model, whereas a thinning which favors small grains does not. Our most interesting finding concerns the case where only disjoint grains are retained, which corresponds to the well-known Matérn type I thinning. In the resulting germ-grain model, typical grains have exponentially small sizes, but rather surprisingly, the long-range dependence property is still present. As a byproduct, we obtain new mechanisms for generating homogeneous and isotropic random point configurations having a power-law correlation decay.

preprint2012arXiv

Juggler's exclusion process

Juggler's exclusion process describes a system of particles on the positive integers where particles drift down to zero at unit speed. After a particle hits zero, it jumps into a randomly chosen unoccupied site. We model the system as a set-valued Markov process and show that the process is ergodic if the family of jump height distributions is uniformly integrable. In a special case where the particles jump according to a set-avoiding memoryless distribution, the process reaches its equilibrium in finite nonrandom time, and the equilibrium distribution can be represented as a Gibbs measure conforming to a linear gravitational potential.

preprint2011arXiv

A local limit theorem for a transient chaotic walk in a frozen environment

This paper studies particle propagation in a one-dimensional inhomogeneous medium where the laws of motion are generated by chaotic and deterministic local maps. Assuming that the particle's initial location is random and uniformly distributed, this dynamical system can be reduced to a random walk in a one-dimensional inhomogeneous environment with a forbidden direction. Our main result is a local limit theorem which explains in detail why, in the long run, the random walk's probability mass function does not converge to a Gaussian density, although the corresponding limiting distribution over a coarser diffusive space scale is Gaussian.

preprint2011arXiv

Stochastic order characterization of uniform integrability and tightness

We show that a family of random variables is uniformly integrable if and only if it is stochastically bounded in the increasing convex order by an integrable random variable. This result is complemented by proving analogous statements for the strong stochastic order and for power-integrable dominating random variables. Especially, we show that whenever a family of random variables is stochastically bounded by a p-integrable random variable for some p>1, there is no distinction between the strong order and the increasing convex order. These results also yield new characterizations of relative compactness in Wasserstein and Prohorov metrics.

preprint2010arXiv

A Dilution Test for the Convergence of Subseries of a Monotone Series

Cauchy's condensation test allows to determine the convergence of a monotone series by looking at a weighted subseries that only involves terms of the original series indexed by the powers of two. It is natural to ask whether the converse is also true: Is it possible to determine the convergence of an arbitrary subseries of a monotone series by looking at a suitably weighted version of the original series? In this note we show that the answer is affirmative and introduce a new convergence test particularly designed for this purpose.

preprint2010arXiv

Stability of a spatial polling system with greedy myopic service

This paper studies a spatial queueing system on a circle, polled at random locations by a myopic server that can only observe customers in a bounded neighborhood. The server operates according to a greedy policy, always serving the nearest customer in its neighborhood, and leaving the system unchanged at polling instants where the neighborhood is empty. This system is modeled as a measure-valued random process, which is shown to be positive recurrent under a natural stability condition that does not depend on the server's scan range. When the interpolling times are light-tailed, the stable system is shown to be geometrically ergodic. The steady-state behavior of the system is briefly discussed using numerical simulations and a heuristic light-traffic approximation.

preprint2010arXiv

Stability of parallel queueing systems with coupled service rates

This paper considers a parallel system of queues fed by independent arrival streams, where the service rate of each queue depends on the number of customers in all of the queues. Necessary and sufficient conditions for the stability of the system are derived, based on stochastic monotonicity and marginal drift properties of multiclass birth and death processes. These conditions yield a sharp characterization of stability for systems, where the service rate of each queue is decreasing in the number of customers in other queues, and has uniform limits as the queue lengths tend to infinity. The results are illustrated with applications where the stability region may be nonconvex.

preprint2007arXiv

Scaling limits for random fields with long-range dependence

This paper studies the limits of a spatial random field generated by uniformly scattered random sets, as the density $λ$ of the sets grows to infinity and the mean volume $ρ$ of the sets tends to zero. Assuming that the volume distribution has a regularly varying tail with infinite variance, we show that the centered and renormalized random field can have three different limits, depending on the relative speed at which $λ$ and $ρ$ are scaled. If $λ$ grows much faster than $ρ$ shrinks, the limit is Gaussian with long-range dependence, while in the opposite case, the limit is independently scattered with infinite second moments. In a special intermediate scaling regime, there exists a nontrivial limiting random field that is not stable.