CN109637672B - 一种基于药物效用度的方剂相似度计算方法 - Google Patents

一种基于药物效用度的方剂相似度计算方法 Download PDF

Info

Publication number
CN109637672B
CN109637672B CN201811436546.6A CN201811436546A CN109637672B CN 109637672 B CN109637672 B CN 109637672B CN 201811436546 A CN201811436546 A CN 201811436546A CN 109637672 B CN109637672 B CN 109637672B
Authority
CN
China
Prior art keywords
drug
similarity
prescription
calculate
subsets
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811436546.6A
Other languages
English (en)
Other versions
CN109637672A (zh
Inventor
李建强
马泽锐
赵亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Lingxiao Smart Cloud (Beijing) Technology Co.,Ltd.
Original Assignee
Beijing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201811436546.6A priority Critical patent/CN109637672B/zh
Publication of CN109637672A publication Critical patent/CN109637672A/zh
Application granted granted Critical
Publication of CN109637672B publication Critical patent/CN109637672B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H70/00—ICT specially adapted for the handling or processing of medical references
    • G16H70/40—ICT specially adapted for the handling or processing of medical references relating to drugs, e.g. their side effects or intended usage
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/22—Matching criteria, e.g. proximity measures

Landscapes

  • Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Chemical & Material Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Medicinal Chemistry (AREA)
  • Pharmacology & Pharmacy (AREA)
  • Toxicology (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • Medical Treatment And Welfare Office Work (AREA)

Abstract

本发明公开了一种基于药物效用度的方剂相似度计算方法，包括：将方剂的药物分为相同药物子集和不同药物子集，对每味药物赋予不同权重即药物效用度，分别计算两类药物子集对方剂相似度的贡献值，同时最终得到方剂之间的相似度。本发明提出的基于药物效用度的方剂相似度计算方法，能够更加准确地计算方剂相似度，该算法应用于中医方剂处理时的方剂相似度计算。

Description

一种基于药物效用度的方剂相似度计算方法
技术领域
本发明涉及一种方剂相似度计算方法，特别涉及应用在中医领域的一种基于药物效用度的方剂相似度计算方法。
背景技术
中医药起源于中国，有着数千年的悠久历史，底蕴十分丰富为中华民族的繁衍昌盛和文明进步曾作出过巨大的贡献，对世界医学的发展和文化交流也产生了一定的影响。近年来，国家对于中医药的研究和开发越来越重视，出台了一系列中医药相关的政策法规，我国的中医药现代化正面临快速发展的时期。
方剂是根据配伍原则，总结临床经验，以若干药物配合组成的药方。在进行方剂分析时，往往需要比较方剂之间的相似性，来判断两个方剂之间功效的差异程度。
现有的方剂相似度计算方法存在以下不足：在计算相似度时只是单纯地把药物名和剂量作为方剂的特征，而没有把药物的性味归经等属性作为特征进行计算；或者是把所有药物的属性按照相同的权重进行相加计算，没有考虑到每一味药物对于方剂的重要程度是不相同的。
发明内容
本发明的目的在于提供一种用于计算方剂相似度的方法，其可以克服现有方法和装置的缺陷，计算药物子集的相似度检验方剂相似度。
为实现上述目的，本发明采用如下的技术方案：
一种基于药物效用度的方剂相似度计算方法，包括：
步骤1：对方剂数据进行处理，对方剂数据中的方剂中药物组成和特征进行标准化；
步骤2：基于方剂中的药物频数，计算药物的效用度；
步骤3：计算相同药物子集间的相似度；
步骤4：计算不同药物子集间的相似度；
步骤5：基于步骤3和步骤4两个计算，采用加权平均的方法计算方剂相似度。
步骤1：对方剂的药物组成进行处理，去除药物组成中关于药物的剂量描述，保留药物组成中药物名称的描述。对药物名称进行标准化处理，参照相关中药名称标准规范，把中药名称中出现别名、通假字或其他使用不规范的名称统一修改为标准名称。
步骤2：计算药物的效用度。效用度用来度量某个药物对于所属类别方剂的重要程度MI。计算方法如下：
式(1)中，count Mi为药物Mi在类别CX中出现的频数，P_set为所有方剂的集合，||P_set||为所有方剂的集合长度，count P∈P_set为包含药物Mi的方剂频数。
步骤3：计算相同药物子集间的相似度J。
式(2)中，MI(Mi，CX)为药物效用度系数，LS(P1，P2)为两个方剂相同药物的个数，P1、P2表示作为比较的两个方剂。
步骤4：计算不同药物子集间的相似度。
式(3)中：LD(P1)，LD(P2)为方剂内不同药物的个数，Similarity(Mi，Mj)为药物相似度，为相似药物对的平均效用度。
步骤5：计算总的方剂相似度。基于步骤3和步骤4，分别计算相同药物子集和不同药物子集的相似度，对这两个计算结果进行综合。
SF(P1，P2)＝wj*J(P1，P2)+wd*D(P1，P2) (4)
其中，wj，wd分别为相同药物子集和不同药物子集的相似度权重，且wj>wd。
与已有技术相比，本发明有益效果在于：
1、本发明使用药物效用度，对不同药物赋予不同的权重，在方剂相似度计算时更加准确。
2、本发明在计算方剂差异度时不仅考虑到了相同药物，还对不同的药物建立了相似药物对，全面衡量了两个方剂间的相似度。
附图说明
结合附图，从下面对本发明实施案例的详细描述，将更好的理解本发明，其中：
图1示出了根据本发明的计算方剂相似度的方法的流程图。
具体实施方式
本发明所述的方剂相似度的计算方法流程图如图1所示。
本实施实例中，如图1所示，基于药物效用度的方剂差异度计算方法，包括以下步骤：
步骤1：对方剂的药物组成进行处理，去除药物组成中关于药物的剂量描述，保留药物组成中药物名称的描述。对药物名称进行标准化处理，参照相关中药名称标准规范，把中药名称中出现别名、通假字或其他使用不规范的名称统一修改为标准名称。
步骤2：计算药物的效用度。效用度用来度量某个药物对于所属类别方剂的重要程度。计算方法如下：
式(1)中，Mi为药物，CX为所属方剂的类别，count Mi为药物Mi在类别CX中出现的频数，为类别CX中所有药物的频数。
式(1)中，P_set为所有方剂的集合，||P_set||为所有方剂的集合长度，count P∈P_set为P_set中包含药物Mi的方剂频数。
步骤3：计算相同药物子集间的相似度。相同药物子集是两个方剂药物交集部分，计算公式如下：
式(2)中，|P1∩P2|为P1、P2两个方剂的药物交集个数，P1∪P2为P1、P2两个方剂的药物并集个数，MI(Mi，CX)为药物效用度系数，LS(P1，P2)为两个方剂相同药物的个数。
步骤4：计算不同药物子集间的相似度，不同药物子集是两个方剂药物差集部分，计算公式如下：
式(3)中，LD(P1)，LD(P2)为方剂内不同药物的个数，Similarity(Mi，Mj)为药物相似度，为相似药物对的平均效用度。
式(3)中Similarity(Mi，Mj)为药物相似度，可以使用已有的药物相似度计算方法，这里列出一些常见的计算方法，如：
使用列联表计算药物的相似度,将药物的功效、主治等属性，构成one-hot向量，计算公式为：
式(5)中a为药物Mi，Mj的相同属性数量，b为药物Mi中与Mj不同的属性数量，c为药物Mj中与Mi不同的属性数量，d为药物Mi，Mj都不包含的属性数量，n＝a+b+c+d。
式(3)中选择相似药物对的步骤是，对于方剂中的每个药物，基于Similarity(Mi，Mj)，在另一个方剂中找到与其相似度最高的一个药物，组成一组相似药物对。与此类似，分别找到方剂中其他药物所对应的相似药物对，最终能够组成Min{LD(P1)，LD(P2)组相似药物对。
式(3)中为相似药物对的平均效用度，由于相似药物对有两味药物，所以使用它们药物效用度的平均值作为药物的权重系数。
式(3)中使用做归一化处理，使式(2)能够和式(3)做相加计算。
步骤5：计算总的方剂相似度。基于步骤3和步骤4，分别计算了相同药物子集和不同药物子集的相似度，对这两个计算结果进行加权相加。
SF(P1，P2)＝wj*J(P1，P2)+wd*D(P1，P2) (4)
式(4)中，wj，wd分别为相同药物子集和不同药物子集的相似度权重，且wj>wd。J(P1，P2)为相同药物子集间的相似度，D(P1，P2)为不同药物子集间的相似度。
本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本公开的其它实施方案。本申请旨在涵盖本公开的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本公开的一般性原理并包括本公开未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本公开的真正范围和精神由下面的权利要求指出。
应当理解的是，本公开并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本公开的范围仅由所附的权利要求来限制。

Claims (1)

1.一种基于药物效用度的方剂相似度计算方法，其特征在于：该方法包括：
步骤1：对方剂数据进行处理，对方剂数据中的方剂中药物组成和特征进行标准化；
步骤2：基于方剂中的药物频数，计算药物的效用度；
步骤3：计算相同药物子集间的相似度；
步骤4：计算不同药物子集间的相似度；
步骤5：基于步骤3和步骤4两个计算，采用加权平均的方法计算方剂相似度；
步骤1：对方剂的药物组成进行处理，去除药物组成中关于药物的剂量描述，保留药物组成中药物名称的描述；对药物名称进行标准化处理，参照相关中药名称标准规范，把中药名称中出现别名、通假字或其他使用不规范的名称统一修改为标准名称；
步骤2：计算药物的效用度；效用度用来度量某个药物对于所属类别方剂的重要程度MI；计算方法如下：
式(1)中，count Mi为药物Mi在类别CX中出现的频数，P_set为所有方剂的集合，||P_set||为所有方剂的集合长度，count P∈P_set为包含药物Mi的方剂频数；
步骤3：计算相同药物子集间的相似度J；
式(2)中，MI(Mi,CX)为药物效用度系数，LS(P1,P2)为两个方剂相同药物的个数，P1、P2表示作为比较的两个方剂；
步骤4：计算不同药物子集间的相似度；
式(3)中：LD(P1),LD(P2)为方剂内不同药物的个数，Similarity(Mi,Mj)为药物相似度，为相似药物对的平均效用度；
步骤5：计算总的方剂相似度；基于步骤3和步骤4，分别计算相同药物子集和不同药物子集的相似度，对这两个计算结果进行综合；
SF(P1,P2)＝wj*J(P1,P2)+wd*D(P1,P2) (4)
其中，wj,wd分别为相同药物子集和不同药物子集的相似度权重，且wj>wd。
CN201811436546.6A 2018-11-28 2018-11-28 一种基于药物效用度的方剂相似度计算方法 Active CN109637672B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811436546.6A CN109637672B (zh) 2018-11-28 2018-11-28 一种基于药物效用度的方剂相似度计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811436546.6A CN109637672B (zh) 2018-11-28 2018-11-28 一种基于药物效用度的方剂相似度计算方法

Publications (2)

Publication Number Publication Date
CN109637672A CN109637672A (zh) 2019-04-16
CN109637672B true CN109637672B (zh) 2023-11-07

Family

ID=66070032

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811436546.6A Active CN109637672B (zh) 2018-11-28 2018-11-28 一种基于药物效用度的方剂相似度计算方法

Country Status (1)

Country Link
CN (1) CN109637672B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN204349099U (zh) * 2015-02-02 2015-05-20 浙江创佳电子有限公司 前锁式贴片式fpc连接器
GB201507501D0 (en) * 2015-04-30 2015-06-17 Fujitsu Ltd A similarity-computation apparatus, a side effect determining apparatus and a system for calculating similarities between drugs and using the similarities
CN105550526A (zh) * 2016-01-04 2016-05-04 北京科技大学 一种基于中药属性的中医方剂功效定量方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101436534B1 (ko) * 2011-06-01 2014-09-01 가부시키가이샤 쓰무라 집합데이타의 유사성 평가방법, 유사성 평가프로그램 및 유사성 평가장치
JP5997120B2 (ja) * 2013-09-26 2016-09-28 富士フイルム株式会社 薬剤照合装置、薬剤照合システム、及び薬剤照合方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN204349099U (zh) * 2015-02-02 2015-05-20 浙江创佳电子有限公司 前锁式贴片式fpc连接器
GB201507501D0 (en) * 2015-04-30 2015-06-17 Fujitsu Ltd A similarity-computation apparatus, a side effect determining apparatus and a system for calculating similarities between drugs and using the similarities
CN105550526A (zh) * 2016-01-04 2016-05-04 北京科技大学 一种基于中药属性的中医方剂功效定量方法及系统

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
中医药方剂相似度模型;操牡丹等;《计算机工程》;20090820(第16期);第275-277页 *
利用效用度挖掘核心药物及配伍规律;周伟;王峰;王崇骏;谢俊元;;计算机科学与探索(第11期);第994-1000页 *
基于SimRank的中药"效-效"相似关系挖掘;田玲 等;《计算机工程》;20080620(第12期);第242-244页 *
基于中药功效的聚类分析;何前锋等;《中国中医药信息杂志》;20040615(第06期);全文 *
基于文本分类技术计算中医方剂相似度;顾铮;《微计算机信息》;20100425(第12期);第199-201页 *
相似度算法和凝聚法对中药功效及性味的分类研究;曾萍等;《福建电脑》;20100525(第05期);全文 *

Also Published As

Publication number Publication date
CN109637672A (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
CN105005589B (zh) 一种文本分类的方法和装置
Silva et al. Adverse psychosocial working conditions and poor quality of life among financial service employees in Brazil
CN105468742A (zh) 恶意订单识别方法及装置
CN102708168A (zh) 一种教学资源搜索结果排序的系统及方法
JP2016181036A (ja) 作業配分支援装置、作業配分支援プログラムおよび作業配分支援方法
CN103164416A (zh) 一种用户关系的识别方法及设备
Crespo et al. The effect on citation inequality of differences in citation practices at the web of science subject category level
CN115311478A (zh) 一种基于图深度聚类的联邦影像分类方法、存储介质
Su Balancing multiple baseline characteristics in randomized clinical trials
CN108628967A (zh) 一种基于学习生成网络相似度的网络学习群组划分方法
CN109637672B (zh) 一种基于药物效用度的方剂相似度计算方法
WO2023178979A1 (zh) 问题标注方法、装置、电子设备及存储介质
CN110084464A (zh) 一种基于物联网的智慧城市火灾危险性实时评估方法
Heng et al. Estimation and Hypothesis Testing of Strain‐Specific Vaccine Efficacy With Missing Strain Types With Application to a COVID‐19 Vaccine Trial
CN112348356A (zh) 一种作品质量确定方法、装置、计算机设备及可读存储介质
CN104376064B (zh) 一种挖掘用户年龄样本的方法和装置
CN113486219A (zh) 一种处理数据的方法、装置、电子设备以及存储介质
CN109544296B (zh) 一种用于交易平台的商品挂单调整系统
CN106485404A (zh) 一种基于triz理论在模具制造中的应用评价方法
Ching ISI perceptions and hard facts: An empirical study from Taiwan
Huang et al. A study on the allocation of special education resources in China
CN106529212B (zh) 基于序列依赖频率矩阵的生物序列进化信息提取方法
CN104090950B (zh) 融入簇存在强度的数据流聚类方法
Mosteller et al. Combination of results of stated precision: I. The optimistic case
CN109215793A (zh) 一种金融行业用户体验健康度评估系统及其使用方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20241212

Address after: 518000 1002, Building A, Zhiyun Industrial Park, No. 13, Huaxing Road, Henglang Community, Longhua District, Shenzhen, Guangdong Province

Patentee after: Shenzhen Wanzhida Technology Co.,Ltd.

Country or region after: China

Address before: 100124 No. 100 Chaoyang District Ping Tian Park, Beijing

Patentee before: Beijing University of Technology

Country or region before: China

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20250124

Address after: 100000, 08B03, 7th Floor, No. 2, East Third Ring North Road, Chaoyang District, Beijing

Patentee after: Lingxiao Smart Cloud (Beijing) Technology Co.,Ltd.

Country or region after: China

Address before: 518000 1002, Building A, Zhiyun Industrial Park, No. 13, Huaxing Road, Henglang Community, Longhua District, Shenzhen, Guangdong Province

Patentee before: Shenzhen Wanzhida Technology Co.,Ltd.

Country or region before: China