数据分享|逻辑回归、随机森林、SVM支持向量机预测心脏病风险数据和模型诊断可视化...

原文链接:http://tecdat.cn/?p=24973 

简介

世界卫生组织估计全世界每年有 1200 万人死于心脏病。在美国和其他发达国家,一半的死亡是由于心血管疾病。心血管疾病的早期预后可以帮助决定改变高危患者的生活方式,从而减少并发症。本研究旨在查明心脏病最相关/风险因素,并使用机器学习预测总体风险。

数据准备 

来源

该数据集查看文末了解数据获取方式来自对居民正在进行的心血管研究。分类目标是预测患者未来是否有 10 年患冠心病 (CHD) 的风险。数据集提供了患者的信息。它包括超过 4,000 条记录和 15 个属性。

变量

每个属性都是一个潜在的风险因素。有人口、行为和医疗风险因素。

人口统计:
• 性别:男性或女性(标量)
• 年龄:患者年龄;(连续 - 尽管记录的年龄已被截断为整数,但年龄的概念是连续的)
行为
• 当前吸烟者:患者是否是当前吸烟者(标量)
• 每天吸烟数:此人一天内平均吸烟的香烟数量。(可以认为是连续的,因为一个人可以拥有任意数量的香烟,甚至半支香烟。)
• BP Meds:患者是否服用降压药(标量)
•中风:患者之前是否有中风(标量)
•  Hyp:患者是否患有高血压(标量)
• 糖尿病:患者是否患有糖尿病(标量)
• Tot Chol:总胆固醇水平(连续)
• Sys BP:收缩压(连续)
• Dia BP:舒张压(连续)
• BMI:体重指数(连续)
• 心率:心率(连续 - 在医学研究中,心率等变量虽然实际上是离散的,但由于存在大量可能值而被认为是连续的。)
• 葡萄糖:葡萄糖水平(连续)
预测变量(目标)
• 10 年患冠心病 CHD 的风险(二元:“1”表示“是”,“0”表示“否”)

心脏病预测

# 获取数据
rdaa <- read.csv(路径)
# 这边可以考虑增加变量收缩压与舒张压之差、描述收缩压、舒张压与高血压等级的变量

# 看数据结构
str(ata)

c0f6ce35902919a647bd440f1abb7a86.png

# 考虑增加变量bplevel
raw_data <- sqldf

# 对变量类别进行区分

ra_da <- map
str(ra_da )

d6c117eeb08ef456e8c2189c120cdb2f.png

数据预处理

查看和处理缺失值

# 这里我们使用mice包进行缺失值处理
aggr

1bb018c5fbb0f309a592b8a4d09f0866.jpeg

matplot

d9de05601d9de2dcf9352ead44693c9c.png


点击标题查阅往期内容

0d37eba4f562e73b36434d7b2dff3e1d.jpeg

R语言逻辑回归、Naive Bayes贝叶斯、决策树、随机森林算法预测心脏病

outside_default.png

左右滑动查看更多

outside_default.png

01

06e2c80c8b67da45a0a6719142de3b58.png

02

0ec56ed43bf1653e0d3f77ee90efe2a0.png

03

94b75a8f8579a0f89a36d57990a00b81.png

04

49a3da6811648aaa60ebc8df0aad679f.png

由上图可以看出,除了glucose变量,其它变量的缺失比例都低于5%,而glucose变量缺失率超过了10%。对此的处理策略是保留glucose变量的缺失值,直接删除其它变量的缺失值。现在处理glucose的缺失值,

# 处理glucose列
lee_a <- subset & !is.na & !is.na & !is.na & !is.na & !is.na
# 查看glce与其它变量的线性相关性确定mice的填充策略
gcog = glm(lcse ~ .)
smry(glseg)

e201b29dcdd5de729d852ae924d4b7bb.png

c2b444dcc61fd8e2ba12f3c4dce4456e.png

填充,排除不重要的变量。至于为什么不选diaBP,主要是后面的相关性分析中,这两个变量会造成多重共线性。

mice%in%  m=5,  "pmm", mai = 50, sd=2333, pint= FALSE)
#查看填充结果
smr(mc_od)

a919e47161a57bcbe1c34d074fc431dd.png

# 查看原始数据和插补后的数据分布情况
epot(mi_md)

550734de580ae213909419d5a590f952.jpeg

sipt(mcod, pch=12)

2a35052676a1971e350bea8e4445dcc3.jpeg

# 填充数据
mi_t <- complete
fir_aa$loe <- miout$guose
sum(is.na(flda))

e91777707b16cc9126c980161bf62d01.png

删除重复行

