1 Introduction
Many natural phenomena decompose into latent features. For example, visual scenes can be decomposed into objects; genetic regulatory networks can be decomposed into transcription factors; music can be decomposed into spectral components. In these examples, multiple latent features can be simultaneously active, and each can influence the observed data. Dimensionality reduction methods, such as principal component analysis, factor analysis, and probabilistic matrix factorization, provide a statistical approach to inferring latent features
[4]. These methods characterize a small set of dimensions, or features, and model each data point as a weighted combination of these features. Dimensionality reduction can improve predictions and identify hidden structures in observed data.Dimensionality reduction methods typically require that the number of latent features (i.e., the number of dimensions) be fixed in advance. Researchers have recently proposed a more flexible approach based on Bayesian nonparametric models, where the number of features is inferred from the data through a posterior distribution. These models are usually based on the Indian buffet process (IBP; [16, 17]), a prior over binary matrices with a finite number of rows (corresponding to data points) and an infinite number of columns (corresponding to latent features). Using the IBP as a building block, Bayesian nonparametric latent feature models have been applied to several statistical problems (e.g., [19, 20, 22, 23]). Since the number of features is effectively unbounded, these models are sometimes known as “infinite” latent feature models.
The IBP assumes that data are exchangeable
: permuting the order of rows leaves the probability of a feature matrix unchanged. This assumption may be appropriate for some data sets, but for many others we expect dependencies between data points and, consequently, between their latent representations. As examples, the latent features describing human motion are autocorrelated over time; the latent features describing environmental risk factors are autocorrelated over space. In this paper, we present a generalization of the IBP—the
distance dependent IBP (ddIBP)—that addresses this limitation. The ddIBP allows infinite latent feature models to capture nonexchangeable structure.The problem of adapting nonparametric models to nonexchangeable data has been studied extensively in the mixturemodeling literature. In particular, variants of the Dirichlet process mixture model allow dependencies between data points (e.g., [9, 15, 7, 12, 24, 1]). These dependencies may be spatial, temporal or more generally covariatedependent; the effect of such dependencies is to induce sharing of features between nearby data points.
Among these methods is the distance dependent Chinese restaurant process (ddCRP; [6]). The ddCRP is a nonexchangeable generalization of the Chinese restaurant process (CRP), the prior over partitions of data that emerges in Bayesian nonparametric mixture modeling [5, 14, 25]. The ddCRP models nonexchangeability by using using distances between data points—nearby data points (e.g., in time or space) are more likely to be assigned to the same mixture component. The ddIBP extends these ideas to infinite latent feature models, where distances between data points influence featuresharing, and nearby data points are more likely to share latent features.
We review the IBP in Section 2.1 and develop the ddIBP in Section 2.2. Like the ddCRP, the ddIBP lacks marginal invariance, which means that removing one observation changes the distribution over the other observations. We discuss this property further in Section 2.3. Although many Bayesian nonparametric models have this property, we view it as a particular modeling choice that may be appropriate for some problems but not for others.
Several other infinite latent feature models have been developed to capture dependencies between data in different ways, for example using phylogenetic trees [21] or latent Gaussian processes [29]. Of particular relevance to this work is the model of Zhou et al. [31], which uses a hierarchical beta process to couple data. These and other related models are discussed further in Section 3. In Section 4, we characterize the featuresharing properties of the ddIBP and compare it to those of the model proposed by Zhou et al. [31]. We find that the different models capture qualitatively distinct dependency structures.
Exact posterior inference in the ddIBP is intractable. We present an approximate inference algorithm based on Markov chain Monte Carlo (MCMC; [27]) in Section 5, and we apply this algorithm in Section 6 to infer the latent features in a linearGaussian model. The experimental results presented in Section 7 suggest that the ddIBP is an effective tool for modeling latent structure in data with dependencies between observations.
2 The distance dependent Indian buffet process
We first review the definition of the IBP and its role in defining infinite latent feature models. We then introduce the ddIBP.
2.1 The Indian buffet process
The IBP is a prior over binary matrices with an infinite number of columns [16, 17]. In the Indian buffet metaphor, rows of correspond to customers and columns correspond to dishes. In data analysis, the customers represent data points and the dishes represent features. Let denote the entry of at row and column . Whether customer has decided to sample dish (that is, whether ) corresponds to whether data point possesses feature .
The IBP is defined as a sequential process. The first customer enters the restaurant and samples the first
number of dishes, where the hyperparameter
is a scalar. In the binary matrix, this corresponds to the first row being a contiguous block of ones, whose length is the number of dishes sampled (), followed by an infinite block of zeros.Subsequent customers enter, sampling each previously sampled dish according to its popularity,
(1) 
where is the number of customers that sampled dish prior to customer . (We emphasize that Eq. 1 applies only to dishes that were previously sampled, i.e., for which .) Then, each customer samples new dishes. Again these are represented as a contiguous block of ones in the columns beyond the last dish sampled by a previous customer.
Though described sequentially, Griffiths and Ghahramani [16] showed that the resulting rows of the binary matrix are exchangeable (up to a permutation of the columns). This means that the order of the customers does not affect the probability of the resulting binary matrix. This is seen in the BetaBernoulli perspective, which we review in Section 3. In the next section, we develop a generalization of the IBP that relaxes this assumption.
2.2 The distance dependent Indian buffet process
Like the IBP, the ddIBP is a distribution over binary latent feature matrices with a finite number of rows and an infinite number of columns. Each pair of customers has an associated distance, e.g., distance in time or space, or based on a covariate. Two customers that are close together in this distance will be more likely to share the same dishes (that is, features) than two customers that are far apart.
The ddIBP can be understood in terms of the following sequential construction. First, each customer selects a Poissondistributed number of dishes (feature columns). The dishes selected by a customer in this phase of the construction are said to be “owned” by this customer. A dish is either unowned, or is owned by exactly one customer. This step is akin to the selection of new dishes in the IBP.
Then, for each owned dish, customers connect to one another. The probability that one customer connects to another decreases in the distance between them. Note that customers do not sample each dish, as in the IBP, but rather connect to other customers.
Finally, dish inheritance is computed: A customer inherits a dish if its owner (from the first step) is reachable in the connectivity graph for that dish. This inheritance is computed deterministically from the connections generated in the previous step.^{1}^{1}1If one insists upon a complete gastronomical metaphor, customer connectivity can be thought of as “I’ll have what he’s having.” The dishes that each customer samples are those that he inherits or owns. Thus, similarity of sampled dishes between nearby customers is induced via distancedependent connection probabilities.
We now more formally describe the probabilistic generative process of the binary matrix Z. First, we introduce some notation and terminology.

