声明:代码的运行环境为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