Source author record

Jacek Wesolowski

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

ResearcherUnclaimed source record

math.PR math.ST Statistics Theory

Catalog footprint

What is connected

7works

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.

preprint2015arXiv

An eigenproblem approach to optimal equal-precision sample allocation in subpopulations

Allocation of samples in stratified and/or multistage sampling is one of the central issues of sampling theory. In a survey of a population often the constraints for precision of estimators of subpopulations parameters have to be taken care of during the allocation of the sample. Such issues are often solved with mathematical programming procedures. In many situations it is desirable to allocate the sample, in a way which forces the precision of estimates at the subpopulations level to be both: optimal and identical, while the constraints of the total (expected) size of the sample (or samples, in two-stage sampling) are imposed. Here our main concern is related to two-stage sampling schemes. We show that such problem in a wide class of sampling plans has an elegant mathematical and computational solution. This is done due to a suitable definition of the optimization problem, which enables to solve it through a linear algebra setting involving eigenvalues and eigenvectors of matrices defined in terms of some population quantities. As a final result we obtain a very simple and relatively universal method for calculating the subpopulation optimal and equal-precision allocation which is based on one of the most standard algorithms of linear algebra (available e.g. in R software). Theoretical solutions are illustrated through a numerical example based on the Labour Force Survey. Finally, we would like to stress that the method we describe, allows to accommodate quite automatically for different levels of precision priority for subpopulations.

preprint2015arXiv

Regression version of the Matsumoto-Yor type characterization of the gamma and Kummer distributions

In this paper we study a Matsumoto-Yor type property for the gamma and Kummer inde- pendent variables discovered in Koudou and Vallois (2012). We prove that constancy of regressions of U = (1 + 1/(X + Y ))=(1 + 1/X) given V = X + Y and of 1/U given V , where X and Y are indepen- dent and positive random variables, characterizes the gamma and Kummer distributions. This result completes characterizations by independence of U and V obtained, under smoothness assumptions for densities, in Koudou and Vallois (2011, 2012). Since we work with differential equations for the Laplace transforms, no density assumptions are needed.

preprint2014arXiv

A new prior for the discrete DAG models with a restricted set of directions

In this paper, we first develop a new family of conjugate prior distributions for the cell parameters of discrete graphical models Markov with respect to a set P of moral directed acyclic graphs with skeleton a given decomposable graph G. Such families arise when the set of conditional independences between discrete variables is given and can be represented by a decomposable graph and additionally, the direction of certain edges is imposed by the practitioner. This family, which we call the P-Dirichlet, is a generalization of the hyper Dirichlet given in Dawid and Lauritzen (1993): it keeps the strong directed hyper Markov property for every DAG in P but increases the flexibility in the choice of its parameters, i.e. the hyper parameters. Our second contribution is a characterization of the P-Dirichlet, which yields, as a corollary, a characterization of the hyper Dirichlet and a characterization of the Dirichlet also. Like that given by Geiger and Heckerman (1997), our characterization of the Dirichlet is based on local and global independence of the probability parameters but we need not make the assumption of the existence of a positive density function. We use the method of moments for our proofs.

preprint2014arXiv

Exploring recursion for optimal estimators under cascade rotation

We are concerned with optimal linear estimation of means on subsequent occasions under sample rotation where evolution of samples in time is designed through a cascade pattern. It has been known since the seminal paper of Patterson (1950) that when the units are not allowed to return to the sample after leaving it for certain period (there are no gaps in the rotation pattern), one step recursion for optimal estimator holds. However, in some important real surveys, e.g. Current Population Survey in the US or Labour Force Survey in many countries in Europe, units return to the sample after being absent in the sample for several occasions (there are gaps in rotation patterns). In such situations difficulty of the question of the form of the recurrence for optimal estimator increases drastically. This issue has not been resolved yet. Instead alternative sub-optimal approaches were developed, as K-composite estimation (see e.g. Hansen et al. (1955)), AK-composite estimation (see e.g. Gurney and Daly (1965) or time series approach (see e.g. Binder and Hidiroglou (1988)). In the present paper we overcome this long-standing difficulty, that is, we present analytical recursion formulas for the optimal linear estimator of the mean for schemes with gaps in rotation patterns. It is achieved under some technical conditions: ASSUMPTION I and ASSUMPTION II (numerical experiments suggest that these assumptions might be universally satisfied). To attain the goal we develop an algebraic operator approach which allows to reduce the problem of recursion for the optimal linear estimator to two issues: (1) localization of roots (possibly complex) of a polynomial Q_p defined in terms of the rotation pattern (Q_p happens to be conveniently expressed through Chebyshev polynomials of the first kind), (2) rank of a matrix S defined in terms of the rotation pattern and the roots of the polynomial Q_p.

preprint2013arXiv

Asymptotic normality through factorial cumulants and partitions identities

In the paper we develop an approach to asymptotic normality through factorial cumulants. Factorial cumulants arise in the same manner from factorial moments, as do (ordinary) cumulants from (ordinary) moments. Another tool we exploit is a new identity for "moments" of partitions of numbers. The general limiting result is then used to (re-)derive asymptotic normality for several models including classical discrete distributions, occupancy problems in some generalized allocation schemes and two models related to negative multinomial distribution.

preprint2012arXiv

The Lukacs theorem and the Olkin-Baker equation

The Olkin-Baker functional equation is closely related to the celebrated Lukacs characterization of the gamma distribution. Its deeper understanding is essential to settle a challenging question of multivariate extensions of the Lukacs theorem. In this paper, first, we provide a new approach to the additive Olkin-Baker equation which holds almost everywhere on (0,\infinity)^2 (with respect to the Lebesgue measure on R^2) under measurability assumption. Second, this new approach is adapted to the case when unknown functions are allowed to be non-measurable and the complete solution is given in such a general case. Third, the Olkin-Baker equation holding outside of a set from proper linearly invariant ideal of subsets of R^2 is considered.

preprint2011arXiv

Stitching pairs of Levy processes into harnesses

We consider natural exponential families of Levy processes with randomized parameter. Such processes are Markov, and under suitable assumptions, pairs of such processes with shared randomization can be stitched together into a single harness. The stitching consists of deterministic reparametrization of the time for both processes, so that they run on adjacent time intervals, and of the choice of the appropriate law at the boundary. Processes in the Levy-Meixner class have an additional property that they are quadratic harnesses, and in this case stitching constructions produce quadratic harnesses.