Latent semantic analysis
'''Latent semantic analysis (LSA)''' is a technique in [[natural language processing]], in particular in [[vectorial semantics]], of analyzing relationships between a set of documents and the terms they contain by producing a set of concepts related to the documents and terms.  

LSA was patented in [[1988]] ([http://patft.uspto.gov/netacgi/nph-Parser?patentnumber=4839853 US Patent 4,839,853]) by [[Scott Deerwester]], [[Susan Dumais]], [[George Furnas]], [[Richard Harshman]], [[Thomas Landauer]], [[Karen Lochbaum]] and [[Lynn Streeter]]. In the context of its application to [[information retrieval]], it is sometimes called '''latent semantic indexing (LSI)'''.

== Occurrence matrix ==
LSA can use a [[term-document matrix]] which describes the occurrences of terms in documents; it is a [[sparse matrix]] whose rows correspond to [[terminology|terms]] and whose columns correspond to documents, typically [[stemming|stemmed]] words that appear in the documents. A typical example of the weighting of the elements of the matrix is [[tf-idf]] (term frequency–inverse document frequency): the element of the matrix is proportional to the number of times the terms appear in each document, where rare terms are upweighted to reflect their relative importance.

This matrix is also common to standard semantic models, though it is not necessarily explicitly expressed as a matrix, since the mathematical properties of matrices are not always used.

LSA transforms the occurrence matrix into a relation between the terms and some ''concepts'', and a relation between those concepts and the documents. Thus the terms and documents are now indirectly related through the concepts.

== Applications ==

The new concept space typically can be used to:
* Compare the documents in the concept space ([[data clustering]], [[document classification]])......
* Find similar documents across languages, after analyzing a base set of translated documents ([[cross language retrieval]]).
* Find relations between terms ([[synonymy]] and [[polysemy]]).
* Given a query of terms, translate it into the concept space, and find matching documents ([[information retrieval]]).

Synonymy and polysemy are fundamental problems in [[natural language processing]]: 
* Synonymy is the phenomenon where different words describe the same idea. Thus, a query in a search engine may fail to retrieve a relevant document that does not contain the words which appeared in the query. For example, a search for "doctors" may not return a document containing the word "physicians", even though the words have the same meaning.
* Polysemy is the phenomenon where the same word has multiple meanings. So a search may retrieve irrelevant documents containing the desired words in the wrong meaning. For example, a botanist and a computer scientist looking for the word "tree" probably desire different sets of documents.

== Rank lowering ==
After the construction of the occurrence matrix, LSA finds a low-[[rank (matrix theory)|rank]] approximation to the [[term-document matrix]]. There could be various reasons for these approximations:

* The original term-document matrix is presumed too large for the computing resources; in this case, the approximated low rank  matrix is interpreted as an ''approximation'' (a "least and necessary evil").
* The original term-document matrix is presumed ''noisy'': for example, anecdotal instances of terms are to be eliminated. From this point of view, the approximated matrix is interpreted as a ''de-noisified matrix'' (a better matrix than the original).
* The original term-document matrix is presumed overly [[Sparse matrix|sparse]] relative to the "true" term-document matrix.  That is, the original matrix lists only the words actually ''in'' each document, whereas we might be interested in all words ''related to'' each document--generally a much larger set due to [[synonymy]].

The consequence of the rank lowering is that some dimensions are combined and depend on more than one term:

:: {(car), (truck), (flower)} -->  {(1.3452 * car + 0.2828 * truck), (flower)}

This mitigates synonymy, as the rank lowering is expected to merge the dimensions associated with terms that have similar meanings. It also mitigates polysemy, since components of polysemous words that point in the "right" direction are added to the components of words that share a similar meaning. Conversely, components that point in other directions tend to either simply cancel out, or, at worst, to be smaller than components in the directions corresponding to the intended sense.

== Derivation ==
Let <math>X</math> be a matrix where element <math>(i,j)</math> describes the occurrence of term <math>i</math> in document <math>j</math> (this can be, for example, the frequency). <math>X</math> will look like this:

:<math>
\begin{matrix} 
 & \textbf{d}_j \\
 & \downarrow \\
\textbf{t}_i^T \rightarrow &
\begin{bmatrix} 
x_{1,1} & \dots & x_{1,n} \\
\vdots & \ddots & \vdots \\
x_{m,1} & \dots & x_{m,n} \\
\end{bmatrix}
\end{matrix}
</math>

Now a row in this matrix will be a vector corresponding to a term, giving its relation to each document: 

:<math>\textbf{t}_i^T = \begin{bmatrix} x_{i,1} & \dots & x_{i,n} \end{bmatrix}</math>

Likewise, a column in this matrix will be a vector corresponding to a document, giving its relation to each term:

:<math>\textbf{d}_j = \begin{bmatrix} x_{1,j} \\ \vdots \\ x_{m,j} \end{bmatrix}</math>

Now the [[dot product]] <math>\textbf{t}_i^T \textbf{t}_p</math> between two term vectors gives the [[correlation]] between the terms over the documents. The [[matrix product]] <math>X X^T</math> contains all these dot products. Element <math>(i,p)</math> (which is equal to element <math>(p,i)</math>) contains the dot product <math>\textbf{t}_i^T \textbf{t}_p</math> (<math> = \textbf{t}_p^T \textbf{t}_i</math>). Likewise, the matrix <math>X^T X</math> contains the dot products between all the document vectors, giving their correlation over the terms: <math>\textbf{d}_j^T \textbf{d}_q = \textbf{d}_q^T \textbf{d}_j</math>.

Now assume that there exists a decomposition of <math>X</math> such that <math>U</math> and <math>V</math> are [[orthonormal matrix|orthonormal matrices]] and <math>\Sigma</math> is a [[diagonal matrix]]. This is called a [[singular value decomposition]] (SVD):

:<math>
X = U \Sigma V^T
</math>

The matrix products giving us the term and document correlations then become

:<math>
\begin{matrix}
X X^T &=& (U \Sigma V^T) (U \Sigma V^T)^T = (U \Sigma V^T) (V^{T^T} \Sigma^T U^T) = U \Sigma V^T V \Sigma^T U^T = U \Sigma \Sigma^T U^T \\
X^T X &=& (U \Sigma V^T)^T (U \Sigma V^T) = (V^{T^T} \Sigma^T U^T) (U \Sigma V^T) = V \Sigma U^T U \Sigma V^T = V \Sigma^T \Sigma V^T
\end{matrix}
</math>

Since <math>\Sigma \Sigma^T</math> and <math>\Sigma^T \Sigma</math> are diagonal we see that <math>U</math> must contain the [[eigenvector]]s of <math>X X^T</math>, while <math>V</math> must be the eigenvectors of <math>X^T X</math>. Both products have the same non-zero eigenvalues, given by the non-zero entries of <math>\Sigma \Sigma^T</math>, or equally, by the non-zero entries of <math>\Sigma^T\Sigma</math>. Now the decomposition looks like this:

:<math>
\begin{matrix} 
 & X & & & U & & \Sigma & & V^T \\
 & (\textbf{d}_j) & & & & & & & (\hat \textbf{d}_j) \\
 & \downarrow & & & & & & & \downarrow \\
(\textbf{t}_i^T) \rightarrow 
&
\begin{bmatrix} 
x_{1,1} & \dots & x_{1,n} \\
\\
\vdots & \ddots & \vdots \\
\\
x_{m,1} & \dots & x_{m,n} \\
\end{bmatrix}
&
=
&
(\hat \textbf{t}_i^T) \rightarrow
&
\begin{bmatrix} 
\begin{bmatrix} \, \\ \, \\ \textbf{u}_1 \\ \, \\ \,\end{bmatrix} 
\dots
\begin{bmatrix} \, \\ \, \\ \textbf{u}_l \\ \, \\ \, \end{bmatrix}
\end{bmatrix}
&
\cdot
&
\begin{bmatrix} 
\sigma_1 & \dots & 0 \\
\vdots & \ddots & \vdots \\
0 & \dots & \sigma_l \\
\end{bmatrix}
&
\cdot
&
\begin{bmatrix} 
\begin{bmatrix} & & \textbf{v}_1 & & \end{bmatrix} \\
\vdots \\
\begin{bmatrix} & & \textbf{v}_l & & \end{bmatrix}
\end{bmatrix}
\end{matrix}
</math>

The values <math>\sigma_1, \dots, \sigma_l</math> are called the singular values, and <math>u_1, \dots, u_l</math> and <math>v_1, \dots, v_l</math> the left and right singular vectors.
Notice how the only part of <math>U</math> that contributes to <math>\textbf{t}_i</math> is the <math>i\textrm{'th}</math> row. 
Let this row vector be called <math>\hat \textrm{t}_i</math>. 
Likewise, the only part of <math>V^T</math> that contributes to <math>\textbf{d}_j</math> is the <math>j\textrm{'th}</math> column, <math>\hat \textrm{d}_j</math>. 
These are ''not'' the eigenvectors, but ''depend'' on ''all'' the eigenvectors. 

It turns out that when you select the <math>k</math> largest singular values, and their corresponding singular vectors from <math>U</math> and <math>V</math>, you get the rank <math>k</math> approximation to X with the smallest error ([[Frobenius norm]]). The amazing thing about this approximation is that not only does it have a minimal error, but it translates the term and document vectors into a concept space. The vector <math>\hat \textbf{t}_i</math> then has <math>k</math> entries, each giving the occurrence of term <math>i</math> in one of the <math>k</math> concepts. Likewise, the vector <math>\hat \textbf{d}_j</math> gives the relation between document <math>j</math> and each concept. We write this approximation as 

:<math>X_k = U_k \Sigma_k V_k^T</math>

You can now do the following:
* See how related documents <math>j</math> and <math>q</math> are in the concept space by comparing the vectors <math>\hat \textbf{d}_j</math> and <math>\hat \textbf{d}_q</math> (typically by [[vector space model|cosine similarity]]). This gives you a clustering of the documents.
* Comparing terms <math>i</math> and <math>p</math> by comparing the vectors <math>\hat \textbf{t}_i</math> and <math>\hat \textbf{t}_p</math>, giving you a clustering of the terms in the concept space.
* Given a query, view this as a mini document, and compare it to your documents in the concept space.

To do the latter, you must first translate your query into the concept space. It is then intuitive that you must use the same transformation that you use on your documents:

:<math>\textbf{d}_j = U_k \Sigma_k \hat \textbf{d}_j</math>

:<math>\hat \textbf{d}_j = \Sigma_k^{-1} U_k^T \textbf{d}_j</math>

This means that if you have a query vector <math>q</math>, you must do the translation <math>\hat \textbf{q} = \Sigma_k^{-1} U_k^T \textbf{q}</math> before you compare it with the document vectors in the concept space. You can do the same for pseudo term vectors:

:<math>\textbf{t}_i^T = \hat \textbf{t}_i^T \Sigma_k V_k^T</math>

:<math>\hat \textbf{t}_i^T = \textbf{t}_i^T V_k^{-T} \Sigma_k^{-1} = \textbf{t}_i^T V_k \Sigma_k^{-1}</math>

:<math>\hat \textbf{t}_i = \Sigma_k^{-1}  V_k^T \textbf{t}_i</math>

== Implementation ==

The [[Singular Value Decomposition|SVD]] is typically computed using large matrix methods (for example, [[Lanczos method]]s) but may also be computed incrementally and with greatly reduced resources via a [[neural network]]-like approach, which does not require the large, full-rank matrix to be held in memory ([http://www.dcs.shef.ac.uk/~genevieve/gorrell_webb.pdf Gorrell and Webb, 2005]).

A fast, incremental, low-memory, large-matrix SVD algorithm has recently been developed ([http://www.merl.com/publications/TR2006-059/ Brand, 2006]). Unlike Gorrell and Webb's (2005) stochastic approximation, Brand's (2006) algorithm provides an exact solution.

== Limitations ==
LSA has two drawbacks:

* The resulting dimensions might be difficult to interpret. For instance, in 
:: {(car), (truck), (flower)} -->  {(1.3452 * car + 0.2828 * truck), (flower)}
:the (1.3452 * car + 0.2828 * truck) component could be interpreted as "vehicle". However, it is very likely that cases close to 
:: {(car), (bottle), (flower)} -->  {(1.3452 * car + 0.2828 * bottle), (flower)}
:will occur. This leads to results which can be justified on the mathematical level, but have no interpretable meaning in natural language.

* The [[probabilistic model]] of LSA does not match observed data: LSA assumes that words and documents form a joint [[normal distribution|Gaussian]] model ([[ergodic hypothesis]]), while a [[Poisson distribution]] has been observed.  Thus, a newer alternative is [[probabilistic latent semantic analysis]], based on a [[multinomial distribution|multinomial]] model, which is reported to give better results than standard LSA {{Fact|date=April 2008}}.

== See also ==
* [[Vectorial semantics]]
* [[DSIR model]]
* [[Latent Dirichlet allocation]]
* [[Spamdexing]]
* An example of the application of [http://blog.targetwoman.com/latent-semantic-analysis/] Latent Semantic Analysis in Natural language Processing 
* [[Probabilistic latent semantic analysis]]
* [[Latent semantic mapping]]
* [[Latent Semantic Structure Indexing]]
* [[Principal components analysis]]
* [[Compound term processing]]

==External links==
* [http://www.seobook.com/lsi/lsa_definition.htm Latent Semantic Indexing], a non mathematical introduction and explanation of LSI

* [http://knowledgesearch.org/ The Semantic Indexing Project], an open source program for latent semantic indexing

* [http://senseclusters.sourceforge.net SenseClusters], an open source package for Latent Semantic Analysis and other methods for clustering similar contexts

== References ==
* {{cite web
 | url=http://lsa.colorado.edu/
 | title=The Latent Semantic Indexing home page
}}
* {{cite journal
 | url=http://www.merl.com/publications/TR2006-059/
 | title=Fast Low-Rank Modifications of the Thin Singular Value Decomposition
 | author=Matthew Brand
 | journal=Linear Algebra and Its Applications
 | volume=415
 | pages=20–30
 | date=2006
 | doi=10.1016/j.laa.2005.07.021
}} -- a [http://www.eecs.umich.edu/~wingated/resources.html MATLAB implementation of Brand's algorithm] is available
* {{cite journal
 | url=http://lsa.colorado.edu/papers/dp1.LSAintro.pdf
 | title=Introduction to Latent Semantic Analysis
 | author=[[Thomas Landauer]], P. W. Foltz, & D. Laham
 | journal=Discourse Processes
 | volume=25
 | pages=259–284
 |date=1998
}}
* {{cite journal
 | url=http://lsi.research.telcordia.com/lsi/papers/JASIS90.pdf 
 | title=Indexing by Latent Semantic Analysis
 | author=[[Scott Deerwester|S. Deerwester]], [[Susan Dumais]], [[George Furnas|G. W. Furnas]], [[Thomas Landauer|T. K. Landauer]], [[Richard Harshman|R. Harshman]]
 | journal=Journal of the American Society for Information Science
 | volume=41
 | issue=6
 | pages=391–407
 |date=1990 
 | doi=10.1002/(SICI)1097-4571(199009)41:6<391::AID-ASI1>3.0.CO;2-9
}} Original article where the model was first exposed.
* {{cite journal
 | url=http://citeseer.ist.psu.edu/berry95using.html
 | title=Using Linear Algebra for Intelligent Information Retrieval
 | author=Michael Berry, [[Susan Dumais|S.T. Dumais]], G.W. O'Brien
 |date=1995
}} [http://lsirwww.epfl.ch/courses/dis/2003ws/papers/ut-cs-94-270.pdf PDF]. Illustration of the application of LSA to document retrieval.
* {{cite web
 | url=http://iv.slis.indiana.edu/sw/lsa.html
 | title=Latent Semantic Analysis
 | publisher=InfoVis
}}
* {{cite conference
 | url=http://www.cs.brown.edu/people/th/papers/Hofmann-UAI99.pdf
 | title=Probabilistic Latent Semantic Analysis
 | author=T. Hofmann
 | booktitle=Uncertainty in Artificial Intelligence
 |date=1999
}}
* {{cite conference
 | url=http://www.dcs.shef.ac.uk/~genevieve/gorrell_webb.pdf
 | title=Generalized Hebbian Algorithm for Latent Semantic Analysis
 | author=G. Gorrell and B. Webb
 | booktitle=Interspeech
 |date=2005
}}
* {{cite web
 | url=http://cran.at.r-project.org/web/packages/lsa/index.html
 | title=An Open Source LSA Package for R
 | publisher=CRAN
 | author=Fridolin Wild
 | date=[[November 23]] [[2005]]
 | accessdate=2006-11-20
}}
* {{ cite web
 | url=http://www.welchco.com/02/14/01/60/96/02/2901.HTM
 | title=A Solution to Plato's Problem: The Latent Semantic Analysis Theory of Acquisition, Induction, and Representation of Knowledge
 | author=[[Thomas Landauer]]
 | accessdate=2007-07-02
}}
* {{cite web
 | url=http://scgroup.hpclab.ceid.upatras.gr/scgroup/Projects/TMG/
 | title=A MATLAB Toolbox for generating term-document matrices from text collections
 | author=Dimitrios Zeimpekis and E. Gallopoulos
 | date=September 11, 2005
 | accessdate=2006-11-20
}}

[[Category:Information retrieval]]
[[Category:Natural language processing]]

[[de:Latent Semantic Indexing]]
[[fr:Analyse sémantique latente]]
[[ja:潜在意味解析]]
[[zh:潜在语义学]]
[[ru:Латентно-семантический анализ]]