Source author record

Jacek Wesołowski

Jacek Wesołowski 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

9works
6topics
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

9 published item(s)

preprint2022arXiv

About an extension of the Matsumoto-Yor property

If $α,β>0$ are distinct and if $A$ and $B$ 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 $A$ and $B$ 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 $(R_n)$ defined by the recurrence $R_n=Q_n+M_nR_{n-1}$, $n\ge1$, where $R_0$ is arbitrary and $(Q_n,M_n)$, $n\ge1$, are i.i.d. copies of a two-dimensional random vector $(Q,M)$, and $(Q_n,M_n)$ is independent of $R_{n-1}$. It is well known that if $E{\ln}|M|<0$ and $E{\ln^+}|Q|<\infty$, then the sequence $(R_n)$ converges in distribution to a random variable $R$ given by $R\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 $(R_n)$ itself does not converge. We assume that $E{\ln}|M|$ exists but that it is non-negative and we ask if in this situation the sequence $(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 $R$ which are solutions of the distributional equation $R\stackrel{d}{=}Q+MR$, where $(Q,M)$ is independent of $R$ and $|M|\le 1$. Goldie and Grübel showed that the tails of $R$ are no heavier than exponential and that if $Q$ is bounded and $M$ resembles near 1 the uniform distribution, then the tails of $R$ are Poissonian. In this paper, we further investigate the connection between the tails of $R$ and the behavior of $M$ near 1. We focus on the special case when $Q$ is constant and $M$ is nonnegative.

preprint2010arXiv

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

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