基于数据挖掘技术探讨奔豚气病用药规律*

薛皓方，刘 锐，黄文珊，古展鑫

1 广西中医药大学研究生学院，广西 南宁 530001； 2 广西中医药大学附属瑞康医院，广西 南宁 530001

奔豚，因其发作时胸腹如有小豚奔闯而得名，是一种发作时觉气自下腹上冲心胸，直达咽喉的疾病，可伴有寒热往来、吐脓、腹部绞痛、胸闷气急、头昏目眩、心悸，烦躁不安等症状。从证候表现看，奔豚相当于西医的胃肠神经官能症（肠道积气和蠕动亢进或痉挛状态）及冠心病、心血管神经症［1］等。奔豚病名最早见于《灵枢·邪气藏府病形》篇［2］，《难经·五十六难》谓：“肾之积名曰奔豚，发于少腹，上至心下，若豚状，或上下无时”［3］，将其归属于五积之肾积。本研究运用数据挖掘方法，分析《中医方剂大词典》中治疗奔豚气病处方潜在的共性，探索本病病因病机和治疗原则。以期为临床治疗本病提供可供参考的临床药对、药组及新处方。

1 资料与方法

1.1 资料来源与筛选以《中医方剂大辞典》［4］为方剂来源，以“奔豚”“奔豚气”为主要检索词进行检索，筛选出主治奔豚气病的方剂。

1.2 数据处理

1.2.1 数据规范化处理 由于中医古籍中药物别称、错称、生僻名使用频繁，根据《中药学》［5］和2015年版《中华人民共和国药典》［6］一部（简称《中国药典》）规范中药名称，如“小香”统一为“小茴香”。

1.2.2 量化数据 将数据库中的中药字段采用二分类量化处理，如药物出现则记为“1”，没有出现则记为“0”。参照《中药学》［5］对药物的四气五味性味归经进行统一规范，如药性为温者在温项下赋值为“1”，余各项赋值为“0”。以此类推。

1.3 统计学方法采用Microsoft Excel 2016建立数据库并进行频数统计；运用SPSS Modeler 18.0软件的Apriori算法进行关联规则分析并利用Web 结点构建药物关联规则网状图；应用SPSS23.0 对频数较高的药物进行因子分析提取载荷较高的公因子。

2 结果

共筛选出符合条件的外用方剂76 首，包含135 味中药。方源来自《圣济总录》《圣惠方》《普济方》《医方类聚》《证治宝鉴》《朱氏集验方》等39部中医著作及验方汇编。

2.1 奔豚气病方剂药物频数分析在治疗奔豚气病的76首方剂中，使用频次≥10次的药物共计22味，频数排名前五的药物是肉桂、小茴香、木香、附子、沉香。见表1。

表1 奔豚气病方剂药物频数分析 次

2.2 高频药物性味频数分析根据《中药学》［5］和《中国药典》2015 年版分类标准，对高频药物（频次≥10）进行性味、归经统计分析。如果一味中药有多个性味，则全部统计在内。76 首方剂中使用频次≥10 次的药物共计22 味，累计出现药物四性46 次，五味37 次，归经69 次。可以看出药物四性以凉（26 次，56.5%）为最；五味以辛（16 次，43.2%）为最，苦（11 次，29.7%）；归经以归脾（16次，23.2%）、肝（13 次，18.8%）和胃（12 次，17.4%）为主。见图1—2。

图1 高频药物四气五味频数统计

图2 高频药物归经和功效频数统计

2.3 高频药物的关联规则分析[7]运用SPSS Modeler 18.0 软件Apriori 算法对频数≥10 次位中药进行关联规则分析，设置支持度≥10%，置信度为100%，最大前项数5，提升度＞1，得到13项强关联规则。见表2。

表2 高频药物关联规则分析（置信度为100%）

2.4 高频药物的因子分析运用SPPS 23.0 对频数≥10的22味高频药物进行因子分析，共得出6个共因子。得出KMO值为0.712，Bartlett球形度检验P=0.000，表明数据相关度较强，可以进行因子分析。采用主成分分析提取初始值特征值＞1的6个公因子，见表3。累计方差贡献率63.516%，表示涵盖了大部分信息，通过具有Kaiser 标准化的正交旋转法得到旋转空间中的成分图，见表3。

表3 高频药物因子分析

2.5 关联规则图将22 味高频中药按照彼此之间的链接进行分析，对链接≥45 的进行描绘，颜色越深代表关联性越强，如图3所示。

图3 常用药物网络图

3 讨论

奔豚气病因病机复杂，历代医家各有不同的认识，但总体可归纳为以下两个方面：肝郁化热，冲气上逆；阳虚致寒气或水饮上冲心胸［3，8-9］。本研究运用关联规则和因子分析对历代医家治疗奔豚气病的处方进行深入挖掘，分析处方背后潜在的共性，总结用药规律。

药物频数分析发现，排名前三的为肉桂、小茴香和木香；药物功效分类，历代医家治疗奔豚气病以理气药（木香、沉香、李根皮、川楝子、半夏、青皮、陈皮）和温里药（肉桂、小茴香、附子、吴茱萸）并重，其次为补虚药（甘草、人参、当归、白芍）；用药以药性凉、味辛、苦，归脾、胃、肝经为主。由此得知，历代医家认识到肝郁火逆、阳虚致寒气或水饮上逆均可发为奔豚气。中医认为，脾胃为气机升降之枢纽，脾主生，胃主降，如中焦气机不畅当以理气为先；辛能散能行，苦可降泄，有辛开苦降、调理气机升降之意；阳气虚弱，阳不制阴，寒气或冲气上逆可发为奔豚气，正如张志聪《灵枢集注》云：“肾为生气之源，正气虚寒，则为沉厥；虚气反逆，故为奔豚。”历代医家多选用温里药和补虚药治疗该药，体现了“虚则补之”的治疗理念。

