近年来,关联规则、聚类分析、贝叶斯网络、神经网络、复杂网络分析等数据挖掘技术为名老中医学术思想传承工作提供了新的方法学基础[。但对许多中医背景的传承人来说,数据挖掘作为多学科交叉的新产物难以熟练运用[。

针对这一现状,为探寻更加实用、高效的数据挖掘方法,本研究组基于复杂网络分析和关联规则分析创建“组方配伍网络分析法”,对天津中医药大学张伯礼教授治疗胸痹的处方进行可视化网络重建及用药规律分析。总结出的重点药物、治法及学术特点,与相关文献[结论基本一致,得到了张教授本人认可。

1 总体思路

相关研究

为体现中医思维,在挖掘组方思路时,借助复杂网络分析软件记录药物的录入顺序,模拟传统师承的“抄方学习模式”,以体现专家遣方用药思路的连续性、系统性;在挖掘配伍思路时,借助关联规则分析法,以获得不同药物间具体的配伍概率。最终,整合两类数据构建“组方配伍网络”,既能从纵向体现组方思路,又能从横向展示配伍关系。

2 研究方法

2.1 数据来源

张伯礼教授的处方数据调取自天津中医药大学附属保康医院门诊系统(2016年),具体内容包括:数据入库顺序、就诊时间、处方编号、患者编号、患者姓名、录入药物、录入顺序、药物剂量、用药周期、ICD诊断名、诊断备注、药物金额及医师姓名。

2.2 处方纳入标准

纳入标准定为:1)ICD诊断名为“胸痹”。2)无第二诊断。3)无其他诊断备注。

按上述标准,共369个胸痹处方被纳入。

2.3 处方排除标准

为保证研究结论能够体现组方思路的连续性,排除标准如下:1)处方药物录入起始序号不为1。2)处方药物录入顺序不连续(缺项)。

最终,13个不合要求的处方被剔除,保留合格处方356个。

2.4 数据编辑与处理

2.4.1 建立处方库

使用Excel 2016(以下简称Excel)从原始资料中提取处方编号、录入药物、录入顺序进行建库。

2.4.2 数据处理

网络的构建及分析采用Gephi 0.9.2[(以下简称Gephi)完成,由介数中心性、K-壳分解法

1)组方数据挖掘:组方数据挖掘借助复杂网络分析软件完成,并以组方网络图的形式呈现。

将处方库中的数据转制为Gephi可调用的CSV格式,所有处方的“首尾”加入“开始”和“结束”标记,即“A-B-C-D-E”转为“开始-A-B-C-D-E-结束”,以避免两个不同处方“首尾”的药物节点生成连接。

数据导入Gephi生成组方网络后:①使用内置工具统计介数中心性;②使用内置模块化工具完成分类;③以不同颜色区分不同模块;④以介数中心性降序决定节点大小和节点名称汉字大小;⑤以K-壳分解法提取核心群。

2)配伍数据挖掘:配伍数据挖掘借助R的关联规则工具箱和复杂网络分析软件共同完成,以配伍网络图的形式呈现。

为收集所有的“两两药物组合”情况,本次研究将支持度和置信度的值均设为0,获得32 400个组合。剔除17 760个频数为0的组合及7 230个镜像关联情况,最终获得7 230个共现频次至少为1的“两两药物组合”,用于构建配伍网络。

将筛选完毕的关联数据转制为Gephi可调用的CSV格式,导入Gephi生成配伍网络后,网络参数调节同“组方网络”。

3)建立组方配伍网络:组方配伍网络的数据可以由组方、配伍网络的CSV文件直接拼接获得,导入Gephi生成组方配伍网络后,网络参数调节同“配伍网络”。

2.5 比较与验证

2.5.1 模型比较

参考频次统计结果,从网络概况、节点、社团模块、链路4个方面,将组方网络、配伍网络、组方配伍网络及各自核心群进行对比,评价组方配伍网络及其核心群能否如实反映处方信息,体现组方思路和配伍关系。

2.5.2 文献验证

以“张伯礼”、“胸痹”、“学术思想”为关键字,分别在知网、万方、维普三大检索平台进行期刊文献检索,经去重和剔除无关文献后,保留继承总结张伯礼教授学术思想的相关文献7篇。结合文献内容,评价组方配伍网络及其核心群的链路、社团模块、节点能否反映张伯礼教授治疗胸痹的学术思想、治则治法、重点药物。