Dishes (columns of ) are identified with the natural numbers . The set of dishes owned by customer is . The cardinality of this set is . These sets are disjoint, so for . The total number of owned dishes is . The set of dishes owned by customers excluding is .

Each dish is associated with a set of customertocustomer assignments, specified by the connectivity matrix , where indicates that customer connects to customer for dish . Given , the customers form a set of (possibly cyclic) directed graphs, one for each dish. The
ownership vector
is , where indicates the customer who owns dish , so . 
The distance matrix between customers is D, where the distance between customers and is . A customer’s selfdistance is 0: . We call the distance matrix sequential when for . In this special case, customers can only connect to previous customers.

The decay function maps distance to a quantity, which we call proximity, that controls the probabilities of customer links. We require that and . We obtain the normalized proximity matrix by applying the decay function to each customer pair and normalizing by customer. That is, , where .
Using this notation, we generate the feature indicator matrix as follows:

Assign dish ownership. For each customer , allocate unowned dishes to the customer’s set of owned dishes, . For each , set the ownership .

Assign customer connections. For each customer and dish , draw a customer assignment according to . Note that customers can connect to themselves. In this case, they do not inherit a dish unless they own it (see the next step).

Compute dish inheritance. We say that customer inherits dish if there exists a path along the directed graph for dish from customer to the dish’s owner (i.e., if is reachable from ), where the directed graph is defined by column of . The owner of a dish automatically inherits it.^{2}^{2}2Although customer can link to other customers for dish even if , these connections are ignored in determining dish inheritance when . We encode reachability with . If customer is reachable from customer for dish then . Otherwise .

Compute the feature indicator matrix. For each customer and dish we set if inherits , otherwise .
An example of customer assignments sampled from the ddIBP is shown in Figure 2. In this example, customer 1 owns dish 1; customers 24 all reach customer 1, either directly or through a chain, and thereby inherit the dish (indicated by gray shading). Consequently, feature 1 is active for customers 14. Dish 2 is owned by customer 2; only customer 1 reaches customer 2, and hence feature 2 is active for customers 1 and 2. Dish 3 is owned by customer 2, but no other customers reach customer 2, and hence feature 3 is active only for that customer.
The generative process of the ddIBP defines the following joint distribution of the ownership vector and connectivity matrix,
(2) 
Consider the first term. Recall that the set of dishes each customer owns and the total number of owned dishes are both functions of the ownership vector . Thus, the probability of the ownership vector is
(3) 
where is a deterministic function of .
Consider the second term. The conditional distribution of the connectivity matrix depends on the total number of owned dishes and the normalized proximity matrix (derived from the distances and decay function),
(4) 
The dependence on comes from , which is determined by .
Random feature models (and the traditional IBP) operate with a random binary matrix Z. In the ddIBP, Z
is a (deterministic) manytoone function of the random variables
and , which we denote by . We compute the probability of a binary matrix by marginalizing out the appropriate configurations of these variables(5) 
The ddIBP reduces to the standard IBP in the special case when for all and the distance matrix is sequential. (Recall: is sequential if for .) To see this, consider the probability that the th dish is sampled by the th customer (that is, ). This probability is the proportion of previous customers that already reach because the probabiity of connecting to each customer is proportional to one. This probability is , which is the same as in the IBP. This is akin to the relationship between the ddCRP and the traditional CRP under the same condition.
Many different decay functions are possible within this framework. Figure 3 shows samples of using four decay functions and a sequential distance defined by absolute temporal distance ( for and for ).

The constant, . This is the standard IBP.

The exponential, .

The logistic, .

