Supervised Topic Models

David M. BleiJon D. McAuliffe

article2007NeurIPS1,828 citations

Proposes supervised latent Dirichlet allocation (sLDA) to jointly model text content and document labels, enabling the discovery of topic representations directly optimized for predictive accuracy over traditional two-stage regression approaches.

Listen

Organizations increasingly need to extract actionable intelligence and predictions from large collections of electronic text, such as public policy records, consumer reviews, and internal communications. Traditional topic modeling techniques uncover the broad themes in a corpus without supervision, but these unsupervised themes often capture irrelevant general patterns rather than the specific features needed to forecast an external outcome, such as a sentiment rating, financial metric, or legislative vote.

The article introduces and evaluates supervised latent Dirichlet allocation, a statistical modeling framework that jointly discovers underlying text themes while simultaneously predicting an external response variable for each document.

To evaluate this framework, the authors tested the model across two distinct real-world datasets: a collection of 5,006 movie reviews used to predict numeric ratings, and legislative records from the 109th and 110th U.S. Senates (288 and 213 amendments, respectively) used to predict ideological voting tone based on amendment text. Using five-fold cross-validation, the model was benchmarked against two standard industry baselines: standard regularized linear regression (the lasso) and a two-stage pipeline where unsupervised topics are extracted first and then fed into a separate regression model.

The empirical findings show that jointly modeling text and outcomes delivers superior predictive performance. First, the supervised topic approach consistently outperformed two-stage unsupervised topic regression across both datasets and across varying numbers of topics. Second, on the legislative amendment data, the supervised topic model substantially outperformed regularized regression, achieving an 80% relative improvement in predictive explanatory power for the 109th Senate (a predictive R-squared of 0.27 versus 0.15) and a 43% relative improvement for the 110th Senate (0.23 versus 0.16). Third, on the movie review data, the model slightly outperformed the regularized baseline with an explanatory power of 0.432 compared to 0.426. Finally, unlike standard regression methods that function as black-box predictors, the proposed model successfully identified intuitive, interpretable thematic word clusters directly linked to positive or negative outcome values.

These results demonstrate that when the primary objective is outcome prediction, statistical dimensionality reduction must be guided directly by the target metric. Relying on disconnected two-stage text analysis pipelines degrades predictive accuracy and increases analytical risk. Furthermore, because this framework accommodates various outcome data typessuch as continuous values, counts, and binary classificationsit provides decision-makers with a flexible tool that combines quantitative forecasting with qualitative, interpretable insights.

For organizations managing complex text repositories, the article supports adopting supervised topic models over disjointed unsupervised workflows when labeled outcome data exists. Practitioners should explore extensions of the framework, including semi-supervised implementations where only a subset of documents have associated outcome labels, or integrating additional metadata covariates alongside document text. Decision-makers should note that model performance is constrained by text quality and vocabulary preprocessing choices, and legislative outcomes remain inherently difficult to predict solely from text because external political factors also drive voting behavior. Overall, confidence in the model's predictive advantage over conventional two-stage topic regression remains high.

  • Paper: Stochastic variational inference, Matt Hoffman et al. (2012). It scales variational inference for latent Dirichlet allocation and related topic models to massive datasets via stochastic optimization.
  • Paper: Learning Word Vectors for Sentiment Analysis, Andrew L. Maas et al. (2011). It extends joint probabilistic topic and rating modeling to continuous word vectors for sentiment and polarity prediction.
  • Paper: Semi-supervised Learning with Deep Generative Models, Diederik P. Kingma et al. (2014). It advances beyond linear supervised topic models by unifying generative latent variable modeling and classification via deep neural variational inference.
Cover for Supervised Topic Models

Abstract

We introduce supervised latent Dirichlet allocation (sLDA), a statistical model of labelled documents. The model accommodates a variety of response types. We derive an approximate maximum-likelihood procedure for parameter estimation, which relies on variational methods to handle intractable posterior expectations. Prediction problems motivate this research: we use the fitted model to predict response values for new documents. We test sLDA on two real-world problems: movie ratings predicted from reviews, and the political tone of amendments in the U.S. Senate based on the amendment text. We illustrate the benefits of sLDA versus modern regularized regression, as well as versus an unsupervised LDA analysis followed by a separate regression.

