用R进行epiinfo 数据的比对

来源:互联网 发布:js img.value 编辑:程序博客网 时间:2024/06/06 00:30
library(RODBC)library(compareDF)mydata <- function(m,n){  path<- paste0("C:/Users/Administrator/Desktop/",m)  setwd(path)  channel<-odbcConnectAccess("dongtai6.mdb")  data1<-sqlFetch(channel, "dongtai6")  data2<-sqlFetch(channel, "dongtai61")  data3<-sqlFetch(channel, "dongtai63")  data4<-sqlFetch(channel, "dongtai64")  data5<-sqlFetch(channel, "dongtai65")  data_1 <- cbind(data2,data3,data4,data5)  #去掉GlobalRecordId变量  data_1$GlobalRecordId <- NULL  data_1$GlobalRecordId <- NULL  data_1$GlobalRecordId <- NULL  data_1$GlobalRecordId <- NULL  clerkcdc <- rep(n,length(data1$UniqueKey))  data<- data.frame(data_1,clerkcdc)}#项目A_firstA1sun <- mydata("epiinfo数据比对在R中的实现/项目A/项目Afirst/项目A[孙]","Asun")A1qian <- mydata("epiinfo数据比对在R中的实现/项目A/项目Afirst/项目Af[钱]","Aqian")A1zhao<- mydata("epiinfo数据比对在R中的实现/项目A/项目Afirst/项目Af[赵]","Azhao")Afirst <- rbind(A1sun,A1qian,A1zhao)#项目A_second   C:\Users\Administrator\Desktop\epiinfo数据比对在R中的实现\项目A\项目Bsecond\项目As[赵]A2qian <- mydata("epiinfo数据比对在R中的实现/项目A/项目Bsecond/项目As[钱]","A2qian")A2sun <- mydata("epiinfo数据比对在R中的实现/项目A/项目Bsecond/项目As[孙]","A2sun")A2zhao <- mydata("epiinfo数据比对在R中的实现/项目A/项目Bsecond/项目As[赵]","A2zhao")Asecond <- rbind(A2qian,A2sun,A2zhao)#B_first  B1qian <- mydata("epiinfo数据比对在R中的实现/项目B/第一次录入/项目Bf[钱]","B1qian")B1zhao<- mydata("epiinfo数据比对在R中的实现/项目B/第一次录入/项目Bf[赵]","B1zhao")Bfirst <- rbind(B1qian,B1zhao)#B_second   B2qian <- mydata("epiinfo数据比对在R中的实现/项目B/第二次录入/项目Bs[钱 ]","B2qian")B2zhao<- mydata("epiinfo数据比对在R中的实现/项目B/第二次录入/项目Bs[赵]","B2zhao")Bsecond <- rbind(B2qian,B2zhao)allfirst <- rbind(Afirst,Bfirst)allsecond <- rbind(Asecond,Bsecond)##  去掉三个变量var_dto6f <- names(allfirst) %in% c("entrydate","clerkcdc","relationship")allfirstNew <-allfirst[!var_dto6f]var_dto6s <- names(allsecond) %in% c("entrydate","clerkcdc","relationship")allsecondNew <-allsecond[!var_dto6s]com_df <- compare_df(allfirstNew,allsecondNew,c("childnumber"),exclude = NULL, limit_html = 800)com_df$html_output

这里写图片描述