Search papers, labs, and topics across Lattice.
This paper introduces a novel method for estimating the expected output of wide, randomly initialized MLPs over Gaussian inputs without relying on Monte Carlo sampling. The approach propagates approximate distributions of activations through the network using cumulants and Hermite expansions. Results demonstrate that this estimator achieves a target mean squared error with significantly fewer FLOPs than sampling, particularly for estimating the probabilities of rare events.
Skip the sampling: accurately predict the behavior of wide, random MLPs with a fraction of the compute, especially when assessing rare, high-stakes outcomes.
By far the most common way to estimate an expected loss in machine learning is to draw samples, compute the loss on each one, and take the empirical average. However, sampling is not necessarily optimal. Given an MLP at initialization, we show how to estimate its expected output over Gaussian inputs without running samples through the network at all. Instead, we produce approximate representations of the distributions of activations at each layer, leveraging tools such as cumulants and Hermite expansions. We show both theoretically and empirically that for sufficiently wide networks, our estimator achieves a target mean squared error using substantially fewer FLOPs than Monte Carlo sampling. We find moreover that our methods perform particularly well at estimating the probabilities of rare events, and additionally demonstrate how they can be used for model training. Together, these findings suggest a path to producing models with a greatly reduced probability of catastrophic tail risks.