2.5.3 专家确认

经相关文献验证后,将组方配伍网络以图片输出,并对各社团模块、节点、链路的内容及其意义向张伯礼教授做说明,确认分析结果是否能得到专家本人认可。

3 研究结果

3.1 频数统计

借助R内置的关联规则分析包,统计不同药物的使用频次、药物组合的共现频次。截取部分频数统计结果制成

表 1 不同药物使用频次前20名

表 2 药物组合使用频次前20名

3.2 网络模型对比

3.2.1 概况对比

本次研究共构建了6个网络,分别是组方网络、配伍网络、组方配伍网络及其各自的核心群。具体参数见

表 3 网络概况对比

组方网络节点连接相对稀疏,核心群最小,模块数最多;配伍网络和组方配伍网络节点连接较紧密,核心群较大,模块数方面组方配伍网络比配伍网络多。

3.2.2 节点对比

介数中心性能够反映所评价节点连接其他节点的能力,即所评价药物与其他药物的配伍关系是否密切。各网络介数中心性前10的药物如下:1)组方网络及其核心:柴葛根、夏枯草、杜仲、知母、黄芩、茯苓、生龙齿、玉竹、桑寄生、干姜。2)配伍网络和组方配伍网络为:丹参、郁金、茯苓、黄连、麦门冬、生龙齿、生牡蛎、首乌藤、半夏、延胡索。3)配伍网络核心和组方配伍网络核心为:郁金、丹参、茯苓、延胡索、柴葛根、降香、黄连、麦门冬、生龙齿、首乌藤。

参考频次统计情况,组方网络(及其核心群)重要节点与频次统计结论的一致性较差,而配伍网络、组方配伍网络(及其核心群)较好。

3.2.3 模块对比

社区探测算法可以依据节点连接的疏密情况,将整个网络划分为若干网络模块。基于方剂学理论,中医处方至少应包括“君臣佐使”4个基本模块。结合

从连接情况来看:1)组方网络的各模块内部连接较紧密,而模块之间连接较稀疏。2)配伍网络中只有1个模块内部连接紧密,其余模块及模块间连接较为稀疏。3)组方配伍网络各模块内部、模块间的连接都较为紧密。

3.2.4 链路对比

药物节点相连构成的链路反映了不同药物间存在的配伍关系。从链路的节点构成来看:1)组方网络及其核心群的链路可体现某一用药目的,如丹参、郁金、元胡、降香体现了活血化瘀、行气止痛的用药目的,但不同链路间是否存在联系难以确定。2)配伍网络及其核心群的链路反映使用频次较高的药物组合,但各链路内部药物重复率高。3)组方配伍网络及其核心群的链路节点重复较少,不同链路间部分节点相互连接,既可以反映不同药物组合的用药目的,又可以观察到不同药物组合间的协同关系。

3.3 比较与验证

通过不同网络模型之间的比较,经文献验证和与专家本人确认,研究认为组方配伍网络能够如实反映处方信息,张伯礼教授治疗胸痹的学术思想、治则治法、常用药物可从组方配伍网中的链路、社团模块和节点中得到反映。具体内容详见

图 1 2016年度张伯礼教授治疗胸痹的组方配伍网络图

4 讨论

4.1 结果分析

在对相关文献[研究中发现,张伯礼院士治疗胸痹的特点可概括为:1)理论上以“痰瘀学说”为基础,主张“痰瘀并治”,重视痰瘀在疾病发展各阶段的变化。2)治法上依患者具体情况,采取活血化瘀(丹参、郁金、延胡索、降香、三七粉、五灵脂)、益气养阴法(生地、沙参、麦门冬、党参、黄精)、健脾祛湿(党参、茯苓、白术)、通经活络(鸡血藤、桑枝)等方法。3)用药上根据不同病情条件下的痰瘀特点,选取相应的对(队)药予以治疗(湿邪重用萆薢、苍术、蚕砂,湿邪化热用茵陈蒿、知母、黄连、大黄,痰火盛用黄芩、浙贝母,痰湿蕴肺用细辛、干姜、半夏)。4)兼顾气、血、阴、阳及诸脏功能的平衡(女贞子、旱莲草、杜仲、牛膝、当归、何首乌),注重调养心神(酸枣仁、柏子仁、夜交藤、生龙齿、生牡蛎)。