The window, .
Each decay function encourages the sharing of features across nearby rows in a different way.
When combined with an observation model (which specifies how the latent features give rise to observed data), the ddIBP functions as a prior over latent feature representations of a data set. In section 6, we consider a specific example of how the ddIBP can be used to analyze data.
2.3 Marginal invariance and exchangeability
Unlike the traditional IBP, the ddIBP is not (in general) marginally invariant, the property that removing a customer leaves the distribution over latent features for the remaining customers unchanged. (The ddIBP builds on the ddCRP, which is not marginally invariant either.) In some circumstances, marginal invariance is desirable for computational reasons. For example, the conditional distributions over missing data for models lacking marginal invariance require computing ratios of normalization constants. In contrast, marginally invariant models, due to their factorized structure, require less computation for conditional distributions over missing data. In other circumstances, such as exploratory analysis of fully observed datasets, this computational concern is less important. Beyond computational considerations, while marginal invariance may be an appropriate modeling assumption in some data sets, it may be inappropriate in others.
Also unlike the traditional IBP, the ddIBP is not exchangeable in general. To state this formally, let be a permutation of the integers , and for a given binary matrix , let be the matrix created by permuting its rows according to . Let be drawn from the ddIBP with distance matrix , mass parameter and decay function . Then, except in certain special cases (such as when recovers the traditional IBP),
Permuting the data changes its distribution, and so the ddIBP is not exchangeable in general.
Although the ddIBP is not exchangeable, it does have a related symmetry. Let be the matrix with both its rows and its columns permuted according to , and let be drawn from the ddIBP with distance matrix rather than . (We retain the same values for and .) Then, in general,
Thus, if we permute both the data and the distance matrix, probabilities remain unchanged. Permuting both the data and the distance matrix is like first relabeling the data, and then explicitly altering the probability distribution to account for this relabeling. If the ddIBP were exchangeable, one would not need to alter the probability distribution to account for relabeling.
3 Related work
In this section we describe related work on infinite latent feature models that capture external dependence between the data. We focus on the most closely related model, which is the dependent hierarchical beta process (dHBP; [31]). As a prelude to describing the dHBP, we review the connection between the IBP and the beta process.
3.1 The beta process
Recall that the IBP is exchangeable. Consequently, by de Finetti’s theorem [3], the rows of , considered as binary vectors , are conditionally independent,
(6) 
In this marginal distribution, is a random measure on the feature space and is the de Finetti mixing distribution (see [3]). Thibaux and Jordan [28] showed that the de Finetti mixing distribution underlying the IBP is the beta process (BP), parameterized by a positive concentration parameter and a base measure on . A draw is defined by a countably infinite collection of weighted atoms,
(7) 
where is a probability distribution that places a single atom at , and the are independent random variables whose distribution is described as follows. If is continuous, then the atoms and their weights are drawn from a nonhomogeneous Poisson process defined on the space with rate measure
(8) 
If is discrete and of the form , , then has atoms at the same locations as , with . Following Thibaux and Jordan [28], we define the mass parameter as . Note that is not necessarily a probability measure, and hence can take on nonnegative values different from 1.
Conditional on a draw from the beta process, the feature representation of data point is generated by drawing from the Bernoulli process (BeP) with base measure : . If is discrete, then , where . In other words, feature is activated with probability independently for all data points. Sampling from the compound betaBernoulli process is equivalent to sampling directly from the IBP when and [28].
3.2 Dependent hierarchical beta processes
The dHBP [31] builds external dependence between data points using the above BP construction. The dependencies are induced by mixing independent BP random measures, weighted by their proximities .
The dHBP is based on the following generative process,
(9) 
This is equivalent to drawing from a Bernoulli process whose base measure is a linear combination of BP random measures,
(10) 
Dependencies between data points are captured in the dHBP by the proximity matrix , as in the ddIBP.^{3}^{3}3Zhou et al. [31] formalize dependencies in an equivalent manner using a normalized kernel function defined over pairs of covariates associated with the data points. This allows proximal data points (e.g., in time or space) to share more latent features than distant ones.
In Section 4
, we compare the featuresharing properties of the dHBP and ddIBP. Using an asymptotic analysis, we show that the ddIBP offers more flexibility in modeling the proportion of features shared between data points, but less flexibility in modeling uncertainty about these proportions.
3.3 Other nonexchangeable variants
Although still a nascent area of research, several other nonexchangeable priors for infinite latent feature models have been proposed. Williamson, Orbanz and Ghahramani [29] used a hierarchical Gaussian process to couple the latent features of data in a covariatedependent manner. They named this model the dependent Indian buffet process (dIBP). Their framework is flexible: It can couple columns of in addition to rows, while the ddIBP cannot. However, this flexibility comes at a computational cost during inference: Their algorithm requires sampling an extra layer of variables.
Miller, Griffiths and Jordan [21] proposed a “phylogenetic IBP” that encodes treestructured dependencies between data. DoshiVelez and Ghahramani [11] proposed a “correlated IBP” that couples data points and features through a set of latent clusters. Both of these models relax exchangeability, but they do not allow dependencies to be specified directly in terms of distances between data. Furthermore, inference for these models requires more intensive computation than does the standard IBP. The MCMC algorithm presented by Miller et al. [21] for the phylogenetic IBP involves both dynamic programming and auxiliary variable sampling. Similarly, the MCMC algorithm for the correlated IBP involves sampling latent clusters in addition to latent features. Our model also incurs extra computational cost relative to the traditional IBP due to the computation of reachability (quadratic in the number of observations); however, it permits a richer specification of the dependency structure between observations than either the phylogenetic or the correlated IBP.
Recently, Ren et al. [26] presented a novel way of introducing dependency into latent feature models based on the beta process. Instead of defining distances between customers, each dish is associated with a latent covariate vector, and distances are defined between each customer’s (observed) covariates and the dishspecific covariates. Customers then choose dishes with probability proportional to the customerdish proximity. This construction comes with a significant computational advantage for data sets where the time complexity is tied predominantly to the number of observations. The downside of this construction is that the MCMC algorithm used for inference must sample a separate covariate vector for each dish, which may scale poorly if the covariate dimensionality is high.
4 Characterizing featuresharing
In this section, we compare the featuresharing properties of the dHBP and ddIBP. Two data points share a feature if that feature is active for both (i.e., for and a given feature ). This analysis is useful for understanding the types of dependencies induced by the different models, and can help guide the choice of model and hyperparameter settings for particular data analysis problems. We consider an asymptotic regime in which the mass parameter is large ( for the ddIBP and for the dHBP), which simplifies featuresharing properties. Proofs of all propositions in this section may be found in the Appendix.
4.1 Featuresharing in the ddIBP
We first characterize the limiting distributional properties of featuresharing in the ddIBP as . We drop the feature index in the reachability indicator , writing it . We do this because features (that is, columns of the binary matrix ) are exchangeable under the ddIBP and, consequently, the distribution of the random vector is invariant across .
Let denote the number of features held by data point , and let denote the number of features shared by data points and , where .
Proposition 1.
Under the ddIBP,
(11)  
(12) 
The probabilities and depend strongly on the distribution of the connectivity matrix , but do not depend on the ownership vector , since is independent of dish ownership.
We derive the limiting properties of and from properties of the Poisson distribution. In this and following results, indicates convergence in distribution.
Corollary 1.
Let . and converge in distribution under the ddIBP to the following constants as :
(13)  
(14)  
(15) 
This corollary shows that the limiting fraction of shared features in the ddIBP is a constant that may be different for each pair of data points and . In contrast, we show below that the same limiting fraction under the dHBP is random, and takes one of two values. These two values are fixed, and do not depend upon the data points and .
4.2 Featuresharing in the dHBP
Here we characterize the limiting distributional properties of feature sharing in the dHBP as becomes infinitely concentrated (i.e., , analogous to ). In this limit, featuresharing is primarily attributable to dependency induced by the proximity matrix .
Proposition 2.
If is continuous, then under the dHBP,
(16)  
(17) 
We derive the limiting properties of and from properties of the Poisson distribution.
Corollary 2.
Let . Conditional on , and converge in distribution under the dHBP to the following constants as :
(18)  
(19)  
(20) 
Thus, the expected fraction of object ’s features shared with object , , is a factor of bigger when . As , this fraction goes to . As , it goes to 1. We can obtain the unconditional fraction by marginalizing over and :
Corollary 3.
Let . converges in distribution under the dHBP as to a random variable defined by
(21) 
where .
This corollary shows that as grows large, the fraction of shared features becomes one of two values (determined by and ), with a mixing probability determined by the dependency structure. Thus, the dHBP affords substantial flexibility in specifying the mixing probability (via ), but is constrained to two possible values of the limiting fraction.
4.3 Featuresharing in the IBP
For comparison, we briefly describe the featuresharing properties under the traditional IBP.
Under the traditional IBP, by exchangeability, and are equal in distribution to and . The first customer draws a number of dishes. The second customer then chooses whether to sample each of these dishes independently and with probability . Thus, the number of dishes sampled by both the first and second customers is .
This shows that, under the traditional IBP, as with ,
(22) 
4.4 Discussion
Using an asymptotic analysis, the preceding theoretical results show that the ddIBP and dHBP provide different forms of flexibility in specifying the way in which features are shared between data points. This asymptotic analysis takes the limit as the mass parameters and become large. This limit is taken for theoretical tractability, and removes much of the uncertainty that is otherwise present in these models. While such limiting ddIBP and dHBP models are not intended for practical use, their simplicity provides insight into behavior in nonasymptotic regimes.
Under the ddIBP, Corollary 1 shows that the modeler is allowed a great deal of flexibility in specifying the proportions of features shared by data points. Given a matrix specifying the proportion of features that are believed to be shared by pairs of data points, one can (if this matrix is sufficiently wellbehaved) design a distance matrix that causes the ddIBP to concentrate on the desired proportions. While the ddIBP cannot model an arbitrary modelerspecified matrix of proportions, the set of matrices that can be modeled is very large.
In contrast, under the dHBP, Corollary 3 shows that the modeler has less flexibility in specifying the proportions of features shared. Under the dHBP, the modeler chooses two values, and , and the proportion of features shared by any pair of data points in the asymptotic regime must be one of these two values.
Section 4.3 shows that the traditional IBP has the least flexibility. In the asymptotic regime, the proportion of features shared by each pair of data points is a constant.
While the ddIBP has more flexibility in specifying values of the featuresharingproportions than the dHBP, it has less flexibility (at least in this asymptotic regime) in modeling uncertainty about these featuresharing proportions. Under the ddIBP, the proportion of features shared by a pair of data points in the asymptotic regime is a deterministic quantity. Under the dHBP, the proportion of features shared is a random quantity, even in the asymptotic regime. A modeler using the dHBP has full flexibility in choosing the joint probability distribution governing these proportions. One could extend the ddIBP to allow uncertainty about the featuresharingproportions by specifying a hyperprior over distance matrices, but we do not consider this extension further.
Figure 4 illustrates the difference in asymptotic featuresharing behavior between the dHBP and ddIBP. Subfigures in the upper row are draws from the dHBP, and subfigures in the bottom row are draws from the ddIBP. Within a single subfigure, the shade in the cell is the fraction . (The diagonals have been set to to bring out other aspects of the matrix.) Each of the four columns represents a pair of independent draws. To approximate the asymptotic regime considered by the theory, the mass parameters for the two models are set to large values of . The figure shows that, in draws from the dHBP, offdiagonal cells have one of two shades, corresponding to the two possible limiting values for . In the different columns, correspoindng to different independent draws, the patterns are different, showing that remains random under the dHBP, even in the asymptotic regime. In contrast, in draws from the ddIBP, offdiagonal cells take a variety of different values, but remain unchanged across independent draws.
Figure 5 illustrates nonasymptotic featuresharing behavior in a simple setting with only two data points. The figure shows the featuresharing behavior of the dHBP (top) and ddIBP (bottom) at two values for the mass parameter: (top row) and (bottom row). Each subfigure shows the probability mass function as a function the proximity , where for the ddIBP. Because there are only two data points, with and , specifying is sufficient for specifying the full proximity matrix . For the dHBP, we set and . Also facilitating comparison, is the same between both models (when ).
Figure 5 shows that as the proximity increases to , the number of shared features tends to increase under both models. More precisely, concentrates on larger values of as increases. However, the way in which the probability mass functions change with differs between the two models. In the ddIBP, the most likely value of increases smoothly, while under the dHBP it remains roughly constant and then jumps. As one varies across its full range, the set of most likely values for under the ddIBP spans its full range from to , while under the dHBP the most likely value for takes only a few values. Instead, varying under the dHBP allows a variety of bimodal distributions centered near the values from the asymptotic analysis.
This difference in nonasymptotic behaviors mirrors the difference between the two models in the asymptotic regime, where the ddIBP allows featuresharingproportions to be specified almost arbitrarily but allows little flexibility in modeling uncertainty about them, and the dHBP limits the number of possible values for the featuresharing proportions, but allows uncertainty over these values.
5 Inference using Markov chain Monte Carlo sampling
Given a dataset and a latent feature model with parameter , the goal of inference is to compute the joint posterior over the customer assignment matrix , the ddIBP hyperparameter , and likelihood parameter , as given by Bayes’ rule:
(23) 
where the first term is the likelihood (recall that is a deterministic function of and ), the second term is the prior over parameters, the third term is the ddIBP prior over the connectivity matrix , the fourth term is the prior over the ownership vector , and the last term is the prior over . In what follows, we assume that is conditionally independent of and for given and .
Exact inference in this model is computationally intractable. We therefore use MCMC sampling [27] to approximate the posterior with samples. The algorithm can be adapted to different datasets by choosing an appropriate likelihood function. In the next section, we show how to adapt this algorithm to a simple linearGaussian model.^{4}^{4}4Matlab software implementing this algorithm is available at the first author’s homepage: www.princeton.edu/~sjgershm
Our algorithm combines Gibbs and Metropolis updates. For Gibbs updates, we sample a variable from its conditional distribution given the current states of all the other variables. Conjugacy allows simple Gibbs updates for and . Because the ddIBP prior is not conjugate to the likelihood, we use the Metropolis algorithm to sample and . We generate proposals for and , and then accept or reject them based on the likelihood ratio. We further divide these updates into two cases: updates for “owned” (active) dishes and updates of dish ownership.
Sampling . To sample the likelihood parameter , we draw from the following conditional distribution:
(24) 
where the prior and likelihood are problemspecific. To obtain a closedform expression for this conditional distribution, the prior and likelihood must be conjugate. For nonconjugate priors, one can use alternative updates, such as MetropolisHastings or slice sampling
[27]. Generally, updates for will be decomposed into separate updates for each component of . In some cases, can be marginalized analytically; an example is presented in the next section.Sampling . To sample the hyperparameter , we draw from the following conditional distribution:
(25) 
where is determined by and the prior on
is a Gamma distribution with shape
and inverse scale . Using the conjugacy of the Gamma and Poisson distributions, the conditional distribution over is given by:(26) 
Sampling assignments for owned dishes. We update customer assignments for owned dishes (corresponding to “active” features) using Gibbs sampling. For , , and , we draw a sample from the conditional distribution over given the current state of all the other variables:
(27) 
where is the th row of , is the th row of , and is excluding row .^{5}^{5}5We rely on several conditional independencies in this expression; for example, is conditionally independent of given , and . The first factor in Eq. 27 is the likelihood,^{6}^{6}6In calculating the likelihood, we only include the active columns of (i.e., those for which ). and the second factor is the prior, given by . In considering possible assignments of , one of two scenarios will occur: Either data point reaches the owner of (in which case feature becomes active for as well as for all other data points that reach ), or it does not (in which case feature becomes inactive for as well as for all other data points that reach ). This means we only need to consider two different likelihoods when updating .
Sampling dish ownership. We update dish ownership and customer assignments for newly owned dishes (corresponding to features going from inactive to active in the sampling step) using Metropolis sampling. Both a new connectivity matrix and ownership vector are proposed by drawing from the prior, and then accepted or rejected according to a likelihood ratio. In more detail, the update proceeds as follows.

