Researcher profile

Gerard Letac

Gerard Letac contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

preprint2013arXiv

Dirichlet random walks

This article provides tools for the study of the Dirichlet random walk in $\mathbb{R}^d$. By this we mean the random variable $W=X_1Θ_1+\cdots+X_nΘ_n$ where $X=(X_1,\ldots,X_n) \sim \mathcal{D}(q_1,\ldots,q_n)$ is Dirichlet distributed and where $Θ_1,\ldots Θ_n$ are iid, uniformly distributed on the unit sphere of $\mathbb{R}^d$ and independent of $X.$ In particular we compute explicitely in a number of cases the distribution of $W.$ Some of our results appear already in the literature, in particular in the papers by Gérard Le Caër (2010, 2011). In these cases, our proofs are much simpler from the original ones, since we use a kind of Stieltjes transform of $W$ instead of the Laplace transform: as a consequence the hypergeometric functions replace the Bessel functions. A crucial ingredient is a particular case of the classical and non trivial identity, true for $0\leq u\leq 1/2$:$$_2F_1(2a,2b;a+b+\frac{1}{2};u)= \_2F_1(a,b;a+b+\frac{1}{2};4u-4u^2).$$ We extend these results to a study of the limits of the Dirichlet random walks when the number of added terms goes to infinity, interpreting the results in terms of an integral by a Dirichlet process. We introduce the ideas of Dirichlet semigroups and of Dirichlet infinite divisibility and characterize these infinite divisible distributions in the sense of Dirichlet when they are concentrated on the unit ball of $\mathbb{R}^d.$ {4mm}\noindent \textsc{Keywords:} Dirichlet processes, Stieltjes transforms, random flight, distributions in a ball, hyperuniformity, infinite divisibility in the sense of Dirichlet. {4mm}\noindent \textsc{AMS classification}: 60D99, 60F99.

preprint2012arXiv

Perpetuity property of the Dirichlet distribution

Let $X$, $B$ and $Y$ be three Dirichlet, Bernoulli and beta independent random variables such that $X\sim \mathcal{D}(a_0,...,a_d),$ such that $\Pr(B=(0,...,0,1,0,...,0))=a_i/a$ with $a=\sum_{i=0}^da_i$ and such that $Y\sim β(1,a).$ We prove that $X\sim X(1-Y)+BY.$ This gives the stationary distribution of a simple Markov chain on a tetrahedron. We also extend this result to the case when $B$ follows a quasi Bernoulli distribution $\mathcal{B}_k(a_0,...,a_d)$ on the tetrahedron and when $Y\sim β(k,a)$. We extend it even more generally to the case where $X$ is a Dirichlet process and $B$ is a quasi Bernoulli random probability. Finally the case where the integer $k$ is replaced by a positive number $c$ is considered when $a_0=...=a_d=1.$ \textsc{Keywords} \textit{Perpetuities, Dirichlet process, Ewens distribution, quasi Bernoulli laws, probabilities on a tetrahedron, $T_c$ transform, stationary distribution.} AMS classification 60J05, 60E99.

preprint2011arXiv

Bayes factors and the geometry of discrete hierarchical loglinear models

A standard tool for model selection in a Bayesian framework is the Bayes factor which compares the marginal likelihood of the data under two given different models. In this paper, we consider the class of hierarchical loglinear models for discrete data given under the form of a contingency table with multinomial sampling. We assume that the Diaconis-Ylvisaker conjugate prior is the prior distribution on the loglinear parameters and the uniform is the prior distribution on the space of models. Under these conditions, the Bayes factor between two models is a function of their prior and posterior normalizing constants. These constants are functions of the hyperparameters $(m,α)$ which can be interpreted respectively as marginal counts and the total count of a fictive contingency table. We study the behaviour of the Bayes factor when $α$ tends to zero. In this study two mathematical objects play a most important role. They are, first, the interior $C$ of the convex hull $\bar{C}$ of the support of the multinomial distribution for a given hierarchical loglinear model together with its faces and second, the characteristic function $\mathbb{J}_C$ of this convex set $C$. We show that, when $α$ tends to 0, if the data lies on a face $F_i$ of $\bar{C_i},i=1,2$ of dimension $k_i$, the Bayes factor behaves like $α^{k_1-k_2}$. This implies in particular that when the data is in $C_1$ and in $C_2$, i.e. when $k_i$ equals the dimension of model $J_i$, the sparser model is favored, thus confirming the idea of Bayesian regularization.

preprint2010arXiv

Increasing hazard rate of mixtures for natural exponential families

Hazard rates play an important role in various areas, e.g., reliability theory, survival analysis, biostatistics, queueing theory and actuarial studies. Mixtures of distributions are also of a great preeminence in such areas as most populations of components are indeed heterogeneous. In this study we present a sufficient condition for mixtures of two elements\ of the same natural exponential family (NEF) to have an increasing hazard rate. We then apply this condition to some classical NEF's having either quadratic, or cubic variance functions (VF) and others as well. A particular attention is devoted to the hyperbolic cosine NEF having a quadratic VF, the Ressel NEF having a cubic VF and to the Kummer distributions of type 2 NEF. The application of such a sufficient condition is quite intricate and cumbersome, in particular when applied to the latter three NEF's. Various lemmas and propositions are needed then to verify this condition for these NEF's.

preprint2010arXiv

The limiting behavior of some infinitely divisible exponential dispersion models

Consider an exponential dispersion model (EDM) generated by a probability $ μ$ on $[0,\infty )$ which is infinitely divisible with an unbounded Lévy measure $ν$. The Jorgensen set (i.e., the dispersion parameter space) is then $\mathbb{R}^{+}$, in which case the EDM is characterized by two parameters: $θ_{0}$ the natural parameter of the associated natural exponential family and the Jorgensen (or dispersion) parameter $t$. Denote by $EDM(θ_{0},t)$ the corresponding distribution and let $Y_{t}$ is a r.v. with distribution $EDM(θ_0,t)$. Then if $ν((x,\infty ))\sim -\ell \log x$ around zero we prove that the limiting law $F_0$ of $ Y_{t}^{-t}$ as $t\rightarrow 0$ is of a Pareto type (not depending on $ θ_0$) with the form $F_0(u)=0$ for $u<1$ and $1-u^{-\ell }$ for $ u\geq 1$. Such a result enables an approximation of the distribution of $ Y_{t}$ for relatively small values of the dispersion parameter of the corresponding EDM. Illustrative examples are provided.