Hierarchical Dirichlet scaling process

We present the hierarchical Dirichlet scaling process (HDSP), a Bayesian nonparametric mixed membership model. The HDSP generalizes the hierarchical Dirichlet process to model the correlation structure between metadata in the corpus and mixture components. We construct the HDSP based on the normalized gamma representation of the Dirichlet process, and this construction allows incorporating a scaling function that controls the membership probabilities of the mixture components. We develop two scaling methods to demonstrate that different modeling assumptions can be expressed in the HDSP. We also derive the corresponding approximate posterior inference algorithms using variational Bayes. Through experiments on datasets of newswire, medical journal articles, conference proceedings, and product reviews, we show that the HDSP results in a better predictive performance than labeled LDA, partially labeled LDA, and author topic model and a better negative review classification performance than the supervised topic model and SVM.


Introduction
The hierarchical Dirichlet process (HDP) is an important nonparametric Bayesian prior for mixed membership models, and the HDP topic model is useful for a wide variety of tasks involving unstructured text (Teh et al., 2006).To extend the HDP topic model, there has been active research in dependent random probability measures as priors for modeling the underlying association between the latent semantic structure and covariates, such as time stamps and spatial coordinates (Ahmed and Xing, 2010;Ren et al., 2011).
A large body of this research is rooted in the dependent Dirichlet process (DDP) (MacEachern, 1999) where the probabilistic random measure is defined as a function of covariates.Most DDP approaches rely on the generalization of Sethuraman's stick breaking representation of DP (Sethuraman, 1991), incorporating the time difference between two or more data points, the spatial difference among observed data, or the ordering of the data points into the predictor dependent stick breaking process (Duan et al., 2007;Dunson and Park, 2008;Griffin and Steel, 2006).Some of these priors can be integrated into the hierarchical construction of DP (Srebro and Roweis, 2005), resulting in topic models where temporally-or spatially-proximate data are more likely to be clustered.These existing DP approaches, however, cannot model datasets with various types of covariates, including categorical and numerical labels.One reason is that categorical labels cannot be used to directly define the similarity between two documents, unlike temporal or spatial information.Also, labels and documents do not have a one-to-one correspondence, as there may be zero, one, or more labels per document.Furthermore, existing DP approaches cannot be applied to datasets with more than one type of covariates, for example numerical and categorical labels.
We suggest the hierarchical Dirichlet scaling process (HDSP) as a new way of modeling a corpus with various types of covariates such as categories, authors, and numerical ratings.The HDSP models the relationship between topics and covariates by generating dependent random measures in a hierarchy, where the first level is a Dirichlet process, and the second level is a Dirichlet scaling process (DSP).The first level DP is constructed in the traditional way of a stick breaking process, and the second level DSP with a normalized gamma process.With the normalized gamma process, each topic proportion of a document is independently drawn from a gamma distribution and then normalized.Unlike the stick breaking process, the normalized gamma process keeps the same order of the atoms as the first level measure, which allows the topic proportions in the random measure to be controlled.The DSP then uses that controllability to guide the topic proportions of a document by replacing the rate parameter of the gamma distribution with a scaling function that defines the correlation structure between topics and labels.The choice of the scaling function reflects the characteristics of the corpus.We show two scaling functions, the first one for a corpus with categorical labels, and the second for a corpus with both categorical and numerical labels.
The HDSP models the topic proportions of a document as a dependent variable of observable side information.This modeling approach differs from the traditional definition of a generative process where the observable variables are generated from a latent variable or parameter.For example, Zhu et al. (2009) and Mcauliffe and Blei (2007) propose generative processes where the observable labels are generated from a topic proportion of a document.However, a more natural model of the human writing process is to decide what to write about (e.g., categories) before writing the content of a document.This same approach is also successfully demonstrated in Mimno and McCallum (2012).
The outline of this paper is as follows.In Section 2, we describe related work and position our work within the topic modeling literature.In Section 3, we describe the gamma process construction of the HDP and how scale parameters are used to develop the HDSP with two different scaling functions.In Section 4, we derive a variational inference for the latent variables.In Section 5, we verify our approach on a synthetic dataset and demonstrate the improved predictive power on real world corpora.In Section 6, we discuss our conclusions and possible directions for future work.

