《实战》基于电商领域的词性提取及其决策树模型建模

最新推荐文章于 2024-11-08 03:03:20 发布

驭风少年君

最新推荐文章于 2024-11-08 03:03:20 发布

阅读量1.3k

点赞数 1

分类专栏：深度学习与学术学习实战文章标签：决策树机器学习 python

本文链接：https://blog.csdn.net/qq_44951759/article/details/126001215

版权

深度学习与学术学习实战专栏收录该内容

5 篇文章 2 订阅

订阅专栏

——基于电商领域的文本情感分析及主题挖掘——

一、数据预处理
二、词性提取
三、决策树分类模型

一、数据预处理

本次的实现目标
1、利用文本挖掘技术，对碎片化、非结构化的电商网站评论数据进行清洗与处理，转化为结构化数据。
2、参考知网发布的情感分析用词语集，统计评论数据的正负情感指数，然后进行情感分析,通过词云图直观查看正负评论的关键词。
3、比较“机器挖掘的正负情感”与“人工打标签的正负情感”。
4、采用LDA主题模型提取评论关键信息，以了解用户的需求、意见、购买原因、产品的优缺点等。

导入所需要的工具包

import os
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns
from matplotlib.pylab import style #自定义图表风格
style.use('ggplot')
from IPython.core.interactiveshell import InteractiveShell
InteractiveShell.ast_node_interactivity = "all"
plt.rcParams['font.sans-serif'] = ['Simhei'] # 解决中文乱码问题

import re
import jieba.posseg as psg
import itertools
#conda install -c anaconda gensim
from gensim import corpora,models #主题挖掘，提取关键信息

# pip install wordcloud
from wordcloud import WordCloud,ImageColorGenerator
from collections import Counter

from sklearn import tree
from sklearn.model_selection import train_test_split
from sklearn.feature_extraction.text import CountVectorizer
from sklearn.metrics import classification_report
from sklearn.metrics import accuracy_score
import graphviz

导入数据

raw_data=pd.read_csv('./reviews.csv')
raw_data.head()

主要有评论文本、时间、用户名、商品名称、用户情感倾向
在这里插入图片描述
查看数据大致分布，有无缺失值

raw_data.info()

<class ‘pandas.core.frame.DataFrame’>
RangeIndex: 2000 entries, 0 to 1999
Data columns (total 5 columns):
content 2000 non-null object
creationTime 2000 non-null object
nickname 2000 non-null object
referenceName 2000 non-null object
content_type 2000 non-null object
dtypes: object(5)
memory usage: 78.2+ KB

查看数据列名

raw_data.columns

Index([‘content’, ‘creationTime’, ‘nickname’, ‘referenceName’, ‘content_type’], dtype=‘object’)

常看每列下，样本统计

#取值分布
for cate in ['creationTime', 'nickname', 'referenceName', 'content_type']:
    raw_data[cate].value_counts()
#取值分布

2016-06-24 22:35:42 2
2015-06-25 17:36:36 2
2016-06-24 17:42:26 2
2016-06-20 16:03:00 2
2017-10-06 13:12:19 1
…
2017-10-22 08:08:35 1
2017-03-25 15:59:46 1
2017-10-01 10:13:38 1
2017-07-03 13:10:14 1
2017-03-19 10:19:17 1
Name: creationTime, Length: 1996, dtype: int64
j1 34
j6 30
j2 25
jb 25
jf 23
…
a1 1
切药 1
无主 1
小热 1
费能 1
Name: nickname, Length: 1190, dtype: int64
美的（Midea）60升预约洗浴无线遥控电热水器 F60-15WB5(Y) 2000
Name: referenceName, dtype: int64
neg 1000
pos 1000
Name: content_type, dtype: int64

二、词性提取

2.1 去重

先将数据复制一份，再对评论与情感倾向做去重

reviews=raw_data.copy()
reviews=reviews[['content', 'content_type']]
print('去重之前：',reviews.shape[0])
reviews=reviews.drop_duplicates()
print('去重之后：',reviews.shape[0])

去重之前： 2000
去重之后： 1974

2.2 数据清洗

# 清洗之前
content=reviews['content']
for i in range(5,10):
    print(content[i])
    print('-----------')