Table of Contents

  • 1 Introduction
  • 2 Supervised latent Dirichlet allocation
  • 3 Computation with supervised LDA
  • 3.1 Posterior inference
  • 3.2 Parameter estimation
  • 3.3 Prediction
  • 3.4 Examples
  • 4 Empirical study
  • 5 Discussion
  • References

Knowls

  1. Knowl 1 — Supervised Latent Dirichlet Allocation Generative Model

    model/method

    Supervised latent Dirichlet allocation (sLDA) is a generative model that jointly models a document consisting of NN words w1:N=(w1,,wN)w_{1:N} = (w_1, \dots, w_N) from a vocabulary of size VV and an associated continuous or discrete response variable yRy \in \mathbb{R}.

    The model parameters comprise:

    • KK topic multinomial distributions β1:K\beta_{1:K}, where each βk=(βk,1,,βk,V)ΔV1\beta_k = (\beta_{k,1}, \dots, \beta_{k,V}) \in \Delta^{V-1} represents term probabilities for topic kk.
    • A Dirichlet prior parameter vector αR+K\alpha \in \mathbb{R}_+^K.
    • Generalized linear model (GLM) regression coefficients ηRK\eta \in \mathbb{R}^K.
    • A GLM dispersion parameter δR+\delta \in \mathbb{R}_+.

    The generative process for each document-response pair is as follows:

    1. Draw topic proportions θDirichlet(α)\theta \sim \text{Dirichlet}(\alpha).
    2. For each word index n{1,,N}n \in \{1, \dots, N\}: a. Draw topic assignment indicator vector znMultinomial(θ,1)z_n \sim \text{Multinomial}(\theta, 1), where zn{0,1}Kz_n \in \{0, 1\}^K and k=1Kzn,k=1\sum_{k=1}^K z_{n,k} = 1. b. Draw word wnMultinomial(βzn,1)w_n \sim \text{Multinomial}(\beta_{z_n}, 1).
    3. Define the document's empirical topic frequency vector as: zˉ=1Nn=1Nzn\bar{z} = \frac{1}{N} \sum_{n=1}^N z_n
    4. Draw the response variable yy from an exponential dispersion family conditioned on zˉ\bar{z}: p(yz1:N,η,δ)=h(y,δ)exp((ηzˉ)yA(ηzˉ)δ)p(y \mid z_{1:N}, \eta, \delta) = h(y, \delta) \exp\left( \frac{(\eta^\top \bar{z})y - A(\eta^\top \bar{z})}{\delta} \right) where h(y,δ)h(y, \delta) is the base measure, ηzˉ\eta^\top \bar{z} is the natural parameter acting as the linear predictor, and A(ηzˉ)A(\eta^\top \bar{z}) is the log-normalizer (cumulant generating function).
  2. Knowl 2 — Variational Inference for Supervised Latent Dirichlet Allocation

    model/method

    The document-level posterior distribution of the latent topic proportions θ\theta and topic assignment vectors z1:N=(z1,,zN)z_{1:N} = (z_1, \dots, z_N) given words w1:Nw_{1:N} and response yy is approximated by a fully factorized mean-field variational family: q(θ,z1:Nγ,ϕ1:N)=q(θγ)n=1Nq(znϕn)q(\theta, z_{1:N} \mid \gamma, \phi_{1:N}) = q(\theta \mid \gamma) \prod_{n=1}^N q(z_n \mid \phi_n) where γR+K\gamma \in \mathbb{R}_+^K is a variational Dirichlet parameter vector, and each ϕnΔK1\phi_n \in \Delta^{K-1} parameterizes a categorical distribution such that Eq[zn]=ϕn\mathbb{E}_q[z_n] = \phi_n.

    The variational distribution is optimized to maximize the Evidence Lower Bound (ELBO): L(w1:N,yα,β1:K,η,δ)=Eq[logp(θα)]+n=1NEq[logp(znθ)]+n=1NEq[logp(wnzn,β1:K)]+Eq[logp(yz1:N,η,δ)]+H(q)\mathcal{L}(w_{1:N}, y \mid \alpha, \beta_{1:K}, \eta, \delta) = \mathbb{E}_q[\log p(\theta \mid \alpha)] + \sum_{n=1}^N \mathbb{E}_q[\log p(z_n \mid \theta)] + \sum_{n=1}^N \mathbb{E}_q[\log p(w_n \mid z_n, \beta_{1:K})] + \mathbb{E}_q[\log p(y \mid z_{1:N}, \eta, \delta)] + H(q) where H(q)=Eq[logq(θ,z1:N)]H(q) = -\mathbb{E}_q[\log q(\theta, z_{1:N})] is the variational entropy.

    Block coordinate-ascent updates optimize this bound:

    • The Dirichlet update is identical to unsupervised LDA: γnewα+n=1Nϕn\gamma^{\text{new}} \leftarrow \alpha + \sum_{n=1}^N \phi_n
    • The gradient of the ELBO with respect to the variational multinomial ϕn\phi_n accounts for both word likelihood and response likelihood: Lϕn=Eq[logθ]+logβ,wnlogϕn+1+(yNδ)η1δϕnEq[A(ηZˉ)]\frac{\partial \mathcal{L}}{\partial \phi_n} = \mathbb{E}_q[\log \theta] + \log \beta_{\cdot, w_n} - \log \phi_n + \mathbf{1} + \left( \frac{y}{N\delta} \right) \eta - \frac{1}{\delta} \frac{\partial}{\partial \phi_n} \mathbb{E}_q[A(\eta^\top \bar{Z})] where Zˉ=1Nn=1Nzn\bar{Z} = \frac{1}{N} \sum_{n=1}^N z_n, β,wn=(β1,wn,,βK,wn)\beta_{\cdot, w_n} = (\beta_{1, w_n}, \dots, \beta_{K, w_n})^\top, 1\mathbf{1} is the vector of ones, and Eq[logθi]=Ψ(γi)Ψ(j=1Kγj)\mathbb{E}_q[\log \theta_i] = \Psi(\gamma_i) - \Psi\left( \sum_{j=1}^K \gamma_j \right) with Ψ()\Psi(\cdot) denoting the digamma function.
  3. Knowl 3 — Variational Expectation-Maximization for sLDA Parameter Estimation

    algorithm

    Variational Expectation-Maximization (variational EM) estimates the corpus-wide parameters α\alpha, β1:K\beta_{1:K}, η\eta, and δ\delta across a training corpus D={(wd,1:Nd,yd)}d=1D\mathcal{D} = \{(w_{d,1:N_d}, y_d)\}_{d=1}^D by optimizing the corpus-level lower bound L=d=1DLd\mathcal{L} = \sum_{d=1}^D \mathcal{L}_d.

    Input: Corpus of document-response pairs (wd,1:Nd,yd)d=1D(w_{d,1:N_d}, y_d)_{d=1}^D, number of topics KK, convergence threshold ϵ\epsilon
    Output: Fitted parameters β1:K\beta_{1:K}, η\eta, δ\delta, α\alpha
    Initialize β1:K\beta_{1:K} randomly, η\eta uniformly on [1,1][-1, 1], δ\delta to empirical variance of yy, and α=1/K\alpha = 1/K
    repeat
        // E-step: Document-level variational inference
        for d=1d = 1 to DD do
            repeat
                Update γdα+n=1Ndϕd,n\gamma_d \leftarrow \alpha + \sum_{n=1}^{N_d} \phi_{d,n}
                for n=1n = 1 to NdN_d do
                    Update ϕd,n\phi_{d,n} by coordinate ascent or gradient step on Ldϕd,n\frac{\partial \mathcal{L}_d}{\partial \phi_{d,n}}
                    Normalize ϕd,n\phi_{d,n} to sum to 1
                end for
            until per-document ELBO Ld\mathcal{L}_d converges to within ϵ\epsilon
        end for
        // M-step: Corpus-level parameter updates
        for k=1k = 1 to KK and word w=1w = 1 to VV do
            β^k,wd=1Dn=1Nd1(wd,n=w)ϕd,n,k\hat{\beta}_{k,w} \propto \sum_{d=1}^D \sum_{n=1}^{N_d} 1(w_{d,n} = w) \phi_{d,n,k}
        end for
        Update η\eta and δ\delta using the corpus ELBO gradients:
        Lη=1δ(d=1Dϕˉdydd=1DEqd[μ(ηZˉd)Zˉd])\frac{\partial \mathcal{L}}{\partial \eta} = \frac{1}{\delta}\left( \sum_{d=1}^D \bar{\phi}_d y_d - \sum_{d=1}^D \mathbb{E}_{q_d}[\mu(\eta^\top \bar{Z}_d)\bar{Z}_d] \right) where ϕˉd=1Ndn=1Ndϕd,n\bar{\phi}_d = \frac{1}{N_d}\sum_{n=1}^{N_d} \phi_{d,n}
        Update α\alpha using Dirichlet maximum-likelihood on expected topic proportions
    until corpus-level ELBO converges to within ϵ\epsilon
    return β1:K,η,δ,α\beta_{1:K}, \eta, \delta, \alpha
  4. Knowl 4 — Response Prediction for Unlabeled Documents in sLDA

    algorithm

    Given a fitted sLDA model with parameters {α,β1:K,η,δ}\{\alpha, \beta_{1:K}, \eta, \delta\} and a newly observed unlabeled document w1:N=(w1,,wN)w_{1:N} = (w_1, \dots, w_N), predicting the response yy corresponds to approximating the posterior expectation E[Yw1:N,α,β1:K,η,δ]=E[μ(ηZˉ)w1:N,α,β1:K]\mathbb{E}[Y \mid w_{1:N}, \alpha, \beta_{1:K}, \eta, \delta] = \mathbb{E}[\mu(\eta^\top \bar{Z}) \mid w_{1:N}, \alpha, \beta_{1:K}], where μ()=EGLM[Y]\mu(\cdot) = \mathbb{E}_{\text{GLM}}[Y \mid \cdot].

    Input: New document words w1:Nw_{1:N}, fitted model {α,β1:K,η,δ}\{\alpha, \beta_{1:K}, \eta, \delta\}, convergence threshold ϵ\epsilon
    Output: Predicted scalar response y^\hat{y}
    Initialize γα+N/K\gamma \leftarrow \alpha + N/K and ϕn1/K\phi_n \leftarrow 1/K for all n=1,,Nn = 1, \dots, N
    repeat
        γα+n=1Nϕn\gamma \leftarrow \alpha + \sum_{n=1}^N \phi_n
        for n=1n = 1 to NN do
            for k=1k = 1 to KK do
                ϕn,kβk,wnexp(Ψ(γk)Ψ(j=1Kγj))\phi_{n,k} \propto \beta_{k, w_n} \exp\left( \Psi(\gamma_k) - \Psi\left(\sum_{j=1}^K \gamma_j\right) \right)
            end for
            Normalize ϕn\phi_n such that k=1Kϕn,k=1\sum_{k=1}^K \phi_{n,k} = 1
        end for
    until variational parameters γ,ϕ1:N\gamma, \phi_{1:N} converge to within ϵ\epsilon
    Compute ϕˉ=1Nn=1Nϕn\bar{\phi} = \frac{1}{N} \sum_{n=1}^N \phi_n
    Compute prediction y^Eq[μ(ηZˉ)]\hat{y} \approx \mathbb{E}_q[\mu(\eta^\top \bar{Z})]
    // For Gaussian response: y^=ηϕˉ\hat{y} = \eta^\top \bar{\phi}
    // For Poisson response: y^=n=1N(K1+i=1Kϕn,iexp(ηi/N))\hat{y} = \prod_{n=1}^N \left( K - 1 + \sum_{i=1}^K \phi_{n,i} \exp(\eta_i / N) \right)
    return y^\hat{y}
  5. Knowl 5 — Exact Variational Updates and Estimation for Gaussian Response sLDA

    model/method

    When the response distribution is Gaussian with mean μ=ηzˉ\mu = \eta^\top \bar{z} and variance σ2=δ\sigma^2 = \delta, the base measure is h(y,δ)=12πσ2exp(y2/(2σ2))h(y, \delta) = \frac{1}{\sqrt{2\pi\sigma^2}}\exp(-y^2 / (2\sigma^2)) and the log-normalizer is A(ηzˉ)=12(ηzˉ)2A(\eta^\top \bar{z}) = \frac{1}{2}(\eta^\top \bar{z})^2.

    1. Second Moment of Topic Frequencies: Eq[ZˉZˉ]=1N2(n=1Nmnϕnϕm+n=1Ndiag(ϕn))\mathbb{E}_q[\bar{Z}\bar{Z}^\top] = \frac{1}{N^2}\left( \sum_{n=1}^N \sum_{m \neq n} \phi_n \phi_m^\top + \sum_{n=1}^N \text{diag}(\phi_n) \right)

    2. Coordinate Ascent Update for ϕj\phi_j: ϕjnewexp(Eq[logθ]+logβ,wj+(yNσ2)η12N2σ2[2(ηϕj)η+(ηη)])\phi_j^{\text{new}} \propto \exp\left( \mathbb{E}_q[\log \theta] + \log \beta_{\cdot, w_j} + \left(\frac{y}{N\sigma^2}\right)\eta - \frac{1}{2N^2\sigma^2}\left[ 2(\eta^\top \phi_{-j})\eta + (\eta \circ \eta) \right] \right) where ϕj=njϕn\phi_{-j} = \sum_{n \neq j} \phi_n, ηη\eta \circ \eta is the elementwise square, and β,wj=(β1,wj,,βK,wj)\beta_{\cdot, w_j} = (\beta_{1,w_j}, \dots, \beta_{K,w_j})^\top.

    3. Closed-Form M-Step for GLM Parameters: Let y=(y1,,yD)RDy = (y_1, \dots, y_D)^\top \in \mathbb{R}^D and XRD×KX \in \mathbb{R}^{D \times K} be the design matrix whose dd-th row is Zˉd\bar{Z}_d. The updates are: η^new=(Eq[XX])1Eq[X]y=(d=1DEqd[ZˉdZˉd])1(d=1Dϕˉdyd)\hat{\eta}^{\text{new}} = \left( \mathbb{E}_q[X^\top X] \right)^{-1} \mathbb{E}_q[X]^\top y = \left( \sum_{d=1}^D \mathbb{E}_{q_d}[\bar{Z}_d \bar{Z}_d^\top] \right)^{-1} \left( \sum_{d=1}^D \bar{\phi}_d y_d \right) σ^2,new=1D(yyyEq[X](Eq[XX])1Eq[X]y)\hat{\sigma}^{2,\text{new}} = \frac{1}{D} \left( y^\top y - y^\top \mathbb{E}_q[X] \left( \mathbb{E}_q[X^\top X] \right)^{-1} \mathbb{E}_q[X]^\top y \right) This update differs from ordinary least squares of yy on Eq[X]\mathbb{E}_q[X] because Eq[XX]Eq[X]Eq[X]\mathbb{E}_q[X^\top X] \neq \mathbb{E}_q[X]^\top \mathbb{E}_q[X].

  6. Knowl 6 — Variational Inference and Estimation for Overdispersed Poisson Response sLDA

    model/method

    When the response variable yN0y \in \mathbb{N}_0 is count data modeled by an overdispersed Poisson distribution with mean parameter λ=exp(ηzˉ)\lambda = \exp(\eta^\top \bar{z}) and dispersion δ\delta, the density is: p(yλ,δ)=1y!exp(ylogλλδ)p(y \mid \lambda, \delta) = \frac{1}{y!} \exp\left( \frac{y\log\lambda - \lambda}{\delta} \right) Here, A(ηzˉ)=μ(ηzˉ)=exp(ηzˉ)A(\eta^\top \bar{z}) = \mu(\eta^\top \bar{z}) = \exp(\eta^\top \bar{z}).

    1. Expected Log-Normalizer and Derivative: Eq[A(ηZˉ)]=n=1NEq[exp(1NηZn)]=n=1N(K1+i=1Kϕn,iexp(ηi/N))=:C\mathbb{E}_q[A(\eta^\top \bar{Z})] = \prod_{n=1}^N \mathbb{E}_q\left[ \exp\left( \frac{1}{N}\eta^\top Z_n \right) \right] = \prod_{n=1}^N \left( K - 1 + \sum_{i=1}^K \phi_{n,i} \exp(\eta_i / N) \right) =: C ϕnEq[A(ηZˉ)]=Cnexp(η/N)\frac{\partial}{\partial \phi_n} \mathbb{E}_q[A(\eta^\top \bar{Z})] = C_{-n} \exp(\eta / N) where Cn=mn(K1+i=1Kϕm,iexp(ηi/N))C_{-n} = \prod_{m \neq n} \left( K - 1 + \sum_{i=1}^K \phi_{m,i} \exp(\eta_i / N) \right).

    2. Expected Mean-Feature Product: Eq[μ(ηZˉ)Zˉ]=exp(η/N)Nn=1NCnϕn\mathbb{E}_q[\mu(\eta^\top \bar{Z})\bar{Z}] = \frac{\exp(\eta/N)}{N} \sum_{n=1}^N C_{-n}\phi_n

    3. M-Step Parameter Optimization:

    • The gradient of the corpus ELBO w.r.t. η\eta is optimized numerically via convex optimization: Lη=1δ(d=1Dϕˉdydd=1Dexp(η/Nd)Ndn=1NdCd,nϕd,n)\frac{\partial \mathcal{L}}{\partial \eta} = \frac{1}{\delta}\left( \sum_{d=1}^D \bar{\phi}_d y_d - \sum_{d=1}^D \frac{\exp(\eta/N_d)}{N_d} \sum_{n=1}^{N_d} C_{d,-n}\phi_{d,n} \right)
    • The dispersion parameter update has an exact closed form: δ^new=d=1Dη^newϕˉdydd=1DCd(η^new)\hat{\delta}^{\text{new}} = \frac{\sum_{d=1}^D \hat{\eta}_{\text{new}}^\top \bar{\phi}_d y_d}{\sum_{d=1}^D C_d(\hat{\eta}_{\text{new}})}
  7. Knowl 7 — Empirical Predictive Performance of sLDA on Movie Ratings and Senate Ideology

    empirical result

    Supervised LDA was evaluated on two real-world continuous-response prediction tasks using 5-fold cross-validation. Performance was measured by Pearson correlation and predictive R2R^2 (pR2=1(yy^)2(yyˉ)2pR^2 = 1 - \frac{\sum (y - \hat{y})^2}{\sum (y - \bar{y})^2}):

    1. Movie Review Rating Prediction:
    • Corpus: 5,006 movie reviews paired with star ratings (log-transformed), with a 2,180-word vocabulary.
    • Results: sLDA achieved a peak pR2pR^2 of 0.4320.432, improving upon regularized linear regression via lasso (pR2=0.426pR^2 = 0.426) and substantially outperforming unsupervised LDA followed by linear regression (which plateaued at pR20.24pR^2 \approx 0.24).
    1. U.S. Senate Amendment Political Tone Prediction:
    • Task: Predict amendment ideal-point discrimination parameters βi\beta_i (derived from roll call voting records) based on amendment text.
    • 109th Senate (288 amendments, 2,084 vocabulary words): sLDA achieved a best pR2pR^2 of 0.270.27, an 80% relative improvement over lasso (pR2=0.15pR^2 = 0.15) and outperforming unsupervised LDA + regression (pR20.22pR^2 \approx 0.22).
    • 110th Senate (213 amendments, 1,653 vocabulary words): sLDA achieved a best pR2pR^2 of 0.230.23, a 43% relative improvement over lasso (pR2=0.16pR^2 = 0.16) and outperforming unsupervised LDA + regression (pR20.18pR^2 \approx 0.18).

    Across all settings, jointly estimating topics conditioned on the response variable produced topic representations with higher predictive accuracy than two-stage unsupervised LDA feature extraction.