Propose for each data point .

Set . Then populate or depopulate it by performing, for each ,

If , insert new dishes into . Then, for all and , sample according to .

If , remove randomly selected dishes from .
This reallocation of dishes induces a new ownership vector .


Compute the acceptance ratio . Because the prior (conditional on the current state of the Markov chain) is being used as the proposal distribution, the acceptance ratio reduces to a likelihood ratio (the prior and proposal terms cancel out):
(28) 
Draw . Set and if , otherwise leave and unchanged.
Iteratively applying these updates, the sampler will (after a burnin period) draw samples from a distribution that approaches the posterior (Eq. 23) as the burnin period grows large. The time complexity of this algorithm is dominated by the reachability computation, , and the likelihood computation, which is if coded naively (see [16] for a more efficient implementation using rankone updates).
6 A linearGaussian model
As an example of how the ddIBP can be used in data analysis, we incorporate it into a linearGaussian latent feature model (Figure 6). This model was originally studied for the IBP by Griffiths and Ghahramani [16, 17]. The observed data consist of objects, each of which is a dimensional vector of realvalued object properties. We model as a linear combination of binary latent features corrupted by Gaussian noise:
(29) 
where is a matrix of realvalued weights, and is a
matrix of independent, zeromean Gaussian noise terms with standard deviation
. We place a zeromean Gaussian prior on with covariance . Intuitively, the weights capture how the latent features interact to produce the observed data. For example, if each latent feature corresponds to a person in an image, then the weight captures the contribution of person to pixel .Within the algorithm of the previous section, . As a consequence of our Gaussian assumptions, can be marginalized analytically, yielding the likelihood:
(30) 
where tr() is the matrix trace and . In calculating the likelihood, we only include the “active” columns of (i.e., those for which ), and is the number of active columns.
7 Experimental results
In this section we report experimental investigations of the ddIBP and comparisons with alternative models. We first show how the ddIBP can be used as a dimensionality reduction preprocessing technique for classification tasks when the data points are nonexchangeable. We then show how the ddIBP can be applied to missing data problems.
7.1 Dimensionality reduction for classification
The performance of supervised learning algorithms is often enhanced by preprocessing the data to reduce its dimensionality
[4]. Classical techniques for dimensionality reduction, such as principal components analysis and factor analysis, assume exchangeability, as does the infinite latent feature model based on the IBP [16]. For this reason, these techniques may not work as well for preprocessing nonexchangeable data, and this may adversely affect their performance on supervised learning tasks.We investigated this hypothesis using a magnetic resonance imaging (MRI) data set collected from 27 patients with Alzheimer’s disease and 35 healthy controls [8].^{7}^{7}7Available at: http://wiki.stat.ucla.edu/socr/index.php/SOCR_Data_July2009_ID_NI. The observed features consist of 4 structural summary statistics measured in 56 brain regions of interest: (1) surface area; (2) shape index; (3) curvedness; (4) fractal dimension. The classification task is to sort individuals into Alzheimer’s or control classes based on their observed features.
Agerelated changes in brain structure produce natural declines in cognitive function that make diagnosis of Alzheimer’s disease difficult [13]. Thus, it is important to take age into account when designing predictive models. For the ddIBP and dHBP, age is naturally incorporated as a covariate over which we constructed a distance matrix. Specifically, we defined as the absolute age difference between subjects and , with for (i.e., the distance matrix is sequential). This induces a prior belief that individuals with similar ages tend to share more latent features. In the MRI data set, ages ranged from 60 to 90 (median: ).
In detail, we ran 1500 iterations of MCMC sampling on the entire data set using the linearGaussian observation model, and then selected the latent features of the maximum a posteriori sample as input to a supervised learning algorithm (
regularized logistic regression, with the regularization constant set to
). Training was performed on half of the data, and testing on the other half.^{8}^{8}8A few randomly chosen individuals were removed from the test set to make it balanced. The noise hyperparameters of the ddIBP and dHBP ( and ) were updated using MetropolisHastings proposals.We monitored the log of the joint distribution . Visual inspection of the log joint probability traces suggested that the sampler reaches a local maximum within 400500 iterations (Figure 7). This process was repeated for a range of decay parameter () values, using the exponential decay function. The same proximity matrix, , was used for both the ddIBP and dHBP. We performed 5 random restarts of the sampler and recomputed the classification measure for each restart, averaging the resulting measures to reduce sampling variability. For comparison, we also made predictions using the standard IBP, the dIBP [29], and the raw observed features (i.e., no preprocessing). The dIBP was fit using the MCMC algorithm described in Williamson et al. [29], which adaptively samples the parameters controlling dependencies between observations (thus the results do not depend on ).
Classification results are shown in Figure 8
(left), where performance is measured as the area under the receiver operating characteristic curve (AUC). Chance performance corresponds to an AUC of 0.5, perfect performance to an AUC of 1. For a range of
values, the ddIBP produces superior classification performance to the alternative models, with performance increasing as a function of . The dHBP performs worse relative to the raw data for lowvalues. The magnitude of the standard error (across random restarts) is roughly
that of the means.We also ran the ddIBP sampler with (in which case the ddIBP and IBP are equivalent) and found no significant difference between it and the standard IBP sampler with respect to performance on the Alzheimer’s classification and the EEG reconstruction (see next section).
7.2 Reconstructing missing data
As an example of a missing data problem, we use latent feature models to reconstruct missing observations in electroencephalography (EEG) time series. The EEG data^{9}^{9}9Available at: http://mmspl.epfl.ch/page33712.html are from a visual detection experiment in which human subjects were asked to count how many times a particular image appeared on the screen [18]. The data were collected as part of a larger effort to design braincomputer interfaces to assist physically disabled subjects.
Distance between data points was defined using the absolute timedifference. Data were zscored prior to analysis. For 10 of the data points, we removed 2 of the observed features at random. We then ran the MCMC sampler for 1500 iterations, adding Gibbs updates for the missing data by sampling from the observation distribution (Eq.
30) conditional on the current values of the latent features and hyperparameters. We then used the MAP sample for reconstruction. We measured performance by the squared reconstruction error on the missing data. Figure 9 shows the reconstruction results, demonstrating that the ddIBP is effective for reconstructing missing data in this dataset and achieves lower reconstruction error than the alternative models we consider.8 Conclusions
By relaxing the exchangeability assumption for infinite latent feature models, the ddIBP extends their applicability to a richer class of data. We have shown empirically that this innovation fares better than the standard IBP on nonexchangeable data (e.g., timeseries).
We note that the ddIBP is not a standard Bayesian nonparametric distribution, in the sense of arising from a de Finetti mixing distribution. For the standard IBP, the de Finetti mixing distribution has been identified as the beta process [28], but this result does not generalize to the ddIBP due to its nonexchangeability (a consequence of de Finetti’s theorem). Nonetheless, this does not detract from our model’s ability to let the data infer the number of latent features, a property that it shares with other infinite latent feature models.
A number of future directions are possible. First, it may be possible to exploit distance dependence to derive more efficient samplers. In particular, DoshiVelez and Ghahramani [10] have shown that partitioning the data into subsets enables faster Gibbs sampling for the traditional IBP; the window decay function imposes a natural partition of the data into conditionally independent subsets.
Second, application of the ddIBP to other likelihood functions is straightforward. For example, it could be applied to relational data [20, 22] or text data [28]. As pointed out by Miller et al. [22], covariates like age or location often play an important role in link prediction. Whereas Miller et al. [22] incorporated covariates into the likelihood function, one could instead incorporate them into the prior by defining covariatebased distances between data points (e.g., the age difference between two people). A distinction of the latter approach is that it would allow one to model dependencies in terms of latent features. For instance, two people close in age or geographic location may be more likely to share latent interests, a pattern naturally captured by the ddIBP.
Third, modeling shared dependency structure across groups is important for several applications. In fMRI and EEG studies, for example, similar spatial and temporal dependencies are frequently observed across subjects. Modeling shared structure without sacrificing intersubject variability has been addressed with hierarchical models [2, 30]. One way to extend the ddIBP hierarchically would be to allow the parameters of the decay function to vary across individuals while being coupled together by higherlevel variables.
Acknowledgements
SJG was supported by a NSF graduate research fellowship. PIF acknowledges support from NSF Award . We thank Matt Hoffman, Chong Wang, Gungor Polatkan, Sean Gerrish and John Paisley for helpful discussions. We are also grateful to Sinead Williamson and Mingyuan Zhou for sharing their code.
References
 [1] A. Ahmed and E.P. Xing. Timeline: A dynamic hierarchical dirichlet process model for recovering birth/death and evolution of topics in text stream. In Proc. of UAI, 2010.
 [2] C.F. Beckmann, M. Jenkinson, and S.M. Smith. General multilevel linear modeling for group analysis in FMRI. NeuroImage, 20(2):1052–1063, 2003.
 [3] JM Bernardo and A.F.M. Smith. Bayesian Theory. John Wiley & Son Ltd, 1994.
 [4] C.M. Bishop. Pattern Recognition and Machine Learning. Springer, 2006.
 [5] D. Blackwell and J. MacQueen. Ferguson distributions via Pólya urn schemes. The Annals of Statistics, 1(2):353–355, 1973.
 [6] D. Blei and P.I. Frazier. Distance dependent Chinese restaurant processes. International Conference on Machine Learning, 2010.

