Latent Semantic Analysis & Sentiment Classification with Python | by Susan Li | Towards Data Science

Latent Semantic Analysis (LSA) is a theory and method for extracting and representing the contextual-usage meaning of words by statistical computations applied to a large corpus of text.
LSA is an information retrieval technique which analyzes and identifies the pattern in unstructured collection of text and the relationship between them.
LSA itself is an unsupervised way of uncovering synonyms in a collection of documents.
To start, we take a look how Latent Semantic Analysis is used in Natural Language Processing to analyze relationships between a set of documents and the terms that they contain. Then we go steps further to analyze and classify sentiment. We will review Chi Squared for feature selection along the way. Let’s get started!
The data set consists of over 500,000 reviews of fine foods from Amazon that can be downloaded from Kaggle.
import pandas as pddf = pd.read_csv('Reviews.csv')df.head()
TF-IDF is an information retrieval technique that weighs a term’s frequency (TF) and its inverse document frequency (IDF). Each word has its respective TF and IDF score. The product of the TF and IDF scores of a word is called the TFIDF weight of that word.
Put simply, the higher the TFIDF score (weight), the rarer the word and vice versa.
from sklearn.feature_extraction.text import TfidfVectorizertfidf = TfidfVectorizer()tfidf.fit(df['Text'])
X = tfidf.transform(df['Text'])df['Text'][1]
In reference to the above sentence, we can check out tf-idf scores for a few words within this sentence.
print([X[1, tfidf.vocabulary_['peanuts']]])
[0.37995462060339136]
print([X[1, tfidf.vocabulary_['jumbo']]])
[0.530965343023095]
print([X[1, tfidf.vocabulary_['error']]])
[0.2302711360436964]
Among the three words, “peanut”, “jumbo” and “error”, tf-idf gives the highest weight to “jumbo”. Why? This indicates that “jumbo” is a much rarer word than “peanut” and “error”. This is how to use the tf-idf to indicate the importance of words or terms inside a collection of documents.
To classify sentiment, we remove neutral score 3, then group score 4 and 5 to positive (1), and score 1 and 2 to negative (0). After simple cleaning up, this is the data we are going to work with.
import numpy as npdf.dropna(inplace=True)df[df['Score'] != 3]df['Positivity'] = np.where(df['Score'] > 3, 1, 0)cols = ['Id', 'ProductId', 'UserId', 'ProfileName', 'HelpfulnessNumerator', 'HelpfulnessDenominator', 'Score', 'Time', 'Summary']df.drop(cols, axis=1, inplace=True)df.head()
from sklearn.model_selection import train_test_splitX = df.Texty = df.PositivityX_train, X_test, y_train, y_test = train_test_split(X, y, random_state = 0)print("Train set has total {0} entries with {1:.2f}% negative, {2:.2f}% positive".format(len(X_train),                                                                             (len(X_train[y_train == 0]) / (len(X_train)*1.))*100,                                                                            (len(X_train[y_train == 1]) / (len(X_train)*1.))*100))
Train set has total 426308 entries with 21.91% negative, 78.09% positive
print("Test set has total {0} entries with {1:.2f}% negative, {2:.2f}% positive".format(len(X_test),                                                                             (len(X_test[y_test == 0]) / (len(X_test)*1.))*100,                                                                            (len(X_test[y_test == 1]) / (len(X_test)*1.))*100))
Test set has total 142103 entries with 21.99% negative, 78.01% positive
You may have noticed that our classes are imbalanced, and the ratio of negative to positive instances is 22:78.
One of the tactics of combating imbalanced classes is using Decision Tree algorithms, so, we are using Random Forest classifier to learn imbalanced data and set class_weight=balanced .
First, define a function to print out the accuracy score.
from sklearn.feature_extraction.text import CountVectorizerfrom sklearn.ensemble import RandomForestClassifierfrom sklearn.pipeline import Pipelinefrom sklearn.metrics import accuracy_scoredef accuracy_summary(pipeline, X_train, y_train, X_test, y_test):    sentiment_fit = pipeline.fit(X_train, y_train)    y_pred = sentiment_fit.predict(X_test)    accuracy = accuracy_score(y_test, y_pred)    print("accuracy score: {0:.2f}%".format(accuracy*100))    return accuracy
To have efficient sentiment analysis or solving any NLP problem, we need a lot of features. Its not easy to figure out the exact number of features are needed. So we are going to try, 10,000 to 30,000. And print out accuracy scores associate with the number of features.
cv = CountVectorizer()rf = RandomForestClassifier(class_weight="balanced")n_features = np.arange(10000,30001,10000)def nfeature_accuracy_checker(vectorizer=cv, n_features=n_features, stop_words=None, ngram_range=(1, 1), classifier=rf):    result = []    print(classifier)    print("\n")    for n in n_features:        vectorizer.set_params(stop_words=stop_words, max_features=n, ngram_range=ngram_range)        checker_pipeline = Pipeline([            ('vectorizer', vectorizer),            ('classifier', classifier)        ])        print("Test result for {} features".format(n))        nfeature_accuracy = accuracy_summary(checker_pipeline, X_train, y_train, X_test, y_test)        result.append((n,nfeature_accuracy))    return resulttfidf = TfidfVectorizer()print("Result for trigram with stop words (Tfidf)\n")feature_result_tgt = nfeature_accuracy_checker(vectorizer=tfidf,ngram_range=(1, 3))
Nice!
Before we are done here, we should check the classification report.
from sklearn.metrics import classification_reportcv = CountVectorizer(max_features=30000,ngram_range=(1, 3))pipeline = Pipeline([        ('vectorizer', cv),        ('classifier', rf)    ])sentiment_fit = pipeline.fit(X_train, y_train)y_pred = sentiment_fit.predict(X_test)print(classification_report(y_test, y_pred, target_names=['negative','positive']))
Feature selection is an important problem in Machine learning. I will show you how straightforward it is to conduct Chi square test based feature selection on our large scale data set.
We will calculate the Chi square scores for all the features and visualize the top 20, here terms or words or N-grams are features, and positive and negative are two classes. given a feature X, we can use Chi square test to evaluate its importance to distinguish the class.
from sklearn.feature_selection import chi2import matplotlib.pyplot as plt%matplotlib inlineplt.figure(figsize=(12,8))scores = list(zip(tfidf.get_feature_names(), chi2score))chi2 = sorted(scores, key=lambda x:x[1])topchi2 = list(zip(*chi2[-20:]))x = range(len(topchi2[1]))labels = topchi2[0]plt.barh(x,topchi2[1], align='center', alpha=0.5)plt.plot(topchi2[1], x, '-o', markersize=5, alpha=0.8)plt.yticks(x, labels)plt.xlabel('$\chi^2$')plt.show();
We can observe that the features with a high χ2 can be considered relevant for the sentiment classes we are analyzing.
For example, the top 5 most useful feature selected by Chi-square test are “not”, “disappointed”, “very disappointed”, “not buy” and “worst”. I assume they are mostly from negative reviews. The next most useful feature selected by Chi-square test is “great”, I assume it is from mostly the positive reviews.
That’s it for today. Source code can be found on Github. I am happy to hear any questions or feedback.