Coverage note — None was omitted; the extracted knowls fully cover the sLDA model formulation, variational inference, variational EM parameter estimation, out-of-fold prediction, exact Gaussian and Poisson solutions, and empirical cross-validation results.

References

  1. 1.Bickel, P. and Doksum, K. (2007). Mathematical Statistics: Basic Ideas and Selected Topics, volume 1. Pearson Prentice Hall, Upper Saddle River, NJ, 2nd edition.
  2. 2.Bishop, C. (2006). Pattern Recognition and Machine Learning. Springer New York.
  3. 3.Blei, D. and Jordan, M. (2003). Modeling annotated data. In Proceedings of the 26th annual International ACM SIGIR Conference on Research and Development in Information Retrieval, pages 127–134. ACM Press.
  4. 4.Blei, D. and Lafferty, J. (2009). Topic models. In Srivastava, A. and Sahami, M., editors, Text Mining: Theory and Applications. Taylor and Francis.
  5. 5.Blei, D., Ng, A., and Jordan, M. (2003). Latent Dirichlet allocation. Journal of Machine Learning Research, 3:993–1022.
  6. 6.Braun, M. and McAuliffe, J. (2010). Variational inference for large-scale models of discrete choice. Journal of the American Statistical Association.
  7. 7.Brown, L. (1986). Fundamentals of Statistical Exponential Families. Institute of Mathematical Statistics, Hayward, CA.
  8. 8.Chang, J. and Blei, D. (2010). Hierarchical relational models for document networks. Annals of Applied Statistics.
  9. 9.Chang, J., Boyd-Graber, J., Wang, C., Gerrish, S., and Blei, D. (2009). Reading tea leaves: How humans interpret topic models. In Neural Information Processing Systems.
  10. 10.Clinton, J., Jackman, S., and Rivers, D. (2004). The statistical analysis of roll call data. American Political Science Review, 98(2):355–370.
  11. 11.Dempster, A., Laird, N., and Rubin, D. (1977). Maximum likelihood from incomplete data via the EM algorithm. Journal of the Royal Statistical Society, Series B, 39:1–38.
  12. 12.Erosheva, E. (2002). Grade of membership and latent structure models with application to disability survey data. PhD thesis, Carnegie Mellon University, Department of Statistics.
  13. 13.Erosheva, E., Fienberg, S., and Lafferty, J. (2004). Mixed-membership models of scientific publications. Proceedings of the National Academy of Science, 97(22):11885–11892.
  14. 14.Fei-Fei, L. and Perona, P. (2005). A Bayesian hierarchical model for learning natural scene categories. IEEE Computer Vision and Pattern Recognition, pages 524–531.
  15. 15.Flaherty, P., Giaever, G., Kumm, J., Jordan, M., and Arkin, A. (2005). A latent variable model for chemogenomic profiling. Bioinformatics, 21(15):3286–3293.
  16. 16.Fukumizu, K., Bach, F., and Jordan, M. (2004). Dimensionality reduction for supervised learning with reproducing kernel hilbert spaces. Journal of Machine Learning Research, 5:73–99.
  17. 17.Griffiths, T., Steyvers, M., Blei, D., and Tenenbaum, J. (2005). Integrating topics and syntax. In Saul, L. K., Weiss, Y., and Bottou, L., editors, Advances in Neural Information Processing Systems 17, pages 537–544, Cambridge, MA. MIT Press.
  18. 18.Hastie, T., Tibshirani, R., and Friedman, J. (2001). The Elements of Statistical Learning. Springer.
  19. 19.Jordan, M., Ghahramani, Z., Jaakkola, T., and Saul, L. (1999). Introduction to variational methods for graphical models. Machine Learning, 37:183–233.
  20. 20.Marlin, B. (2003). Modeling user rating profiles for collaborative filtering. In Neural Information Processing Systems.
  21. 21.McCallum, A., Pal, C., Druck, G., and Wang, X. (2006). Multi-conditional learning: Generative/discriminative training for clustering and classification. In AAAI.
  22. 22.McCullagh, P. and Nelder, J. A. (1989). Generalized Linear Models. London: Chapman and Hall.
  23. 23.Mimno, D. and McCallum, A. (2007). Organizing the OCA: Learning faceted subjects from a library of digital books. In Joint Conference on Digital Libraries.
  24. 24.Pang, B. and Lee, L. (2005). Seeing stars: Exploiting class relationships for sentiment categorization with respect to rating scales. In Proceedings of the Association of Computational Linguistics.
  25. 25.Pritchard, J., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics, 155:945–959.
  26. 26.Quelhas, P., Monay, F., Odobez, J., Gatica-Perez, D., Tuyelaars, T., and Van Gool, L. (2005). Modeling scenes with local descriptors and latent aspects. In ICCV.
  27. 27.Ramage, D., Hall, D., Nallapati, R., and Manning, C. (2009). Labeled LDA: A supervised topic model for credit attribution in muli-labeled corpora. In Empirical Methods in Natural Language Processing.
  28. 28.Ronning, G. (1989). Maximum likelihood estimation of Dirichlet distributions. Journal of Statistcal Computation and Simulation, 34(4):215–221.
  29. 29.Rosen-Zvi, M., Griffiths, T., Steyvers, M., and Smith, P. (2004). The author-topic model for authors and documents. In Proceedings of the 20th Conference on Uncertainty in Artificial Intelligence, pages 487–494. AUAI Press.
  30. 30.Steyvers, M. and Griffiths, T. (2006). Probabilistic topic models. In Landauer, T., McNamara, D., Dennis, S., and Kintsch, W., editors, Latent Semantic Analysis: A Road to Meaning. Laurence Erlbaum.
  31. 31.Tibshirani, R. (1996). Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society, Series B (Methodological), 58(1):267–288.
  32. 32.Wainwright, M. and Jordan, M. (2008). Graphical models, exponential families, and variational inference. Foundations and Trends in Machine Learning, 1(1–2):1–305.
  33. 33.Wallach, H., Mimno, D., and McCallum, A. (2009). Rethinking lda: Why priors matter. In Bengio, Y., Schuurmans, D., Lafferty, J., Williams, C. K. I., and Culotta, A., editors, Advances in Neural Information Processing Systems 22, pages 1973–1981.
  34. 34.Wang, C., Blei, D., and Li, F. (2009). Simultaneous image classification and annotation. In Computer Vision and Pattern Recognition.
  35. 35.Wei, X. and Croft, B. (2006). LDA-based document models for ad-hoc retrieval. In SIGIR.

