ebnm: An R Package for Solving the Empirical Bayes Normal Means Problem Using a Variety of Prior Families

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Willwerscheid, Jason, Carbonetto, Peter, Stephens, Matthew
Format: Preprint
Published: 2021
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866911819634835456
author Willwerscheid, Jason
Carbonetto, Peter
Stephens, Matthew
author_facet Willwerscheid, Jason
Carbonetto, Peter
Stephens, Matthew
contents The empirical Bayes normal means (EBNM) model is important to many areas of statistics, including (but not limited to) multiple testing, wavelet denoising, and gene expression analysis. There are several existing software packages that can fit EBNM models under different prior assumptions and using different algorithms; however, the differences across interfaces complicate direct comparisons. Further, a number of important prior assumptions do not yet have implementations. Motivated by these issues, we developed the R package ebnm, which provides a unified interface for efficiently fitting EBNM models using a variety of prior assumptions, including nonparametric approaches. In some cases, we incorporated existing implementations into ebnm; in others, we implemented new fitting procedures with a focus on speed and numerical stability. We illustrate the use of ebnm in a detailed analysis of baseball statistics. By providing a unified and easily extensible interface, the ebnm package can facilitate development of new methods in statistics, genetics, and other areas; as an example, we briefly discuss the R package flashier, which harnesses methods in ebnm to provide a flexible and robust approach to matrix factorization.
format Preprint
id arxiv_https___arxiv_org_abs_2110_00152
institution arXiv
publishDate 2021
record_format arxiv
spellingShingle ebnm: An R Package for Solving the Empirical Bayes Normal Means Problem Using a Variety of Prior Families
Willwerscheid, Jason
Carbonetto, Peter
Stephens, Matthew
Computation
Methodology
The empirical Bayes normal means (EBNM) model is important to many areas of statistics, including (but not limited to) multiple testing, wavelet denoising, and gene expression analysis. There are several existing software packages that can fit EBNM models under different prior assumptions and using different algorithms; however, the differences across interfaces complicate direct comparisons. Further, a number of important prior assumptions do not yet have implementations. Motivated by these issues, we developed the R package ebnm, which provides a unified interface for efficiently fitting EBNM models using a variety of prior assumptions, including nonparametric approaches. In some cases, we incorporated existing implementations into ebnm; in others, we implemented new fitting procedures with a focus on speed and numerical stability. We illustrate the use of ebnm in a detailed analysis of baseball statistics. By providing a unified and easily extensible interface, the ebnm package can facilitate development of new methods in statistics, genetics, and other areas; as an example, we briefly discuss the R package flashier, which harnesses methods in ebnm to provide a flexible and robust approach to matrix factorization.
title ebnm: An R Package for Solving the Empirical Bayes Normal Means Problem Using a Variety of Prior Families
topic Computation
Methodology
url https://arxiv.org/abs/2110.00152