赞
踩
importxlrdimportjiebaimportsysimportimportlibimport os #python内置的包,用于进行文件目录操作,我们将会用到os.listdir函数
import pickle #导入cPickle包并且取一个别名pickle #持久化类
importrandomimportnumpy as npimportmatplotlib.pyplot as pltfrom mpl_toolkits.mplot3d importAxes3Dfrom pylab importmplfrom sklearn.naive_bayes import MultinomialNB #导入多项式贝叶斯算法包
from sklearn importsvmfrom sklearn importmetricsfrom sklearn.datasets.base importBunchfrom sklearn.feature_extraction.text importTfidfVectorizer
importlib.reload(sys)#把内容和类别转化成一个向量的形式
trainContentdatasave=[] #存储所有训练和测试数据的分词
testContentdatasave=[]
trainContentdata=[]
testContentdata=[]
trainlabeldata=[]
testlabeldata=[]#导入文本描述的训练和测试数据
defimportTrainContentdata():
file= '20180716_train.xls'wb=xlrd.open_workbook(file)
ws= wb.sheet_by_name("Sheet1")for r inrange(ws.nrows):
col=[]for c in range(1):
col.append(ws.cell(r, c).value)
trainContentdata.append(col)defimportTestContentdata():
file= '20180716_test.xls'wb=xlrd.open_workbook(file)
ws= wb.sheet_by_name("Sheet1")for r inrange(ws.nrows):
col=[]for c in range(1):
col.append(ws.cell(r, c).value)
testContentdata.append(col)#导入类别的训练和测试数据
defimportTrainlabeldata():
file= '20180716_train_label.xls'wb=xlrd.open_workbook(file)
ws= wb.sheet_by_name("Sheet1")for r inrange(ws.nrows):
col=[]for c in range(1):
col.append(ws.cell(r, c).value)
trainlabeldata.append(col)defimportTestlabeldata():
file= '20180716_test_label.xls'wb=xlrd.open_workbook(file)
ws= wb.sheet_by_name("Sheet1")for r inrange(ws.nrows):
col=[]for c in range(1):
col.append(ws.cell(r, c).value)
testlabeldata.append(col)"""def importClassSet():
file = 'ClassSet.xls'
wb = xlrd.open_workbook(file)
ws = wb.sheet_by_name("Sheet1")
for r in range(ws.nrows):
col = []
for c in range(ws.ncols):
col.append(ws.cell(r, c).value)
ClassSet.append(col)"""
defbuildtrainbunch(bunch_path):
bunch= Bunch(label=[],contents=[])for item1 intrainlabeldata:
bunch.label.append(item1)for item2 intrainContentdata:
item2=str(item2)
item2= item2.replace("\r\n", "")
item2= item2.replace(" ", "")
content_seg=jieba.cut(item2)
save2=''
for item3 incontent_seg:if len(item3) > 1 and item3!='\r\n':
trainContentdatasave.append(item3)
save2=save2+","+item3
bunch.contents.append(save2)
with open(bunch_path,"wb") as file_obj:
pickle.dump(bunch, file_obj)print("构建训练数据文本对象结束!!!")defbuildtestbunch(bunch_path):
bunch= Bunch(label=[],contents=[])for item1 intestlabeldata:
bunch.label.append(item1)for item2 intestContentdata:
item2=str(item2)
item2= item2.replace("\r\n", "")
item2= item2.replace(" ", "")
content_seg=jieba.cut(item2)
save2=''
for item3 incontent_seg:if len(item3) > 1 and item3!='\r\n':
testContentdatasave.append(item3)
save2=save2+","+item3
bunch.contents.append(save2)
with open(bunch_path,"wb") as file_obj:
pickle.dump(bunch, file_obj)print("构建测试数据文本对象结束!!!")#读取停用词
def_readfile(path):
with open(path,"rb") as fp:
content=fp.read()returncontent#读取bunch对象
def_readbunchobj(path):
with open(path,"rb") as file_obj:
bunch=pickle.load(file_obj)returnbunch#写入bunch对象
def_writebunchobj(path, bunchobj):
with open(path,"wb") as file_obj:
pickle.dump(bunchobj, file_obj)defvector_space(stopword_path,bunch_path,space_path):
stpwrdlst= _readfile(stopword_path).splitlines()#读取停用词
bunch = _readbunchobj(bunch_path)#导入分词后的词向量bunch对象
#构建tf-idf词向量空间对象
tfidfspace = Bunch(label=bunch.label,tdm=[], vocabulary={})'''权重矩阵tdm,其中,权重矩阵是一个二维矩阵,tdm[i][j]表示,第j个词(即词典中的序号)在第i个类别中的IF-IDF值'''
#使用TfidVectorizer初始化向量空间模型
vectorizer = TfidfVectorizer(stop_words=stpwrdlst, sublinear_tf=True, max_df=0.5, min_df=0.0001,use_idf=False,max_features=10000)#print(vectorizer)
#文本转为词频矩阵,单独保存字典文件
tfidfspace.tdm =vectorizer.fit_transform(bunch.contents)
tfidfspace.vocabulary=vectorizer.vocabulary_#创建词袋的持久化
_writebunchobj(space_path, tfidfspace)print("if-idf词向量空间实例创建成功!!!")deftestvector_space(stopword_path,bunch_path,space_path,train_tfidf_path):
stpwrdlst= _readfile(stopword_path).splitlines()#把停用词变成列表
bunch =_readbunchobj(bunch_path)
tfidfspace= Bunch(label=bunch.label,tdm=[], vocabulary={})'''tdm存放的是计算后得到的TF-IDF权重矩阵.
vocabulary是词向量空间的索引,例如,如果我们定义的词向量空间是(我,喜欢,相国大人),那么vocabulary就是这样一个索引字典
vocabulary={"我":0,"喜欢":1,"相国大人":2},你可以简单的理解为:vocabulary就是词向量空间的坐标轴,索引值相当于表明了第几个维度。'''
#导入训练集的TF-IDF词向量空间 ★★
trainbunch =_readbunchobj(train_tfidf_path)
tfidfspace.vocabulary=trainbunch.vocabulary'''关于参数,你只需要了解这么几个就可以了:
stop_words:
传入停用词,以后我们获得vocabulary_的时候,就会根据文本信息去掉停用词得到
vocabulary:
之前说过,不再解释。
sublinear_tf:
计算tf值采用亚线性策略。比如,我们以前算tf是词频,现在用1+log(tf)来充当词频。
smooth_idf:
计算idf的时候log(分子/分母)分母有可能是0,smooth_idf会采用log(分子/(1+分母))的方式解决。默认已经开启,无需关心。
norm:
归一化,我们计算TF-IDF的时候,是用TF*IDF,TF可以是归一化的,也可以是没有归一化的,一般都是采用归一化的方法,默认开启.
max_df:
有些词,他们的文档频率太高了(一个词如果每篇文档都出现,那还有必要用它来区分文本类别吗?当然不用了呀),所以,我们可以
设定一个阈值,比如float类型0.5(取值范围[0.0,1.0]),表示这个词如果在整个数据集中超过50%的文本都出现了,那么我们也把它列
为临时停用词。当然你也可以设定为int型,例如max_df=10,表示这个词如果在整个数据集中超过10的文本都出现了,那么我们也把它列
为临时停用词。
min_df:
与max_df相反,虽然文档频率越低,似乎越能区分文本,可是如果太低,例如10000篇文本中只有1篇文本出现过这个词,仅仅因为这1篇
文本,就增加了词向量空间的维度,太不划算。
当然,max_df和min_df在给定vocabulary参数时,就失效了。'''vectorizer= TfidfVectorizer(stop_words=stpwrdlst, sublinear_tf=True, max_df=0.7, vocabulary=trainbunch.vocabulary, min_df=0.001)#print(vectorizer)
tfidfspace.tdm=vectorizer.fit_transform(bunch.contents)
_writebunchobj(space_path, tfidfspace)print("if-idf词向量空间实例创建成功!!!")def metrics_result(actual, predict): #metrics.f1_score(y_test, y_pred, average='weighted', labels=np.unique(y_pred))
print('精度:{0:.3f}'.format(metrics.precision_score(actual, predict,average='weighted', labels=np.unique(predict))))print('召回:{0:0.3f}'.format(metrics.recall_score(actual, predict,average='weighted', labels=np.unique(predict))))print('f1-score:{0:.3f}'.format(metrics.f1_score(actual, predict, average='weighted', labels=np.unique(predict))))#准确率和召回率是相互影响的,理想情况下是二者都高,但是一般情况下准确率高,召回率就低;召回率高,准确率就低
if __name__=="__main__":
importTrainContentdata()
importTestContentdata()
importTrainlabeldata()
importTestlabeldata()#导入分词后的词向量bunch对象
train_bunch_path ="F:/goverment/ArticleMining/trainbunch.bat"#Bunch保存路径
test_bunch_path ="F:/goverment/ArticleMining/testbunch.bat"stopword_path="F:/goverment/ArticleMining/hlt_stop_words.txt"train_space_path= "F:/goverment/ArticleMining/traintfdifspace.dat"test_space_path= "F:/goverment/ArticleMining/testtfdifspace.dat"
#对训练和测试集进行bunch操作
buildtrainbunch(train_bunch_path)
buildtestbunch(test_bunch_path)
vector_space(stopword_path,train_bunch_path,train_space_path)
testvector_space(stopword_path,test_bunch_path,test_space_path,train_space_path)#导入训练和测试数据集
train_set=_readbunchobj(train_space_path)
test_set=_readbunchobj(test_space_path)print(train_set.tdm)'''mm=0
ii=0
jj=0
for i in range(3142):
for j in range(3142):
if train_set.tdm[i][j] >mm:
mm=train_set.tdm[i][j]
ii=i
jj=j
print(ii)
print(jj)'''
#test_set.tdm
#train_set.label
#训练分类器:输入词袋向量和分类标签,alpha:0.001 alpha越小,迭代次数越多,精度越高
#低召回、F1: 0.75 rbf:0.59 0.8 rbf 0.578
#c0.75 poly 66.5 精度:0.665 gamma=10 召回:0.330 f1-score:0.416
#C=0.7, kernel='poly', gamma=10 召回:0.331 f1-score:0.417
#alpha:0.001 alpha 越小,迭代次数越多,精度越高
'''clf = MultinomialNB(alpha=0.052).fit(train_set.tdm, train_set.label)
#clf = svm.SVC(C=0.7, kernel='poly', gamma=10, decision_function_shape='ovr')
clf.fit(train_set.tdm, train_set.label)
predicted=clf.predict(test_set.tdm)
tv = TfidfVectorizer()
train_data = tv.fit_transform(X_train)
test_data = tv.transform(X_test)
lr = LogisticRegression(C=3)
lr.fit(train_set.tdm, train_set.label)
predicted=lr.predict(test_set.tdm)
print(lr.score(test_set.tdm, test_set.label))
#print(test_set.tdm)'''clf= SVC(C=1500)
clf.fit(train_set.tdm, train_set.label)
predicted=clf.predict(test_set.tdm)print(clf.score(test_set.tdm, test_set.label))'''from sklearn.neighbors import KNeighborsClassifier
knnclf = KNeighborsClassifier(n_neighbors=9)#default with k=5
knnclf.fit(train_set.tdm,train_set.label)
predicted = knnclf.predict(test_set.tdm)'''a=[]
b=[]for i inrange(len(predicted)):
b.append((int)(float(predicted[i])))
a.append(int(test_set.label[i][0]))
f=open('F:/goverment/ArticleMining/predict.txt', 'w')for i inrange(len(predicted)):
f.write(str(b[i]))
f.write('\n')
f.write("写好了")
f.close()#for i in range(len(predicted)):
#print(b[i])
metrics_result(a, b)
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。