聲明:代碼的運作環境為Python3。Python3與Python2在一些細節上會有所不同,希望廣大讀者注意。本部落格以代碼為主,代碼中會有詳細的注釋。相關文章将會釋出在我的個人部落格專欄《Python自然語言處理》,歡迎大家關注。
1、首先來看一個使用樸素貝葉斯分類器對性别進行分類鑒定的例子。
# 構造特征提取器
def gender_features(word): # 提取出字元串的最後一個字母
return {'last_letter': word[-1]}
names_set = ([(name, 'male') for name in names.words('male.txt')] +
[(name, 'female') for name in names.words('female.txt')]) # 擷取語料庫中人名的名稱及性别
print(names_set[:10])
random.shuffle(names_set) # 随機打亂names_set中的名字
print(names_set[:10])
featuresets = [(gender_features(n), g) for (n, g) in names_set] # 選取特征集合
train_set, test_set = featuresets[500:], featuresets[:500] # 抽取訓練資料和測試資料
classifier = nltk.NaiveBayesClassifier.train(train_set) # 樸素貝葉斯分類器
# 測試
classifier.classify(gender_features('Neo'))
classifier.classify(gender_features('Trinity'))
檢視測試結果:
'male'
'female'
2、計算分類的正确率以及最有效的特征
# 分類正确率判斷
print(nltk.classify.accuracy(classifier, test_set))
0.768
# 最有效的特征
classifier.show_most_informative_features(5) # 輸出5個最有效的特征
Most Informative Features
last_letter = 'a' female : male = 40.2 : 1.0
last_letter = 'k' male : female = 30.8 : 1.0
last_letter = 'f' male : female = 16.6 : 1.0
last_letter = 'p' male : female = 11.9 : 1.0
last_letter = 'v' male : female = 9.8 : 1.0
3、當資料量很大時可以用如下方法進行資料集的劃分。
# 大型資料時的資料集劃分
from nltk.classify import apply_features
train_set = apply_features(gender_features, names_set[500:])
test_set = apply_features(gender_features, names_set[:500])
4、貝葉斯公式及代碼實作。
貝葉斯公式如下圖所示:
貝葉斯分類器為:
用代碼的形式實作貝葉斯分類器:
# 手動計算貝葉斯分類器
# 計算P(特征|類别)
def f_c(data, fea, cla):
cfd = nltk.ConditionalFreqDist((classes, features) for (features, classes) in data)
return cfd[cla].freq(fea)
# 計算P(特征)
def p_feature(data, fea):
fd = nltk.FreqDist(fea for (fea, cla) in data)
return fd.freq(fea)
# 計算P(類别)
def p_class(data, cla):
fd = nltk.FreqDist(cla for (fea, cla) in data)
return fd.freq(cla)
# 計算P(類别│特征)
def res(data, fea, cla):
return f_c(data, fea, cla) * p_class(data, cla) / p_feature(data, fea)
測試貝葉斯公式:
# 構造輸入資料集
data = ([(name[-1], 'male') for name in names.words('male.txt')] +
[(name[-1], 'female') for name in names.words('female.txt')])
random.shuffle(data)
train, test = data[500:], data[:500]
# 計算Neo的為男性的機率
res(train, 'k', 'male')
res(train, 'a', 'female')
測試結果為:
0.955223880597015
0.9829612220916567
5、選擇正确的特征。
(1)過度拟合
# 過度拟合
def gender_features2(name):
features = {}
features["firstletter"] = name[0].lower()
features["lastletter"] = name[-1].lower()
for letter in 'abcdefghijklmnopqrstuvwxyz':
features["count(%s)" % letter] = name.lower().count(letter)
features["has(%s)" % letter] = (letter in name.lower())
return features
featuresets = [(gender_features2(n), g) for (n, g) in names_set]
train_set, test_set = featuresets[500:], featuresets[:500]
classifier = nltk.NaiveBayesClassifier.train(train_set)
print(nltk.classify.accuracy(classifier, test_set))
結果:
0.76
(2)劃分資料集,重新測試
# 資料劃分為訓練集、開發測試集、測試集
train_names = names_set[1500:]
devtest_names = names_set[500:1500]
test_names = names_set[:500]
# 重新訓練模型
train_set = [(gender_features(n), g) for (n, g) in train_names]
devtest_set = [(gender_features(n), g) for (n, g) in devtest_names]
test_set = [(gender_features(n), g) for (n, g) in test_names]
classifier = nltk.NaiveBayesClassifier.train(train_set)
print(nltk.classify.accuracy(classifier, devtest_set))
測試結果:
0.761
(3)将測試出錯的元素列印出來
# 列印錯誤清單
errors = []
for (name, tag) in devtest_names:
guess = classifier.classify(gender_features(name))
if guess != tag:
errors.append((tag, guess, name))
for (tag, guess, name) in sorted(errors): # doctest: +ELLIPSIS +NORMALIZE_WHITESPACE
print('correct=%-8s guess=%-8s name=%-30s' % (tag, guess, name))
結果:
correct=female guess=male name=Adelind
correct=female guess=male name=Aeriel
correct=female guess=male name=Aeriell
correct=female guess=male name=Ag
correct=female guess=male name=Aidan
correct=female guess=male name=Allsun
correct=female guess=male name=Anabel
correct=female guess=male name=Ardelis
correct=female guess=male name=Aryn
correct=female guess=male name=Betteann
correct=female guess=male name=Bill
correct=female guess=male name=Blondell
(4)重構特征,進行訓練
# 重新建構特征
def gender_features(word):
return {'suffix1': word[-1:],
'suffix2': word[-2:]}
# 重新訓練模型
train_set = [(gender_features(n), g) for (n, g) in train_names]
devtest_set = [(gender_features(n), g) for (n, g) in devtest_names]
classifier = nltk.NaiveBayesClassifier.train(train_set)
print(nltk.classify.accuracy(classifier, devtest_set))
測試結果:
0.791
6、文檔分類
from nltk.corpus import movie_reviews
documents = [(list(movie_reviews.words(fileid)), category)
for category in movie_reviews.categories()
for fileid in movie_reviews.fileids(category)]
random.shuffle(documents)
# 文檔分類特征提取器
all_words = nltk.FreqDist(w.lower() for w in movie_reviews.words())
word_features = all_words.most_common()[:2000]
def document_features(document):
document_words = set(document)
features = {}
for (word, freq) in word_features:
features['contains(%s)' % word] = (word in document_words)
return features
print(document_features(movie_reviews.words('pos/cv957_8737.txt')))
# 構造分類器
featuresets = [(document_features(d), c) for (d, c) in documents]
train_set, test_set = featuresets[100:], featuresets[:100]
classifier = nltk.NaiveBayesClassifier.train(train_set)
print(nltk.classify.accuracy(classifier, test_set))
print(classifier.show_most_informative_features(5))
測試結果:
0.86
Most Informative Features
contains(outstanding) = True pos : neg = 10.4 : 1.0
contains(seagal) = True neg : pos = 8.7 : 1.0
contains(mulan) = True pos : neg = 8.1 : 1.0
contains(wonderfully) = True pos : neg = 6.3 : 1.0
contains(damon) = True pos : neg = 5.7 : 1.0
None
7、詞性标注
# 詞性标注
from nltk.corpus import brown
suffix_fdist = nltk.FreqDist()
for word in brown.words():
word = word.lower()
suffix_fdist[word[-1:]] += 1
suffix_fdist[word[-2:]] += 1
suffix_fdist[word[-3:]] += 1
common_suffixes = suffix_fdist.most_common()[:100]
print(common_suffixes)
# 定義特征提取器
def pos_features(word):
features = {}
for (suffix, freq) in common_suffixes:
features['endswith(%s)' % suffix] = word.lower().endswith(suffix)
return features
# 訓練分類器
tagged_words = brown.tagged_words(categories='news')
featuresets = [(pos_features(n), g) for (n, g) in tagged_words]
size = int(len(featuresets) * 0.1)
train_set, test_set = featuresets[:1000], featuresets[2000:3000]
classifier = nltk.DecisionTreeClassifier.train(train_set)
print(nltk.classify.accuracy(classifier, test_set))
print(classifier.classify(pos_features('cats')))
結果:
0.611
NNS
決策樹輸出:
# 決策樹輸出
print(classifier.pseudocode(depth=4))
輸出結果:
if endswith(he) == False:
if endswith(s) == False:
if endswith(.) == False:
if endswith(of) == False: return 'CD'
if endswith(of) == True: return 'IN'
if endswith(.) == True: return '.'
if endswith(s) == True:
if endswith(as) == False:
if endswith('s) == False: return 'NPS'
if endswith('s) == True: return 'NN$'
if endswith(as) == True:
if endswith(was) == False: return 'CS'
if endswith(was) == True: return 'BEDZ'
if endswith(he) == True:
if endswith(the) == False: return 'PPS'
if endswith(the) == True: return 'AT'
8、根據上下文構造特征提取器
# 根據上下文構造特征提取器
def pos_features(sentence, i):
features = {"suffix(1)": sentence[i][-1:],
"suffix(2)": sentence[i][-2:],
"suffix(3)": sentence[i][-3:]}
if i == 0:
features["prev-word"] = "<START>"
else:
features["prev-word"] = sentence[i - 1]
return features
pos_features(brown.sents()[0], 8)
tagged_sents = brown.tagged_sents(categories='news')
featuresets = []
for tagged_sent in tagged_sents:
untagged_sent = nltk.tag.untag(tagged_sent)
for i, (word, tag) in enumerate(tagged_sent):
featuresets.append((pos_features(untagged_sent, i), tag))
size = int(len(featuresets) * 0.1)
train_set, test_set = featuresets[size:], featuresets[:size]
classifier = nltk.NaiveBayesClassifier.train(train_set)
nltk.classify.accuracy(classifier, test_set)
結果:
0.7891596220785678
9、序列分類
# 序列分類
# 定義特征提取器
def pos_features(sentence, i, history):
features = {"suffix(1)": sentence[i][-1:],
"suffix(2)": sentence[i][-2:],
"suffix(3)": sentence[i][-3:]}
if i == 0:
features["prev-word"] = "<START>"
features["prev-tag"] = "<START>"
else:
features["prev-word"] = sentence[i - 1]
features["prev-tag"] = history[i - 1]
return features
# 建構序列分類器
class ConsecutivePosTagger(nltk.TaggerI):
def __init__(self, train_sents):
train_set = []
for tagged_sent in train_sents:
untagged_sent = nltk.tag.untag(tagged_sent)
history = []
for i, (word, tag) in enumerate(tagged_sent):
featureset = pos_features(untagged_sent, i, history)
train_set.append((featureset, tag))
history.append(tag)
self.classifier = nltk.NaiveBayesClassifier.train(train_set)
def tag(self, sentence):
history = []
for i, word in enumerate(sentence):
featureset = pos_features(sentence, i, history)
tag = self.classifier.classify(featureset)
history.append(tag)
return zip(sentence, history)
tagged_sents = brown.tagged_sents(categories='news')
size = int(len(tagged_sents) * 0.1)
train_sents, test_sents = tagged_sents[size:], tagged_sents[:size]
tagger = ConsecutivePosTagger(train_sents)
print(tagger.evaluate(test_sents))
結果:
0.7980528511821975
10、句子分割
# 句子分割
# 擷取已分割的句子資料
sents = nltk.corpus.treebank_raw.sents()
tokens = []
boundaries = set()
offset = 0
for sent in nltk.corpus.treebank_raw.sents():
tokens.extend(sent)
offset += len(sent)
boundaries.add(offset - 1)
# 定義特征提取器
def punct_features(tokens, i):
return {'next-word-capitalized': tokens[i + 1][0].isupper(),
'prevword': tokens[i - 1].lower(),
'punct': tokens[i],
'prev-word-is-one-char': len(tokens[i - 1]) == 1}
# 定義标注
featuresets = [(punct_features(tokens, i), (i in boundaries))
for i in range(1, len(tokens) - 1)
if tokens[i] in '.?!']
# 建構分類器
size = int(len(featuresets) * 0.1)
train_set, test_set = featuresets[size:], featuresets[:size]
classifier = nltk.NaiveBayesClassifier.train(train_set)
nltk.classify.accuracy(classifier, test_set)
結果:
0.936026936026936
基于分類的斷句器:
# 基于分類的斷句器
def segment_sentences(words):
start = 0
sents = []
for i, word in words:
if word in '.?!' and classifier.classify(words, i) == True:
sents.append(words[start:i + 1])
start = i + 1
if start < len(words):
sents.append(words[start:])
11、識别對話行為類型
# 識别對話行為類型
posts = nltk.corpus.nps_chat.xml_posts()[:10000]
# 定義特征提取器
def dialogue_act_features(post):
features = {}
for word in nltk.word_tokenize(post):
features['contains(%s)' % word.lower()] = True
return features
# 訓練分類器
featuresets = [(dialogue_act_features(post.text), post.get('class'))
for post in posts]
size = int(len(featuresets) * 0.1)
train_set, test_set = featuresets[size:], featuresets[:size]
classifier = nltk.NaiveBayesClassifier.train(train_set)
print(nltk.classify.accuracy(classifier, test_set))
結果:
0.668
12、評估
(1)準确度
####評估####
# 建立訓練集與測試集
import random
from nltk.corpus import brown
tagged_sents = list(brown.tagged_sents(categories='news'))
random.shuffle(tagged_sents)
size = int(len(tagged_sents) * 0.1)
train_set, test_set = tagged_sents[size:], tagged_sents[:size]
# 使用同類型檔案
file_ids = brown.fileids(categories='news')
size = int(len(file_ids) * 0.1)
train_set = brown.tagged_sents(file_ids[size:])
test_set = brown.tagged_sents(file_ids[:size])
# 使用不同類型檔案
train_set = brown.tagged_sents(categories='news')
test_set = brown.tagged_sents(categories='fiction')
##準确度##
names_set = ([(name, 'male') for name in names.words('male.txt')] +
[(name, 'female') for name in names.words('female.txt')])
random.shuffle(names_set)
featuresets = [(gender_features(n), g) for (n, g) in names_set]
train_set, test_set = featuresets[500:], featuresets[:500]
classifier = nltk.NaiveBayesClassifier.train(train_set)
print('Accuracy: %4.2f' % nltk.classify.accuracy(classifier, test_set))
結果:
Accuracy: 0.78
(2)精确度與召回率
##精準度與召回率##
from sklearn.metrics import classification_report
test_set_fea = [features for (features, gender) in test_set]
test_set_gen = [gender for (features, gender) in test_set]
pre = classifier.classify_many(test_set_fea)
print(classification_report(test_set_gen, pre))
結果:
precision recall f1-score support
female 0.83 0.82 0.83 316
male 0.70 0.71 0.70 184
avg / total 0.78 0.78 0.78 500
(3)混淆矩陣
##混淆矩陣##
cm = nltk.ConfusionMatrix(test_set_gen, pre)
print(cm)
結果:
| f |
| e |
| m m |
| a a |
| l l |
| e e |
-------+---------+
female |<260> 56 |
male | 54<130>|
-------+---------+
(row = reference; col = test)
(4)決策樹中的熵和資訊增益
# 熵和資訊增益
import math
def entropy(labels):
freqdist = nltk.FreqDist(labels)
probs = [freqdist.freq(l) for l in nltk.FreqDist(labels)]
return -sum([p * math.log(p, 2) for p in probs])
print(entropy(['male', 'male', 'male', 'male']))
-0.0
print(entropy(['male', 'female', 'male', 'male']))
0.8112781244591328
print(entropy(['female', 'male', 'female', 'male']))
1.0
print(entropy(['female', 'female', 'male', 'female']))
0.8112781244591328
print(entropy(['female', 'female', 'female', 'female']))
-0.0