物美价廉啊，特别划算的，而且加热速度快。家里用着不错特别方便
-----------
价格合理，配置挺高，物美价值
-----------
老师按装是快的，装修中。没试。希望是正常的！
-----------
五分是习惯
送的快，装的也快，很好，不过装修进行中，为了吊顶只装了热水器，其它后话吧。
我也是醉了，热水器遥控器没电池，还是安装哥车里有一个备用的给装上了。
要赠送的电动牙刷，空气净化器等待配送中&hellip;&hellip;。
一天时间冰箱降100、洗衣机降40、热水器降100、小厨宝降30、烟灶套降100。不过客服说上报解决，电话通知等待中&hellip;&hellip;
-----------
安装的小哥非常好，工作很尽心，我们家是老房子，安装比一般家要费力。安装的非常不错。因为整栋楼要用热水，裙子只能这样装。给小哥点个赞。热水器已经买二个了，一如既往地好用。美的大品牌，质量非常好
-----------

将文本中的品牌名称与数字字母，用正则清洗掉

#清洗之后，将数字、字母、京东美的电热水器字样都删除
info=re.compile('[0-9a-zA-Z]|京东|美的|电热水器|热水器|')
content=content.apply(lambda x: info.sub('',x))  #替换所有匹配项

for i in range(5,10):
    print(content[i])
    print('-----------')

物美价廉啊，特别划算的，而且加热速度快。家里用着不错特别方便
-----------
价格合理，配置挺高，物美价值
-----------
老师按装是快的，装修中。没试。希望是正常的！
-----------
五分是习惯
送的快，装的也快，很好，不过装修进行中，为了吊顶只装了，其它后话吧。
我也是醉了，遥控器没电池，还是安装哥车里有一个备用的给装上了。
要赠送的电动牙刷，空气净化器等待配送中&;&;。
一天时间冰箱降、洗衣机降、降、小厨宝降、烟灶套降。不过客服说上报解决，电话通知等待中&;&;
-----------
安装的小哥非常好，工作很尽心，我们家是老房子，安装比一般家要费力。安装的非常不错。因为整栋楼要用热水，裙子只能这样装。给小哥点个赞。已经买二个了，一如既往地好用。大品牌，质量非常好
-----------

清洗后，文本中的数字、字母都已经清洗掉了

清洗干净的文本保存到content变量中

print(len(content))

1974

2.3 分词、词性标注、去除停用词、词云图

2.3.1 提取出每条文本的词及其词性

输入：
- content、content_type
- 共有1974条评论句子
输出：
- 构造DF，包含：分词、对应词性、分词所在原句子的id、分词所在原句子的content_type
- 共有6万多行

非结构化数据——>结构化数据

在这里插入图片描述
将每条文本进行分词，并且去除每条文本的词，和对应词的词性。主要通过import jieba.posseg as psg 实现

案例

def cut_word(s):
    for x in psg.cut(s):
        print((x.word,x.flag))

cut_word('物美价廉啊，特别划算的，而且加热速度快。家里用着不错特别方便')

输出结果

('物美价廉', 'l')
('啊', 'zg')
('，', 'x')
('特别', 'd')
('划算', 'v')
('的', 'uj')
('，', 'x')
('而且', 'c')
('加热', 'v')
('速度', 'n')
('快', 'a')
('。', 'x')
('家里', 's')
('用', 'p')
('着', 'uz')
('不错', 'a')
('特别', 'd')
('方便', 'a')

实际操作用匿名函数实现

#分词，由元组组成的list
seg_content=content.apply( lambda s:  [(x.word,x.flag) for x in psg.cut(s)] )  

seg_content.shape
len(seg_content)
print(seg_content[5])

(1974,)
1974
[(‘物美价廉’, ‘l’), (‘啊’, ‘zg’), (‘，’, ‘x’), (‘特别’, ‘d’), (‘划算’, ‘v’), (‘的’, ‘uj’), (‘，’, ‘x’), (‘而且’, ‘c’), (‘加热’, ‘v’), (‘速度’, ‘n’), (‘快’, ‘a’), (‘。’, ‘x’), (‘家里’, ‘s’), (‘用’, ‘p’), (‘着’, ‘uz’), (‘不错’, ‘a’), (‘特别’, ‘d’), (‘方便’, ‘a’)]

统计词评论数

统计出每条评论中，词+符号出现的总数，保存到n_word中

#统计评论词数
n_word=seg_content.apply(lambda s: len(s))

len(n_word)
n_word.head(6)

1974
0 32
1 11
2 6
3 39
4 44
5 18
Name: content, dtype: int64

构建每个词，所对应的评论

示例
seg_content.index评论序列

