Preprint arXiv Research Articles

Smoothed alternatives of the two-sample median and Wilcoxon's rank sum tests

ABSTRACTWe discuss smoothed rank statistics for testing the location shift parameter of the two-sample problem. They are based on discrete test statistics – the median and Wilcoxon's rank sum tests. For the one-sample problem, Maesono et al. [Smoothed nonparametric tests and their properties. arXiv preprint. 2016; ArXiv:1610.02145] reported that some nonparametric discrete tests have a problem with their p-values because of their discreteness. The p-values of Wilcoxon's test are frequently smaller than those of the median test in the tail area. This leads to an arbitrary choice of the median and Wilcoxon's rank sum tests. To overcome this problem, we propose smoothed versions of those tests. The smoothed tests inherit the good properties of the original tests and are asymptotically equivalent to them. We study the significance probabilities and local asymptotic powers of the proposed tests.

An upper diameter bound for compact Ricci solitons with application to the Hitchin–Thorpe inequality. II

In this addendum to the author’s paper [H. Tadano, J. Math. Phys. 58, 023503 (2017)], by using an upper diameter bound for compact shrinking Ricci solitons by Wu [preprint arXiv:1706.07897 (2017)], we shall give a new sufficient condition for compact four-dimensional shrinking Ricci solitons to satisfy the Hitchin–Thorpe inequality. Our new condition sharpens the validity of the Hitchin–Thorpe inequality given by Tadano [J. Math. Phys. 58, 023503 (2017)].

Chapter 4 – Migration research and analysis: Growth, reach and recent contributions

Chapter 4 – Migration research and analysis: Growth, reach and recent contributions

An exposition of the false confidence theorem

A recent paper (Balch et al. (), ‘Satellite conjunction analysis and the false confidence theorem’, arXiv 5 preprint arXiv:1706.08565.) presents the “false confidence theorem,” which has potentially broad implications for statistical inference using Bayesian posterior uncertainty. This theorem says that with an arbitrarily large (sampling/frequentist) probability, there exists a set that does not contain the true parameter value but has an arbitrarily large posterior probability. As the use of Bayesian methods has become increasingly popular in applications of science, engineering and business, it is critically important to understand when Bayesian procedures lead to problematic statistical inferences or interpretations. In this paper, we consider a number of examples demonstrating the paradoxical nature of false confidence to begin to understand the contexts in which the false confidence theorem does (and does not) play a meaningful role in statistical inference. Our examples illustrate that models involving marginalization to non‐linear, not one‐to‐one functions of multiple parameters play a key role in more extreme manifestations of false confidence. © 2018 John Wiley & Sons, Ltd.

Commutativity in the Algorithmic Lovász Local Lemma

We consider the recent formulation of the algorithmic Lovász Local Lemma [N. Harvey and J. Vondrák, in Proceedings of FOCS, 2015, pp. 1327--1345; D. Achlioptas and F. Iliopoulos, in Proceedings of SODA, 2016, pp. 2024--2038; D. Achlioptas, F. Iliopoulos, and V. Kolmogorov, A Local Lemma for Focused Stochastic Algorithms, arXiv preprint, 2018] for finding objects that avoid “bad features,” or “flaws.” It extends the Moser--Tardos resampling algorithm [R. A. Moser and G. Tardos, J. ACM, 57 (2010), 11] to more general discrete spaces. At each step the method picks a flaw present in the current state and goes to a new state according to some prespecified probability distribution (which depends on the current state and the selected flaw). However, the recent formulation is less flexible than the Moser--Tardos method since it requires a specific flaw selection rule, whereas the algorithm of Moser and Tardos allows an arbitrary rule (and thus can potentially be implemented more efficiently). We formulate a new “commutativity” condition and prove that it is sufficient for an arbitrary rule to work. It also enables an efficient parallelization under an additional assumption. We then show that existing resampling oracles for perfect matchings and permutations do satisfy this condition.

A Sparse Estimate for Multisublinear Forms Involving Vector-valued Maximal Functions

We prove a sparse bound for the m -sublinear form associated to vector-valued maximal functions of Fefferman-Stein type. As a consequence, we show that the sparse bounds of multisublinear operators are preserved via l r -valued extension. This observation is in turn used to deduce vector-valued, multilinear weighted norm inequalities for multisublinear operators obeying sparse bounds, which are out of reach for the extrapolation theory developed by Cruz-Uribe and Martell in Limited range multilinear extrapolation with applications to the bilinear Hilbert transform , preprint arXiv:1704.06833 (2017). As an example, vector-valued multilinear weighted inequalities for bilinear Hilbert transforms are deduced from the scalar sparse domination theorem of Domination of multilinear singular integrals by positive sparse forms , preprint arXiv:1603.05317.

