Source author record

Matthew Malloy

Matthew Malloy 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
9topics
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

Geometry of the Minimum Volume Confidence Sets

Computation of confidence sets is central to data science and machine learning, serving as the workhorse of A/B testing and underpinning the operation and analysis of reinforcement learning algorithms. This paper studies the geometry of the minimum-volume confidence sets for the multinomial parameter. When used in place of more standard confidence sets and intervals based on bounds and asymptotic approximation, learning algorithms can exhibit improved sample complexity. Prior work showed the minimum-volume confidence sets are the level-sets of a discontinuous function defined by an exact p-value. While the confidence sets are optimal in that they have minimum average volume, computation of membership of a single point in the set is challenging for problems of modest size. Since the confidence sets are level-sets of discontinuous functions, little is apparent about their geometry. This paper studies the geometry of the minimum volume confidence sets by enumerating and covering the continuous regions of the exact p-value function. This addresses a fundamental question in A/B testing: given two multinomial outcomes, how can one determine if their corresponding minimum volume confidence sets are disjoint? We answer this question in a restricted setting.

preprint2020arXiv

Digital Contact Tracing Using IP Colocation

The spread of an infectious disease through a population can be modeled using a network or a graph. In digital advertising, internet device graphs are graph data sets that organize identifiers produced by mobile phones, PCs, TVs, and tablets as they access media on the internet. Characterized by immense scale, they have become ubiquitous as they enable targeted advertising, content customization and tracking. This paper posits that internet device graphs, in particular those based on IP colocation, can provide significant utility in predicting and modeling the spread of infectious disease. Starting the week of March 16th, 2020, in the United States, many individuals began to `shelter-in-place' as schools and workplaces across the nation closed because of the COVID-19 pandemic. This paper quantifies the effect of the shelter-in-place orders on a large scale internet device graph with more than a billion nodes by studying the graph before and after orders went into effect. The effects are clearly visible. The structure of the graph suggests behavior least conducive to transmission of infection occurred in the US between April 12th and 19th, 2020. This paper also discusses the utility of device graphs for i) contact tracing, ii) prediction of `hot spots', iii) simulation of infectious disease spread, and iv) delivery of advertisement-based warnings to potentially exposed individuals. The paper also posits an overarching question: can systems and datasets amassed by entities in the digital ad ecosystem aid in the fight against COVID-19?

preprint2014arXiv

Adding Context to JWST Surveys with Current and Future 21cm Radio Observations

Infrared and radio observations of the Epoch of Reionization promise to revolutionize our understanding of the cosmic dawn, and major efforts with the JWST, MWA and HERA are underway. While measurements of the ionizing sources with infrared telescopes and the effect of these sources on the intergalactic medium with radio telescopes \emph{should} be complementary, to date the wildly disparate angular resolutions and survey speeds have made connecting proposed observations difficult. In this paper we develop a method to bridge the gap between radio and infrared studies. While the radio images may not have the sensitivity and resolution to identify individual bubbles with high fidelity, by leveraging knowledge of the measured power spectrum we are able to separate regions that are likely ionized from largely neutral, providing context for the JWST observations of galaxy counts and properties in each. By providing the ionization context for infrared galaxy observations, this method can significantly enhance the science returns of JWST and other infrared observations.

preprint2014arXiv

How to Search for Islands of Neutral Hydrogen in the $z \sim 5.5$ IGM

Observations of the Lyman-alpha (Ly-$α$) forest may allow reionization to complete as late as $z \sim 5.5$, provided the ionization state of the intergalactic medium (IGM) is sufficiently inhomogeneous at these redshifts. In this case, significantly neutral islands may remain amongst highly ionized gas with the ionized regions allowing some transmission through the Ly-$α$ forest. This possibility has the important virtue that it is eminently testable with existing Ly-$α$ forest data. In particular, we describe three observable signatures of significantly neutral gas in the $z \sim 5.5$ IGM. We use mock quasar spectra produced from numerical simulations of reionization to develop these tests. First, we quantify how the abundance and length of absorbed regions in the forest increase with the volume-averaged neutral fraction in our reionization model. Second, we consider stacking the transmission profile around highly absorbed regions in the forest. If and only if there is significantly neutral gas in the IGM, absorption in the damping wing of the Ly-$α$ line will cause the transmission to recover slowly as one moves from absorbed to transmitted portions of the spectrum. Third, the deuterium Ly-$β$ line should imprint a small but distinctive absorption feature slightly blueward of absorbed neutral regions in the Ly-$β$ forest. We show that these tests can be carried out with existing Keck HIRES spectra at $z \sim 5.5$, with the damping wing being observable for $< x_{\text{HI}} >\gtrsim 0.05$ and the deuterium feature observable with additional high-resolution spectra for $< x_{\text{HI}} >\gtrsim 0.2$.