Int64Index([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
…
1990, 1991, 1992, 1993, 1994, 1995, 1996, 1997, 1998, 1999],
dtype=‘int64’, length=1974)

for x,y in zip(list(seg_content.index),list(n_word)):
    print([x+1]*y)

[1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
[2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2]
[3, 3, 3, 3, 3, 3]
[4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4]
[5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5]
[6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6]
[7, 7, 7, 7, 7, 7, 7, 7]
[8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]
[9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9]
[10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10]
[11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11]
[12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12]

sum([[2,2],[3,3,3]],[])

脱去列表

[2, 2, 3, 3, 3]

实际应用

#得到各分词在第几条评论
n_content=[ [x+1]*y for x,y in zip(list(seg_content.index),list(n_word))] #[x+1]*y,表示复制y份，由list组成的list
print('-----------得到词对应评论序号---------')
print(n_content[:2])
index_content_long=sum(n_content,[]) #表示去掉[]，拉平,返回list
len(index_content_long)
print('----------去掉[]，将序号列表拉伸-----------')
print(index_content_long[:43])

-----------得到词对应评论序号---------
[[1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1], [2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2]]
63794
----------去掉[]，将序号列表拉伸-----------
[1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2]

分词及词性，去掉[]，拉平

将分词性后的评论也去掉[],合并

#分词及词性，去掉[]，拉平
print('原始数据：','\n',seg_content.head())

seg_content_long=sum(seg_content,[])

type(seg_content_long)
print('去掉[],数据：','\n',seg_content_long[:5])
len(seg_content_long)

原始数据：
0 [(东西, ns), (收到, v), (这么久, r), (，, x), (都, d), …
1 [(安装, v), (师傅, nr), (很, d), (给, p), (力, n), (，…
2 [(还, d), (没, v), (安装, v), (，, x), (基本, n), (满意…
3 [(收到, v), (了, ul), (，, x), (自营, vn), (商品, n), …
4 [(用, p), (了, ul), (几次, m), (才, d), (来, v), (评价…
Name: content, dtype: object
list
去掉[],数据：
[(‘东西’, ‘ns’), (‘收到’, ‘v’), (‘这么久’, ‘r’), (‘，’, ‘x’), (‘都’, ‘d’)]
63794

分别取得词和其词性

#得到加长版的分词、词性
word_long=[x[0] for x in seg_content_long]
nature_long=[x[1] for x in seg_content_long]

print(word_long[:5])
print(nature_long[:5])
len(word_long)
len(nature_long)

[‘东西’, ‘收到’, ‘这么久’, ‘，’, ‘都’]
[‘ns’, ‘v’, ‘r’, ‘x’, ‘d’]
63794
63794

content_type拉长

#content_type拉长
n_content_type=[ [x]*y for x,y in zip(list(reviews['content_type']),list(n_word))] #[x+1]*y,表示复制y份
content_type_long=sum(n_content_type,[]) #表示去掉[]，拉平

len(content_type_long)
print('情感倾向：',content_type_long[:5])

63794
情感倾向： [‘pos’, ‘pos’, ‘pos’, ‘pos’, ‘pos’]

保存到Dataframe中

review_long=pd.DataFrame({'index_content':index_content_long,
                        'word':word_long,
                        'nature':nature_long,
                        'content_type':content_type_long})
review_long.shape
review_long.head()

(63794, 4)

2.3.2 去除标点符号、去除停用词

review_long['nature'].unique()

array([‘ns’, ‘v’, ‘r’, ‘x’, ‘d’, ‘ul’, ‘a’, ‘n’, ‘u’, ‘nr’, ‘p’, ‘y’,
‘vn’, ‘t’, ‘c’, ‘m’, ‘l’, ‘b’, ‘i’, ‘uj’, ‘zg’, ‘s’, ‘uz’, ‘nz’,
‘f’, ‘uv’, ‘ad’, ‘q’, ‘j’, ‘g’, ‘ud’, ‘an’, ‘nrt’, ‘vg’, ‘ng’, ‘k’,
‘o’, ‘mq’, ‘df’, ‘e’, ‘vd’, ‘z’, ‘nt’, ‘tg’, ‘rz’, ‘ug’, ‘yg’, ‘h’,
‘vq’, ‘ag’, ‘rr’], dtype=object)

其中x表示标点符号

#去除标点符号
review_long_clean=review_long[review_long['nature']!='x'] #x表示标点符合
review_long_clean.shape

(51436, 4)

原本有 63794，去除后只有51436行。

#导入停用词
stop_path=open('./stoplist.txt','r',encoding='UTF-8')
stop_words=stop_path.readlines()

len(stop_words)
stop_words[0:5]

#停用词，预处理,去除换行符号
stop_words=[word.strip('\n') for word in stop_words]
stop_words[0:5]

5748
[‘\ufeff \n’, ‘说\n’, ‘人\n’, ‘元\n’, ‘hellip\n’]
['\ufeff ', ‘说’, ‘人’, ‘元’, ‘hellip’]

#得到不含停用词的分词表
word_long_clean=list(set(word_long)-set(stop_words))
len(word_long_clean)

review_long_clean=review_long_clean[review_long_clean['word'].isin(word_long_clean)]
review_long_clean.shape

去除评论中的停用词

#得到不含停用词的分词表
word_long_clean=list(set(word_long)-set(stop_words))
print('文本中全部词一共：',len(set(word_long)))
print('停用词全部词一共：',len(set(stop_words)))
print('去除停用词后，无重复一共：',len(set(word_long_clean)))

review_long_clean=review_long_clean[review_long_clean['word'].isin(word_long_clean)]
print('去除停用词后，剩下：',len(review_long_clean))

文本中全部词一共： 5135
停用词全部词一共： 1905
去除停用词后，无重复一共： 4455
去除停用词后，剩下： 25172

2.3.3 在原df中，再增加一列，该分词在本条评论的位置

```python
#再次统计每条评论的分词数量
n_word=review_long_clean.groupby('index_content').count()['word']

print('去除停用词后，每条文本的分词数量：','\n',n_word)

index_word=[ list(np.arange(1,x+1)) for x in list(n_word)] 
index_word_long=sum(index_word,[]) #表示去掉[]，拉平

len(index_word_long)
print('去除停用词后，每条文本的从1开始编号：','\n',index_word_long[:20])

去除停用词后，每条文本的分词数量：
index_content
1 14
2 4
3 2
4 20
5 23
…
1996 6
1997 20
1998 3
1999 10
2000 4
Name: word, Length: 1964, dtype: int64
25172
去除停用词后，每条文本的从1开始编号：
[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 1, 2, 3, 4, 1, 2]

添加到df中进去

review_long_clean['index_word']=index_word_long
review_long_clean.head()

index_content表示评论编号，index_word表示评论编号

在这里插入图片描述

review_long_clean.to_csv('./1_review_long_clean.csv')   #保存为csv格式
review_long_clean.to_excel('./1_review_long_clean.xlsx')#保存为xlsx格式

2.3.4 提取名词

提取出文本中的名词

n_review_long_clean=review_long_clean[[ 'n' in nat for nat in review_long_clean.nature]]
n_review_long_clean.shape
n_review_long_clean.head()

(10189, 5)

统计，并且保存

n_review_long_clean.nature.value_counts()
n_review_long_clean.to_csv('./1_n_review_long_clean.csv')

2.3.5 绘制词云

全部词的词云图

Counter(review_long_clean.word.values)

Counter({‘东西’: 224,
‘收到’: 85,
‘这么久’: 4,
‘忘’: 5,
‘好评’: 81,
‘品牌’: 120,
‘值得’: 110,
‘信赖’: 100,
‘整体’: 14,
‘个人感觉’: 1,
‘不错’: 309,
‘拥有’: 13,
‘安装’: 1626,
‘师傅’: 502,
‘力’: 68,
‘感谢’: 14,
‘满意’: 185,

font=r"C:\Windows\Fonts\msyh.ttc"

background_image=plt.imread('./pl.jpg')
wordcloud = WordCloud(font_path=font, 
                      max_words = 100,
                      mask=background_image,
                      background_color='white'
                      ) #width=1600,height=1200, mode='RGBA'
wordcloud.generate_from_frequencies(Counter(review_long_clean.word.values))
wordcloud.to_file('1_分词后的词云图.png')

plt.figure(figsize=(20,10))
plt.imshow(wordcloud)
plt.axis('off')
plt.show()

名词的词云

Counter(n_review_long_clean.word.values)

Counter({‘东西’: 224,
‘品牌’: 120,
‘信赖’: 100,
‘整体’: 14,
‘个人感觉’: 1,
‘师傅’: 502,
‘力’: 68,
‘自营’: 23,
‘商品’: 47,
‘发货’: 33,
‘速度’: 129,
‘品质’: 14,
‘效果’: 36,
‘宝贝’: 15,
‘小时’: 85,
‘热水’: 83,
‘评价’: 48,
‘产品’: 118,

font=r"C:\Windows\Fonts\msyh.ttc"

background_image=plt.imread('./pl.jpg')
wordcloud = WordCloud(font_path=font,
                      max_words = 100, 
                      mode='RGBA' ,
                      background_color='white',
                      mask=background_image) #width=1600,height=1200
wordcloud.generate_from_frequencies(Counter(n_review_long_clean.word.values))
wordcloud.to_file('1_分词后的词云图(名词）.png')

plt.figure(figsize=(20,10))
plt.imshow(wordcloud)
plt.axis('off')
plt.show()

三、决策树分类模型

3.1 构造特征空间和标签

将提纯后每条评论的词，再合并起来

#第一步：构造特征空间和标签

Y=[]
for ind in review_long_clean.index_content.unique():
    y=[ word for word in review_long_clean.content_type[review_long_clean.index_content==ind].unique() ]
    Y.append(y)
len(Y)

X=[]
for ind in review_long_clean.index_content.unique():
    term=[ word for word in review_long_clean.word[review_long_clean.index_content==ind].values ]
    X.append(' '.join(term))
len(X)

X
Y

[‘东西收到这么久忘好评品牌值得信赖东西整体个人感觉不错值得拥有’,
‘安装师傅力感谢’,
‘安装满意’,
‘收到自营商品发货速度品质保障安装效果宝贝喜欢冬天小时热水自营值得信赖值得推荐自营’,
‘几次评价产品满意加热保温时间长售后服务特别主动打电话询问送货情况帮安装非常感谢售后大姐服务下次购买’,
‘物美价廉特别划算加热速度家里不错特别’,
‘价格合理配置高物美价值’,
‘老师装装修中试希望’,
‘五分习惯送装装修中吊顶装话醉遥控器电池安装哥车里备用给装赠送电动牙刷空气净化器等待配送中时间冰箱降洗衣机降降小厨宝降烟灶套降客服上报解决电话通知等待中’,
‘安装小哥工作尽心家老房子安装家费力安装不错整栋楼热水裙子只能装小哥点赞二个一如既往好用品牌质量’,
‘前天下单买用上发现街上卖便宜店家网上他家贵老板面网上’,
‘冲着以内变频级能效制冷效果不错’,
‘购物电器问型号型号相识家电下线安装师傅挺好加热很快购物满意走用安装费花元喷头槊料花伞安装麻烦’,
‘物流早上送来中午安装安装安装师傅负责效果售后确实’,
‘满意价格产品质量货运安装调试工作人员服务态度售后服务’,
‘收费安装速度很快管道走整齐不错’,
…]
[[‘pos’],
[‘pos’],
[‘pos’],
[‘pos’],
[‘pos’],
[‘pos’],
[‘pos’],
[‘pos’],
[‘pos’],
[‘pos’],
…]

3.2 构建模型

主要用词袋模型构建文本向量，可以换成TFIDF应该会更好

#第二步：训练集、测试集划分
x_train,x_test,y_train,y_test=train_test_split(X,Y,test_size=0.2,random_state=7)

#第三步：词转向量，01矩阵
count_vec=CountVectorizer(binary=True) 
x_train=count_vec.fit_transform(x_train)
x_test=count_vec.transform(x_test)  

#第四步：构建决策树
dtc=tree.DecisionTreeClassifier(max_depth=5)
dtc.fit(x_train,y_train)
print('在训练集上的准确率：%.2f'% accuracy_score(y_train,dtc.predict(x_train)))

y_true=y_test
y_pred=dtc.predict(x_test)
print(classification_report(y_true,y_pred))
print('在测试集上的准确率：%.2f'% accuracy_score(y_true,y_pred))

DecisionTreeClassifier(max_depth=5)
在训练集上的准确率：0.72
precision recall f1-score support
neg 0.63 0.96 0.76 197
pos 0.92 0.44 0.60 196

accuracy 0.70 393
macro avg 0.78 0.70 0.68 393
weighted avg 0.77 0.70 0.68 393

在测试集上的准确率：0.70

3.3绘制决策树

#第五步：画决策树

cwd=os.getcwd()
dot_data=tree.export_graphviz(dtc
                              ,out_file=None
                              ,feature_names=count_vec.get_feature_names())
graph=graphviz.Source(dot_data)
graph.format='svg'
graph.render(cwd+'/tree',view=True)
graph