[7]
F. Caron, M. Davy, and A. Doucet.
Generalized Polya urn for timevarying Dirichlet process mixtures.
In
Uncertainty in Artificial Intelligence
, 2007.  [8] N. Christou and I.D. Dinov. Confidence interval based parameter estimationâ€”a new socr applet and activity. PloS One, 6(5):e19178, 2011.
 [9] Y. Chung and D.B. Dunson. The local Dirichlet process. Annals of the Institute of Statistical Mathematics, 63(1):59–80, 2011.
 [10] F. DoshiVelez and Z. Ghahramani. Accelerated sampling for the Indian buffet process. In International Conference on Machine Learning, pages 273–280, 2009.
 [11] F. DoshiVelez and Z. Ghahramani. Correlated nonparametric latent feature models. Uncertainty in Artificial Intelligence, 2009.
 [12] J.A. Duan, M. Guindani, and A.E. Gelfand. Generalized spatial Dirichlet process models. Biometrika, 94(4):809–825, 2007.
 [13] T. Erkinjuntti, R. Laaksonen, R. Sulkava, R. Syrjäläinen, and J. Palo. Neuropsychological differentiation between normal aging, alzheimer’s disease and vascular dementia. Acta Neurologica Scandinavica, 74(5):393–403, 1986.
 [14] M.D. Escobar and M. West. Bayesian Density Estimation and Inference Using Mixtures. Journal of the American Statistical Association, 90(430):577–588, 1995.
 [15] J.E. Griffin and M.F.J. Steel. Orderbased dependent Dirichlet processes. Journal of the American statistical Association, 101(473):179–194, 2006.
 [16] T.L. Griffiths and Z. Ghahramani. Infinite latent feature models and the Indian buffet process. Advances in Neural Information Processing Systems, 18, 2005.
 [17] T.L. Griffiths and Z. Ghahramani. The Indian Buffet Process: An Introduction and Review. Journal of Machine Learning Research, 12:1185–1224, 2011.
 [18] U. Hoffmann, J.M. Vesin, T. Ebrahimi, and K. Diserens. An efficient P300based braincomputer interface for disabled subjects. Journal of Neuroscience Methods, 167(1):115–125, 2008.