Data and Decisions about Scholarly Knowledge

Data and Decisions about Scholarly Knowledge Daniel Goroff (bio) and Josh Greenberg (bio) current discussions about scholarly knowledge do not necessarily concern great discoveries, disciplinary advances, or epistemological puzzles. Words that often come up instead are "measurement," "autonomy," and "accountability." Sometimes explicitly and sometimes implicitly, academics worry about potentially punitive policies or other possible threats to the scholarly enterprise as we know it. While it can be easy and engaging enough to speculate about institutions and incentives this way, what ultimately matters are not discussions but decisions. How, after all, should we allocate scarce resources? Scholars are remarkably free to make their own choices, of course, and routinely trade off what others value against their own values. But what about funders and other functionaries? This note presents an idiosyncratic sample of how two program directors at a foundation that supports scientific research see the role that measures of scholarly knowledge play in making decisions. (Note, however, that opinions expressed here do not necessarily represent those of the Alfred P. Sloan Foundation.) Of course, decisions these days are supposed to be made based on data. The revolution in information technology is ultimately responsible for much disruption and anxiety in the world at large as well as in academia. What good is all this data? Perhaps there is some irony here. When trying to publish an empirical paper, academics are often cavalier about overlooking the extent to which correlation masquerades [End Page 733] as causation in contemporary data science. But when data scientists turn their bibliometric analyses on academia, suddenly we hear scholarly cautions against placing too much weight on correlations, measurements, and rankings that could not possibly tell the whole story about what it means to be a scholar. Here are four examples of resource allocation questions that involve the measurement of scholarly knowledge. ALLOCATION OF RESEARCHERS' ATTENTION TO ONE SCIENTIFIC RESULT OR RESOURCE AS OPPOSED TO ANOTHER Everyone knows that finding what you are looking for among all the information available is a great and growing challenge. This is especially true as we move from a system of filtering-then-disseminating (think of traditional journals) to one of disseminating-then-filtering (think of the arXiv and similar preprint servers). Just as search engines help with navigating the web generally, bibliometrics can be quite helpful to scholars in particular. Data analysis can surface and explore hidden relationships, connections, or patterns among various references and resources—including data-sets, models, code, and other artifacts as well as papers. The Alfred P. Sloan Foundation has therefore funded projects ranging from Wikipedia and the Digital Public Library of America to international data citation standards and planning for a global library of mathematics. How to search for mathematical results is a particularly intriguing problem, for example, that for-profit companies have little incentive to address but that could greatly enhance the usefulness of the scientific literature. ALLOCATION OF CURATORIAL STEWARDSHIP TO ONE SCIENTIFIC PRODUCT AS OPPOSED TO ANOTHER We really cannot and should not keep everything, so some triage is necessary. It sometimes seems otherwise because the marginal cost of one more call on a resource that is already online may be negligible. [End Page 734] But editors, publishers, librarians, and archivists still perform important and costly services. Their work requires judgment and not just bibliometrics. In fact, curation by humans is arguably becoming even more valuable as computers and computing begin to dominate so much else. Given all the unfiltered dissemination of everything from books to bits, what is worth curating? In the sciences, this question is tied up with the reproducibility of research. There are calls for openness and transparency, particularly concerning research data, which push towards preserving everything a scholar ever thinks or does. But reproducibility is not a monitoring problem as much as it is a methodological one. Exploratory work with data may be fun, creative, and all the rage, but perhaps all that matters at that stage are the plausible hypotheses generated. Confirmatory research to test those hypotheses then deserves careful methods, fresh data, scholarly scrutiny, and detailed dissemination of both positive as well as null results. Science simply cannot progress efficiently if we continue conflating, in publications and in...

Hexagonal Smoothness-Increasing Accuracy-Conserving Filtering

