- setwd("D:/")
- library(stringr)
- da<-read.csv("E:/data/R/name1.csv")
- da<-da[which(da[,3]=="A"),]
- n<-dim(da)[1]
- da1<-da
- #把数据中的“ID=”去掉
- da1[,9]<-sapply(da[,9],function(x){substr(x,6,7)})
- ID=rep("NA",dim(da)[1])
- da<-data.frame(da,ID)
- ##第四列
- for (i in 1:n)
- {
- dis<-abs(da[,4]-da[i,4])
- re<-da[which(dis<=100),]
- re[,10]<-da1[i,9]
- write.table(re,file="start.csv",row.names = F,col.names = F,append=T,sep=",")
- #新增加第10行,表示了该行是和哪个ID 相邻
- }
- #同理可得到第5列



雷达卡

京公网安备 11010802022788号