Related Work
For model construction, the model most closely related to HDSP is the discrete infinite logistic normal (DILN) model (Paisley et al., 2012) in which the correlations among topics are modeled through the normalized gamma construction.DILN allocates a latent location for each topic in the first level, and then draws the second level random measures from the normalized gamma construction of the DP.Those random measures are then scaled by an exponentiated Gaussian process defined on the latent locations.DILN is a nonparametric counterpart of the correlated topic model (Blei and Lafferty, 2007) in which the logistic normal prior is used to model the correlations between topics.The HDSP is also constructed through the normalized gamma distribution with an informative scaling parameter, but our goal in HDSP is to model the correlations between topics and labels.
The Dirichlet-multinomial regression topic model (DMR-TM) (Mimno and McCallum, 2012) also models the label dependent topic proportions of documents, but it is a parametric model.The DMR-TM places a log-linear prior on the parameter of the Dirichlet distribution to incorporate arbitrary types of observed labels.The DMR-TM takes the "upstream" approach in which the latent variable or latent topics are conditionally generated from the observed label information.The author-topic model (Rosen-Zvi et al., 2004) also takes the same approach, but it is a specialized model for authors of documents.Unlike the "downstream" generative approach used in the supervised topic model (Mcauliffe and Blei, 2007), the maximum margin topic model (Zhu et al., 2009), and the relational topic model (Chang and Blei, 2009), the upstream approach does not require specifying the probability distribution over all possible values of observed labels.
The HDSP is a new way of constructing a dependent random measure in a hierarchy.In the field of Bayesian nonparametrics, the introduction of DDP (Sethuraman, 1991) has led to increased attention in constructing dependent random measures.Most such approaches develop priors to allow covariate dependent variation in the atoms of the random measure (Gelfand et al., 2005;Rao and Teh, 2009) or in the weights of atoms (Griffin and Steel, 2006;Duan et al., 2007;Dunson and Park, 2008).These priors replace the first level of the HDP to incorporate a document-specific covariate for generating a dependent topic proportion.These approaches focus on the spatial distances or the ordering of the covariate, so they cannot be generalized for arbitrary types of label information.The HDSP, on the other hand, can model any types of labels.The HDSP allows covariate dependent variation in the weights of atoms, where the variation is controlled by the scaling function that defines the correlation between atoms and labels.A proper definition of the scaling function gives the flexibility to model various types of labels.
Several topic models for labeled documents use the credit attribution approach where each observed word token is assigned to one of the observed labels.Labeled LDA (L-LDA) allocates one dimension of the topic simplex per label and generates words from only the topics that correspond to the labels in each document (Ramage et al., 2009).An extension of this model, partially labeled LDA (PLDA), adds more flexibility by allocating a predefined number of topics per label and including a background label to handle documents with no labels (Ramage et al., 2011).The Dirichlet process with mixed random measures (DP-MRM) is a nonparametric topic model which generates an unbounded number of topics per label but still excludes topics from labels that are not observed in the document (Kim et al., 2012).

Hierarchical Dirichlet Scaling Process
In this section, we describe the hierarchical Dirichlet scaling process (HDSP).First we review the HDP with an alternative construction using the normalized gamma process construction for the second level DP.We then present the HDSP where the second level DP is replaced by Dirichlet scaling process (DSP).Finally, we describe two scaling functions for the DSP to incorporate categorical and numerical labels.

The normalized gamma process construction of HDP
The HDP1 consists of two levels of the DP where the random measure drawn from the upper level DP is the base distribution of the lower level DP.The formal definition of the hierarchical representation is as follows: where H is a base distribution, α, and β are concentration parameters for each level respectively, and index m represents multiple draws from the second level DP.For the mixed membership model, x mn , observation n in group m, can be drawn from where f (•) is a data distribution parameterized by θ.In the context of topic models, the base distribution H is usually a Dirichlet distribution over the vocabulary, so the atoms of the first level random measure G 0 are an infinite set of topics drawn from H. The second level random measure G m is distributed based on the first level random measure G 0 , so the second level shares the same set of topics, the atoms of the first level random measure.
The constructive definition of the DP can be represented as a stick breaking process (Sethuraman, 1991), and in the HDP inference algorithm based on stick breaking, the first level DP is given by the following conditional distributions: where V k defines a corpus level topic distribution for topic φ k .The second level random measures are conditionally distributed on the first level discrete random measure G 0 : where the second level atom θ ml corresponds to one of the first level atoms φ k .This stick breaking construction is the most widely used method for the hierarchical construction (Wang et al., 2011;Teh et al., 2006).
An alternative construction of the HDP is based on the normalized gamma process (Paisley et al., 2012).While the first level construction remains the same, the gamma process changes the second level construction from Eq. 4 to where Gamma(x; a, b) = b a x (a−1) e −bx /Γ(a).Unlike the stick breaking construction, the atom of the π mk of the gamma process is the same as the atom of the kth stick of the first level.Therefore, during inference, the model does not need to keep track of which second level atoms correspond to which first level atoms.Furthermore, by placing a proper random variable on the rate parameter of the gamma distribution, the model can infer the correlations among the topics (Paisley et al., 2012) through the Gaussian process (Rasmussen and Williams, 2005).The normalized gamma process itself is not an appropriate construction method for the approximate posterior inference algorithm based on the variational truncation method (Blei and Jordan, 2006) because, unlike the stick breaking process, the probability mass of a random measure constructed by the normalized gamma process is not limited to the first few number of atoms.But once the base distribution of second level DP is constructed by the stick breaking process of first level DP, the total mass of the second level base distribution G 0 is limited to the first few number of atoms, and then the truncation based posterior inference algorithm approximates the true posterior of the normalized gamma construction.

