CN105868296B - 一种基于快速剪枝策略的高效用序列模式的用药ddd值数据分析方法 - Google Patents

一种基于快速剪枝策略的高效用序列模式的用药ddd值数据分析方法 Download PDF

Info

Publication number
CN105868296B
CN105868296B CN201610172655.6A CN201610172655A CN105868296B CN 105868296 B CN105868296 B CN 105868296B CN 201610172655 A CN201610172655 A CN 201610172655A CN 105868296 B CN105868296 B CN 105868296B
Authority
CN
China
Prior art keywords
sequence
value
utility
item
swu
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
CN201610172655.6A
Other languages
English (en)
Other versions
CN105868296A (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.)
Yinjiang Technology Co.,Ltd.
Original Assignee
Enjoyor Co Ltd
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 Enjoyor Co Ltd filed Critical Enjoyor Co Ltd
Priority to CN201610172655.6A priority Critical patent/CN105868296B/zh
Publication of CN105868296A publication Critical patent/CN105868296A/zh
Application granted granted Critical
Publication of CN105868296B publication Critical patent/CN105868296B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24Querying
    • G06F16/245Query processing
    • G06F16/2458Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2465Query processing support for facilitating data mining operations in structured databases
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/22Indexing; Data structures therefor; Storage structures
    • G06F16/2282Tablespace storage structures; Management thereof

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Computational Linguistics (AREA)
  • Mathematical Physics (AREA)
  • Fuzzy Systems (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种基于快速剪枝策略的高效用序列模式的用药DDD值数据分析方法,包括以下步骤:步骤1、对整理出用药记录表和药品DDD值表,对用药记录表中每条的用药记录的数据进行预处理;步骤2、计算序列数据库中q‑项的序列权重效用值SWU;步骤3、去除序列权重效用值SWU小于阈值ξ的q‑项,得到新的s‑序列数据库;步骤4、建立q‑项集内连接效用矩阵WICUM;步骤5、建立q‑项集间连接效用矩阵BICUM;步骤6、建立数据库中每个项的效用矩阵;步骤7、建立序列模式树,得到全部的高效用序列模式,这些高效用序列模式代表不同时间段服用各种药物,各种药物间组合后累计DDD值是超出阈值的,为高累计DDD值药物组合。本发明简化了剪枝操作的复杂性并减少了数据库的扫描次数,从而提高了挖掘高效用序列模式的效率。

Description

一种基于快速剪枝策略的高效用序列模式的用药DDD值数据 分析方法
技术领域
本发明涉及一种用药DDD(defined daily dose,平均日剂量)值数据分析方法。
背景技术
WHO在1969年制定了解剖-治疗-化学的药物分类系统,确定了将限定日剂量(defined daily dose DDD)作为用药频度分析的单位。并给其下定义为:用于主要治疗目的的成人的药物平均日剂量。
目前国内滥用抗生素的情况非常普遍,滥用抗生素会诱发细菌耐药,病原微生物为躲避药物在不断变异,耐药菌株也随之产生。抗生素在杀菌同时,也会造成人体损害。所以挖掘高累计DDD值得药物组合,及时控制使用抗生素的用量将产生重大意义。
序列是指不同项集的有序排列。例如<(电脑),(打印机)>,这个序列代表大多数顾客会在买电脑以后,过段时间会买打印机。序列模式挖掘和关联规则挖掘不同,关联规则挖掘注重项集内的联系,序列模式挖掘注重项集间的联系,并且会按照一定顺序进行排列,例如前面例子<(电脑),(打印机)>和<(打印机),(电脑)>是两个不同的序列,但在关联规则挖掘中这是两个相同的项集。
高效用序列模式是指在序列数据库中序列的效用值超过预先定义效用值的序列。挖掘高效用序列模式就是找出序列数据库中所有超过最小效用阈值的序列模式。
目前高效用序列模式挖掘算法典型代表有SPAN算法和HuspExt算法,这两种算法都用到基于序列最大剩余效用值得到上界值,但计算序列最大剩余效用值花费的代价太大,仍存在候选集大,挖掘效率低等问题。
发明内容
为了克服现有用药DDD值数据分析方法存在的候选集大、挖掘效率低的不足,本发明提供一种基于快速剪枝策略的高效用序列模式的用药DDD值数据分析方法,简化了剪枝操作的复杂性并减少了数据库的扫描次数,从而提高了挖掘高效用序列模式的效率。
本发明解决其技术问题所采用的技术方案是:
一种基于快速剪枝策略的高效用序列模式的用药DDD值数据分析方法,包括以下步骤:
步骤1、整理出用药记录表和药品DDD值表,对用药记录表中每条的用药记录的数据进行预处理;
步骤2、计算序列数据库中q-项的序列权重效用值SWU;
步骤3、去除序列权重效用值SWU小于阈值ξ的q-项,得到新的s-序列数据库;
步骤4、建立q-项集内连接效用矩阵WICUM;
步骤5、建立q-项集间连接效用矩阵BICUM;
步骤6、建立数据库中每个项的效用矩阵;
步骤7、建立序列模式树,过程如下:
7.1创建一个根结点;
7.2按照字母序建立1-q-序列的子结点;
7.3建立2-q-序列的子结点,按照字母序进行连接,连接分为q-项集内的连接和q-项集间的连接,产生的2-q-序列的子结点分为大小为1的2-q-序列和大小为2的2-q-序列,具体步骤为:
7.3.1 q-项集内的连接产生大小为1的2-q-序列,查找WICUM结构表找出相应的连接项的WICUM值,当WICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.3.3;
7.3.2 q-项集间的连接产生大小为2的2-q-序列,查找BICUM结构表找出相应的连接项的BICUM值。当BICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点,否则生成结点,转到步骤7.3.3;
7.3.3建立连接后新结点的效用矩阵,通过效用矩阵计算新结点的效用值
7.3.4计算新结点的效用值上界,新结点的效用值上界为新结点在效用矩阵中效用值与剩余效用值总和,当效用值上界大于阈值ξ时,则向下生成结点,否则不向下生成结点;
7.4建立3-q-序列的子结点,按照字母序进行连接,连接分为q-项集内的连接和q-项集间的连接,产生的3-q-序列的子结点分为大小为1的3-q-序列、大小为2的3-q-序列和大小为3的3-q-序列,具体步骤为:
7.4.1由大小为1的2-q-序列进行q-项集内的连接产生大小为1的3-q-序列,查找WICUM结构表找出相应的连接项的WICUM值。当WICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.4.5;
7.4.2由大小为1的2-q-序列进行q-项集间的连接产生大小为2的3-q-序列,查找BICUM结构表找出相应的连接项的BICUM值。当BICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.4.5;
7.4.3由大小为2的2-q-序列进行q-项集内的连接产生大小为2的3-q-序列,查找WICUM结构表找出相应的连接项的WICUM值。当WICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点,否则生成结点,转到步骤7.4.5;
7.4.4由大小为2的2-q-序列进行q-项集间的连接产生大小为3的3-q-序列,查找BICUM结构表找出相应的连接项的BICUM值;当BICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.4.5;
7.4.5建立连接后新结点的效用矩阵,通过效用矩阵计算新结点的效用值;
7.4.6计算新结点的效用值上界,新结点的效用值上界为新结点在效用矩阵中效用值与剩余效用值总和,当效用值上界大于阈值ξ时,则向下生成结点,否则不向下生成结点;
7.5循环构建结点,直到序列模式树构建完成,最终得到全部的高效用序列模式,所述高效用序列模式代表不同时间段服用各种药物,当各种药物间组合后累计DDD值超出预设阈值,判定为高累计DDD值药物组合。
本发明中,为了描述有关高效用序列的问题,与序列模式区别开来。将序列模式中的项、项集、序列、序列数据库扩展为q-项、q-项集、q-序列和q-序列数据库,设D={s1,s2,s3,…,sn}是一个包含序列效用信息的q-序列数据库,其中该数据库包含m个不同的项,即I={i1,i2,i3,…,im},sj(j=1,2,3,…,n)是数据库D中第j个q-序列(j被称为q-序列的ID),q-序列s=<l1,l2,l3,…,lh>(h为q-序列中q-项的个数或q-项集的个数),其中l代表q-项集,它的数据形式为是一个q-项并且1≤k≤n′。为了方便起见,当q-项集中只有一个q-项则去掉中括号,q-项与项不同,比项多包含数量信息也称为内部效用值q(i,s),q-项也关联一个外部效用值记为p(i)这个外部效用值可以理解为项的单价或者权重。
定义1 q-项i在q-序列s中的效用值,记为u(i,s),其定义如下
u(i,s)=p(i)*q(i,s)
定义3 q-项集l在q-序列s中的效用值,记为u(l,s),其定义如下
定义4 q-项集l在q-序列数据库中的效用值,记为u(l),其定义如下
定义5一个q-序列的效用值,记为su(s),其定义如下
su(s)=∑i∈1u(i,s)
定义6 q-序列数据库D={s1,s2,s3,…,sn}的效用值,记为u(D),其定义如下:
u(D)=∑s∈Dsu(s)
定义7给定一个序列t=<t1t2…tn>和一个q-序列s′=<(i1,q1)(i2,q2)…(in,qn)>,可以用t来匹配s′,当且仅当tk=ik并且1≤k≤n,定义为t~s′
定义8高效用序列模式。因为一个序列t在q-序列的环境下可能有多个效用值,定义选择最大的效用值作为序列效用值,记为umat(t),其定义如下
当且仅当umat(t)>ξ,则序列t是高效用序列模式,其中ξ是人为设定的最小效用值。
定义9q-序列的长度和大小。q-序列s的长度为q-序列s中包含q-项的个数,记为len(s)。例如表2中len(s1)=5,len(s3)=3,len(s4)=5。q-序列s的长度为q-序列s中含有q-项集的个数,记为size(s)。例如size(s1)=4,size(s3)=2,size(s4)=5。
定义10(SWU序列权重效用值)一个序列t的SWU值,记为SWU(t),其定义如下
例如表2中SWU(<(ac)>)=su(s3)+su(s5)=26+17=43和SWU(<a>)=su(s1)+su(s3)+su(s4)+su(s5)=32+26+32+17=107。
性质1 q-序列的序列权重效用值满足闭包属性,如果SWU(s)小于阈值ξ,则q-序列s以及s的所有超集都不是高效用序列模式。
定义11 q-项集内连接效用矩阵WICUM。WICUM结构保存的存在于q-序列数据库中q-序列里项集内一对q-项的SWU值,如果k-q-序列连接成为(k+1)-q-序列的时候,序列大小不变项集内发生连接时,根据性质1,如果k-q-序列的最后一项ia与要连接的项ib的SWU(<(iaib)>)小于阈值ξ,则序列<(iaib)>以及<(iaib)>的所有超集都不是高效用序列模式。
当q-序列数据库D中有m个q-项,WICUM的结构如下表所示
i<sub>1</sub> i<sub>2</sub> i<sub>m-1</sub>
i<sub>2</sub> SWU(&amp;lt;(i<sub>1</sub>i<sub>2</sub>)&amp;gt;)
i<sub>3</sub> SWU(&amp;lt;(i<sub>1</sub>i<sub>3</sub>)&amp;gt;) SWU(&amp;lt;(i<sub>2</sub>i<sub>3</sub>)&amp;gt;)
i<sub>m</sub> SWU(&amp;lt;(i<sub>1</sub>i<sub>m</sub>)&amp;gt;) SWU(&amp;lt;(i<sub>2</sub>i<sub>m</sub>)&amp;gt;) SWU(&amp;lt;(i<sub>m-1</sub>i<sub>m</sub>)&amp;gt;)
定义12 q-项集间连接效用矩阵BICUM。BICUM结构保存的存在于q-序列数据库中q-序列里项集间一对q-项的SWU值,如果k-q-序列连接成为(k+1)-q-序列的时候,序列大小加1,q-序列项集间发生连接时,根据性质1,如果k-q-序列的最后一项ia与要连接的项ib的SWU(&lt;iaib&gt;)小于阈值ξ,则序列&lt;iaib&gt;以及&lt;iaib&gt;的所有超集都不是高效用序列模式。BICUM与WICUM不同的是,项集间的一对项连接,是会出现相同项的连接,并且连接顺序的不同,所对应的值也会不同,例如在序列中&lt;(aa)&gt;是不存在的,但&lt;aa&gt;是存在的,并且&lt;ab&gt;和&lt;ba&gt;是不同序列。
当q-序列数据库D中有m个q-项,BICUM的结构如下表所示
定义13如果序列t存在q-序列s中,则在q-序列s中排在序列t后面的q-项集合记为s/t。
定义14给定一个只含一个项的序列t的效用矩阵。序列t的效用矩阵包含q-序列的序列号sm(m=1,2,3,…)、q-序列sm的项集号ln(n=1,2,3,…)、序列t的效用值和剩余效用值。序列t在q-序列s,q-项集l中的剩余效用值为排在t(不包含t)的后面的所有项的效用值总和,记为ru(t,l,s),当序列t在序列s中的某项集内不存在则t在s该项集的效用值和剩余效用值都为0。在其定义如下
序列t的效用矩阵如下表所示,其中u(t,l)为序列t在q-项集l范围内的效用值。
l<sub>1</sub> l<sub>2</sub> l<sub>1</sub>
s<sub>1</sub> (u(t,l<sub>1</sub>),ru(t,l<sub>1</sub>,s<sub>1</sub>)) (u(t,l<sub>2</sub>),ru(t,l<sub>2</sub>,s<sub>1</sub>)) (u(t,l<sub>n</sub>),ru(t,l<sub>n</sub>,s<sub>1</sub>))
s<sub>2</sub> (u(t,l<sub>1</sub>),ru(t,l<sub>1</sub>,s<sub>2</sub>)) (u(t,l<sub>2</sub>),ru(t,l<sub>2</sub>,s<sub>2</sub>)) (u(t,l<sub>n</sub>),ru(t,l<sub>n</sub>,s<sub>2</sub>))
s<sub>m</sub> (u(t,l<sub>1</sub>),ru(t,l<sub>1</sub>,s<sub>m</sub>)) (u(t,l<sub>2</sub>),ru(t,l<sub>2</sub>,s<sub>m</sub>)) (u(t,l<sub>n</sub>),ru(t,l<sub>n</sub>,s<sub>m</sub>))
性质2给定一个序列t和q-序列数据库D,t以及t扩展后最大的效用值不会超过
i∈s′∧t~s′∧s′∈s∧s∈D∧l∈s(u(s′)+ru(i,l,s))
称上面这个值为序列t与i的上界效用值。
定义15 2-q-序列的效用矩阵:2-q-序列的效用矩阵的效用矩阵不用再重新扫描数据库生成,只要从1-q-序列效用矩阵合成生成即可,生成过程如3和图4所示。
同理,3-q-序列也是从2-q-序列矩阵生成,所有(k+1)-q-序列都可以从k-q-序列矩阵生成。
本发明的有益效果主要表现在:采用快速剪枝策略,建立WICUM和BICUM,可以在计算序列最大剩余效用值前,对需要扩展的模式进行判断,直接查看WICUM和BICUM上的值来得出是否需要进行计算序列最大剩余效用值,避免过多计算序列剩余效用最大值,提高算法的效率。
附图说明
图1是项a和项b的2-序列的序列模式树的示意图。
图2是基于快速剪枝策略的高效用序列模式的用药DDD值数据分析方法的流程图。
图3是项集内连接的2-q-序列效用矩阵的示意图。
图4是项集间连接的2-q-序列效用矩阵的示意图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1和图2,一种基于快速剪枝策略的高效用序列模式的用药DDD值数据分析方法,包括以下步骤:
步骤1、整理出用药记录表和药品DDD值表,对用药记录表中每条的用药记录的数据进行预处理;
步骤2、计算序列数据库中q-项的序列权重效用值SWU;
步骤3、去除序列权重效用值SWU小于阈值ξ的q-项,得到新的s-序列数据库;
步骤4、建立q-项集内连接效用矩阵WICUM;
步骤5、建立q-项集间连接效用矩阵BICUM;
步骤6、建立数据库中每个项的效用矩阵;
步骤7、建立序列模式树,过程如下:
7.1创建一个根结点;
7.2按照字母序建立1-q-序列的子结点;
7.3建立2-q-序列的子结点,按照字母序进行连接,连接分为q-项集内的连接和q-项集间的连接,产生的2-q-序列的子结点分为大小为1的2-q-序列和大小为2的2-q-序列,具体步骤为:
7.3.1 q-项集内的连接产生大小为1的2-q-序列,查找WICUM结构表找出相应的连接项的WICUM值,当WICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.3.3;
7.3.2 q-项集间的连接产生大小为2的2-q-序列,查找BICUM结构表找出相应的连接项的BICUM值。当BICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点,否则生成结点,转到步骤7.3.3;
7.3.3建立连接后新结点的效用矩阵,通过效用矩阵计算新结点的效用值
7.3.4计算新结点的效用值上界,新结点的效用值上界为新结点在效用矩阵中效用值与剩余效用值总和,当效用值上界大于阈值ξ时,则向下生成结点,否则不向下生成结点;
7.4建立3-q-序列的子结点,按照字母序进行连接,连接分为q-项集内的连接和q-项集间的连接,产生的3-q-序列的子结点分为大小为1的3-q-序列、大小为2的3-q-序列和大小为3的3-q-序列,具体步骤为:
7.4.1由大小为1的2-q-序列进行q-项集内的连接产生大小为1的3-q-序列,查找WICUM结构表找出相应的连接项的WICUM值。当WICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.4.5;
7.4.2由大小为1的2-q-序列进行q-项集间的连接产生大小为2的3-q-序列,查找BICUM结构表找出相应的连接项的BICUM值。当BICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.4.5;
7.4.3由大小为2的2-q-序列进行q-项集内的连接产生大小为2的3-q-序列,查找WICUM结构表找出相应的连接项的WICUM值。当WICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点,否则生成结点,转到步骤7.4.5;
7.4.4由大小为2的2-q-序列进行q-项集间的连接产生大小为3的3-q-序列,查找BICUM结构表找出相应的连接项的BICUM值;当BICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.4.5;
7.4.5建立连接后新结点的效用矩阵,通过效用矩阵计算新结点的效用值;
7.4.6计算新结点的效用值上界,新结点的效用值上界为新结点在效用矩阵中效用值与剩余效用值总和,当效用值上界大于阈值ξ时,则向下生成结点,否则不向下生成结点;
7.5循环构建结点,直到序列模式树构建完成,最终得到全部的高效用序列模式,所述高效用序列模式代表不同时间段服用各种药物,当各种药物间组合后累计DDD值超出预设阈值时,判定为高累计DDD值药物组合。
进一步,所述步骤4中的q-项集内连接效用矩阵WICUM,WICUM结构保存的存在于q-序列数据库中q-序列里项集内一对q-项的SWU值,如果k-q-序列连接成为(k+1)-q-序列的时候,序列大小不变项集内发生连接时,如果k-q-序列的最后一项ia与要连接的项ib的SWU(&lt;(iaib)&gt;)小于阈值ξ,则序列&lt;(iaib)&gt;以及&lt;(iaib)&gt;的所有超集都不是高效用序列模式;
当q-序列数据库D中有m个q-项,WICUM的结构如下表I所示:
i<sub>1</sub> i<sub>2</sub> i<sub>m-1</sub>
i<sub>2</sub> SWU(&amp;lt;(i<sub>1</sub>i<sub>2</sub>)&amp;gt;)
i<sub>3</sub> SWU(&amp;lt;(i<sub>1</sub>i<sub>3</sub>)&amp;gt;) SWU(&amp;lt;(i<sub>2</sub>i<sub>3</sub>)&amp;gt;)
i<sub>m</sub> SWU(&amp;lt;(i<sub>1</sub>i<sub>m</sub>)&amp;gt;) SWU(&amp;lt;(i<sub>2</sub>i<sub>m</sub>)&amp;gt;) SWU(&amp;lt;(i<sub>m-1</sub>i<sub>m</sub>)&amp;gt;)
表I。
再进一步,所述步骤5中的q-项集间连接效用矩阵BICUM,BICUM结构保存的存在于q-序列数据库中q-序列里项集间一对q-项的SWU值,如果k-q-序列连接成为(k+1)-q-序列的时候,序列大小加1,q-序列项集间发生连接时,如果k-q-序列的最后一项ia与要连接的项ib的SWU(&lt;iaib&gt;)小于阈值ξ,则序列&lt;iaib&gt;以及&lt;iaib&gt;的所有超集都不是高效用序列模式;
当q-序列数据库D中有m个q-项,BICUM的结构如下表II所示:
i<sub>1</sub> i<sub>2</sub> i<sub>m</sub>
i<sub>1</sub> SWU(&amp;lt;i<sub>1</sub>i<sub>1</sub>&amp;gt;) SWU(&amp;lt;i<sub>1</sub>i<sub>2</sub>&amp;gt;) SWU(&amp;lt;i<sub>1</sub>i<sub>m</sub>&amp;gt;)
i<sub>2</sub> SWU(&amp;lt;i<sub>2</sub>i<sub>1</sub>&amp;gt;) SWU(&amp;lt;i<sub>2</sub>i<sub>2</sub>&amp;gt;) SWU(&amp;lt;i<sub>2</sub>i<sub>m</sub>&amp;gt;)
i<sub>m</sub> SWU(&amp;lt;i<sub>m</sub>i<sub>1</sub>&amp;gt;) SWU(&amp;lt;i<sub>m</sub>i<sub>2</sub>&amp;gt;) SWU(&amp;lt;i<sub>m</sub>i<sub>m</sub>&amp;gt;)
表II。
更进一步,给定一个只含一个项的序列t的效用矩阵,序列t的效用矩阵包含q-序列的序列号sm(m=1,2,3,…)、q-序列sm的项集号ln(n=1,2,3,…)、序列t的效用值和剩余效用值,序列t在q-序列s,q-项集l中的剩余效用值为排在t的后面的所有项的效用值总和,记为ru(t,l,s),当序列t在序列s中的某项集内不存在则t在s该项集的效用值和剩余效用值都为0,在其定义如下
序列t的效用矩阵如下表所示:
l<sub>1</sub> l<sub>2</sub> l<sub>n</sub>
s<sub>1</sub> (u(t,l<sub>1</sub>),ru(t,l<sub>1</sub>,s<sub>1</sub>)) (u(t,l<sub>2</sub>),ru(t,l<sub>2</sub>,s<sub>1</sub>)) (u(t,l<sub>n</sub>),ru(t,l<sub>n</sub>,s<sub>1</sub>))
s<sub>2</sub> (u(t,l<sub>1</sub>),ru(t,l<sub>1</sub>,s<sub>2</sub>)) (u(t,l<sub>2</sub>),ru(t,l<sub>2</sub>,s<sub>2</sub>)) (u(t,l<sub>n</sub>),ru(t,l<sub>n</sub>,s<sub>2</sub>))
s<sub>m</sub> (u(t,l<sub>1</sub>),ru(t,l<sub>1</sub>,s<sub>m</sub>)) (u(t,l<sub>2</sub>),ru(t,l<sub>2</sub>,s<sub>m</sub>)) (u(t,l<sub>n</sub>),ru(t,l<sub>n</sub>,s<sub>m</sub>))
给定一个序列t和q-序列数据库D,t以及t扩展后最大的效用值不会超过
i∈s'^s′~t∧s′∈s^s∈D^l∈s(u(s′)+ru(i,l,s))
称上面这个值为序列t与i的上界效用值;
2-q-序列的效用矩阵。2-q-序列的效用矩阵的效用矩阵不用再重新扫描数据库生成,只要从1-q-序列效用矩阵合成生成即可;
同理,3-q-序列也是从2-q-序列矩阵生成,所有(k+1)-q-序列都从k-q-序列矩阵生成。
本实施例中,以青霉素类和头孢菌素类用药记录为例,取五种这两类药物的具体药品,分别为双氯西林、头孢拉定、氨曲南、阿莫西林和头孢噻啶,如表1所示。表1为青霉素类和头孢菌素类常用药品名称代码表
药品名称 双氯西林 头孢拉定 头孢雷特 阿莫西林 头孢唑林
药品代码 a b c d e
表1
表2所示的为用药记录表,这表中有6个序列(s1,s2,s3,s4,s5,s6),每个序列代表一个人一段时间的用药记录,其中(d:1)里的数字代表使用该药品的天数。(d:1)表示为该病人服用了一天的阿莫西林。因为每位医生对每位病人开的处方用药量会根据每个人情况不同,用药量也会不同。所以用药量很难进行确定,为了确定每天的用药量,取该药量一天用量的平均值来确定它的用药量也就是它的DDD值,表3所示就是每种药物的DDD值。s1=&lt;(d:1);[(e:3)(c:1)];(b:7);(a:2)&gt;表示为该病人服用四个周期的抗生素,第一个周期服用了一天阿莫西林以后第二周期又同时服用了了3天的头孢唑林和1天的头孢雷特,第三周期又服用了7天的头孢拉定,第四周期服用了2天的双氯西林。
事务号 用药记录
s<sub>1</sub> &lt;(d:1);[(e:3)(c:1)];(b:7);(a:2)&gt;
s<sub>2</sub> &lt;(b:2);[(c:1)(d:3)];[(b:1)(d:2)]&gt;
s<sub>3</sub> &lt;[(c:4)(a:2)];(b:3)&gt;
s<sub>4</sub> &lt;(b:1);(b:6);(a:1);(a:2);(c:3)&gt;
s<sub>5</sub> &lt;[(a:1)(c:2)];(b:2);(d:1);(b:1)&gt;
s<sub>6</sub> &lt;[(c:3)(b:4)];(d:2)&gt;
表2
表3为药品的DDD值
药品代码 a b c d e
DDD(g) 2 2 4 1 3
表3
将用药数据表按照字母序进行排序得到表4,表4为整理后的用药记录表。
事务号 用药记录
s<sub>1</sub> &lt;(d:1);[(c:1)(e:3)];(b:7);(a:2)&gt;
s<sub>2</sub> &lt;(b:2);[(c:1)(d:3)];[(b:1)(d:2)]&gt;
s<sub>3</sub> &lt;[(a:2)(c:4)];(b:3)&gt;
s<sub>4</sub> &lt;(b:1);(b:6);(a:1);(a:2);(c:3)&gt;
s<sub>5</sub> &lt;[(a:1)(c:2)];(b:2);(d:1);(b:1)&gt;
s<sub>6</sub> &lt;[(b:4)(c:3)];(d:2)&gt;
表4
计算每一个序列的序列效用值su。su相当于用药数据表当中一条记录所有药品的累计DDD值。例如su(s1)=u(d)+u([(c,e)])+u(b)+u(a)=1*1+1*4+3*3+7*2+2*2=32,su(s2)=15,su(s3)=26,su(s4)=32,su(s5)=17,su(s6)=22,得到表5,表5表示用药记录表各事务效用值。
表5
计算每一类药品的序列权重效用值SWU,项a的SWU值就是用药记录表中存在所有包含a项记录的累计DDD值总和。SWU(a)为利尿剂的SWU值,利尿剂存在于s1,s3,s4,s5四个记录当中,计算出这四个记录的累计DDD值为SWU(a)=su(s1)+su(s3)+su(s4)+su(s5)=32+26+32+17=107;同样SWU(b)=su(s1)+su(s2)+su(s3)+su(s4)+su(s5)+su(s6)=32+15+26+32+17+22=144;SWU(c)=144;SWU(d)=86;SWU(e)=32,得到表6,表6序列权重效用值。
a b c D e
SWU 107 144 144 86 32
表6
当效用阈值设为33时,则SWU(e)=32&lt;33,因为当全部存在e的记录的累计DDD值都不大于阈值时,则所有与e连接的项的累计DDD值也不可能大于33.所以去除e项,整理数据库得到表7,表7为去除e后的用药记录表。
表7
建立WICUM,WICUM结构保存的存在于q-序列数据库中q-序列里项集内一对q-项的SWU值。建立用药记录的WICUM结构,用药记录的WICUM结构中保存的是同时服用两种药品的SWU值,因[(a,b)]和[(b,a)]是相同的,所以WICUM斜三角形状的。因用药记录表当中没有同时存在ab和ad所以SWU([(ab)])=0,SWU([(ad)])=0。SWU([(ac)])=su(s3)+su(s5)=26+17=43,SWU([(bc)])=su(s6)=22,SWU([(bd)])=su(s2)=15,SWU([(cd)])=su(s2)=15。最终构成了WICUM如表8所示,表8表示WICUM。
表8
建立BICUM,BICUM结构保存的存在于q-序列数据库中q-序列里项集间一对q-项的SWU值,建立用药记录的BICUM结构,用药记录的BICUM结构中保存的是先后服用两种药品的SWU值,因[a,b]和[b,a]不同的,[a,b]表示服用a以后下一个周期再服用b,[b,a]表示先服用b以后下一个周期再服用a。所以BICUM的结构是与WICUM的结构不一样的。因用药记录表中没有[c,c],则SWU([c,c])=0,SWU([a,a])=su(s2)=32,其余见表9所示,表9表示BICUM。
用药代号 a b c d
A 32 43 32 17
B 55 64 47 54
C 23 81 0 54
D 23 55 23 15
表9
建立数据库中每个项的效用矩阵,如表10所示为项a的效用矩阵。l代表周期,表10矩阵的第四行第二列(4,22)前面的4表示药品a出现在事务s3的l1周期服用了DDD值量为4,后面的22代表事务s3中排在药品a后面的累计DDD值为22,表10为项a的效用矩阵
l<sub>1</sub> l<sub>2</sub> l<sub>3</sub> l<sub>4</sub> l<sub>5</sub>
s<sub>1</sub> 0 0 0 (4,0) 0
s<sub>3</sub> (4,22) 0 0 0 0
s<sub>4</sub> 0 0 (2,16) (4,12) 0
s<sub>5</sub> (2,15) 0 0 0 0
s<sub>6</sub> 0 0 0 0 0
表10
表11为项b的效用矩阵:
l<sub>1</sub> l<sub>2</sub> l<sub>3</sub> l<sub>4</sub> l<sub>5</sub>
s<sub>1</sub> 0 0 (14,4) 0 0
s<sub>2</sub> (4,11) 0 (2,2) 0 0
s<sub>3</sub> 0 (6,0) 0 0 0
s<sub>4</sub> (2,30) (12,18) 0 0 0
s<sub>5</sub> 0 (4,3) 0 (2,0) 0
s<sub>6</sub> 0 (8,14) 0 0 0
表11
表12为项c的效用矩阵,其余项的效用矩阵就不一一列出。
l<sub>1</sub> l<sub>2</sub> l<sub>3</sub> l<sub>4</sub> l<sub>5</sub>
s<sub>1</sub> 0 (4,27) 0 0 0
s<sub>2</sub> 0 (4,7) 0 0 0
s<sub>3</sub> (16,6) 0 0 0 0
s<sub>4</sub> 0 0 0 0 (12,0)
s<sub>5</sub> (8,7) 0 0 0 0
s<sub>6</sub> (12,2) 0 0 0 0
表12
建立序列模式树,首先创建一个根结点,然后按照字母序建立1-q-序列的子结点。
建立2-q-序列的子结点,按照字母序进行连接。2-q-序列的子结点分为大小为1的2-q-序列和大小为2的2-q-序列。
建立大小为1的2-q-序列查找WICUM结构表找出相应的连接项的WICUM值。
以a结点为例,生成&lt;(a,b)&gt;大小为1的2-q-序列时,WICUM(&lt;(a,b)&gt;)=0,所以就不生成&lt;(a,b)&gt;结点。
生成&lt;(a,c)&gt;1-q-序列时,WICUM(&lt;(a,c)&gt;)=43&gt;33,超过了阈值。
构建&lt;(a,c)&gt;的效用矩阵。表13为&lt;(a,c)&gt;的效用矩阵:
l<sub>1</sub> l<sub>2</sub> l<sub>3</sub> l<sub>4</sub> l<sub>5</sub>
s<sub>3</sub> (20,6) 0 0 0 0
s<sub>5</sub> (10,7) 0 0 0 0
表13
计算&lt;(a,c)&gt;的真正效用值也就是&lt;(a,c)&gt;的累计DDD值,u(&lt;(a,c)&gt;,l1,s3)+u(&lt;(a,c)&gt;,l1,s5)=20+10=30。
计算&lt;(a,c)&gt;的效用值上界,&lt;(a,c)&gt;的效用值上界为&lt;(a,c)&gt;在效用矩阵中所有的效用值与剩余效用值总和,为u(&lt;(a,c)&gt;,l1,s3)+ru(&lt;(a,c)&gt;,l1,s3)+u(&lt;(a,c)&gt;,l1,s5)+r u(&lt;(a,c)&gt;,l1,s5)=20+6+10+7=43大于阈值33时,则继续向下生成结点
生成&lt;(a,d)&gt;结点时,WICUM(&lt;(a,d)&gt;)=0,所以不生成&lt;(a,d)&gt;结点。
建立长度为2的2-q-序列查找BICUM结构表找出相应的连接项的BICUM值。
生成&lt;a,a&gt;2-q-序列时,查找BICUM结构表,BICUM(&lt;a,a&gt;)=32,因为32小于效用阈值33,不生成&lt;a,a&gt;结点;
生成&lt;a,b&gt;2-q-序列时,查找BICUM结构表,BICUM(&lt;a,b&gt;)=55,超过了阈值。
构建&lt;a,b&gt;的效用矩阵,表14为&lt;a,b&gt;的效用矩阵:
l<sub>1</sub> l<sub>2</sub> l<sub>3</sub> l<sub>4</sub> l<sub>5</sub>
s<sub>3</sub> 0 (10,0) 0 0 0
s<sub>5</sub> 0 (6,3) 0 (4,0) 0
表14
计算&lt;a,b&gt;的真正效用值,umax(&lt;a,b&gt;,s3)+umax(&lt;(a,b)&gt;,s5)=10+6=16。
计算&lt;a,b&gt;的效用值上界,&lt;a,b&gt;的效用值上界为&lt;a,b&gt;在效用矩阵中所有的效用值与剩余效用值总和,为u(&lt;a,b&gt;,l2,s3)+ru(&lt;a,b&gt;,l2,s3)+u(&lt;a,b&gt;,l2,s5)+ru(&lt;a,b&gt;,l2,s5)+u(&lt;a,b&gt;,l4,s5)+ru(&lt;a,b&gt;,l4,s5)=10+6+3+4=23小于阈值33时,则不向下生成结点;
生成&lt;a,c&gt;结点时,查找BICUM表,BICUM(&lt;a,c&gt;)=23&lt;33,不生成&lt;a,c&gt;结点;
生成&lt;a,d&gt;结点时,查找BICUM表,BICUM(&lt;a,d&gt;)=23&lt;33,不生成&lt;a,d&gt;结点。
遍历所有与b结合的项,按照字母序进行连接,最终生成a,b的2-q-序列的全部子结点为&lt;(a,c)&gt;、&lt;a,b&gt;、&lt;b,a&gt;、&lt;b,b&gt;、&lt;b,c&gt;和&lt;b,d&gt;。a,b的序列模式树如图1所示。
建立3-q-序列的子结点,对序列树中的2-q-序列进行从左到右扫描,并且按照字母序进行连接。
从树的最左端&lt;(a,c)&gt;开始,建立大小为1的3-q-序列。按照字母序只能与d进行连接,查找WICUM结构表找出相应的&lt;(c,d)&gt;的WICUM值为15小于阈值33,则不生成&lt;(a,c,d)&gt;结点。
建立大小为2的3-q-序列。&lt;(a,c)&gt;按照字母序分别与a,b,c,d连接,首先与a进行连接,查找BICUM结构表找出相应的&lt;c,a&gt;的BICUM值为23小于阈值33,不生成&lt;(a,c),a&gt;结点;
&lt;(a,c)&gt;与b进行项集间连接,查找BICUM结构表找出相应的&lt;c,b&gt;的BICUM值为81大于阈值33;
构建&lt;(a,c),b&gt;的效用矩阵,表15为&lt;(a,c),b&gt;的效用矩阵:
l<sub>1</sub> l<sub>2</sub> l<sub>3</sub> l<sub>4</sub> l<sub>5</sub>
s<sub>3</sub> 0 (26,0) 0 0 0
s<sub>5</sub> 0 (14,3) 0 0 0
表15
计算&lt;(a,c),b&gt;的真正效用值,u(&lt;(a,c),b&gt;,l2,s3)+u(&lt;(a,c),b&gt;,l5,s3&gt;)=40。
计算&lt;(a,c),b&gt;的上界效用值,&lt;(a,c),b&gt;的效用值上界为&lt;(a,c),b&gt;在效用矩阵中所有的效用值与剩余效用值总和,为u(&lt;(a,c),b&gt;,l2,s3)+u(&lt;(a,c),b&gt;,l5,s3&gt;)+ru(&lt;(a,c),b&gt;,l2,s3)+ru(&lt;(a,c),b&gt;,l2,s5)=26+14+3=43大于阈值33,则继续向下生成结点。
&lt;(a,c)&gt;与c进行项集间连接,查找BICUM结构表找出相应的&lt;c,c&gt;的BICUM值为0小于阈值33,对&lt;(a,c),c&gt;进行剪枝。
&lt;(a,c)&gt;与d进行项集间连接,查找BICUM结构表找出相应的&lt;c,d&gt;的BICUM值为54大于阈值33,
构建&lt;(a,c),d&gt;的效用矩阵,表16为&lt;(a,c),d&gt;的效用矩阵:
l<sub>1</sub> l<sub>2</sub> l<sub>3</sub> l<sub>4</sub> l<sub>5</sub>
s<sub>5</sub> 0 0 0 (11,2) 0
表16
计算&lt;(a,c),d&gt;的真正效用值,u(&lt;(a,c),d&gt;,l4,s5)=11。
计算&lt;(a,c),d&gt;的上界效用值。&lt;(a,c),d&gt;的效用值上界为&lt;(a,c),d&gt;在效用矩阵中所有的效用值与剩余效用值总和,为u(&lt;(a,c),d&gt;,l4,s5)+ru(&lt;(a,c),d&gt;,l4,s5)=11+2=13小于阈值33。则不进行继续向下生成结点。
接着扫描序列树的第二项2-q-序列&lt;a,b&gt;,因前面计算过&lt;a,b&gt;的上界效用值为23小于阈值33,则不进行向下进行扩展。
其余的2-q-序列也按照也同样的步骤进行扩展
序列树的3-q-序列生成完后,继续生成4-q-序列,5-q-序列...,直到最终建立完整的序列模式树,保存所有高效用序列模式,算法结束。
最后得到三个高效用序列模式分别为&lt;b,a&gt;效用值为34,&lt;c,b&gt;效用值为58和&lt;(a,c),b&gt;效用值为40。
说明服用了头孢拉定以后再服用双氯西林累计DDD值为34、服用了头孢雷特再服用头孢拉定累计DDD值为58、服用了双氯西林和头孢雷特再服用了头孢拉定的累计DDD值为40。
说明医生在开出头孢拉定以后接下来的周期开出双氯西林;开出头孢雷特以后接下来的周期开出头孢拉定;开出双氯西林和头孢雷特以后接下来的周期开出头孢拉定,有这样三种用药习惯。
说明这三种用药习惯,用药量太大。这样会降低抗生素的使用效果。挖掘出这样的序列模式将对有关部门发现滥用抗生素等情况有着重大的意义,同时也辅助医生在用药的时候注意过多使用某种抗生素。

Claims (4)

1.一种基于快速剪枝策略的高效用序列模式的用药DDD值数据分析方法,其特征在于:包括以下步骤:
步骤1、整理出用药记录表和药品DDD值表,对用药记录表中每条用药记录的用药数据进行预处理;
步骤2、计算序列数据库中q-项的序列权重效用值SWU,将序列模式中的项、项集、序列、序列数据库扩展为q-项、q-项集、q-序列和q-序列数据库,设D={s1,s2,s3,…,sn}是一个包含序列效用信息的q-序列数据库,其中该数据库包含m个不同的项,即I={i1,i2,i3,…,im},sj是数据库D中第j个q-序列,j=1,2,3,…,n,j被称为q-序列的ID,q-序列s=&lt;l1,l2,l3,…,lh&gt;,h为q-序列中q-项的个数或q-项集的个数,其中l代表q-项集,它的数据形式为(ijk,qk)是一个q-项并且1≤k≤n′,当q-项集中只有一个q-项则去掉中括号,q-项与项不同,比项多包含数量信息也称为内部效用值q(i,s),q-项也关联一个外部效用值记为p(i),这个外部效用值为项的单价或者权重;
步骤3、去除序列权重效用值SWU小于阈值ξ的q-项,得到新的s-序列数据库;
步骤4、建立q-项集内连接效用矩阵WICUM;
步骤5、建立q-项集间连接效用矩阵BICUM;
步骤6、建立数据库中每个项的效用矩阵;
步骤7、建立序列模式树,过程如下:
7.1创建一个根结点;
7.2按照字母序建立1-q-序列的子结点;
7.3建立2-q-序列的子结点,按照字母序进行连接,连接分为q-项集内的连接和q-项集间的连接,产生的2-q-序列的子结点分为大小为1的2-q-序列和大小为2的2-q-序列,具体步骤为:
7.3.1q-项集内的连接产生大小为1的2-q-序列,查找WICUM结构表找出相应的连接项的WICUM值,当WICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.3.3;
7.3.2q-项集间的连接产生大小为2的2-q-序列,查找BICUM结构表找出相应的连接项的BICUM值,当BICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点,否则生成结点,转到步骤7.3.3;
7.3.3建立连接后新结点的效用矩阵,通过效用矩阵计算新结点的效用值7.3.4计算新结点的效用值上界,新结点的效用值上界为新结点在效用矩阵中效用值与剩余效用值总和,当效用值上界大于阈值ξ时,则向下生成结点,否则不向下生成结点;
7.4建立3-q-序列的子结点,按照字母序进行连接,连接分为q-项集内的连接和q-项集间的连接,产生的3-q-序列的子结点分为大小为1的3-q-序列、大小为2的3-q-序列和大小为3的3-q-序列,具体步骤为:
7.4.1由大小为1的2-q-序列进行q-项集内的连接产生大小为1的3-q-序列,查找WICUM结构表找出相应的连接项的WICUM值,当WICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.4.5;
7.4.2由大小为1的2-q-序列进行q-项集间的连接产生大小为2的3-q-序列,查找BICUM结构表找出相应的连接项的BICUM值,当BICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.4.5;
7.4.3由大小为2的2-q-序列进行q-项集内的连接产生大小为2的3-q-序列,查找WICUM结构表找出相应的连接项的WICUM值,当WICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点,否则生成结点,转到步骤7.4.5;
7.4.4由大小为2的2-q-序列进行q-项集间的连接产生大小为3的3-q-序列,查找BICUM结构表找出相应的连接项的BICUM值;当BICUM值小于阈值ξ时,则不进行连接操作,进行剪枝,这个连接项不产生结点;否则生成结点,转到步骤7.4.5;
7.4.5建立连接后新结点的效用矩阵,通过效用矩阵计算新结点的效用值;
7.4.6计算新结点的效用值上界,新结点的效用值上界为新结点在效用矩阵中效用值与剩余效用值总和,当效用值上界大于阈值ξ时,则向下生成结点,否则不向下生成结点;
7.5循环构建结点,直到序列模式树构建完成,最终得到全部的高效用序列模式,所述高效用序列模式代表不同时间段服用各种药物,当各种药物间组合后累计DDD值超出预设阈值,判定为高累计DDD值药物组合。
2.如权利要求1所述的基于快速剪枝策略的高效用序列模式的用药DDD值数据分析方法,其特征在于:所述步骤4中的q-项集内连接效用矩阵WICUM,WICUM结构保存的存在于q-序列数据库中q-序列里项集内一对q-项的SWU值,如果k-q-序列连接成为(k+1)-q-序列的时候,序列大小不变项集内发生连接时,如果k-q-序列的最后一项ia与要连接的项ib的SWU(&lt;(iaib)&gt;)小于阈值ξ,则序列&lt;(iaib)&gt;以及&lt;(iaib)&gt;的所有超集都不是高效用序列模式;
当q-序列数据库D中有m个q-项,WICUM的结构如下表I所示:
i<sub>1</sub> i<sub>2</sub> i<sub>m-1</sub> i<sub>2</sub> SWU(&amp;lt;(i<sub>1</sub>i<sub>2</sub>)&amp;gt;) i<sub>3</sub> SWU(&amp;lt;(i<sub>1</sub>i<sub>3</sub>)&amp;gt;) SWU(&amp;lt;(i<sub>2</sub>i<sub>3</sub>)&amp;gt;) i<sub>m</sub> SWU(&amp;lt;(i<sub>1</sub>i<sub>m</sub>)&amp;gt;) SWU(&amp;lt;(i<sub>2</sub>i<sub>m</sub>)&amp;gt;) SWU(&amp;lt;(i<sub>m-1</sub>i<sub>m</sub>)&amp;gt;)
表I。
3.如权利要求1或2所述的基于快速剪枝策略的高效用序列模式的用药DDD值数据分析方法,其特征在于:所述步骤5中的q-项集间连接效用矩阵BICUM,BICUM结构保存的存在于q-序列数据库中q-序列里项集间一对q-项的SWU值,如果k-q-序列连接成为(k+1)-q-序列的时候,序列大小加1,q-序列项集间发生连接时,如果k-q-序列的最后一项ia与要连接的项ib的SWU(&lt;iaib&gt;)小于阈值ξ,则序列&lt;iaib&gt;以及&lt;iaib&gt;的所有超集都不是高效用序列模式;
当q-序列数据库D中有m个q-项,BICUM的结构如下表II所示:
i<sub>1</sub> i<sub>2</sub> i<sub>m</sub> i<sub>1</sub> SWU(&amp;lt;i<sub>1</sub>i<sub>1</sub>&amp;gt;) SWU(&amp;lt;i<sub>1</sub>i<sub>2</sub>&amp;gt;) SWU(&amp;lt;i<sub>1</sub>i<sub>m</sub>&amp;gt;) i<sub>2</sub> SWU(&amp;lt;i<sub>2</sub>i<sub>1</sub>&amp;gt;) SWU(&amp;lt;i<sub>2</sub>i<sub>2</sub>&amp;gt;) SWU(&amp;lt;i<sub>2</sub>i<sub>m</sub>&amp;gt;) i<sub>m</sub> SWU(&amp;lt;i<sub>m</sub>i<sub>1</sub>&amp;gt;) SWU(&amp;lt;i<sub>m</sub>i<sub>2</sub>&amp;gt;) SWU(&amp;lt;i<sub>m</sub>i<sub>m</sub>&amp;gt;)
表II。
4.如权利要求1或2所述的基于快速剪枝策略的高效用序列模式的用药DDD值数据分析方法,其特征在于:给定一个只含一个项的序列t的效用矩阵,序列t的效用矩阵包含q-序列的序列号sm(m=1,2,3,…)、q-序列sm的项集号ln(n=1,2,3,…)、序列t的效用值和剩余效用值,序列t在q-序列s,q-项集l中的剩余效用值为排在t的后面的所有项的效用值总和,记为ru(t,l,s),当序列t在序列s中的某项集内不存在则t在s该项集的效用值和剩余效用值都为0,在其定义如下
序列t的效用矩阵如下表所示:
给定一个序列t和q-序列数据库D,t以及t扩展后最大的效用值不会超过
i∈s′∧s′~t∧s′∈s∧s∈D∧l∈s(u(s′)+rui,l,s)
称上面这个值为序列t与i的上界效用值;
2-q-序列的效用矩阵的效用矩阵不用再重新扫描数据库生成,只要从1-q-序列效用矩阵合成生成即可;
同理,3-q-序列也是从2-q-序列矩阵生成,所有(k+1)-q-序列都从k-q-序列矩阵生成。
CN201610172655.6A 2016-03-24 2016-03-24 一种基于快速剪枝策略的高效用序列模式的用药ddd值数据分析方法 Active CN105868296B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610172655.6A CN105868296B (zh) 2016-03-24 2016-03-24 一种基于快速剪枝策略的高效用序列模式的用药ddd值数据分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610172655.6A CN105868296B (zh) 2016-03-24 2016-03-24 一种基于快速剪枝策略的高效用序列模式的用药ddd值数据分析方法

Publications (2)

Publication Number Publication Date
CN105868296A CN105868296A (zh) 2016-08-17
CN105868296B true CN105868296B (zh) 2019-02-05

Family

ID=56625504

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610172655.6A Active CN105868296B (zh) 2016-03-24 2016-03-24 一种基于快速剪枝策略的高效用序列模式的用药ddd值数据分析方法

Country Status (1)

Country Link
CN (1) CN105868296B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108733705B (zh) * 2017-04-20 2023-04-04 哈尔滨工业大学深圳研究生院 一种高效用序列模式挖掘方法及装置
CN109460424A (zh) * 2018-10-18 2019-03-12 哈尔滨工业大学(深圳) 高效用序列模式处理方法、装置和计算机设备
CN113377766B (zh) * 2021-05-21 2022-09-13 哈尔滨工业大学(深圳) 基于效用的序列数据库对比挖掘方法、装置及计算机设备
CN113792099B (zh) * 2021-08-12 2023-08-25 上海熙业信息科技有限公司 一种基于历史效用表剪枝的数据流高效用项集挖掘系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1701343A (zh) * 2002-09-20 2005-11-23 德克萨斯大学董事会 用于信息发现以及关联分析的计算机程序产品、系统以及方法
CN102662948A (zh) * 2012-02-23 2012-09-12 浙江工商大学 一种快速发现效用模式的数据挖掘方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8732188B2 (en) * 2007-02-18 2014-05-20 Abbott Diabetes Care Inc. Method and system for providing contextual based medication dosage determination

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1701343A (zh) * 2002-09-20 2005-11-23 德克萨斯大学董事会 用于信息发现以及关联分析的计算机程序产品、系统以及方法
CN102662948A (zh) * 2012-02-23 2012-09-12 浙江工商大学 一种快速发现效用模式的数据挖掘方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"空间高效用Co-location模式挖掘技术初探";杨世晟,王丽珍等;《小型微计算机系统》;20141031;第2302-2307页

Also Published As

Publication number Publication date
CN105868296A (zh) 2016-08-17

Similar Documents

Publication Publication Date Title
CN105868296B (zh) 一种基于快速剪枝策略的高效用序列模式的用药ddd值数据分析方法
CN107656952B (zh) 平行智能病例推荐模型的建模方法
CN106228023B (zh) 一种基于本体和主题模型的临床路径挖掘方法
Lahrouz et al. Deterministic and stochastic stability of a mathematical model of smoking
Bast et al. Qlever: A query engine for efficient sparql+ text search
CN104199855B (zh) 一种针对中医药学信息的检索系统和方法
CN110825862B (zh) 基于药学知识图谱的智能问答方法及装置
CN113707297A (zh) 医疗数据的处理方法、装置、设备及存储介质
CN107391927B (zh) 一种预测药物和疾病对应关系的方法和电子设备
WO2020258483A1 (zh) 一种基于高效的负序列挖掘模式的临床用药行为分析系统及其工作方法
CN110019410A (zh) 针对中医药临床案例信息的大数据挖掘系统
CN107818816A (zh) 基于人工智能学习的医生电子病历自动生成系统
CN103914486B (zh) 文档的搜索及展现的系统
CN113421653B (zh) 医疗信息的推送方法及装置、存储介质、计算机设备
CN111599487A (zh) 一种基于关联分析的中药配伍辅助决策方法
Banek et al. Integrating different grain levels in a medical data warehouse federation
CN107832633A (zh) 关系事务数据发布的隐私保护方法
Kulshrestha et al. Prediction of antimicrobial resistance for disease-causing agents using machine learning
Denham Alyawarra kinship, infant carrying, and alloparenting
Chaturvedi et al. Fuzzy c-means based inference mechanism for association rule mining: a clinical data mining approach
JP2021077382A (ja) 知識管理システム
Bokhari et al. A framework for clustering dental patients' records using unsupervised learning techniques
Shoaib et al. URL ordering based performance evaluation of Web crawler
Schapranow et al. The medical knowledge cockpit: real-time analysis of big medical data enabling precision medicine
Syed et al. Document oriented graphical analysis and prediction

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 310012 1st floor, building 1, 223 Yile Road, Hangzhou City, Zhejiang Province

Patentee after: Yinjiang Technology Co.,Ltd.

Address before: 310012 1st floor, building 1, 223 Yile Road, Hangzhou City, Zhejiang Province

Patentee before: ENJOYOR Co.,Ltd.