Citation

MLA
Blei, D. M., and J. D. McAuliffe. “Supervised Topic Models”. arXiv, 2010, http://arxiv.org/abs/1003.0783v1.
APA
Blei, D. M., & McAuliffe, J. D. (2010). Supervised Topic Models. arXiv. http://arxiv.org/abs/1003.0783v1
Chicago
Blei, D. M., and J. D. McAuliffe. 2010. “Supervised Topic Models”. arXiv. http://arxiv.org/abs/1003.0783v1.
Harvard
Blei, D.M. and McAuliffe, J.D. (2010) “Supervised Topic Models”, arXiv [Preprint]. Available at: http://arxiv.org/abs/1003.0783v1.
Vancouver
1. Blei DM, McAuliffe JD (2010) Supervised Topic Models. arXiv

BibTeX

@article{blei2010supervised,
  title = {Supervised Topic Models},
  author = {Blei, David M. and McAuliffe, Jon D.},
  year = {2010},
  journal = {arXiv},
  url = {http://arxiv.org/abs/1003.0783v1},
  eprint = {1003.0783}
}
Metadata:arXiv

Source Code

This paper has an official code repository available. Click below to access the source code.

View Repository

Access the Paper

This paper is available from its original source. Click below to access the PDF.

Open PDF

License: Authors