Hierarchical Dirichlet scaling process
The HDSP generalizes the HDP by modeling mixture proportions dependent on covariates.As a topic model, the HDSP assumes that topics and labels are correlated, and the topic proportions of a document are proportional to the correlations between the topics and the observed labels of the document.We develop the Dirichlet scaling process (DSP) with the normalized gamma construction of the DP, where the rate parameter of the gamma distribution is replaced by the scaling function.This scaling function serves the central role of defining the correlation structure between a topic and labels.Formally, the HDSP  consists of DP and DSP in a hierarchy: where the first level random measure G 0 is drawn from the DP with concentration parameter α and base distribution H.The second level random measure G m for document m is drawn from the DSP parameterized by the concentration parameter β, base distribution G 0 , observed labels of document r m , and scaling function s(•) with scaling parameter w.
As in the HDP, the first level of HDSP is a DP where the base distribution is the product of two distributions for data distribution and scaling parameter w.Specifically, the base distribution H is Dir(η) ⊗ L w where η is the parameter of the word-topic distribution, and L w is a prior distribution for the scaling parameter w.The form of the resulting random measure is where p k is the stick length for topic k, At the second level construction, w k becomes the parameter to guide the proportion of topic k's for each document.
At the second level of HDSP, label-dependent random measures are drawn from the DSP.First, as in the HDP, draw a random measure G m ∼ DP(β, H) for document m.Second, scale the weights of the atoms based on a scaling function parameterized by w k and the observed labels.Let r mj be the value of observed label j in document m, then G m is scaled as follows: where s w k (•) is the scaling function parameterized by the scaling parameter w k .Topic k is scaled by the scaling weight, s w k (r mj ), and therefore, the topic proportions of a document is proportional to the scaling weights of the observed labels.The scaling function should be carefully chosen to reflect the underlying relationship between topics and labels.We show two concrete examples of scaling functions in Section 3.3.The constructive definition of HDSP is similar to the HDP, but the difference comes from the scaling function.The stick breaking process is used to construct the first level random measure: where the pair {φ k , w k } drawn i.i.d.from two base distributions forms an atom of the resulting measure.
Based on the discrete first level random measure, the second level random measure is constructed by the normalized gamma process.As in the HDP, the weight of atom k is drawn from a gamma distribution with parameter βp k , and then scaled by the scaling weight The scaling weight can be directly incorporated into the second parameter of the gamma distribution because the scaled gamma random variable y = kx ∼ Gamma(a, 1) is equal to y ∼ Gamma(a, k −1 ), Then, the random variables are normalized to form a proper probability random measure For the mixed membership model, nth observation in mth group is drawn as follows: where f is a data distribution parameterized by φ k .For topic modeling, G m and x mn correspond to document m and word n in document m, respectively.

Scaling functions
Now we propose two scaling functions to express the correlation between topics and labels of documents.A scaling method is properly defined by two factors: 1) a proper prior over the scaling parameter w k , 2) a plausible scaling function between topic specific scaling parameter w k and the observed labels of document r m .
Scaling function 1: We design the first scaling function to model categorical side information such as authors, tags, and categories.For a corpus with J unique labels, then w k is a J-dimensional parameter where each dimension matches to a corresponding label.We define the scaling function as the product of scaling parameters that correspond to the observed labels: where r mj is an indicator variable whose value is one when label j is observed in document m and zero otherwise.w kj is a scaling parameter of topic k for label j.We place a inverse gamma prior over the weight variable w kj .With this scaling function, the proportion of topic k for document m is scaled as follows: The scaled gamma distribution is equal to the gamma distribution with the rate parameter of inverse scaling factor, so we can rewrite the above equation as follows: Finally, we normalize these random variables to make a probabilistic random measure summed up to unity for document m: Scaling function 2: The above scaling function models categorical side information, but many datasets, such as product reviews have numerical ratings as well as categorical information.We propose the second scaling function that can model both numerical and categorical information.Again, let w k be J-dimensional scaling parameter where each dimension matches to a corresponding label.The second scaling function is defined as follows: where w kj is the scaling parameter of label j for topic k, and r mj is the observed value of label j of document m.We place a normal prior over the scaling parameter w k .The scaling function is an inverse log-linear to the weighted sum of document's labels.Unlike the previous scaling function which only considers whether a label is observed in a document, this scaling function incorporates the value of the observed label.With this scaling function, the proportion of topic k for document m is scaled as follows Again, we can rewrite this equations as π mk is proportional to the inverse weighted sum of observed labels.Again, we normalize π mk to construct a proper random measure.
The choice of scaling function reflects the modeler's perspective with respect to the underlying relationship between topics and labels.The first scaling function scales each topic by the product of the scaling parameters of the observed labels.This reflects the modeler's assumption that a document with a set of observed labels is likely to exhibit topics that have high correlation with all of the observed labels.With the second scaling function, the scaling weight changes exponentially as the value of label changes.This reflects the modeler's assumption that two documents with the same set of observed labels but with different values are likely to exhibit different topics.