Discontinuous Galerkin (DG) methods are a popular class of numerical techniques to solve partial differential equations due to their higher order of accuracy. However, the inter-element discontinuity of a DG solution hinders its utility in various applications, including visualization and feature extraction. This shortcoming can be alleviated by postprocessing of DG solutions to increase the inter-element smoothness. A class of postprocessing techniques proposed to increase the inter-element smoothness is SIAC filtering. In addition to increasing the inter-element continuity, SIAC filtering also raises the convergence rate from order $$k + 1$$ to order $$2k + 1$$ . Since the introduction of SIAC filtering for univariate hyperbolic equations by Cockburn et al. (Math Comput 72(242):577–606, 2003), many generalizations of SIAC filtering have been proposed. Recently, the idea of dimensionality reduction through rotation has been the focus of studies in which a univariate SIAC kernel has been used to postprocess a two-dimensional DG solution (Docampo-Sanchez et al. in Multi-dimensional filtering: reducing the dimension through rotation, 2016. arXiv preprint arXiv:1610.02317 ). However, the scope of theoretical development of multidimensional SIAC filters has never gone beyond the usage of tensor product multidimensional B-splines or the reduction of the filter dimension. In this paper, we define a new SIAC filter called hexagonal SIAC (HSIAC) that uses a nonseparable class of two-dimensional spline functions called hex splines. In addition to relaxing the separability assumption, the proposed HSIAC filter provides more symmetry to its tensor-product counterpart. We prove that the superconvergence property holds for a specific class of structured triangular meshes using HSIAC filtering and provide numerical results to demonstrate and validate our theoretical results.

Decomposing graphs into paths and trees

In [Bensmail, J., A. Harutyunyan, T.-N. Le and S. Thomassé, Edge-partitioning a graph into paths: beyond the Barát-Thomassen conjecture, arXiv preprint arXiv:1507.08208 (2015)], the authors conjecture that for a fixed tree T, the edge set of any graph G of size divisible by size of T with sufficiently high degree can be decomposed into disjoint copies of T, provided that G is sufficiently highly connected in terms of maximal degree of T. In [Bensmail, J., A. Harutyunyan, T.-N. Le and S. Thomassé, Edge-partitioning a graph into paths: beyond the Barát-Thomassen conjecture, arXiv preprint arXiv:1507.08208 (2015)], the conjecture was proven for trees of maximal degree 2 (i.e., paths). In particular, it was shown that in the case of paths, the conjecture holds for 24-edge-connected graphs. We improve this result showing that 3-edge-connectivity suffices, which is best possible. We disprove the conjecture for trees of maximum degree greater than two and prove a relaxed version of the conjecture that concerns decomposing the edge set of a graph into disjoint copies of two fixed trees of coprime sizes.

HFAG 2016 and PDG 2016 τ lepton averages and |Vus| determination from τ data

The HFAG-Tau sub-group has performed in 2016 two similar global fits on the τ lepton branching fractions measurements. One, unitarity constrained, has been published in the 2016 edition of the Review of Particle Physics [C. Patrignani, et al., Review of particle physics, Chin. Phys. C40 (2016) 100001. doi:10.1088/1674-1137/40/10/100001] and a second one, without the unitarity constraint, is currently preliminary and is expected to be published soon on the HFAG 2016 report [HFAG, Averages of b-hadron, c-hadron, and τ-lepton properties as of Summer 2016, forthcoming as arXiv preprint]. The resulting τ branching fractions are used to test the Standard Model lepton universality predictions and to measure the Cabibbo-Kobayashi-Maskawa (CKM) matrix element |Vus|.

Brownian dynamics of confined suspensions of active microrollers.

