-
Multivariate Quantiles: Geometric and Measure-Transportation-Based Contours
Authors:
Marc Hallin,
Dimitri Konen
Abstract:
Quantiles are a fundamental concept in probability and theoretical statistics and a daily tool in their applications. While the univariate concept of quantiles is quite clear and well understood, its multivariate extension is more problematic. After half a century of continued efforts and many proposals, two concepts, essentially, are emerging: the so-called (relabeled) geometric quantiles, extend…
▽ More
Quantiles are a fundamental concept in probability and theoretical statistics and a daily tool in their applications. While the univariate concept of quantiles is quite clear and well understood, its multivariate extension is more problematic. After half a century of continued efforts and many proposals, two concepts, essentially, are emerging: the so-called (relabeled) geometric quantiles, extending the characterization of univariate quantiles as minimizers of an L1 loss function involving the check functions, and the more recent center-outward quantiles based on measure transportation ideas. These two concepts yield distinct families of quantile regions and quantile contours. Our objective here is to present a comparison of their main theoretical properties and a numerical investigation of their differences.
△ Less
Submitted 4 January, 2024;
originally announced January 2024.
-
Output-sensitive Complexity of Multi-Objective Integer Network Flow Problems
Authors:
David Könen,
Michael Stiglmayr
Abstract:
This paper addresses the output-sensitive complexity for linear multi-objective integer minimum cost flow (MOIMCF) problems and provides insights about the time complexity for enumerating all supported nondominated vectors. The paper shows that there can not exist an output-polynomial time algorithm for the enumeration of all supported nondominated vectors that determine the vectors in an ordered…
▽ More
This paper addresses the output-sensitive complexity for linear multi-objective integer minimum cost flow (MOIMCF) problems and provides insights about the time complexity for enumerating all supported nondominated vectors. The paper shows that there can not exist an output-polynomial time algorithm for the enumeration of all supported nondominated vectors that determine the vectors in an ordered way in the outcome space unless NP = P. Moreover, novel methods for identifying supported nondominated vectors in bi-objective minimum cost flow (BOIMCF) problems are proposed, accompanied by a numerical comparison between decision- and objective-space methods. A novel, equivalent and more compact formulation of the minimum cost flow ILP formulation used in the e-constrained-scalarization approach is introduced, demonstrating enhanced efficiency in the numerical tests
△ Less
Submitted 4 December, 2023;
originally announced December 2023.
-
An output-polynomial time algorithm to determine all supported efficient solutions for multi-objective integer network flow problems
Authors:
David Könen,
Michael Stiglmayr
Abstract:
This paper addresses the problem of enumerating all supported efficient solutions for a linear multi-objective integer minimum cost flow problem (MOIMCF).
First, we highlight an inconsistency in various definitions of supported nondominated vectors for multi-objective integer linear programs (MOILP). Several characterizations for supported nondominated vectors/efficient solutions are used in the…
▽ More
This paper addresses the problem of enumerating all supported efficient solutions for a linear multi-objective integer minimum cost flow problem (MOIMCF).
First, we highlight an inconsistency in various definitions of supported nondominated vectors for multi-objective integer linear programs (MOILP). Several characterizations for supported nondominated vectors/efficient solutions are used in the literature, which are equivalent in the non-integer case. However, they may lead to different sets of supported nondominated vectors/efficient solutions for MOILPs. This motivates us to summarize equivalent definitions and characterizations for supported efficient solutions and to distinguish between supported and weakly supported efficient solutions.
In this paper we derive an output-polynomial time algorithm to determine all supported efficient solutions for MOIMCF problems.
This is the first approach that solves this general problem in output-polynomial time.
Moreover, we prove that the existence of an output-polynomial time algorithm to determine all weakly supported nondominated vectors (or all weakly supported efficient solutions) for a MOIMCF problem with a fixed number of d>3 objectives can be excluded, unless P = NP.
△ Less
Submitted 22 May, 2023;
originally announced May 2023.
-
PDE characterisation of geometric distribution functions and quantiles
Authors:
Dimitri Konen
Abstract:
We show that in any Euclidean space, an arbitrary probability measure can be reconstructed explicitly by its geometric (or spatial) distribution function. The reconstruction takes the form of a (potentially fractional) linear PDE, where the differential operator is given in closed form. This result implies that, contrary to a common belief in the statistical depth community, geometric cdf's in pri…
▽ More
We show that in any Euclidean space, an arbitrary probability measure can be reconstructed explicitly by its geometric (or spatial) distribution function. The reconstruction takes the form of a (potentially fractional) linear PDE, where the differential operator is given in closed form. This result implies that, contrary to a common belief in the statistical depth community, geometric cdf's in principle provide exact control over the probability content of all depth regions. We present a comprehensive study of the regularity of the geometric cdf, and show that a continuous density in general does not give rise to a geometric cdf with enough regularity to reconstruct the density pointwise. Surprisingly, we prove that the reconstruction displays different behaviours in odd and even dimension: it is local in odd dimension and completely nonlocal in even dimension. We investigate this issue and provide a partial counterpart for even dimensions, and establish a general representation formula of the geometric cdf of spherically symmetric probability laws in odd dimensions. We provide explicit examples of the reconstruction of a density from its geometric cdf in dimension 2 and 3.
△ Less
Submitted 10 April, 2024; v1 submitted 24 August, 2022;
originally announced August 2022.
-
Timing jitter characterization of free-running dual-comb laser with sub-attosecond resolution using optical heterodyne detection
Authors:
Sandro L. Camenzind,
Daniel Koenen,
Benjamin Willenberg,
Justinas Pupeikis,
Christopher R. Phillips,
Ursula Keller
Abstract:
Pulse trains emitted from dual-comb systems are designed to have low relative timing jitter, making them useful for many optical measurement techniques such as optical ranging and spectroscopy. However, the characterization of low-jitter dual-comb systems is challenging because it requires measurement techniques with high sensitivity. Motivated by this challenge, we developed a technique based on…
▽ More
Pulse trains emitted from dual-comb systems are designed to have low relative timing jitter, making them useful for many optical measurement techniques such as optical ranging and spectroscopy. However, the characterization of low-jitter dual-comb systems is challenging because it requires measurement techniques with high sensitivity. Motivated by this challenge, we developed a technique based on an optical heterodyne detection approach for measuring the relative timing jitter of two pulse trains. The method is suitable for dual-comb systems with essentially any repetition rate difference. Furthermore, the proposed approach allows for continuous and precise tracking of the sampling rate. To demonstrate the technique, we perform a detailed characterization of a single-mode-diode pumped $\mathrm{Yb:CaF_2}$ dual-comb laser from a free-running polarization-multiplexed cavity. This new laser produces 115 fs pulses at 160 MHz repetition rate, with 130 mW of average power in each comb. The detection noise floor for the relative timing jitter between the two pulse trains reaches $8.0 \times 10^{-7} \, \mathrm{fs}^2/\mathrm{Hz} \; ( \sim 896 \: \mathrm{zs}/\sqrt{\mathrm{Hz}} )$, and the relative root mean square (rms) timing jitter is 13 fs when integrating from 100 Hz to 1 MHz. This performance indicates that the demonstrated laser is highly compatible with practical dual-comb spectroscopy, ranging, and sampling applications. Furthermore, our results show that the relative timing noise measurement technique can characterize dual-comb systems operating in free-running mode or with finite repetition rate differences while providing a sub-attosecond resolution, which was not feasible with any other approach before.
△ Less
Submitted 21 January, 2022; v1 submitted 11 November, 2021;
originally announced November 2021.
-
Finding all minimum cost flows and a faster algorithm for the K best flow problem
Authors:
David Könen,
Daniel R. Schmidt,
Christiane Spisla
Abstract:
This paper addresses the problem of determining all optimal integer solutions of a linear integer network flow problem, which we call the all optimal integer flow (AOF) problem. We derive an O(F (m + n) + mn + M ) time algorithm to determine all F many optimal integer flows in a directed network with n nodes and m arcs, where M is the best time needed to find one minimum cost flow. We remark that…
▽ More
This paper addresses the problem of determining all optimal integer solutions of a linear integer network flow problem, which we call the all optimal integer flow (AOF) problem. We derive an O(F (m + n) + mn + M ) time algorithm to determine all F many optimal integer flows in a directed network with n nodes and m arcs, where M is the best time needed to find one minimum cost flow. We remark that stopping Hamacher's well-known method for the determination of the K best integer flows at the first sub-optimal flow results in an algorithm with a running time of O(F m(n log n + m) + M ) for solving the AOF problem. Our improvement is essentially made possible by replacing the shortest path sub-problem with a more efficient way to determine a so called proper zero cost cycle using a modified depth-first search technique. As a byproduct, our analysis yields an enhanced algorithm to determine the K best integer flows that runs in O(Kn3 + M ). Besides, we give lower and upper bounds for the number of all optimal integer and feasible integer solutions. Our bounds are based on the fact that any optimal solution can be obtained by an initial optimal tree solution plus a conical combination of incidence vectors of all induced cycles with bounded coefficients.
△ Less
Submitted 27 January, 2022; v1 submitted 21 May, 2021;
originally announced May 2021.