HDSP as a dependent Dirichlet process
We can view the HDSP as an alternative construction of the hierarchical dependent Dirichlet process (DDP) via a hierarchy consisting of a stick breaking process and a normalized gamma process.Let us compare the HDSP approach to the general DDP approach for topic modeling.The formal definition of DDP is: where the resulting random measure G 0 is a function of some covariates.Using G 0 as the base distribution of a DP for a document with a covariate, the random measure corresponding to document m is constructed as follows: where G 0 (r m ) is the base distribution for the document with same covariate r m (Srebro and Roweis, 2005).
Similarly, the HDSP constructs a dependent random measure with covariates.However, unlike the DDP-DP approach, G 0 is no longer a function of covariates.The HDSP defines a single global random measure G 0 and then scales G 0 based on the covariates with the scaling function.The advantage of the HDSP is that it only requires a proper, but relatively simple, scaling function that reflects the correlation between covariates and topics, whereas the DDP requires a complex dependent process for different types of covariates (Griffin and Steel, 2006).

Variational Inference for HDSP
The posterior inference for Bayesian nonparametric models is important because it is intractable to compute the posterior over an infinite dimensional space.Approximation algorithms, such as marginalized MCMC (Escobar and West, 1995;Teh et al., 2006) and variational inference (Blei and Jordan, 2006;Teh et al., 2008), have been developed for the Bayesian nonparametric mixture models.We develop a mean field variational inference (Jordan et al., 1999;Wainwright and Jordan, 2008) algorithm for approximate posterior inference of the HDSP topic model.The objective of variational inference is to minimize the KL divergence between a distribution over the hidden variables and the true posterior, which is equivalent to maximizing the lower bound of the marginal log likelihood of observed data.
In this section, we first derive the inference algorithm for the first scaling function with a fully factorized variational family.Variational inference algorithms can be easily modularized with the fully factorized variational family, and the variation in a model only affects the update rules for the modified parts of the model.Therefore, for the second scaling function, we only need to update the part of the inference algorithm related to the new scaling function.

Variational inference for the first scaling function
For the first scaling function, we use a fully factorized variational distribution and perform a mean-field variational inference.There are five latent variables of interest: the corpus level stick proportion V k , the document level stick proportion π mk , the scaling parameter between topic and label w kj , the topic assignment for each word z mn , and the word topic distribution φ k .Thus the variational distribution q(z, π, V, w, φ) can be factorized into where the variational distributions are q(z mn ) = Multinomial(z mn |γ mn ) For the corpus level stick proportion V k , we use the delta function as a variational distribution for simplicity and tractability in inference steps as demonstrated in (Liang et al., 2007).Infinite dimensions over the posterior is a key problem in Bayesian nonparametric models and requires an approximation method.In variational treatment, we truncate the unbounded dimensionality to T by letting V T = 1.Thus the model still keeps the infinite dimensionality while allowing approximation to be carried out under the bounded variational distributions.
Using standard variational theory, we derive the evidence lower bound (ELBO) of the marginal log likelihood of the observed data where H(q) is the entropy for the variational distribution.By taking the derivative of this lower bound, we derive the following coordinate ascent algorithm.Document-level Updates: At the document level, we update the variational distribution for the topic assignment z mn and the document level stick proportion π mk .The update for q(z mn |γ mn ) is Updating q(π mk |a π mk , b π mk ) requires computing the expectation term E[ln T k=1 π mk ].Following Blei and Lafferty (2007), we approximate the lower bound of the expectation by using the first-order Taylor expansion, where the update for Then, the update for π mk is Note again r mj is equal to 1 when jth label is observed in mth document, otherwise 0.
Corpus-level Updates: At the corpus level, we update the variational distribution for the scaling parameter w kj , corpus level stick length V k and word topic distribution η ki .
The optimal form of a variational distribution can be obtained by exponentiating the variational lower bound with all expectations except the parameter of interest (Bishop and Nasrabadi, 2006).For w kj , we can derive the optimal form of variational distribution as follows where m = {m : r mj = 1} and j /j = {j : r mj = 1, j = j}.See Appendix A for the complete derivation.There is no closed form update for V k , instead we use the steepest ascent algorithm to jointly optimize where ψ(•) is a digamma function.Finally, the update for the word topic distribution q(φ k |η k ) is where i is a word index, and 1 is an indicator function (Blei et al., 2003).
The expectations of latent variables under the variational distribution q are

