[Rscript]逻辑回归识别学生群体的R实现

2024-04-16 09:08

本文主要是介绍[Rscript]逻辑回归识别学生群体的R实现,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

library(ggplot2)
library(pROC)###读取数据并查看数据情况
setwd("D:\\student_recognition")
student<-read.csv("student.csv")
str(student) #查看变量
apply(student,2,function(x){mean(x=='NULL'|x=="")
})  #查看缺失值情况###数据预处理
#略###描述性分析
#绘制箱线图,可见https://blog.csdn.net/TOMOCAT/article/details/80559006
#分组查看变量
library(plyr) #可进行类似数据透视表的操作,将数据分割成更小的数据,对分割后的数据进行操作
ddply(student[,c("is_student","student_app")],c("is_student"),summarise,mean_app=mean(student_app))
##   whether_student  mean_app
## 1               0 0.1745353
## 2               1 0.4435253###逻辑回归建模
model_student<-student[,c("is_student",cont_var,class_var)]
#这里cont_var是经过数据清洗和预处理之后的用于建模的连续变量,class_var是分类变量
#计算建模的平均auc
auc<-NULL
for(i in 1:5){ind<-sample(1:dim(model_student)[1],round(dim(model_student)[1]*0.3)) #对数据集按照记录数划分成训练集和测试集train<-model_student[-ind,]test<-model_student[ind,]for (var in cont_var){train[,var]<-(train[,var]-mean(train[,var]))/sd(train[,var])test[,var]<-(test[,var]-mean(test[,var]))/sd(test[,var])} #连续变量标准化for (var in class_var){train[,var]<-as.factor(train[,var])test[,var]<-as.factor(test[,var])} #离散变量转化为因子型lm_res<-glm(is_student~.,data=train,family="binomial")pre<-predict(lm_res,test,type="reponse")auc0<-auc(test$is_student,pre)print(auc0)auc=c(auc,auc0)
}
mean(auc) #相当于五次留出法的平均auc值
#从所有变量建模到用逐步回归筛选变量
lm_res_all<-glm(is_student~.,data=model_student,family="binomial")
step_lm<-step(lm_res_all) #逐步回归筛选变量
#绘制逻辑回归变量系数大小图
coef = as.data.frame(step_lm$coefficients)
coef$var = row.names(coef)
colnames(coef)[1] = "coef"
coef = coef[-1,]
coef$pos = coef$coef>0
ggplot(coef,aes(x = reorder(var,-coef),y = coef,fill = pos))+geom_bar(stat = 'identity',position = 'identity')+scale_x_discrete(labels = c("有无学生类APP","score","学校线路占比","总出行次数","平均出行距离","出行时段-早高峰","工作日出行时间标准差","出行时段-平峰","周末出行次数","周末出行时间标准差","工作日出行时段-平峰","周末出行时段-平峰","工作日出行时段-早高峰","学校线路出行次数","周末出行时段-早高峰","工作日出行时段-无","工作日平均出行距离","周末平均出行距离","出行时间标准差","工作日出行次数","周末出行时段-无"))+theme(axis.text.x = element_text(angle = 50,hjust = 1,vjust = 1))+labs(fill = "系数正负",x = "变量",y="系数")
#可视化结果参照我的R可视化相关博文

 

这篇关于[Rscript]逻辑回归识别学生群体的R实现的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

基于Python实现一个图片拆分工具

《基于Python实现一个图片拆分工具》这篇文章主要为大家详细介绍了如何基于Python实现一个图片拆分工具,可以根据需要的行数和列数进行拆分,感兴趣的小伙伴可以跟随小编一起学习一下... 简单介绍先自己选择输入的图片,默认是输出到项目文件夹中,可以自己选择其他的文件夹,选择需要拆分的行数和列数,可以通过

Python中将嵌套列表扁平化的多种实现方法

《Python中将嵌套列表扁平化的多种实现方法》在Python编程中,我们常常会遇到需要将嵌套列表(即列表中包含列表)转换为一个一维的扁平列表的需求,本文将给大家介绍了多种实现这一目标的方法,需要的朋... 目录python中将嵌套列表扁平化的方法技术背景实现步骤1. 使用嵌套列表推导式2. 使用itert

Python使用pip工具实现包自动更新的多种方法

《Python使用pip工具实现包自动更新的多种方法》本文深入探讨了使用Python的pip工具实现包自动更新的各种方法和技术,我们将从基础概念开始,逐步介绍手动更新方法、自动化脚本编写、结合CI/C... 目录1. 背景介绍1.1 目的和范围1.2 预期读者1.3 文档结构概述1.4 术语表1.4.1 核

在Linux中改变echo输出颜色的实现方法

《在Linux中改变echo输出颜色的实现方法》在Linux系统的命令行环境下,为了使输出信息更加清晰、突出,便于用户快速识别和区分不同类型的信息,常常需要改变echo命令的输出颜色,所以本文给大家介... 目python录在linux中改变echo输出颜色的方法技术背景实现步骤使用ANSI转义码使用tpu

Python使用python-can实现合并BLF文件

《Python使用python-can实现合并BLF文件》python-can库是Python生态中专注于CAN总线通信与数据处理的强大工具,本文将使用python-can为BLF文件合并提供高效灵活... 目录一、python-can 库:CAN 数据处理的利器二、BLF 文件合并核心代码解析1. 基础合

Python使用OpenCV实现获取视频时长的小工具

《Python使用OpenCV实现获取视频时长的小工具》在处理视频数据时,获取视频的时长是一项常见且基础的需求,本文将详细介绍如何使用Python和OpenCV获取视频时长,并对每一行代码进行深入解析... 目录一、代码实现二、代码解析1. 导入 OpenCV 库2. 定义获取视频时长的函数3. 打开视频文

golang版本升级如何实现

《golang版本升级如何实现》:本文主要介绍golang版本升级如何实现问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录golanwww.chinasem.cng版本升级linux上golang版本升级删除golang旧版本安装golang最新版本总结gola

SpringBoot中SM2公钥加密、私钥解密的实现示例详解

《SpringBoot中SM2公钥加密、私钥解密的实现示例详解》本文介绍了如何在SpringBoot项目中实现SM2公钥加密和私钥解密的功能,通过使用Hutool库和BouncyCastle依赖,简化... 目录一、前言1、加密信息(示例)2、加密结果(示例)二、实现代码1、yml文件配置2、创建SM2工具

Mysql实现范围分区表(新增、删除、重组、查看)

《Mysql实现范围分区表(新增、删除、重组、查看)》MySQL分区表的四种类型(范围、哈希、列表、键值),主要介绍了范围分区的创建、查询、添加、删除及重组织操作,具有一定的参考价值,感兴趣的可以了解... 目录一、mysql分区表分类二、范围分区(Range Partitioning1、新建分区表:2、分

MySQL 定时新增分区的实现示例

《MySQL定时新增分区的实现示例》本文主要介绍了通过存储过程和定时任务实现MySQL分区的自动创建,解决大数据量下手动维护的繁琐问题,具有一定的参考价值,感兴趣的可以了解一下... mysql创建好分区之后,有时候会需要自动创建分区。比如,一些表数据量非常大,有些数据是热点数据,按照日期分区MululbU