[19]
D. Knowles and Z. Ghahramani.
Infinite sparse factor analysis and infinite independent components analysis.
Independent Component Analysis and Signal Separation, pages 381–388, 2007.  [20] E. Meeds, Z. Ghahramani, R.M. Neal, and S.T. Roweis. Modeling dyadic data with binary latent factors. Advances in Neural Information Processing Systems, 19, 2007.
 [21] K.T. Miller, T.L. Griffiths, and M.I. Jordan. The phylogenetic indian buffet process: A nonexchangeable nonparametric prior for latent features. 2008.
 [22] K.T. Miller, T.L. Griffiths, and M.I. Jordan. Nonparametric latent feature models for link prediction. Advances in Neural Information Processing Systems, 2009.
 [23] D.J. Navarro and T.L. Griffiths. Latent features in similarity judgments: A nonparametric Bayesian approach. Neural computation, 20(11):2597–2628, 2008.
 [24] V. Rao and Y.W. Teh. Spatial normalized gamma processes. Advances in Neural Information Processing Systems, 22:1554–1562, 2009.

[25]
C.E. Rasmussen.
The infinite Gaussian mixture model.
Advances in Neural Information Processing Systems, 12:554–560, 2000.  [26] L. Ren, Y. Wang, D. Dunson, and L. Carin. The kernel beta process. Advances in Neural Information Processing Systems, 2011.
 [27] C.P. Robert and G. Casella. Monte Carlo statistical methods. Springer Verlag, 2004.
 [28] R. Thibaux and M.I. Jordan. Hierarchical beta processes and the Indian buffet process. In International Conference on Artificial Intelligence and Statistics, volume 11, pages 564–571. Citeseer, 2007.
 [29] S. Williamson, P. Orbanz, and Z. Ghahramani. Dependent Indian buffet processes. International Conference on Artificial Intelligence and Statistics, 2010.