关联规则作为最成熟且应用广泛的数据挖掘方法，可从给定事项中发现两个或多个变量之间满足一定支持度和置信度的某种规律性，发现各事物之间较强的关联组合。在关联规则中，支持度简单的字面理解就是支持的程度，表示前项与后项在一个数据集中同时出现的频率。本文对使用频次较高的药物进行关联规则挖掘分析，旨在发现临床常配伍使用的药对或药团。

本研究关联规则分析共发现13 个强关联的药组，其中支持度最高的为李根皮→半夏-当归（支持度13.16%）和半夏→李根皮-当归（支持度13.16%）。在关联规则分析中，提升度可反应关联的强弱，从关联强度高的药组中可发现：1）半夏、李根皮、当归、川芎、甘草五味药物联系密切，同时这五味药物也是治疗肝郁化热，冲气上逆常用方剂奔豚汤的主要药物组成。其中李根皮又名李根白皮、甘李根白皮，具有清热，下气，解毒之功效，历代本草著作均将其作为治疗奔豚气的专药，可清热平冲，平复奔豚气。《长沙药解》认为李根皮味涩，性寒，入足厥阴肝经，可下肝气之奔冲，清风木之郁热。李根皮作为《金匮要略》奔豚汤中的主药，其泄奔豚气效果极佳。然如今李根皮货品短缺，现常用具有敛肝降逆之桑白皮、代赭石，疏肝之川楝子、枳实、香橼，清热之椿白皮、栀子等代替［10-11］。现代研究表明，奔豚汤治疗神经官能症、焦虑症惊恐发作、抑郁症和神经性头痛等具有较好的疗效，其作用可能是通过下调海马区GAT-1和GAT-3 蛋白表达，进而调节中枢神经系统氨基酸类神经递质的含量［12］。2）小茴香（散寒止痛，理气和胃）、莪术（破血行气，消积止痛）、沉香（行气止痛，温中降逆）、木香（行气止痛）、川楝子（行气止痛，疏肝泄热）五味药物相互之间联系密切，且均具有疏肝泄热、行气止痛的功效，提示在临床治疗肝郁化火型奔豚气时，应注意以上诸药的配伍使用。王松龄教授在奔豚汤中加入川楝子、沉香等药治疗肝郁化火型奔豚气取得良好的临床疗效［13］。

因子分析是一种可在众多变量中提取隐藏但具有代表性因子的统计技术。本研究对高频次变量（药物）进行因子分析，得出的每一个公因子代表一种治法。治法是联系辨证理论和遣药组方的纽带，对每一个公因子分析来推测本病的病因病机。通过因子分析共得出7 个公因子，整体反应了疏肝泄热、行气止痛、温阳散寒和温阳利水的治疗方法。公因子1 由当归、川芎、李根皮、白芍、半夏、甘草组成，为奔豚汤的主要药物组成。结合关联规则分析可发现以上诸味药物之间联系密切，常配伍使用，说明历代医家治疗奔豚气病多以奔豚汤为底方［14］。方中以芍药、甘草奠安中气为君，以半夏开其结气，当归、川芎入血和心气，李根皮性大寒，以折其冲逆之气，此乃将欲降之，以先升之之理。耿琦等［9］通过对现代半个多世纪以来228 例有关奔豚气病的回顾性分析发现，将近半数患者与情志性发病因素有关，且多选用奔豚汤加减治疗。因子2 由青皮（疏肝破气），沉香（行气止痛、温中降逆），木香（行气止痛）、莪术（破血行气）、小茴香（散寒止痛、理气和胃）、槟榔（行气利水）组成，结合因子3 中的全蝎（熄风止痉、攻毒散结、通络止痛）、延胡索（活血、行气、止痛）、川楝子（行气止痛、疏肝泄热）共同体现了疏肝泄热、行气止痛的治疗方法。《诸病源候论·奔豚气候》云：“夫奔豚气者，肾之积气，起于惊、恐、忧思所生。”［15］中医认为惊则伤心，恐则伤肾，忧思伤肝，肝气郁结，血行不畅，或肝郁化火，随冲气上逆可发为奔豚气。关于奔豚气病的治疗，张仲景在《金匮要略方论·奔豚气病脉证并治第八》中曰：“奔豚气上冲胸腹痛，往来寒热，奔豚汤主之。”又曰：“发汗后，烧针令其汗，针处被寒，核起而赤者，必发奔豚。气从小腹上至心，灸其核上各一壮，与桂枝加桂汤主之。”奔豚气病有寒热之分，热者治之以奔豚汤，寒者则用桂枝加桂汤［16］。因子4 包含的肉桂、人参、附子、茯苓与因子5 中的吴茱萸均为温经散寒、温阳利水之常用药物。提示在治疗因外邪致阳虚寒逆引起气冲或误汗阳气受损，水饮上冲，所致的奔豚气病时，可使用肉桂、附子、人参、茯苓、吴茱萸的配伍使用。这与心阳虚衰，寒水上冲，欲作奔豚［17-18］的病因病机相吻合。

综上所述，本研究通过数据挖掘对《中医方剂大辞典》中收录的历代医家治疗奔豚气病处方进行分析，发现：奔豚气病的病因病机主要为肾脏寒气上冲、肝脏气火上逆、肝郁化火，阳虚寒气或水气上逆贯穿奔豚气病发生和发展的全过程，当以疏肝泄热、行气止痛和温经散寒、温阳化饮为主要治法，共得出3 首临床备选处方。由于本病临床不常见，现代文献关于本病的报道较少，有待进一步实践验证。