We develop efficient numerical methods for performing many-body Brownian dynamics simulations of a recently observed fingering instability in an active suspension of colloidal rollers sedimented above a wall [M. Driscoll, B. Delmotte, M. Youssef, S. Sacanna, A. Donev, and P. Chaikin, Nat. Phys. (2016), preprint arXiv:1609.08673. We present a stochastic Adams-Bashforth integrator for the equations of Brownian dynamics, which has the same cost but is more accurate than the widely used Euler-Maruyama scheme, and use a random finite difference to capture the stochastic drift proportional to the divergence of the configuration-dependent mobility matrix. We generate the Brownian increments using a Krylov method and show that for particles confined to remain in the vicinity of a no-slip wall by gravity or active flows, the number of iterations is independent of the number of particles. Our numerical experiments with active rollers show that the thermal fluctuations set the characteristic height of the colloids above the wall, both in the initial condition and the subsequent evolution dominated by active flows. The characteristic height in turn controls the time scale and wavelength for the development of the fingering instability.

A novel Monte Carlo approach to hybrid local volatility models

We present in a Monte Carlo simulation framework, a novel approach for the evaluation of hybrid local volatility [Risk, 1994, 7, 18–20], [Int. J. Theor. Appl. Finance, 1998, 1, 61–110] models. In particular, we consider the stochastic local volatility model—see e.g. Lipton et al. [Quant. Finance, 2014, 14, 1899–1922], Piterbarg [Risk, 2007, April, 84–89], Tataru and Fisher [Quantitative Development Group, Bloomberg Version 1, 2010], Lipton [Risk, 2002, 15, 61–66]—and the local volatility model incorporating stochastic interest rates—see e.g. Atlan [ArXiV preprint math/0604316, 2006], Piterbarg [Risk, 2006, 19, 66–71], Deelstra and Rayée [Appl. Math. Finance, 2012, 1–23], Ren et al. [Risk, 2007, 20, 138–143]. For both model classes a particular (conditional) expectation needs to be evaluated which cannot be extracted from the market and is expensive to compute. We establish accurate and ‘cheap to evaluate’ approximations for the expectations by means of the stochastic collocation method [SIAM J. Numer. Anal., 2007, 45, 1005–1034], [SIAM J. Sci. Comput., 2005, 27, 1118–1139], [Math. Models Methods Appl. Sci., 2012, 22, 1–33], [SIAM J. Numer. Anal., 2008, 46, 2309–2345], [J. Biomech. Eng., 2011, 133, 031001], which was recently applied in the financial context [Available at SSRN 2529691, 2014], [J. Comput. Finance, 2016, 20, 1–19], combined with standard regression techniques. Monte Carlo pricing experiments confirm that our method is highly accurate and fast.

On the Location of Eigenvalues of Real Matrices

The research in this paper is motivated by a recent work of I. Barany and J. Solymosi [I. Barany and J. Solymosi. Gershgorin disks for multiple eigenvalues of non-negative matrices. Preprint arXiv no. 1609.07439, 2016.] about the location of eigenvalues of nonnegative matrices with geometric multiplicity higher than one. In particular, an answer to a question posed by Barany and Solymosi, about how the location of the eigenvalues can be improved in terms of their geometric multiplicities is obtained. New inclusion sets for the eigenvalues of a real square matrix, called Ger\v{s}gorin discs of the second type, are introduced. It is proved that under some conditions, an eigenvalue of a real matrix is in a Ger\v{s}gorin disc of the second type. Some relationships between the geometric multiplicities of eigenvalues and these new inclusion sets are established. Some other related results, consequences, and examples are presented. The results presented here apply not only to nonnegative matrices, but extend to all real matrices, and some of them do not depend on the geometric multiplicity.

Multiparametric Characterization of Single, Unlabeled Proteins in Solution

With large diversity in structure and function, as well as clinical relevance, proteins represent an important target for biophysical characterization. Resistive pulse-based nanopore sensing is a compelling platform for this task, as it provides information about individual proteins during their translocation through the zeptoliter sensing volume inside of a nanopore. Previous work in our group used lipid-coated synthetic nanopores to extract five distinct protein descriptors by analyzing modulations in ionic current, ΔI, resulting from the translocation and rotation of individual lipid-tethered proteins [1]. However, while coating of nanopore walls with fluid lipid bilayers and tethering proteins with lipid anchors is useful for preventing non-specific protein adsorption, extending translocation times, and increasing specificity, it is technically demanding and its success depends on nanopore diameter, geometry, and surface chemistry. This work explores the possibility of using the detergent Tween-20 for nanopore surface coating [2] combined with high bandwidth recording electronics to characterize freely translocating, untethered proteins on a single molecule level. Here, we utilize the dependence of ΔI on the orientation of non-spherical proteins transiting a nanopore to determine their intrinsic shapes, volumes, and dipole moments in solution. The ability to thoroughly examine unlabeled, natively-folded proteins in an aqueous sample on a single molecule level signifies an important step toward the use of nanopores for proteomic and diagnostic applications.[1] Yusko, E.C., et al., Real-time shape approximation and 5-D fingerprinting of single proteins. arXiv preprint arXiv:1510.01935, 2015.[2] Hu, Rui, et al. Intrinsic and membrane-facilitated α-synuclein oligomerization revealed by label-free detection through solid-state nanopores. Scientific reports 6, 2016.

CM values of regularized theta lifts and harmonic weak Maaß forms of weight 1

We study special values of regularized theta lifts at complex multiplication (CM) points. In particular, we show that CM values of Borcherds products can be expressed in terms of finitely many Fourier coefficients of certain harmonic weak Maa{\ss} forms of weight one. As it turns out, these coefficients are logarithms of algebraic integers whose prime ideal factorization is determined by special cycles on an arithmetic curve. Our results imply a conjecture of Duke and Li and give a new proof of the modularity of a certain arithmetic generating series of weight one studied by Kudla, Rapoport and Yang. The results of the paper are much improved in comparison to the 2012 preprint arxiv:1208.2386 which contained partial results in the same direction. Moreover, they are also an improvement of the main result of the authors thesis (CM values of regularized theta lifts, TU Darmstadt, 2013).

The maximal excess charge for a family of density-matrix-functional theories including Hartree-Fock and Müller theories

We will give a proof that the maximal excess charge for an atom described by a family of density-matrix-functionals, which includes Hartree-Fock and Müller theories, is bounded by a universal constant. We will use the new technique introduced by Frank et al. [preprint arXiv:1608.05625 (2016)].

Cherednik and Hecke Algebras of Varieties with a Finite Group Action

This paper is an expanded and updated version of the preprint arXiv:math/0406499. It includes a more detailed description of the basics of the theory of Cherednik and Hecke algebras of varieties started in arXiv:math/0406499, as well as a new Section 4, which reviews the developments in this theory since 2004 with references to the relevant literature. Let $G$ be a finite group of linear transformations of a finite dimensional complex vector space $V$. To this data one can attach a family of algebras $H_{t,c}(V,G)$, parametrized by complex numbers $t$ and conjugation invariant functions $c$ on the set of complex reflections in $G$, which are called rational Cherednik algebras. These algebras have been studied for over 15 years and revealed a rich structure and deep connections with algebraic geometry, representation theory, and combinatorics. In this paper, we define global analogs of Cherednik algebras, attached to any smooth algebraic or analytic variety $X$ with a finite group $G$ of automorphisms. We show that many interesting properties of Cherednik algebras (such as the PBW theorem, universal deformation property, relation to Calogero-Moser spaces, action on quasiinvariants) still hold in the global case, and give several interesting examples. Then we define the KZ functor for global Cherednik algebras, and use it to define (in the case $\pi_2(X)\otimes \Bbb Q=0$) a flat deformation of the orbifold fundamental group of the orbifold $X/G$, which we call the Hecke algebra of $X/G$. This includes usual, affine, and double affine Hecke algebras for Weyl groups, Hecke algebras of complex reflection groups, as well as many new examples.

The accessible lasso models

ABSTRACTA new line of research [Lee et al. Exact post-selection inference with the lasso. arXiv preprint arXiv:1311.6238 (2013)] on the lasso [Tibshirani. Regression shrinkage and selection via the lasso. J R Stat Soc Ser B (Methodological). 1996;267–288] exploits the beautiful geometric fact that the lasso fit is the residual from projecting the response vector y onto a certain convex polytope [Tibshirani et al. Degrees of freedom in lasso problems. Ann Stat. 2012;40:1198–1232]. This geometric picture also allows an exact geometric description of the set of accessible lasso models for a given design matrix, that is, which configurations of the signs of the coefficients it is possible to realize with some choice of y. In particular, the accessible lasso models are those that correspond to a face of the convex hull of all the feature vectors together with their negations. This convex hull representation then permits the enumeration and bounding of the number of accessible lasso models, which in turn provides a direct proof of model selection inconsistency when the size of the true model is greater than half the number of observations.

One-species Vlasov-Poisson-Landau system for soft potentials in ℝ3

We consider the global classical solution near a global Maxwellian to the one-species Vlasov-Poisson-Landau system in the whole space R x 3 . It is shown that our global solvability result is obtained under the weaker smallness condition on the initial perturbation than that of Duan et al., [preprint arXiv:1112.3261 (2011)] and Lei et al., [Kinet. Relat. Models 7(3), 551–590 (2014)].

Arithmetic mirror symmetry for genus 1 curves with n marked points

We establish a \({\mathbb {Z}}[[t_1,\ldots , t_n]]\)-linear derived equivalence between the relative Fukaya category of the 2-torus with n distinct marked points and the derived category of perfect complexes on the n-Tate curve. Specialising to \(t_1= \cdots =t_n=0\) gives a \({\mathbb {Z}}\)-linear derived equivalence between the Fukaya category of the n-punctured torus and the derived category of perfect complexes on the standard (Neron) n-gon. We prove that this equivalence extends to a \({\mathbb {Z}}\)-linear derived equivalence between the wrapped Fukaya category of the n-punctured torus and the derived category of coherent sheaves on the standard n-gon. The corresponding results for \(n=1\) were established in Lekili and Perutz (Arithmetic mirror symmetry for the 2-torus (preprint) arXiv:1211.4632, 2012).