R语言Apriori算法关联规则对中药用药复方配伍规律药方挖掘可视化

本文主要是介绍R语言Apriori算法关联规则对中药用药复方配伍规律药方挖掘可视化,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

全文链接:http://tecdat.cn/?p=32316

我们常说的中药挖掘,一般是用药挖掘,还有穴位的挖掘,主要是想找出一些用药的规律点击文末“阅读原文”获取完整代码数据)。

相关视频

在中医挖掘中,数据的来源比较广泛,有的是通过临床收集用药处方,比如,一个著名老中医针对某一疾病的用药情况;有的是通过古籍,古代流传下来的药方;还有一种情况是在论文数据框里查找专门治疗某一疾病的文献,从中找到处方,用来分析。

Apriori算法是一种最有影响的挖掘关联规则频繁项集的算法。其核心是基于两阶段频集思想的递推算法。该关联规则在分类上属于单维、单层、布尔关联规则,Apriori 算法采用了逐层搜索的迭代的方法,算法简单明了,没有复杂的理论推导,也易于实现。

由于Apriori算法的特性,十分适合中药处方、膏方、方剂的挖掘,甚至于穴位的挖掘。

本文帮助客户得出不同处方的药物组合和频率,挖掘出药方内在的规律。

中药处方数据

e3d63ef9852efe9a0c459d0c44ddc069.png

读取数据

a_df3=read.xlsx("挖掘用.xlsx",startRow=0, colNames = F)

7693bacbbb4509b204380ca24d9888ad.png

转换数据结构