[30]
M.W. Woolrich, T.E.J. Behrens, C.F. Beckmann, M. Jenkinson, and S.M. Smith.
Multilevel linear modelling for FMRI group analysis using Bayesian inference.
NeuroImage, 21(4):1732–1747, 2004. 
[31]
M. Zhou, H. Yang, G. Sapiro, D. Dunson, and L. Carin.
Dependent hierarchical beta process for image interpolation and denoising.
International Conference on Artificial Intelligence and Statistics, 2011.
Appendix: proofs
Recall that is the number of features held by data point , and is the number of features shared by data points and , where .
Proposition 1
Under the ddIBP,
(31)  
(32) 
Proof.
For each feature, there is some probability that it is turned on for data point , and some probability that it is shared by and (note that features are exchangeable in the ddIBP). The total number of features across all data points is distributed according to , where . Since each feature is turned on independently, the total number of active features for a single data point is distributed according to . Similarly, the total number of features shared by data points and is . The activation and coactivation probabilities are given by:
(33)  
(34) 
where we have used that . ∎
Proposition 2
If is continuous, then under the dHBP,
(35)  
(36) 
Proof.
We write the random measures and in the generative model defining the dHBP in Section 3.2 as the following mixtures over point masses.
(37)  
(38) 
Recall that where .
Let be the random variable that is if the Bernoulli process draw has atom , and if not. We have . Because is continuous, for and the random variables and satisfy
(39) 
We first show that is Poisson distributed with mean .
Let denote the probability that has atom conditioned on (this value does not depend on ). That is, .
For a given , the density of conditioned on is:
(40) 
We can use this density to calculate the success probability :
(41) 
where we have used the tower property of conditional expectation in the second and third equalities.
For a given , let denote the number of atoms in with . This number is Poissondistributed with mean .
Let be the number of such atoms that are also in . Because is the sum of independent Bernoulli trials that each have success probability , it follows that and
(42) 
Because , it follows that , where
Comments
There are no comments yet.