Source author record

Jacek Wesołowski

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

9works
6topics
4close collaborators

Published work

9 publications

preprint2022arXiv

About an extension of the Matsumoto-Yor property

If α,β>0α,β>0 are distinct and if AA and BB are independent non-degenerate positive random variables such that $$S=\tfrac{1}{B}\,\tfrac{βA+B}{αA+B}\quad \mbox{and}\quad T=\tfrac{1}{A}\,\tfrac{βA+B}{αA+B} $$ are independent, we prove that this happens if and only if the AA and BB have generalized inverse Gaussian distributions with suitable parameters. Essentially, this has already been proved in Bao and Noack (2021) with supplementary hypothesis on existence of smooth densities. The sources of these questions are an observation about independence properties of the exponential Brownian motion due to Matsumoto and Yor (2001) and a recent work of Croydon and Sasada (2000) on random recursion models rooted in the discrete Korteweg - de Vries equation, where the above result was conjectured. We also extend the direct result to random matrices proving that a matrix variate analogue of the above independence property is satisfied by independent matrix-variate GIG variables. The question of characterization of GIG random matrices through this independence property remains open.

preprint2016arXiv

Linearity of regression for weak records, revisited

Since many years characterization of distribution by linearity of regression of non-adjacent weak records E(W_{i+s}|W_i) = β_1 W_i+β_0 for discrete observations has been known to be a difficult question. Lopez- Blazquez (2004) proposed an interesting idea of reducing it to the adjacent case and claimed to have the characterization problem completely solved. We will explain that, unfortunately, there is a major aw in the proof given in that paper. This aw is related to fact that in some situations the operator responsible for reduction of the non-adjacent case to the adjacent one is not injective. The operator is trivially injective when 0<β_1<1. We show that when β_1>=1 the operator is injective when s = 2, 3, 4. Therefore in these cases the method proposed by Lopez-Blazquez is valid. We also show that the operator is not injective when β_1 >=1 and s >= 5. Consequently, in this case the reduction methodology does not work and thus the characterization problem remains open.

preprint2015arXiv

Limit Theorems for Empirical Rényi Entropy and Divergence with Applications to Molecular Diversity Analysis

Quantitative methods for studying biodiversity have been traditionally rooted in the classical theory of finite frequency tables analysis. However, with the help of modern experimental tools, like high throughput sequencing, we now begin to unlock the outstanding diversity of genomic data in plants and animals reflective of the long evolutionary history of our planet. This molecular data often defies the classical frequency/contingency tables assumptions and seems to require sparse tables with very large number of categories and highly unbalanced cell counts, e.g., following heavy tailed distributions (for instance, power laws). Motivated by the molecular diversity studies, we propose here a frequency-based framework for biodiversity analysis in the asymptotic regime where the number of categories grows with sample size (an infinite contingency table). Our approach is rooted in information theory and based on the Gaussian limit results for the effective number of species (the Hill numbers) and the empirical Renyi entropy and divergence. We argue that when applied to molecular biodiversity analysis our methods can properly account for the complicated data frequency patterns on one hand and the practical sample size limitations on the other. We illustrate this principle with two specific RNA sequencing examples: a comparative study of T-cell receptor populations and a validation of some preselected molecular hepatocellular carcinoma (HCC) markers.

preprint2011arXiv

Renorming divergent perpetuities

We consider a sequence of random variables (Rn)(R_n) defined by the recurrence Rn=Qn+MnRn−1R_n=Q_n+M_nR_{n-1}, n≥1n\ge1, where R0R_0 is arbitrary and (Qn,Mn)(Q_n,M_n), n≥1n\ge1, are i.i.d. copies of a two-dimensional random vector (Q,M)(Q,M), and (Qn,Mn)(Q_n,M_n) is independent of Rn−1R_{n-1}. It is well known that if Eln⁡∣M∣<0E{\ln}|M|<0 and Eln⁡+∣Q∣<∞E{\ln^+}|Q|<\infty, then the sequence (Rn)(R_n) converges in distribution to a random variable RR given by R=d∑k=1∞Qk∏j=1k−1MjR\stackrel{d}{=}\sum_{k=1}^{\infty}Q_k\prod_{j=1}^{k-1}M_j, and usually referred to as perpetuity. In this paper we consider a situation in which the sequence (Rn)(R_n) itself does not converge. We assume that Eln⁡∣M∣E{\ln}|M| exists but that it is non-negative and we ask if in this situation the sequence (Rn)(R_n), after suitable normalization, converges in distribution to a non-degenerate limit.

preprint2010arXiv

Free Quadratic Harness

Free quadratic harness is a Markov process from the class of quadratic harnesses, i.e. processes with linear regressions and quadratic conditional variances. The process has recently been constructed for a restricted range of parameters in the paper "Askey-Wilson polynomials, quadratic harnesses and martingales" by W. Bryc and J. Wesołowski using Askey--Wilson polynomials. Here we provide a self-contained construction of the free quadratic harness for all values of parameters.

preprint2010arXiv

Perpetuities with thin tails revisited

We consider the tail behavior of random variables RR which are solutions of the distributional equation R=dQ+MRR\stackrel{d}{=}Q+MR, where (Q,M)(Q,M) is independent of RR and ∣M∣≤1|M|\le 1. Goldie and Grübel showed that the tails of RR are no heavier than exponential and that if QQ is bounded and MM resembles near 1 the uniform distribution, then the tails of RR are Poissonian. In this paper, we further investigate the connection between the tails of RR and the behavior of MM near 1. We focus on the special case when QQ is constant and MM is nonnegative.

preprint2010arXiv

Why Jordan algebras are natural in statistics:quadratic regression implies Wishart distributions

If the space Q\mathcal{Q} of quadratic forms in Rn\mathbb{R}^n is splitted in a direct sum Q1⊕...⊕Qk\mathcal{Q}_1\oplus...\oplus \mathcal{Q}_k and if XX and YY are independent random variables of Rn\mathbb{R}^n, assume that there exist a real number aa such that E(X∣X+Y)=a(X+Y)E(X|X+Y)=a(X+Y) and real distinct numbers b1,...,bkb_1,...,b_k such that E(q(X)∣X+Y)=biq(X+Y)E(q(X)|X+Y)=b_iq(X+Y) for any qq in Qi.\mathcal{Q}_i. We prove that this happens only when k=2k=2, when Rn\mathbb{R}^n can be structured in a Euclidean Jordan algebra and when XX and YY have Wishart distributions corresponding to this structure.

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.