Abstract:We define a family of probability distributions for random count matrices with a potentially unbounded number of rows and columns. The three distributions we consider are derived from the gamma-Poisson, gamma-negative binomial, and beta-negative binomial processes, which we refer to generically as a family of negative-binomial processes. Because the models lead to closed-form update equations within the context of a Gibbs sampler, they are natural candidates for nonparametric Bayesian priors over count matrices. A key aspect of our analysis is the recognition that although the random count matrices within the family are defined by a row-wise construction, their columns can be shown to be independent and identically distributed (iid). This fact is used to derive explicit formulas for drawing all the columns at once. Moreover, by analyzing these matrices’ combinatorial structure, we describe how to sequentially construct a column-iid random count matrix one row at a time, and derive the predictive distribution of a new row count vector with previously unseen features. We describe the similarities and differences between the three priors, and argue that the greater flexibility of the gamma- and beta-negative binomial processes—especially their ability to model over-dispersed, heavy-tailed count data—makes these well suited to a wide variety of real-world applications. As an example of our framework, we construct a naive-Bayes text classifier to categorize a count vector to one of several existing random count matrices of different categories. The classifier supports an unbounded number of features and, unlike most existing methods, it does not require a predefined finite vocabulary to be shared by all the categories, and needs neither feature selection nor parameter tuning. Both the gamma- and beta-negative binomial processes are shown to significantly outperform the gamma-Poisson process when applied to document categorization, with comparable performance to other state-of-the-art supervised text classification algorithms. Supplementary materials for this article are available online.

The Babe Ruth Algorithm: a fast, unbiased procedure to randomize presence-absence data matrices with fixed row and column totals

Probabilistic Analysis for the Random Assignment Problem

Assessing Data Mining Results on Matrices with Randomization

Randomization or Condensation?

Priors for Random Count Matrices Derived from a Family of Negative Binomial Processes

Randomization Tests for Adaptively Collected Data

A Fast Random Sampling Algorithm for Sparsifying Matrices

The Randomness Recycler: A new technique for perfect sampling

Biased versus unbiased numerical methods for stochastic simulations

The necessity and power of random, under-sampled experiments in biology

Algorithmic approach for an unique definition of the next generation matrix

Pattern detection in bipartite networks: a review of terminology, applications and methods

Improved estimates for the number of non-negative integer matrices with given row and column sums

Faster Randomized Methods for Orthogonality Constrained Problems

Missing Not at Random in Matrix Completion: The Effectiveness of Estimating Missingness Probabilities Under a Low Nuclear Norm Assumption

Fast randomized algorithms for low-rank matrix approximations with applications in global comparative analysis of a class of data sets

Randomized Approximation of the Gram Matrix: Exact Computation and Probabilistic Bounds

A Representation of Binary Matrices

A more practical approach for the Benjamini-Hochberg FDR controlling procedure for huge-scale testing problems

How often does the best team win? A unified approach to understanding randomness in North American sport