# 查看有无重复行并删除重复行
sum(duplicated

27e6cf3cadaa928764691580be26e2d2.png

comd_ata <- comdta\[!duplicated(), \]

查看离群点

#查看异常值
gplot(coedta)+geom_boxplot(ae(ftr(1),age))

ae6741150be302fc922ebef94deabb90.png

ggplot(copd\_dta)+geom\_boxplot(aes(factor(1cigDy))

b8105dd3f5aaa524d8e60442523d855a.png

ggplot(coea)+geom_boxplot(aes(factor(1),ttl))

f04a3fc08041b687f10b316a912865ce.png

ggplot(colt\_ta)+geom\_boxplot(aes(factor(1),syBP))

75bb6aa125e77187027c63e10c6cb73b.png

ggplot(comeaa)+geom_boxplot(aes(factor(1),daP))

cbe66b0f14a4a8183ef3ac9ac21aa32d.png

ggplot()+gem_boxplot(aes(factor(1),BMI))

f0ff9ecf8c48351a19451ca1d2f5f8e3.png

bc45e910dc91ffee7e4bf931d46f5acf.png

f9c7bb29975a35631231a7c8adf32e34.png

# 查看cigsPerDay
cigs\_sub <- comled\_dta
# 查看totChol,删除异常点
# 查看sysBP, 删除异常点
# 查看BMI

totChol: 总胆固醇水平大于240mg/dl已属于非常高,故删去水平值为600mg/dl的记录。sysBP: 去掉收缩压为295mg/dl的记录

# 删除各变量离群点
competedata
# 分类型变量列联分析
ggplot+geom_boxplot

01e9b92624b7c83c94da564b98f6a9a3.png

ggplot+geom_boxplot(aes,totChol,fill=TenYerCHD))

48da04348ec696f998b7474348f01502.png

cometddata %>% fitr %>% 
ggplot

5d5a53cd2124fa28476822bf95e7790c.png

1f86eda4f5b13f45e24158326ac8c8e5.png

a3061a98c2f3adca3b4a53a74a76db35.png

1f06ab54475099230879b29264dbfd4f.png

97dba4e4349edc0d2264fe6bf3547245.png

由图像知,glucose和hearRate变量有不显着的风险

table1=table
chisq.test

3f7f84f68901a747559dc4fa7b0960e4.png

table1

1325895b4fdb73dc354e5b4132c2b538.png

table2=table
chisq.test

da3ac30ff7af951156952889def7deaa.png

table3=table
chisq.test

54d0f69ecefd7baff6bdcfbf1b850ef8.png

chisq.test

72c203bdb822663ce90a113e92bc7d08.png

9adad3007cac37ec4e620d58acde5712.png

2a7cbbac0a4a55f43d95f0cf64f14856.png

9345c7b201259010b5d8bfa7fc505a2f.png

5d18c52ebcd263beb7243f1890a5166e.png

ggpairs

ad2be13bbff794637640124356495f89.jpeg

diaBP和sysBP有多重共线性的问题。 

currentSmoker变量可能不显着,下面进入模型部分。

模型

# 划分数据集

split = sample.split

train = subset

逻辑回归

# 逻辑回归模型 - 使用所有变量
fultaog = glm
summary(fulog)

8a1ab2653de1606aa16d45682fae41ab.png

a392e1c5d0280cba4c8cb404a98a2b27.png

fldaog = glm
summary(fuatLg)

3c0b44d101dd4550f74223e42f32cded.png

prdts = predict
glm_le <- table
ACCU

84980aaa81730d893df112db23638e52.png

随机森林

rfoel <- randomForest
# 获得重要性
imprace

cb00fcd340a6cb3cf262b6bf827be772.png

相关视频:Boosting原理与R语言提升回归树BRT预测短鳍鳗分布

# 选择重要的因素
rfmdel <- randomForest
# 误差
plot

7778d09145a1f083a7409ead3a657a3d.jpeg

# 获取重要性
ggplot +
   geom_bar
   geom_text

a8848e9f5ad4baf2d799d67389a8dcd8.png

这里有患病风险的误差不降反升,需要探究其中原因

# 绘制分类图像
pred<-predict
pdou_1<-predict  #输出概率
table <- table
sum(diag/sum #预测准确率

d2503a8720ff1f251ae5365e90773f06.png

plot(margin

3ef6aba200df1473347f122c1ebcfe3d.png

SVM支持向量机

# 先进行模型调优
tud <- tune.svm
summary(tud )

0c07b4d82f57bae1b72c09db923035ee.png

f157e30a64c390a09a0e4e8a8bd81271.png

fd849bb551b718a265ff7a74d73527bc.png

cfaaeb558f7b651b92ee1c97f6415d45.png

# 使用turning函数得到最佳参数设置支持向量机
mel.nd <- svm
cost=tuned$
summary(modted)

69a5cb9b3f31bdcfad11f7b9e75d5957.png

# 调用predict函数基于刚配置好的SVM模型进行类标号的预测:
sm.ne.ed <- predict
sv.tuedtble <- table
sm.ue.tbe

3365852ca0c1ae10018cdfa05e58ecf5.png

acy.s.vm <- sum(diag)/sum

01f5b3debd6da70ef188aee3ade927fa.png

模型诊断

根据上面三个模型的结果,可以看出预测结果的类别数量分布非常不均衡

sum

0e57aa6bc6fc9e92f21e6df725c5fd41.png

sum(TeYaHD == 0)

6950a2e72e360b21cbc2c95651f65196.png

针对这一现象,需要采取方法平衡数据集。

数据获取

在下面公众号后台回复“心脏病风险数据”,可获取完整数据。


316c4be2abc4fd19009d1b1efceed933.png

本文摘选《R语言逻辑回归、随机森林、SVM支持向量机预测FRAMINGHAM心脏病风险和模型诊断可视化》,点击“阅读原文”获取全文完整资料。


点击标题查阅往期内容

高维数据惩罚回归方法:主成分回归PCR、岭回归、lasso、弹性网络elastic net分析基因数

分类回归决策树交互式修剪和更美观地可视化分析细胞图像分割数据集

实现广义相加模型GAM和普通最小二乘(OLS)回归

贝叶斯线性回归和多元线性回归构建工资预测模型

Python支持向量回归SVR拟合、预测回归数据和可视化准确性检查实例

逻辑回归(对数几率回归,Logistic)分析研究生录取数据实例

广义线性模型glm泊松回归的lasso、弹性网络分类预测学生考试成绩数据和交叉验证

非线性回归nls探索分析河流阶段性流量数据和评级曲线、流量预测可视化

特征选择方法——最佳子集回归、逐步回归

线性回归和时间序列分析北京房价影响因素可视化案例

贝叶斯分位数回归、lasso和自适应lasso贝叶斯分位数回归分析免疫球蛋白、前列腺癌数据

用加性多元线性回归、随机森林、弹性网络模型预测鲍鱼年龄和可视化

PYTHON用户流失数据挖掘:建立逻辑回归、XGBOOST、随机森林、决策树、支持向量机、朴素贝叶斯和KMEANS聚类用户画像

PYTHON集成机器学习:用ADABOOST、决策树、逻辑回归集成模型分类和回归和网格搜索超参数优化

R语言集成模型:提升树boosting、随机森林、约束最小二乘法加权平均模型融合分析时间序列数据

Python对商店数据进行lstm和xgboost销售量时间序列建模预测分析

R语言用主成分PCA、 逻辑回归、决策树、随机森林分析心脏病数据并高维可视化

R语言基于树的方法:决策树,随机森林,Bagging,增强树

R语言用逻辑回归、决策树和随机森林对信贷数据集进行分类预测

spss modeler用决策树神经网络预测ST的股票

R语言中使用线性模型、回归决策树自动组合特征因子水平

R语言中自编基尼系数的CART回归决策树的实现

R语言用rle,svm和rpart决策树进行时间序列预测

python在Scikit-learn中用决策树和随机森林预测NBA获胜者

python中使用scikit-learn和pandas决策树进行iris鸢尾花数据分类建模和交叉验证

R语言里的非线性模型:多项式回归、局部样条、平滑样条、 广义相加模型GAM分析

R语言用标准最小二乘OLS,广义相加模型GAM ,样条函数进行逻辑回归LOGISTIC分类

R语言ISLR工资数据进行多项式回归和样条回归分析

R语言中的多项式回归、局部回归、核平滑和平滑样条回归模型

R语言用泊松Poisson回归、GAM样条曲线模型预测骑自行车者的数量

R语言分位数回归、GAM样条曲线、指数平滑和SARIMA对电力负荷时间序列预测

R语言样条曲线、决策树、Adaboost、梯度提升(GBM)算法进行回归、分类和动态可视化

如何用R语言在机器学习中建立集成模型?

R语言ARMA-EGARCH模型、集成预测算法对SPX实际波动率进行预测

在python 深度学习Keras中计算神经网络集成模型

R语言ARIMA集成模型预测时间序列分析

R语言基于Bagging分类的逻辑回归(Logistic Regression)、决策树、森林分析心脏病患者

R语言基于树的方法:决策树,随机森林,Bagging,增强树

R语言基于Bootstrap的线性回归预测置信区间估计方法

R语言使用bootstrap和增量法计算广义线性模型(GLM)预测置信区间

R语言样条曲线、决策树、Adaboost、梯度提升(GBM)算法进行回归、分类和动态可视化

Python对商店数据进行lstm和xgboost销售量时间序列建模预测分析

R语言随机森林RandomForest、逻辑回归Logisitc预测心脏病数据和可视化分析

R语言用主成分PCA、 逻辑回归、决策树、随机森林分析心脏病数据并高维可视化

Matlab建立SVM,KNN和朴素贝叶斯模型分类绘制ROC曲线

matlab使用分位数随机森林(QRF)回归树检测异常值

49c177ade6dee4bba76baa6daf9cbb39.png

c6511bfc39753a9181ba635618bea03c.jpeg

af57002a091d409e348765c0211f7c8a.png

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值