Variational inference for the second scaling function
Introducing a new scaling function requires a new approximation method.We first choose the part of ELBO which requires new treatment as the scaling function changes.From Equation 25, we take the terms that are related to the scaling function s: To update the scaling parameters, we need a proper prior and variational distribution.
For the second scaling function, the normal distribution with zero mean and variance σ is used as a prior of w kj , and the delta function is used as the variational distribution of w kj .Newton-Raphson optimization method are used to update the weight parameters.The Newton-Raphson optimization finds a stationary point of a function by iterating: where H(w k ) and ∂ ∂w k are the Hessian matrix and gradient at the point w old k .The lower bound with respect to the parameter w kj is, Then, the gradient and Hessian matrix of w kj are Because the w kj is depends on both label j and topic k, we iteratively update w kj until converged.
The update rules for the variational parameter of π mk and V k need to be changed to accommodate the change of scaling function.The variational parameters of π mk are approximated by using the first-order Taylor expansion, To update V k , we take the same approach in which the variational distribution is a delta function of current V k .Again, we use the steepest ascent algorithm to jointly optimize V k , and the gradient of The update rules for π mk and V k only requires the expectation of the log scaling function.
The update rules for the other parameters remain the same as the previous section.
Introducing a new scaling function requires a new inference algorithm, and this can be cumbersome.Once one defines a tractable expectation of the log of a scaling function, a recently suggested Black-box method (Ranganath et al., 2014) can be an alternative to update the function-specific parameters instead of deriving function-specific approximation algorithms.

Experiments
In this section, we describe how the HDSP performs with real and synthetic data.We fit the HDSP topic model with three different types of data and compare the results with several comparison models.First, we test the model with synthetic data to verify the approximate inference.Second, we train the model with categorical data whose label information is represented by binary values.Third, we train the model with mixed-type of data whose label information has both numerical and categorical values.

Synthetic data
There is no naturally-occurring dataset with the observable weights between topics and labels, so we synthesize data based on the model assumptions to verify our model and the approximate inference.First, we check the difference between the original topics and the inferred topics via simple visualization.Then, we focus on the differences between the inferred and synthetic weights.For all experiments with synthetic data, the datasets are generated by following the model assumptions with the first scaling function, and the posterior inferences are done with the first scaling function.We set the truncation level T at twice the number of topics.We terminate the variational inference when the fractional change of the lower bound falls below 10 −3 , and we average all results over 10 individual runs with different initializations.
With the first experiment, we show that HDSP correctly recovers the underlying topics and scaling parameter between topics and labels.For the dataset, we generate 500 documents using the following steps.We define five topics over ten terms shown in Figure 2 ).As the volume of space for locations increases, the mean absolute error also increases (left).However, the model preserves the relative weights between topics and labels, shown by the high and stabilized correlation between the original ordering and the recovered ordering of label-topic pairs in terms of the weights between the two (right).This is a key characteristic of the HDSP model which scales the mixture components according to the inverse of the weights.
and the scaling parameter of five topics and four labels shown in Figure 2(d).For each document, we randomly draw N m from the Poisson distribution and r mj from the Bernoulli distribution.The average length of a document is 20, and the average number of labels per document is 2. We generate topic proportions of corpus and documents by using Equations 10 and 13.For each word in a document, we draw the topic and the word by using Equation 14.We set both α and β to 1.
Figure 2 shows the results of the HDP and the HDSP on the synthetic dataset.Figure 2(b) and Figure 2(c) are the heat maps of topics inferred from each model.We match the inferred topics to the original topics using KL divergence between the two sets of topic distributions.There are no significant differences between the inferred topics of HDSP and HDP.In addition to the topics, HDSP infers the scaling parameters between topics and labels, which are shown in Figure 2(e).The results show that the relative differences between original scaling parameters are preserved in the inferred parameters through the variational inference.
With the second experiment, we show that the inferred parameters preserve the relative differences between labels and topics in the dataset.For this experiment, we generate 1,000 documents with ten randomly drawn topics from Dirichlet(0.1) with the vocabulary size of 20.To generate the weights between topics and labels, we randomly place the topics and labels into three dimensional euclidean space, and use the distance between a topic and label as a scaling parameter.The location of topics and labels are uniformly drawn from three dimensional euclidean space, so the total volume is x 3 , then we vary the x value from 1 to 20 for each experiment.As the volume of space increases, the potential scaling parameter between a topic and label increases, and the scaling effect on a topic proportion of document also increases.We compute the mean absolute error (MAE) and the spearman's rank correlation coefficient ρ between the original parameters and the inferred parameters.The spearman's ρ is designed to measure the ranking correlation of two lists.Figure 3 shows the results.The MAE increases as the volume of the space increases.However, spearman's ρ stabilizes, indicating that the relative differences are preserved even when the MAE increases.Since there are an infinite number of configurations of scaling parameters that generate the same expectation E[p(π m |βp, w j )] given π m and βp, preserving the relative differences verifies our model's capability of capturing the underlying structure of topics and labels.