参照张伯礼教授治疗胸痹的特点[,结合组方配伍网络的节点、链路、社团模块情况,本研究组认为:1)社团模块及其内部介数中心性高的节点,能够反映特定的用药目的,如:模块1体现了健脾祛湿法(党参、茯苓、白术)。模块2以牛蒡子、射干、橘红、浙贝母等化痰药为主,体现了治痰的用药变化。模块3体现了活血化瘀法(丹参、郁金、延胡索、降香、三七粉、五灵脂)。模块4以柴葛根、杜仲、夏枯草、鸡血藤等药物为主,包含药物最多,是兼顾气、血、阴、阳及诸脏功能的临证加减用药模块。模块5以首乌藤、生龙齿、柏子仁、酸枣仁等为主,体现了注重调养心神。2)不同模块间重点药物的连接,能够反映专家治疗某一疾病的学术特点,如:本研究中,模块1与模块3的重点药物及其构成的链路,体现了张伯礼教授治疗胸痹的“痰瘀学说”。

4.2 组方配伍网络的优势

基于药物录入顺序构建的组方网络及其核心,有利于体现组方思路,但缺乏关联数据支撑;基于关联分析构建的配伍网络及其核心,有利于挖掘核心配伍组合,但会使处方信息碎片化,难以揭示组方思路。

作为组方网络和配伍网络的有机结合体,组方配伍网络既可以直观展示组方思路,又可以在关联数据的基础上明确具体的药物配伍关系。以链路、社团模块反映理论和治法,基于介数中心性筛选重点药物,符合“理、法、方、药”的中医临床思维。

4.3 问题与展望

名老中医经验的存在形式既有显性的,也有隐性的。单纯照搬某种方法的分析结果,很难得到名老中医本人的认可。因此,需要在中医理论指导下综合运用多种方法,对名老中医经验方进行系统、全面地挖掘和总结

作为一种多分析方法综合应用的尝试,“组方配伍网络”在参考传统师承“抄方”的基础上,综合运用关联分析、复杂网络分析两种方法,为总结名老中医治疗某一病证(疾病)的“理、法、方、药”提供了一套较为完备的解决方案。

但完备并不等于完美,它也存在一些问题需要深入探讨,例如:如何确定社团模块划分的最优解析度。析其原因,复杂网络分析法虽已应用于中医药数据挖掘领域,但具体操作中相关参数设置尚无统一标准可以借鉴。因此,需要在大量数据实验的基础上,由相关研究人员共同寻找最优参数,制定统一标准。