preprint2014arXiv

On Modeling and Measuring the Temperature of the z~5 IGM

The temperature of the low-density intergalactic medium (IGM) at high redshift is sensitive to the timing and nature of hydrogen and HeII reionization, and can be measured from Lyman-alpha forest absorption spectra. Since the memory of intergalactic gas to heating during reionization gradually fades, measurements as close as possible to reionization are desirable. In addition, measuring the IGM temperature at sufficiently high redshifts should help to isolate the effects of hydrogen reionization since HeII reionization starts later, at lower redshift. Motivated by this, we model the IGM temperature at z>5 using semi-numeric models of patchy reionization. We construct mock Lyman-alpha forest spectra from these models and consider their observable implications. We find that the small-scale structure in the Lyman-alpha forest is sensitive to the temperature of the IGM even at redshifts where the average absorption in the forest is as high as 90%. We forecast the accuracy at which the z~5 IGM temperature can be measured using existing samples of high resolution quasar spectra, and find that interesting constraints are possible. For example, an early reionization model in which reionization ends at z~10 should be distinguishable -- at high statistical significance -- from a lower redshift model where reionization completes at z~6. We discuss improvements to our modeling that may be required to robustly interpret future measurements.

preprint2013arXiv

lil' UCB : An Optimal Exploration Algorithm for Multi-Armed Bandits

The paper proposes a novel upper confidence bound (UCB) procedure for identifying the arm with the largest mean in a multi-armed bandit game in the fixed confidence setting using a small number of total samples. The procedure cannot be improved in the sense that the number of samples required to identify the best arm is within a constant factor of a lower bound based on the law of the iterated logarithm (LIL). Inspired by the LIL, we construct our confidence bounds to explicitly account for the infinite time horizon of the algorithm. In addition, by using a novel stopping time for the algorithm we avoid a union bound over the arms that has been observed in other UCB-type algorithms. We prove that the algorithm is optimal up to constants and also show through simulations that it provides superior performance with respect to the state-of-the-art.

preprint2013arXiv

On Finding the Largest Mean Among Many

Sampling from distributions to find the one with the largest mean arises in a broad range of applications, and it can be mathematically modeled as a multi-armed bandit problem in which each distribution is associated with an arm. This paper studies the sample complexity of identifying the best arm (largest mean) in a multi-armed bandit problem. Motivated by large-scale applications, we are especially interested in identifying situations where the total number of samples that are necessary and sufficient to find the best arm scale linearly with the number of arms. We present a single-parameter multi-armed bandit model that spans the range from linear to superlinear sample complexity. We also give a new algorithm for best arm identification, called PRISM, with linear sample complexity for a wide range of mean distributions. The algorithm, like most exploration procedures for multi-armed bandits, is adaptive in the sense that the next arms to sample are selected based on previous samples. We compare the sample complexity of adaptive procedures with simpler non-adaptive procedures using new lower bounds. For many problem instances, the increased sample complexity required by non-adaptive procedures is a polynomial factor of the number of arms.

preprint2011arXiv

On the Limits of Sequential Testing in High Dimensions

This paper presents results pertaining to sequential methods for support recovery of sparse signals in noise. Specifically, we show that any sequential measurement procedure fails provided the average number of measurements per dimension grows slower then log s / D(f0||f1) where s is the level of sparsity, and D(f0||f1) the Kullback-Leibler divergence between the underlying distributions. For comparison, we show any non-sequential procedure fails provided the number of measurements grows at a rate less than log n / D(f1||f0), where n is the total dimension of the problem. Lastly, we show that a simple procedure termed sequential thresholding guarantees exact support recovery provided the average number of measurements per dimension grows faster than (log s + log log n) / D(f0||f1), a mere additive factor more than the lower bound.

preprint2011arXiv

Sequential Analysis in High Dimensional Multiple Testing and Sparse Recovery

This paper studies the problem of high-dimensional multiple testing and sparse recovery from the perspective of sequential analysis. In this setting, the probability of error is a function of the dimension of the problem. A simple sequential testing procedure is proposed. We derive necessary conditions for reliable recovery in the non-sequential setting and contrast them with sufficient conditions for reliable recovery using the proposed sequential testing procedure. Applications of the main results to several commonly encountered models show that sequential testing can be exponentially more sensitive to the difference between the null and alternative distributions (in terms of the dependence on dimension), implying that subtle cases can be much more reliably determined using sequential methods.