Categorical data
We evaluate the performance of HDSP and compare it with the HDP, labeled LDA (L-LDA), partially labeled LDA (PLDA), and author-topic model (ATM).For the HDSP, we use both scaling functions and denote the model with the second scaling function as wHDSP.We use three multi-labeled corpora: RCV2 (newswire from Reuter's), OHSUMED3 (a subset of the Medline journal articles), and NIPS (proceedings of NIPS conference).For RCV and OHSUMED, we use multi-category information of documents as labels, and for NIPS, we use authors of papers as labels.The average number of labels per article is 3.2 for RCV, 5.2 for OHSUMED, and 2.4 for NIPS.Table 1 contains the details of the datasets.

Experimental settings
For the HDP and HDSP, we initialize the word-topic distribution with three iterations of LDA for fast convergence to the posterior while preventing the posterior from falling into a local mode of LDA and then reorder these topics by the size of the posterior word count.For all experiments, we set the truncation level T to 200.We terminate variational inference when the fractional change of the lower bound falls below 10 −3 , and we optimize all hyperparameters during inference except η.For the L-LDA and PLDA, we implement the collapsed Gibbs sampling algorithm.For each model, we run 5,000 iterations, the first 3,000 as burn-in and then using the samples thereafter with gaps of 100 iterations.For PLDA, we set the number of topics for each label to two and five (PLDA-2, PLDA-5).For For HDSP (first scaling function), wHDSP (second scaling function), L-LDA, ATM, and PLDA, the perplexity is measured given documents and observed labels.For HDP, the model only uses the words of the documents.The HDSP which, instead of excluding topics from unobserved labels, scales all topics according to observed labels, shows the best heldout perplexity.
the ATM, we set the number of topics to 50, 100, and 150.We try five different values for the topic Dirichlet parameter η: η = 0.1, 0.25, 0.5, 0.75, 1.0.Finally all results are averaged over 20 runs with different random initialization.We do not report the standard errors because they are small enough to ignore.

Evaluation metric
The goal of our model is to construct the dependent random probability measure given multiple labels.Therefore, our interest is to see the increments of predictive performance when the label information is given.
The predictive probability given label information for held-out documents are approximated by the conditional marginal, where D train = {x train , r train } is the training data, x is the vector of N words of a held-out document, r are the labels of the held-out document, z n is the latent topic of word n, and π k is the kth topic proportion of the held-out document.Since the integral is intractable, we approximate the probability where φk and πk are the variational expectations of φ k and π k given label r .This approximated likelihood is then used to compute the perplexity of the held-out document Lower perplexity indicates better performance.We also take the same approach to compute the perplexity for L-LDA, PLDA and HDP, but HDP does not use the labels of held-out documents.To measure the predictive performance, we leave 20% of the documents for testing and use the remaining 80% to train the models.

Experimental results
Figure 4 shows the predictive performance of our model against the comparison models.
For the OHSUMED and RCV corpora, both HDSP and wHDSP outperform all others.Among these models, L-LDA restricts the modeling flexibility the most; the PLDA relaxes that restriction by adding an additional latent label and allowing multiple topics per label.HDSP and wHDSP further increase the modeling flexibility by allowing all topics to be generated from each label.This is reflected in the results of predictive performance of the three models; L-LDA shows the worst performance, then PLDA, and HDSP and wHDSP show the lowest perplexity.For the NIPS data, we compare HDSP and wHDSP to ATM, and again, HDSP and wHDSP show the lowest perplexity.
To visualize the relationship between topics and labels, we embed the inferred topics and the labels into the two dimensional euclidean space by using multidimensional scaling (Kruskal, 1964) on the inferred parameters of HDSP.In Figure 5, we choose and display a few representative topics and authors from NIPS.For instance, Geoffrey Hinton and Yann LeCun are closely located to the neural network related topics such as 'learning, network error' and 'recognition character, network', and the reinforcement learning researcher Richard Sutton is closely located to the 'state, learning policy' topic.Figure 6 shows the embedded labels and topics from OHSUMED.The labels 'Preschool', 'Pregnancy', and 'Infant' are closely located to one another with similar topics.While the model explicitly models the Figure 5: Relative locations of observed labels (red) and latent topics (blue) inferred by HDSP from the NIPS corpus correlation between topics and labels, embedding them together shows that the correlation among labels, as well as among topics, can also be inferred.
Figure 7 and 8 show the expected topic distributions of HDSP given different sets of labels.When multiple labels are given, the model expects high probabilities for the topics that are similar to all given labels.For example, when 'Market' and 'Sports' labels are given, the model expects high probabilities on sports related topics and relatively high probability on 'Market' related topics based on the weights between topics and two labels.

Modeling data with missing labels
We also test our model with partially labeled data which have not been previously covered in topic modeling.Many real-world data fall into this category where some of the data are labeled, others are incompletely labeled, and the rest are unlabeled.For this experiment, we randomly remove existing labels from the RCV and OHSUMED corpora.To remove observed labels in the training corpus, we use Bernoulli trials with varying parameters to analyze how the proportion of observed labels affects the heldout predictive performance of the model.
Figure 9 shows the predictive perplexity with varying parameters of Bernoulli distribution from 0.1 to 0.9.For both scaling functions, the perplexity decreases as the model observes more labels.Compared to the PLDA (with the parameter setting for optimal performance), the HDSP achieves similar perplexity with only 20% of the labels.One notable phenomenon is that the HDSP outperforms wHDSP on both datasets when the number of observed labels is less than 50% of the total number of labels.

Mixed-Type data
In this section, we present the performance of second scaling function with a corpus of product reviews which has real-valued ratings and category information.
The first scaling function is only applicable to categorical side information, so we use the second scaling function (wHDSP) which can model numerical as well as categorical side information of documents.To evaluate the performance of wHDSP with numerical side information, we train the model with the Amazon review data collected from seven categories of electronic products: air conditioner, canister vacuum, coffee machine, digital SLR, laptop, MP3 player, and space heater.Amazon uses a five-star rating system, so each review contains one numerical rating ranging from one to five.Table 2 shows the number of reviews for each rating and category.Recall that r is a vector whose values denote the observation of the labels.For each review, we set the dimension of r to eight in which the first dimension is a numerical rating of a review, and then the remaining seven dimensions match the seven product categories.We set the value of each dimension to one if the review belongs to the corresponding category, and zero otherwise.
To evaluate the performance of wHDSP, we classify the ratings of the reviews based on a trained model.We use 90% of the corpus to train models and the remaining 10% of the corpus to test the models.To classify the rating of each review in the test set,    we compute the perplexity of the given review with varying ratings from one to five, and choose the rating that shows the lowest perplexity.Generally, computing the perplexity of heldout document requires complex approximation schemes (Wallach et al., 2009), but we compute the perplexity based on the expected topic distribution given category and rating information, which requires a finite number of computations.We compare the wHDSP with the supervised LDA (SLDA), LDA+SVM, as well as classifiers Naive Bayes, SVM, and decision trees (CART).For the LDA+SVM approach, we first train the LDA model and then use the inferred topic proportion and categories as features of the SVM.For the SLDA model, the category information cannot be used because the model is designed to learn and predict the single response variable.For both models, we set the number of topics to 50, 100, and 200.
In many applications, classifying negative feedback of users is more important than classifying positive feedback.From the negative feedback, companies can identify possible problems of their products and services and use the information to design their next product or improve their services.In most online reviews, however, the proportion of negative feedback is smaller than the proportion of positive feedback.For example, in the Amazon data, about 51% of reviews are rated as five-star, and 72% rated as four or five.A classifier trained by such skewed data is likely to be biased toward the majority class.We report the classification results of Amazon dataset in Table 3 and Table 4. Table 3 shows the results for each rating in terms of F1, and Table 4 shows the results in terms of micro and macro F1.The wHDSP outperforms the other models in terms of macro F1 but performs worse than sLDA in terms of micro F1.As we noted earlier, classifying negative reviews may be more important in many applications.Both the SLDA with 100 topics and the wHDSP are comparable in classifying the most negative (one-star) reviews.However, the confusion matrices and Table 3 indicate that the SLDA dichotomously learns the decision boundaries where the most reviews are classified into one-star or five-star.For example, the SLDA with 50 and 100 topics did not classify any two-star and three-star reviews correctly.The wHDSP learns the decision boundaries for classifying subtle differences between fivestar rating reviews.These patterns are shown clearly with the confusion matrices in Figure 10 where the diagonal entries are the numbers of correctly classified reviews.For example, the wHDSP classified only eight one-star reviews as five-star reviews, but the SLDA50 assigned 68 reviews as five-star reviews.
We perform a rating prediction task with and without the category information of reviews to see the effect of using both the category and rating information on the wHDSP and LDA+SVM approaches.The results represented by wHDSP* in Table 4 and Figure 10(b) show the performances of rating prediction with the wHDSP trained without category information.For wHDSP*, the model performs worse than wHDSP, which indicates the model, without category information, cannot distinguish the review ratings which depend  on topical context.The LDA+SVM without categories achieves 0.309 macro F1 and 0.533 micro F1, which are comparable to the LDA+SVM with the category information.Unlike the wHDSP, the decision boundaries of SVM are not improved with the additional category information.The result supports that for learning decision boundaries between ratings over different categories, the approach of including category information to train topics is more effective than using topics and the category information independently.

Discussions
We have presented the hierarchical Dirichlet scaling process (HDSP), a Bayesian nonparametric prior for a mixed membership model that lets us analyze underlying semantics and observable side information.The combination of the stick breaking process with the normalized gamma process in HDSP is a more controllable construction of the hierarchical Dirichlet process because each atom of the second level measure inherits from the first level measure in order.HDSP also allows more flexibility and the capability of modeling side information by the scaling functions that plug into the rate parameter of the gamma distribution.The choice of the scaling function is the most important part of the model in terms of establishing a link between topics and observed labels.We developed two scaling functions but the choice of scaling function depends on the modeler's intention.For example, the well known linking functions from the generalized linear model can be used as scaling functions, or one can use several scaling functions together on purpose.We showed that the application of HDSP to topic modeling correctly recovers the topics and topic-label weights of synthetic data.Experiments with the real dataset show that the first scaling function is more suited for partially labeled data, and the second first scaling function is more suited for a dataset with both numerical and categorical labels.Hierarchical Dirichlet scaling process opens up a number of interesting research questions that should be addressed in future work.First, in the two scaling functions we proposed to model the correlation structure between topics and side information, we simply defined the relationship between topic k and label j through the scaling parameter w kj .However, this approach does not consider the correlation within topics and labels.Taking inspiration from previous work (Blei and Lafferty, 2007;Mimno et al., 2007;Paisley et al., 2012) that showed correlations among topics, we can define a scaling function with a prior over the topics and labels to capture their complex relationships.Second, our posterior inference algorithm based on mean-field variational inference is tested with tens of thousands documents.However, modern data analysis requires inference of massive and/or streaming data.For a fast and efficient posterior inference, we can apply parallel or distributed algorithms based on a stochastic update (Hoffman et al., 2013;Ahn et al., 2014).Furthermore, we fix the number of labels before training but we need to find a way to model the unbounded number of labels for streaming data.(42) where D is the training set of documents and labels.Ψ denotes the set of variational parameters, Θ denotes the set of model parameters.
The optimal form of the variational distribution can be obtained by exponentiating the variational lower bound with all expectations except the parameter of interest (Bishop and Nasrabadi, 2006).For π mk , we can derive the optimal form of variational distribution as follows q(π mk ) ∝ exp E q −π mk [ln p(π mk |z, π −mk , r m , V )] (47) We take the same approach described in (Paisley et al., 2012), and the only difference comes from the product of the inverse distance term.
For w kj , we can derive the optimal form of the variational distribution as follows q(w kj ) ∝ exp E q −w kj Therefore, the optimal form of variational distribution for w kj is where m = {m : r mj = 1} and j /j = {j : r mj = 1, j = j}.

B. Posterior Word Count
Like the HDP and other nonparametric topic models, our model also uses only a few topics even though we set the truncation level to 200.
topic k and label j : k th topic : first and second level topic proportion : topic of n th word in m th document : j th label of m th document {0, 1} : observed n th word in m th document

Figure 1 :
Figure 1: Graphical model of the hierarchical Dirichlet scaling process.

Figure 2 :
Figure 2: Experiments with synthetic data.(a) is the synthetic topic distribution of 5 topics over 10 terms.(b) and (c) are topic distributions inferred by the HDSP and the HDP.Both models recover the original topics.(d) shows the heat map of original scaling parameters between the topics and labels.(e) shows the heat map of the recovered parameters by HDSP.

Figure 3 :
Figure 3: Spearman's correlation coefficient and mean absolute error of the synthetic datawith various volume of space (x 3 ).As the volume of space for locations increases, the mean absolute error also increases (left).However, the model preserves the relative weights between topics and labels, shown by the high and stabilized correlation between the original ordering and the recovered ordering of label-topic pairs in terms of the weights between the two (right).This is a key characteristic of the HDSP model which scales the mixture components according to the inverse of the weights.

Figure 4 :
Figure 4: Perplexity of held-out documents.For HDSP (first scaling function), wHDSP(second scaling function), L-LDA, ATM, and PLDA, the perplexity is measured given documents and observed labels.For HDP, the model only uses the words of the documents.The HDSP which, instead of excluding topics from unobserved labels, scales all topics according to observed labels, shows the best heldout perplexity.

Figure 6 :
Figure 6: Relative locations of observed labels (red) and latent topics (blue) inferred by HDSP from the OHSUMED corpus.

Figure 7 :
Figure 7: Expected topic distributions given labels from RCV and OHSUMED.Topics are sorted by their posterior word counts, and the top 20 topics are displayed with the top 10 words (stemmed).From top to bottom, we compute an expected topic distribution given a randomly selected set of labels.

Figure 8 :Figure 9 :
Figure 8: Expected topic distributions given labels from NIPS.Topics are sorted by their posterior word counts, and the top 20 topics are displayed with the top 10 words.From top to bottom, we compute an expected topic distribution given a set of representative labels.

Figure 10 :
Figure 10: Confusion matrices from classification results of wHDSP and the other models.Diagonal entries indicate the number of correctly classified reviews per rating.Despite the highly skewed data distribution (c), wHDSP achieves relatively better classification results for negative and neutral reviews as shown in (a).Confusion matrices are column-wise normalized.

Figure 11 :
Figure 11: Posterior number of words.We set the truncation level to 200, but only a few topics are used during inference.

Table 1 :
Datasets used for the experiments in 5.2.As the last two columns show, we experiment on datasets with a varied number of unique labels, as well as the average number of labels per document.

Table 2 :
The number of reviews for each rating and category in the Amazon dataset.The dataset contains 24,259 reviews collected from seven different product categories.The description reveals a highly skewed distribution of reviews where 52% of reviews are rated as five-star.The percentage denotes the proportion of reviews per each rating.

Table 3 :
F1 of wHDSP and the other models for the Amazon review corpus.wHDSP and SLDA perform comparably on one-star ratings but wHDSP outperforms SLDA on middle range ratings (two, three, and four stars).

Table 4 :
Macro and micro F1 of the wHDSP and the other models.The left table shows the performance of classification with five-star ratings system.The models with asterisk are trained without category information.Note that the SLDA cannot incorporate two different types of labels together.