R语言公交地铁路线网络图实现数据挖掘实战

/*
* 提示：该行代码过长，系统自动注释不进行高亮。一键复制会移除系统注释 
* install.packages("dplyr") library("dplyr")#读取dplyr包用以排序 ###对数据读取 data=read.table("E:\201501一卡通进出站.txt",stringsAsFactors=F) ##对数据列进行命名 colnames(data)=c("逻辑卡号", "交易日期" , "交易时间", "票种", "交易代码", "交易车站", "上次交易车站") ###对数据进行预处理 for( ii in 20150101:20150131){#每天的数据 data1=data[which(data[,2]==ii),]#筛选出日期为20150101这天的数据 data2=data1[,c(2,3,6,7)]#筛选出"交易日期" ,"交易时间", "交易车站","上次交易车站"的数据 data2#查看数据 data2=data2[order(data2$交易车站),] line1=data2[substr(data2$交易车站,1,1)=="1",]#1号线 line2=data2[substr(data2$交易车站,1,1)=="2",]#2号线 ###筛选出车站为243 bus=unique(data2[,3])####################每个站的数据 for(busi in 1:length(bus)){ index=which(data2[,3]==bus[busi])#筛选出车站为243的数据行号 data3=data2[index,]#获取交易车站为243的数据 ###data3=data2[order(data2$交易车站)，]#如果不筛选车站，直接按交易车站递增排序 data4=arrange(data3,交易日期,交易时间)#对时间排序，先按年份递增排序，然后按照时间递增排序 ###按每十分钟时间分割 for (time in 6:21){ for(i in 1:6){ index=intersect(which(data4[,2]>time*10000+(i-1)*1000),which(data4[,2]<=time*10000+1000*i)) datat=data4[index,] outnum=length(which(datat[,4]!=0)) innum=length(which(datat[,4]==0)) if(i!=6)cat(file=paste("E:\",bus[busi],"车站",ii,"日一卡通进出站时间.txt"),append=TRUE,ii,"日",time,"点",i-1,"0分到",i,"0分的出站人数为",outnum," ","进站人数为",innum,"n") else cat(file=paste("E:\",bus[busi],"车站",ii,"日一卡通进出站时间.txt"),append=TRUE,ii,"日",time,"点",i-1,"0分到",time+1,"点0分的出站人数为",outnum," ","进站人数为",innum,"n") #cat(file="E:\243车站一卡通进出站时间.txt",append=TRUE,time,"点",i-1,"0分到",time+1,"点0分的出站人数为",outnum," ","进站人数为",innum,"n") } } #筛选出出站人数 dataout=data3[which(data3[,4]!=0),]#上次交易车站不为0，为出站人数 datain=data3[which(data3[,4]==0),] ###将数据进行输出 write.table(data4,paste("E:\",ii,"日 ",bus[busi],"车站一卡通进出站整理.txt"))#将数据整理好输出到指定的目录文件名 } } ####################################################################################3 ################1,2号线########## data2=data2[order(data2$交易车站),] line1=data2[substr(data2$交易车站,1,1)=="1",]#1号线 line2=data2[substr(data2$交易车站,1,1)=="2",]#2号线 #########1号线 data4=arrange(line1,交易日期,交易时间)#对时间排序，先按年份递增排序，然后按照时间递增排序 ###按每十分钟时间分割 cat(file="E:\1号线一卡通进出站时间.txt",append=TRUE, " 点", " 分"," 出站人数", " ","进站人数 " ,"n") for (time in 6:21){ for(i in 1:6){ index=intersect(which(data4[,2]>time*10000+(i-1)*1000),which(data4[,2]<=time*10000+1000*i)) datat=data4[index,] outnum=length(which(datat[,4]!=0)) innum=length(which(datat[,4]==0)) if(i!=6)cat(file="E:\1号线一卡通进出站时间.txt",append=TRUE,time," ",i-1,"0 "," ",outnum," "," ",innum,"n")#cat(time,"点",i-1,"0分到",i,"0分的出站人数为",outnum," ","进站人数为",innum,"n") else cat(file="E:\1号线一卡通进出站时间.txt",append=TRUE,time," ",i-1,"0 "," ",outnum," "," ",innum,"n")#cat(time,"点",i-1,"0分到",time+1,"点0分的出站人数为",outnum," ","进站人数为",innum,"n") # #cat(file="E:\20150101日243车站一卡通进出站时间.txt",append=TRUE,time,"点",i-1,"0分到",time+1,"点0分的出站人数为",outnum," ","进站人数为",innum,"n") } } #筛选出出站人数 dataout=data3[which(data3[,4]!=0),]#上次交易车站不为0，为出站人数 datain=data3[which(data3[,4]==0),] numout=dim(dataout)[1]#出站人数总和 numin=dim(datain)[1]#进站人数总和 ###将数据进行输出 write.table(data4,"E:\1号线一卡通进出站整理.txt")#将数据整理好输出到指定的目录文件名 ########2号线 data4=arrange(line2,交易日期,交易时间)#对时间排序，先按年份递增排序，然后按照时间递增排序 ###按每十分钟时间分割 cat(file="E:\2号线一卡通进出站时间.txt",append=TRUE, " 点", " 分"," 出站人数", " ","进站人数 " ,"n") for (time in 6:21){ for(i in 1:6){ index=intersect(which(data4[,2]>time*10000+(i-1)*1000),which(data4[,2]<=time*10000+1000*i)) datat=data4[index,] outnum=length(which(datat[,4]!=0)) innum=length(which(datat[,4]==0)) if(i!=6)cat(file="E:\2号线一卡通进出站时间.txt",append=TRUE,time," ",i-1,"0 "," ",outnum," "," ",innum,"n")#cat(time,"点",i-1,"0分到",i,"0分的出站人数为",outnum," ","进站人数为",innum,"n") else cat(file="E:\2号线一卡通进出站时间.txt",append=TRUE,time," ",i-1,"0 ", " ",outnum," "," ",innum,"n")#cat(time,"点",i-1,"0分到",time+1,"点0分的出站人数为",outnum," ","进站人数为",innum,"n") # #cat(file="E:\TB related\Service\temp\20150101日243车站一卡通进出站时间.txt",append=TRUE,time,"点",i-1,"0分到",time+1,"点0分的出站人数为",outnum," ","进站人数为",innum,"n") } } #筛选出出站人数 dataout=data3[which(data3[,4]!=0),]#上次交易车站不为0，为出站人数 datain=data3[which(data3[,4]==0),] ###将数据进行输出 write.table(data4,"E:\2号线一卡通进出站整理.txt")#将数据整理好输出到指定的目录文件名 #########1,2总和 data4=arrange(line1,交易日期,交易时间)#对时间排序，先按年份递增排序，然后按照时间递增排序 data44=arrange(line2,交易日期,交易时间)#对时间排序，先按年份递增排序，然后按照时间递增排序 cat(file="E:\1,2号线一卡通进出站时间.txt",append=TRUE, " 点", " 分"," 出站人数", " ","进站人数 " ,"n") for (time in 6:21){ for(i in 1:6){ index=intersect(which(data4[,2]>time*10000+(i-1)*1000),which(data4[,2]<=time*10000+1000*i)) index2=intersect(which(data44[,2]>time*10000+(i-1)*1000),which(data44[,2]<=time*10000+1000*i)) datat=data4[index,] datat1=data44[index2,] outnum=length(which(datat[,4]!=0)) outnum1=length(which(datat1[,4]!=0)) innum=length(which(datat[,4]==0)) innum1=length(which(datat1[,4]==0)) if(i!=6)cat(file="E:\1,2号线一卡通进出站时间.txt",append=TRUE,time," ",i-1,"0 "," ",outnum+outnum1," "," ",innum+innum1,"n")#cat(time,"点",i-1,"0分到",i,"0分的出站人数为",outnum," ","进站人数为",innum,"n") else cat(file="E:\1,2号线一卡通进出站时间.txt",append=TRUE,time," ",i-1,"0 ", " ",outnum+outnum1," "," ",innum+innum1,"n")#cat(time,"点",i-1,"0分到",time+1,"点0分的出站人数为",outnum," ","进站人数为",innum,"n") # #cat(file="E:\20150101日243车站一卡通进出站时间.txt",append=TRUE,time,"点",i-1,"0分到",time+1,"点0分的出站人数为",outnum," ","进站人数为",innum,"n") } } }
*/

