Source author record

Anya Katsevich

Anya Katsevich 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 math.FA math.RA

Catalog footprint

What is connected

5works

5topics

3close 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.

preprint2024arXiv

On the Approximation Accuracy of Gaussian Variational Inference

The main computational challenge in Bayesian inference is to compute integrals against a high-dimensional posterior distribution. In the past decades, variational inference (VI) has emerged as a tractable approximation to these integrals, and a viable alternative to the more established paradigm of Markov Chain Monte Carlo. However, little is known about the approximation accuracy of VI. In this work, we bound the TV error and the mean and covariance approximation error of Gaussian VI in terms of dimension and sample size. Our error analysis relies on a Hermite series expansion of the log posterior whose first terms are precisely cancelled out by the first order optimality conditions associated to the Gaussian VI optimization problem.

preprint2022arXiv

From local equilibrium to numerical PDE: Metropolis crystal surface dynamics in the rough scaling limit

We derive the PDE governing the hydrodynamic limit of a Metropolis rate crystal surface height process in the "rough scaling" regime introduced by Marzuola and Weare. The PDE takes the form of a continuity equation, and the expression for the current involves a numerically computed multiplicative correction term similar to a mobility. The correction accounts for the fact that, unusually, the local equilibrium distribution of the process is not a local Gibbs measure even though the global equilibrium distribution is Gibbs. We give definitive numerical evidence of this fact, originally suggested in Gao, et. al., Pure and Applied Analysis (2021). In that paper, an approximate PDE -- our PDE, but without the correction term -- was derived for the limit of the Metropolis rate process under the assumption of a local Gibbs distribution. Our main contribution is to present a numerical method to compute the corrected macroscopic current, which is given by a function of the third spatial derivative of the height profile. Our method exploits properties of the local equilibrium (LE) state of the third order finite difference process. We find that the LE state of this process is not only useful for deriving the PDE; it also enjoys nonstandard properties which are interesting in their own right. Namely, we demonstrate that the LE state is a "rough LE", a novel kind of LE state discovered in our recent work on an Arrhenius rate crystal surface process.

preprint2022arXiv

The Local Equilibrium State of a Crystal Surface Jump Process in the Rough Scaling Regime

We investigate the local equilibrium (LE) distribution of a crystal surface jump process as it approaches its hydrodynamic (continuum) limit in a nonstandard scaling regime introduced by Marzuola and Weare. The atypical scaling leads to a local equilibrium state whose structure is novel, to the best of our knowledge. The distinguishing characteristic of the new, "rough" LE state is that the ensemble average of single lattice site observables do not vary smoothly across lattice sites. We investigate numerically and analytically how the rough LE state affects the convergence mechanism via three key limits, and show that by comparison, more standard, "smooth" LE states satisfy stronger versions of these limits.

preprint2020arXiv

Likelihood Maximization and Moment Matching in Low SNR Gaussian Mixture Models

We derive an asymptotic expansion for the log likelihood of Gaussian mixture models (GMMs) with equal covariance matrices in the low signal-to-noise regime. The expansion reveals an intimate connection between two types of algorithms for parameter estimation: the method of moments and likelihood optimizing algorithms such as Expectation-Maximization (EM). We show that likelihood optimization in the low SNR regime reduces to a sequence of least squares optimization problems that match the moments of the estimate to the ground truth moments one by one. This connection is a stepping stone toward the analysis of EM and maximum likelihood estimation in a wide range of models. A motivating application for the study of low SNR mixture models is cryo-electron microscopy data, which can be modeled as a GMM with algebraic constraints imposed on the mixture centers. We discuss the application of our expansion to algebraically constrained GMMs, among other example models of interest.

preprint2014arXiv

On De Graaf spaces of pseudoquotients

A space of pseudoquotients $\mathcal{B}(X,S)$ is defined as equivalence classes of pairs $(x,f)$, where $x$ is an element of a non-empty set $X$, $f$ is an element of $S$, a commutative semigroup of injective maps from $X$ to $X$, and $(x,f) \sim (y,g)$ if $gx=fy$. In this note we consider a generalization of this construction where the assumption of commutativity of $S$ by Ore type conditions. As in the commutative case, $X$ can be identified with a subset of $\mathcal{B}(X,S)$ and $S$ can be extended to a group $G$ of bijections on $\mathcal{B}(X,S)$. We introduce a natural topology on $\mathcal{B}(X,S)$ and show that all elements of $G$ are homeomorphisms on $\mathcal{B}(X,S)$.