如何用计算机应用于中医方剂,方剂“组方配伍网络”分析方法研究——以名老中医门诊处方数据为例...相关推荐

  1. 白宁超计算机科学院,基于主动学习的传统中医症状本体构建方法研究综述.PDF...

    基于主动学习的传统中医症状本体构建方法研究综述.PDF 计算机技术应用 • the Application of Computer Technology 基于主动学习的传统中医症状本体构建方法研究综述 ...

  2. 计算机应用 网络管理开发,基于XML的iBAC网络管理系统的研究与开发-计算机应用技术专业论文.docx...

    ⅢY ⅢY iii■l 洲8 mmj■I ㈣0 Ⅲ4 6 ㈣2 学位论文数据集 中图分类号TP311.1学科分类号520.3040 论文编号10010200705 12密级 学位授予单位代码10010 ...

  3. python情绪分析的意义_如何用Python和R对故事情节做情绪分析?

    想知道一部没看过的影视剧能否符合自己口味,却又怕被剧透?没关系,我们可以用情绪分析来了解故事情节是否足够跌宕起伏.本文一步步教你如何用Python和R轻松愉快完成文本情绪分析.一起来试试吧. 烦恼 追 ...

  4. 笔记本改装家庭文件服务器,如何用一个废旧的笔记本打造一个家庭网络服务器?...

    原标题:如何用一个废旧的笔记本打造一个家庭网络服务器? 猿 份 本来好好的笔记本的,可偏偏屏幕说烂就烂,那就干脆不要了,改造改造吧! 正文共:1098 字 预计阅读时间:3 分钟 他本来长这样的,如图 ...

  5. 上网本改文件服务器,如何用一个废旧的笔记本打造一个家庭网络服务器?

    原标题:如何用一个废旧的笔记本打造一个家庭网络服务器? 本来好好的笔记本的,可偏偏屏幕说烂就烂,那就干脆不要了,改造改造吧! 他本来长这样的,如图: 什么,你看见了水印-.老哥,别在意这些细节. 简单 ...

  6. 非计算机专业特别是人文社科,中职非计算机专业《计算机应用基础》教学必须重视学情分析...

    摘要:中职非计算机专业的<计算机应用基础>课程的教学过于一统化,学情分析是系统教学设计不可或缺的重要组成部分,要提高该课程的教学有效性,优化教学效果,教师必须在课前.课中.课后关注学情.研 ...

  7. 准确率99.9%!如何用深度学习最快找出放倒的那张X光胸片(代码+数据)

    医学图像数据的质量一直是个老大难题.难以清理的数据制约着许多深度学习的应用. 而实际上,深度学习本身就是清洗医疗数据的好帮手. 今天,我们就来讲一个案例,展示如何用深度学习迅速清洗一个杂乱的医疗图像数 ...

  8. 计算机应用基础怎么教作业,北京语言大学网络教育2017春计算机应用基础作业(1)及答案...

    以下为湖南自考网指导老师为大家整理的,北京语言大学网络教育2017春计算机应用基础作业(1)及答案,请准备报考北京语言大学网络教育的考生,考前多做作业,熟悉考试重点. 北京语言大学网络教育2017春 ...

  9. 如何用Apache POI操作Excel文件-----如何在已有的Excel文件中插入一行新的数据?

    在POI的第一节入门中,我们提供了两个简单的例子,一个是如何用Apache POI新建一个工作薄,另外一个例子是,如果用Apache POI新建一个工作表.那么在这个章节里面,我将会给大家演示一下,如 ...

  10. 武汉理工大学计算机应用基础作业,2017年武汉理工大学网络教育计算机应用基础作业2...

    2017年武汉理工大学网络教育计算机应用基础作业2 (6页) 本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦! 11.90 积分 一.单项选择题(共20道小题,共 ...

最新文章

  1. QT 中使用 c++ 的指针
  2. 畅玩mt3单机游戏服务器维护,【梦幻西游】MT3仿端手工游戏服务端源码[教程+授权物品后台]...
  3. Java多线程相关的常用接口
  4. 利用均差的牛顿插值法(Newton)
  5. redis smembersmap_redis 基本指令以及数据类型
  6. 二分法求数组最大最小_JavaScript之如何求数组的最大值和最小值
  7. 【报告分享】2020巨量引擎营销通案.pdf(附下载链接)
  8. 2011版计算机网络原理,2011年7月自考04741计算机网络原理真题及答案
  9. svn merger的时候 报远程主机强迫_SVN与Git比较(二)版本库与工作区
  10. java实验三多态性_Java实验3 类的多态性和接口
  11. Kafka配置4--Windows下配置Kafka的SSL证书
  12. matlab半波整流怎么做,基于Matlab的单相半波可控整流电路的设计与仿真.doc
  13. 阿里云服务器地域及可用区的选择问题
  14. python字典get用法_详细解析python字典get()实例教程
  15. 2081.09.22 Kuma(非旋treap)
  16. ios-音乐播放器的简单实现
  17. Webpack经典入门
  18. day25组合总和3_电话号码的字母总和
  19. 数字证书颁发及认证原理
  20. 用windows系统下的DOS命令将腾讯视频客户端下载的qlv文件转换成MP4格式(图文详解)

热门文章

  1. vue 判断设备是手机端还是pc端
  2. 简单有效的记录日常收支
  3. 计算机考研408(计算机网络(持续更新))
  4. 【javaEE】——多线程进阶(锁策略:面试相关考点)04
  5. oracle exadata x7发布,Oracle Exadata X8 Hardware DataSheet (配置说明)
  6. LDA模型原理+代码+实操
  7. spotfire Document Property
  8. Spotfire 连接mysql数据库
  9. 2021年全年详细工作日、周末、节假日数据表sql_excel
  10. ldap 统一认证 java_基于LDAP的统一身份认证系统的设计与实现