a_list=list(0)  
for(i in 1:nrow(a_df3)){  ##删除事务中的重复项目a_list[[i]]= unique(strsplit(a_df3[i,],",")[[1]])

476a5e9818fa67a9aa5119bbfbb44c69.png

将数据转换成事务类型

trans2 <- as(a_list, "tran

查看每个商品的出现频率

67f3d8a956cb99b296a36b94b66f659b.png

可以看到每个物品出现的频率,从而判断哪些物品的支持度较高。

关联规则挖掘

药对挖掘

at(dat1,parameter=list(support=0.3,minlen=2,maxle

d653f942c21678a123f59b17f79841f8.png

得到频繁规则挖掘

inspect(frequent

a870adfd2c35f66dfb6459b084d5c8b2.png

查看求得的频繁项集

spect(sort(frequentsets,by="suppo

ca2ac80bfba4228f2deec8eebc13df22.png

根据支持度对求得的频繁项集排序并查看(等价于inspect(sort(frequentsets)[1:10])。

建立模型

apriori(dat1,parame

76855d593db066459915debe73eabcfa.png

设置支持度为0.01,置信度为0.3

summary(rules)#查看规则

b09d535ea925e1b909992167ffca2403.png

查看部分规则

6668fbc45e14ce7244161728e896f905.png

查看置信度 支持度和提升度

5a29edf4cbf4313c9a8fb2306558bb21.png

可视化

绘制不同规则图形来表示支持度,置信度和提升度。

17bedc303a91031ac3e0068e12dabfa3.png

通过该图可以看到规则前项和规则后项分别有哪些物品 以及每个物品的支持度大小,支持度越大则圆圈越大。


点击标题查阅往期内容

ee9233a046830e5720abba1e8a1553ba.jpeg

R语言APRIORI关联规则、K-MEANS均值聚类分析中药专利复方治疗用药规律网络可视化

outside_default.png

左右滑动查看更多

outside_default.png

01

24d4107387d2938b7ef62fc50986ca27.png

02

e931953c6f6285b3c118f084632e40dc.png

03

64de1180d4ddcbf8316a68c672cc30a6.png

04

61c1ad430fecc68adb524623e2081cc7.png

f20b024dc2ae23a388ca942da2cd59bd.png

ules, method = NULL,  measure = "support", shading = "lift", int

2a01e871e77ede4c1e6a9c24ab0ba432.png

从该图可以看到支持度和置信度的关系,置信度越高提升度也越高。

ef314670d58d96bf8b7b032ffa220a42.png

faa0238f6b496ec597ea6e66fe4fcbee.png

从该图可以看到支持度和置信度的关系,提升度越高置信度也越高。

ules, method="matrix3D", measure="lift

59846d3fcee420e701499b1d08fdef89.png

84158b70fbad417403ffa59e7415b16c.png

从上图可以看到不同物品之间的关联关系,图中的点越大说明该物品的支持度越高,颜色越深说明该物品的提升度越高。

plot(rules, method="doubledecker" )

查看最高的支持度样本规则

ules::inspect(head(rules

185f41429cd3e8ef9fb95560b0dd6a4c.png

查看最高置信度样本规则

sort(rules, by="confidenc
nspect(head(rules

15c0903cab174d37d30b4529eb1bd9fc.png

sort(rules, by="lift

aff0a643da643a39c23ade90efab8a49.png

得到有价值规则子集

rules,subset=confidence>0.3 & support>0.2 & lift>=1
summary(x)

d67ac2937ea3d2a934a5321822d191b4.png

按照支持度排序

sort(x,by="support

e5c366829e6a767c987d60e3b6f3a445.png

按照置信度排序

inspect(sort(x,by="confide

5044d9af74570600923cb08aa6d40c1e.png

对有价值的x集合进行数据可视化。

method="grouped")

dac7ae32f7ee0e0ee25a3882d37162f5.png

组合挖掘

at(dat1,parameter=list(support=0.22,minlen=3,maxle

28c64c46106dc7af5667bca643a5436b.png

得到频繁规则挖掘

nspect(frequents

a8bf28e0347939cc2f1022edb299e5ab.png

查看求得的频繁项集

nspect(sort(frequentsets,by="sup

58b78f559ae88a24d5b7ebcacf59ed7d.png

根据支持度对求得的频繁项集排序并查看(等价于inspect(sort(frequentsets)[1:10])

建立模型

apriori(dat1,parameter=list(support=0.24

8863f831f5a5e319d2cec0aac5f942cb.png

设置支持度为0.01,置信度为0.3。

summary(rules)#查看规则

0dd40caac20c51821f8a58533d50ef67.png

查看部分规则

e93ccde7ff783ed61996b1a6ae990f93.png

查看置信度、支持度和提升度

51e46d2fe722a18415fefb47b1b7514f.png

可视化

15438fc7c5b37ad5d613e31c511a28df.png

d2114fd6cd94bb73647e94b22f497d04.png

253ece3a0a76b5ff2898acd6b038f55b.png

275a8e3b1ebbe7b1ac6473a17e153138.png

5988be1f940f24f69832b7423e4e1b2a.png

fbd8b8d866e9ccd28755b0235a98f3b1.png

从该图可以看到支持度和置信度的关系,提升度越高置信度也越高。

57ad48d3072d1ab0a2970ec70bdc782f.png

f0df3d3c6588bca2fb6c7fd4823b78cc.png

c727394b6c00078aeeb33775356f3153.png

查看最高的支持度样本规则

6bc8ef585d55054dae88cb386d07ba75.png

查看最高置信度样本规则

a7fd19604764f414db94f22eb643adb2.png

查看最高提升度样本规则

4d6f720770bd1199c31f830906d83522.png

confidence>0.3 & support>0.3 & lift>=1)    #得到有价值规则子集
summary(x)

23161cb9b10e59aaee4d224064ffcdf7.png

aspect(sort(x,by="support"))    #按照支持度排序##    lhs            rhs    support   confidence lift     
## 45 {川芎,黄芪} => {地龙} 0.3189655 0.7872340  1.602090  
## 43 {地龙,黄芪} => {川芎} 0.3189655 0.9024390  1.586105  
## 44 {川芎,地龙} => {黄芪} 0.3189655 0.8043478  1.481021  
## 42 {川芎,黄芪} => {当归} 0.3103448 0.7659574  1.615474  
## 41 {川芎,当归} => {黄芪} 0.3103448 0.8181818  1.506494  
## 40 {当归,黄芪} => {川芎} 0.3103448 0.8571429  1.506494  
## 37 {当归,地龙} => {川芎} 0.3017241 0.9210526  1.618820  
## 38 {川芎,当归} => {地龙} 0.3017241 0.7954545  1.618820  
## 39 {川芎,地龙} => {当归} 0.3017241 0.7608696  1.604743pect(sort(x,by="confidence"))    #按照置信度排序##    lhs            rhs    support   confidence lift     
## 37 {当归,地龙} => {川芎} 0.3017241 0.9210526  1.618820  
## 43 {地龙,黄芪} => {川芎} 0.3189655 0.9024390  1.586105  
## 40 {当归,黄芪} => {川芎} 0.3103448 0.8571429  1.506494  
## 41 {川芎,当归} => {黄芪} 0.3103448 0.8181818  1.506494  
## 44 {川芎,地龙} => {黄芪} 0.3189655 0.8043478  1.481021  
## 38 {川芎,当归} => {地龙} 0.3017241 0.7954545  1.618820  
## 45 {川芎,黄芪} => {地龙} 0.3189655 0.7872340  1.602090  
## 42 {川芎,黄芪} => {当归} 0.3103448 0.7659574  1.615474  
## 39 {川芎,地龙} => {当归} 0.3017241 0.7608696  1.604743

对有价值的x集合进行数据可视化

aa5130a4a26c93843dcdb3aaf7f1ddd2.png

ee8505f72db3ae8f38decf7f221cfab1.png


8e0b9fa777376373ae882ea72d283231.png

点击文末“阅读原文”

获取全文完整代码数据资料。

本文选自《R语言Apriori算法关联规则对中药用药复方配伍规律药方挖掘可视化》。

点击标题查阅往期内容

非线性混合效应 NLME模型对抗哮喘药物茶碱动力学研究

Python面板时间序列数据预测:格兰杰因果关系检验Granger causality test药品销售实例与可视化

R语言用关联规则和聚类模型挖掘处方数据探索药物配伍中的规律

用SPSS Modeler的Web复杂网络对所有腧穴进行关联规则分析

PYTHON在线零售数据关联规则挖掘APRIORI算法数据可视化

R语言关联规则模型(Apriori算法)挖掘杂货店的交易数据与交互可视化

R语言关联挖掘实例(购物篮分析)

python关联规则学习:FP-Growth算法对药品进行“菜篮子”分析

基于R的FP树fp growth 关联数据挖掘技术在煤矿隐患管理

python关联规则学习:FP-Growth算法对药品进行“菜篮子”分析

通过Python中的Apriori算法进行关联规则挖掘

Python中的Apriori关联算法-市场购物篮分析

R语言用关联规则和聚类模型挖掘处方数据探索药物配伍中的规律

在R语言中轻松创建关联网络

python主题建模可视化LDA和T-SNE交互式可视化

R语言时间序列数据指数平滑法分析交互式动态可视化

用R语言制作交互式图表和地图

如何用r语言制作交互可视化报告图表

K-means和层次聚类分析癌细胞系微阵列数据和树状图可视化比较

KMEANS均值聚类和层次聚类:亚洲国家地区生活幸福质量异同可视化分析和选择最佳聚类数

PYTHON实现谱聚类算法和改变聚类簇数结果可视化比较

有限混合模型聚类FMM、广义线性回归模型GLM混合应用分析威士忌市场和研究专利申请数据

R语言多维数据层次聚类散点图矩阵、配对图、平行坐标图、树状图可视化城市宏观经济指标数据

r语言有限正态混合模型EM算法的分层聚类、分类和密度估计及可视化

Python Monte Carlo K-Means聚类实战研究

R语言k-Shape时间序列聚类方法对股票价格时间序列聚类

R语言对用电负荷时间序列数据进行K-medoids聚类建模和GAM回归

R语言谱聚类、K-MEANS聚类分析非线性环状数据比较

R语言实现k-means聚类优化的分层抽样(Stratified Sampling)分析各市镇的人口

R语言聚类有效性:确定最优聚类数分析IRIS鸢尾花数据和可视化

Python、R对小说进行文本挖掘和层次聚类可视化分析案例

R语言k-means聚类、层次聚类、主成分(PCA)降维及可视化分析鸢尾花iris数据集

R语言有限混合模型(FMM,finite mixture model)EM算法聚类分析间歇泉喷发时间

R语言用温度对城市层次聚类、kmean聚类、主成分分析和Voronoi图可视化

R语言k-Shape时间序列聚类方法对股票价格时间序列聚类

R语言中的SOM(自组织映射神经网络)对NBA球员聚类分析

R语言复杂网络分析:聚类(社区检测)和可视化

R语言中的划分聚类模型

基于模型的聚类和R语言中的高斯混合模型

r语言聚类分析:k-means和层次聚类

SAS用K-Means 聚类最优k值的选取和分析

用R语言进行网站评论文本挖掘聚类

基于LDA主题模型聚类的商品评论文本挖掘

R语言鸢尾花iris数据集的层次聚类分析

R语言对用电负荷时间序列数据进行K-medoids聚类建模和GAM回归

R语言聚类算法的应用实例

7713e12c938468cc37e5de19f9f5cf20.png

a21239e7923f5e9f876e3207448dc73f.jpeg

be8fe5f9fdd727cfa9824e9df8d0e6fc.png

这篇关于R语言Apriori算法关联规则对中药用药复方配伍规律药方挖掘可视化的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/408307

相关文章

不懂推荐算法也能设计推荐系统

本文以商业化应用推荐为例,告诉我们不懂推荐算法的产品,也能从产品侧出发, 设计出一款不错的推荐系统。 相信很多新手产品,看到算法二字,多是懵圈的。 什么排序算法、最短路径等都是相对传统的算法(注:传统是指科班出身的产品都会接触过)。但对于推荐算法,多数产品对着网上搜到的资源,都会无从下手。特别当某些推荐算法 和 “AI”扯上关系后,更是加大了理解的难度。 但,不了解推荐算法,就无法做推荐系

康拓展开(hash算法中会用到)

康拓展开是一个全排列到一个自然数的双射(也就是某个全排列与某个自然数一一对应) 公式: X=a[n]*(n-1)!+a[n-1]*(n-2)!+...+a[i]*(i-1)!+...+a[1]*0! 其中,a[i]为整数,并且0<=a[i]<i,1<=i<=n。(a[i]在不同应用中的含义不同); 典型应用: 计算当前排列在所有由小到大全排列中的顺序,也就是说求当前排列是第

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

综合安防管理平台LntonAIServer视频监控汇聚抖动检测算法优势

LntonAIServer视频质量诊断功能中的抖动检测是一个专门针对视频稳定性进行分析的功能。抖动通常是指视频帧之间的不必要运动,这种运动可能是由于摄像机的移动、传输中的错误或编解码问题导致的。抖动检测对于确保视频内容的平滑性和观看体验至关重要。 优势 1. 提高图像质量 - 清晰度提升:减少抖动,提高图像的清晰度和细节表现力,使得监控画面更加真实可信。 - 细节增强:在低光条件下,抖

【数据结构】——原来排序算法搞懂这些就行,轻松拿捏

前言:快速排序的实现最重要的是找基准值,下面让我们来了解如何实现找基准值 基准值的注释:在快排的过程中,每一次我们要取一个元素作为枢纽值,以这个数字来将序列划分为两部分。 在此我们采用三数取中法,也就是取左端、中间、右端三个数,然后进行排序,将中间数作为枢纽值。 快速排序实现主框架: //快速排序 void QuickSort(int* arr, int left, int rig

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

科研绘图系列:R语言扩展物种堆积图(Extended Stacked Barplot)

介绍 R语言的扩展物种堆积图是一种数据可视化工具,它不仅展示了物种的堆积结果,还整合了不同样本分组之间的差异性分析结果。这种图形表示方法能够直观地比较不同物种在各个分组中的显著性差异,为研究者提供了一种有效的数据解读方式。 加载R包 knitr::opts_chunk$set(warning = F, message = F)library(tidyverse)library(phyl

透彻!驯服大型语言模型(LLMs)的五种方法,及具体方法选择思路

引言 随着时间的发展,大型语言模型不再停留在演示阶段而是逐步面向生产系统的应用,随着人们期望的不断增加,目标也发生了巨大的变化。在短短的几个月的时间里,人们对大模型的认识已经从对其zero-shot能力感到惊讶,转变为考虑改进模型质量、提高模型可用性。 「大语言模型(LLMs)其实就是利用高容量的模型架构(例如Transformer)对海量的、多种多样的数据分布进行建模得到,它包含了大量的先验

秋招最新大模型算法面试,熬夜都要肝完它

💥大家在面试大模型LLM这个板块的时候,不知道面试完会不会复盘、总结,做笔记的习惯,这份大模型算法岗面试八股笔记也帮助不少人拿到过offer ✨对于面试大模型算法工程师会有一定的帮助,都附有完整答案,熬夜也要看完,祝大家一臂之力 这份《大模型算法工程师面试题》已经上传CSDN,还有完整版的大模型 AI 学习资料,朋友们如果需要可以微信扫描下方CSDN官方认证二维码免费领取【保证100%免费

dp算法练习题【8】

不同二叉搜索树 96. 不同的二叉搜索树 给你一个整数 n ,求恰由 n 个节点组成且节点值从 1 到 n 互不相同的 二叉搜索树 有多少种?返回满足题意的二叉搜索树的种数。 示例 1: 输入:n = 3输出:5 示例 2: 输入:n = 1输出:1 class Solution {public int numTrees(int n) {int[] dp = new int