#读取数据 ljhdat1=readLines("E:/ shanghai_1.txt" )ljhdat2=readLines("E:/ shanghai_2.txt") ljhdat3=readLines("E:/ shanghai_3.txt")ljhdat4=readLines("E:/ shanghai_4.txt") ljhdat5=readLines("E:/ shanghai_5.txt")bus=""#建立巴士信息库for(i in 1:length(ljhdat1)){if(ljhdat1[i]=="")bus=c(bus,ljhdat1[i-1])#提取每个巴士的路线信息}for(i in 1:length(ljhdat2)){if(ljhdat2[i]=="")bus=c(bus,ljhdat2[i-1])#提取每个巴士的路线信息}for(i in 1:length(ljhdat3)){if(ljhdat3[i]=="")bus=c(bus,ljhdat3[i-1])#提取每个巴士的路线信息}for(i in 1:length(ljhdat4)){if(ljhdat4[i]=="")bus=c(bus,ljhdat4[i-1])#提取每个巴士的路线信息}for(i in 1:length(ljhdat5)){if(ljhdat5[i]=="")bus=c(bus,ljhdat5[i-1])#提取每个巴士的路线信息}bus;bus=bus[-1]route=list(0)#建立路线信息#######################分割路线得到站点信息 #################################route[[1]]=unlist(strsplit(bus[1],split=" "))[-1]route[[1]]=route[[1]][-which(route[[1]]=="#")]#删除#号n=length(route[[1]])library(igraph)d = data.frame(route[[1]][1:n-1] ,route[[1]][2:n ]#建立邻接矩阵)g = graph.data.frame(d, directed = TRUE)plot(g )################################分割所有路线得到站点信息###########################library(igraph)route1=character(0);