Scripts relacionados con el acceso y análisis en bases de datos Access.
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

328 lines
14 KiB

2 years ago
  1. require(RODBC)
  2. sqlLastDrop<-function(conn, tablename, droplast=1,dbtype=NULL){
  3. if(sqlTables(conn) %>% filter(TABLE_NAME == "UMID") %>% nrow > 0){dbtype<-"UM"}
  4. if(sqlTables(conn) %>% filter(TABLE_NAME == "OVID") %>% nrow > 0){dbtype<-"OV"}
  5. table<-sqlFetch(conn, tablename)
  6. table<-table[1:(nrow(table)-droplast),]
  7. if (dbtype == "OV"){sqlSave(conn, table, tablename = tablename, safer = F)}
  8. if (dbtype == "UM"){
  9. sqlDrop(conn, tablename)
  10. sqlSave(conn, table, tablename = tablename, safer = F)}
  11. }
  12. sqlInitialize<-function(){
  13. library(tidyverse)
  14. library(RODBC)
  15. library(openxlsx)
  16. ## Conexión a la base de datos
  17. source("ruta_database.R", encoding = "UTF-8")
  18. }
  19. sqlBackUp<-function(dbfile=file,conn=dta,bu.dir=NULL){
  20. if(sqlTables(conn) %>% filter(TABLE_NAME == "UMID") %>% nrow > 0){bu.dir<-"BU_UM"}
  21. if(sqlTables(conn) %>% filter(TABLE_NAME == "OVID") %>% nrow > 0){bu.dir<-"BU_OVARIO"}
  22. db=strsplit(dbfile, "/")[[1]]%>% tail(n=1)
  23. bu_path<-gsub(db,bu.dir,dbfile)
  24. if (!dir.exists(bu_path)){
  25. dir.create(bu_path)
  26. print(paste0("Back Up directory ", bu_path, " created"))
  27. }
  28. cp_bu<-paste0(bu_path, "/", format(Sys.time(), format="%Y%m%d"),"-",db)
  29. file.copy(dbfile, cp_bu)
  30. }
  31. sqlShowSamples<-function(conn=dta, nhcs=nhc.test, verb=F, dbtype=NULL){
  32. if(sqlTables(dta) %>% filter(TABLE_NAME == "UMID") %>% nrow > 0){dbtype<-"UM"}
  33. if(sqlTables(dta) %>% filter(TABLE_NAME == "OVID") %>% nrow > 0){dbtype<-"OV"}
  34. if (dbtype == "OV"){
  35. db<-c("dbtables"="SAMPLES", "dbcode"="OVID", "dbsamples"="samples")
  36. query<-paste0("SELECT O.NHC,S.* FROM ",db["dbtables"]," S INNER JOIN ",db["dbcode"]," O ON O.",db["dbcode"],"=S.",db["dbcode"])
  37. }
  38. if (dbtype == "UM"){
  39. db<-c("dbtables"="MUESTRAS", "dbcode"="UMID", "dbsamples"="CODIGO")
  40. query<-paste0("SELECT O.NHC,S.* FROM ",db["dbtables"]," S INNER JOIN ",db["dbcode"]," O ON O.",db["dbcode"],"=S.",db["dbcode"])
  41. }
  42. if (nrow(sqlQuery(conn, query) %>% filter(NHC %in% nhcs)) == 0){
  43. return("No hay muestras de ningún paciente.")
  44. }
  45. if (isFALSE(verb)){
  46. sqlQuery(conn, query) %>% filter(NHC %in% nhcs) %>%
  47. group_by(NHC,UQ(rlang::sym(db["dbcode"]))) %>% summarise(Samples=length(UQ(rlang::sym(db["dbsamples"]))), Names=paste0(UQ(rlang::sym(db["dbsamples"])), collapse = ";")) %>%
  48. merge(data.frame(NHC=nhcs),all=T) %>% mutate(NHC=factor(NHC,levels = unique(nhcs))) %>% arrange(NHC)
  49. }else{
  50. sqlQuery(conn, query) %>% filter(NHC %in% nhcs)
  51. }
  52. }
  53. sqlGenOVID<-function(conn=dta, nhcs=nhc.test, verb=T, sinc=F, dbtype=NULL){
  54. if(sqlTables(dta) %>% filter(TABLE_NAME == "UMID") %>% nrow > 0){dbtype<-"UM"}
  55. if(sqlTables(dta) %>% filter(TABLE_NAME == "OVID") %>% nrow > 0){dbtype<-"OV"}
  56. if (dbtype == "OV"){
  57. db<-c("dbcode"="OVID")
  58. }
  59. if (dbtype == "UM"){
  60. db<-c("dbcode"="UMID")
  61. }
  62. dbid<-sqlFetch(conn,db["dbcode"])
  63. new.nhc<-nhcs[!nhcs %in% dbid$NHC] %>% unique()
  64. if(length(new.nhc) > 0){
  65. next.num<-gsub(db["dbcode"],"",dbid[,db["dbcode"]]) %>% as.numeric %>% max(na.rm=T)+1
  66. last.num<-next.num+(length(new.nhc)-1)
  67. newtab<-data.frame("NHC"=new.nhc, "ID"=sprintf("%s%04d",db["dbcode"],next.num:last.num)) %>% rename(!!db["dbcode"]:="ID")
  68. if(dbtype=="OV"){
  69. dbid<-rbind(dbid,newtab)
  70. }
  71. if(dbtype=="UM"){
  72. dbid<-merge(dbid, newtab, all=T) %>% select(NHC,UMID) %>% arrange(UMID)
  73. # dbid$Id<-as.numeric(rownames(dbid))
  74. dbid$NHC<-as.numeric(dbid$NHC)
  75. }
  76. rownames(dbid)<-as.character(1:nrow(dbid))
  77. dbid<-filter(dbid, NHC %in% new.nhc) %>% mutate(NHC=as.character(NHC))
  78. if (sinc){
  79. ### !! Atención, esto cambia la base de datos:
  80. sqlSave(conn, dbid, tablename=db["dbcode"], append = T, rownames = F)
  81. print("La base ha sido actualizada.")
  82. }
  83. if (verb){
  84. return(dbid)
  85. }
  86. }else{
  87. print("No hay pacientes nuevos.")
  88. }
  89. }
  90. sqlWriteTemp<-function(conn=dta, nhcs=nhc.test, file="queryOV.xlsx", samples.mod=T, clinics.mod=T, dbtype=NULL, today=F){
  91. if(sqlTables(dta) %>% filter(TABLE_NAME == "UMID") %>% nrow > 0){dbtype<-"UM"}
  92. if(sqlTables(dta) %>% filter(TABLE_NAME == "OVID") %>% nrow > 0){dbtype<-"OV"}
  93. if (dbtype=="OV"){
  94. upd.ovid<-sqlFetch(conn, "OVID") %>% filter(NHC %in% nhcs)
  95. if (samples.mod){
  96. ## Generar código para las nuevas muestras
  97. samples<-sqlFetch(conn, "SAMPLES")
  98. if(sum(grepl(paste0("OV",Sys.time() %>% format("%y")), samples$samples)) > 0){
  99. next.samp<-gsub(paste0("OV",Sys.time() %>% format("%y")),"", samples$samples) %>% as.numeric %>% max(na.rm=T)+1
  100. }else{
  101. next.samp<-1
  102. }
  103. last.samp<-next.samp+(length(nhcs)-1)
  104. new.samp<-sprintf("OV%s%02d",Sys.time() %>% format("%y"),next.samp:last.samp)
  105. new.samp.df<-merge(sqlFetch(dta,"OVID") %>% merge(data.frame("NHC"=nhcs)), data.frame("NHC"=nhcs, "samples"=new.samp))
  106. samples.exp<-merge(samples %>% slice(0), new.samp.df %>% select(-NHC), all=T) %>% select(colnames(samples)) %>% arrange(samples)
  107. }
  108. if (clinics.mod){
  109. ## Importar los datos clínicos de pacientes existentes y generar nueva entrada par los nuevos
  110. upd.clinics<-sqlFetch(conn, "CLINICS")
  111. ovid.new<-sqlFetch(conn, "OVID") %>% filter(NHC %in% nhcs)
  112. upd.clinics<-merge(ovid.new,upd.clinics, all.x=T, by="OVID")
  113. upd.clinics$NHC<-as.character(upd.clinics$NHC)
  114. for (i in colnames(upd.clinics)[sapply(upd.clinics, lubridate::is.POSIXct)]){upd.clinics[,i]<-as.Date(upd.clinics[,i])}
  115. }
  116. ## Exportar tablas a la plantilla de entrada para su rellenado
  117. wb <- loadWorkbook(file)
  118. writeData(wb, "NHC", upd.ovid)
  119. if (samples.mod){writeData(wb,"samples",samples.exp)}
  120. if (clinics.mod){writeData(wb,"CLINICS",upd.clinics)}
  121. saveWorkbook(wb,file,overwrite = TRUE)
  122. }
  123. if (dbtype=="UM"){
  124. upd.umid<-sqlFetch(conn, "UMID") %>% filter(NHC %in% nhcs)
  125. if (samples.mod){
  126. ## Generar código para las nuevas muestras
  127. samples<-sqlFetch(conn, "MUESTRAS")
  128. if(sum(grepl(paste0("UM",Sys.time() %>% format("%y")), samples$CODIGO)) > 0){
  129. next.samp<-gsub(paste0("UM",Sys.time() %>% format("%y")),"", samples$CODIGO) %>% as.numeric %>% max(na.rm=T)+1
  130. }else{
  131. next.samp<-1
  132. }
  133. last.samp<-next.samp+(length(nhcs)-1)
  134. new.samp<-sprintf("UM%s%02d",Sys.time() %>% format("%y"),next.samp:last.samp)
  135. new.samp.df<-data.frame("NHC"=nhcs, "CODIGO"=new.samp) %>% merge(sqlFetch(dta,"UMID"), all.x=T) %>% arrange(CODIGO)
  136. samples.exp<-merge(samples %>% slice(0), new.samp.df %>% select(-NHC), all=T) %>% select(colnames(samples)) %>% arrange(CODIGO)
  137. if (today==TRUE){
  138. samples.exp$FECHA_RECEPCION<-format(Sys.Date(), "%d/%m/%y")
  139. }
  140. }
  141. nhc.table<-read.xlsx(file, sheet = "NHC")
  142. if (any(sapply(nhc.table$Samples, function(x) "cnag" %in% strsplit(x,",")[[1]]) == T)){
  143. nhcs.cnag<-nhc.table[sapply(nhc.table$Samples, function(x) "cnag" %in% strsplit(x,",")[[1]]),"NHC"]
  144. umid.cnag<-sqlFetch(conn, "UMID") %>% filter(NHC %in% nhcs.cnag) %>% pull(UMID)
  145. sample.cnag<-samples.exp %>% filter(UMID %in% umid.cnag) %>% pull(CODIGO)
  146. cnag.exp<-merge(data.frame("UMID"=umid.cnag, "CODIGO"=sample.cnag), sqlFetch(dta, "CNAG") %>% slice(0), all=T)
  147. if (today==TRUE){
  148. cnag.exp$FECHA_ENVIO<-format(Sys.Date(), "%d/%m/%y")
  149. }
  150. }else{cnag.exp<-sqlFetch(dta, "CNAG") %>% slice(0)}
  151. if (any(sapply(nhc.table$Samples, function(x) "rna" %in% strsplit(x,",")[[1]]) == T)){
  152. nhcs.rna<-nhc.table[sapply(nhc.table$Samples, function(x) "rna" %in% strsplit(x,",")[[1]]),"NHC"]
  153. umid.rna<-sqlFetch(conn, "UMID") %>% filter(NHC %in% nhcs.rna) %>% pull(UMID)
  154. sample.rna<-samples.exp %>% filter(UMID %in% umid.rna) %>% pull(CODIGO)
  155. rna.exp<-merge(data.frame("UMID"=umid.rna, "CODIGO"=sample.rna), sqlFetch(dta, "RNADNA") %>% slice(0), all=T)
  156. }else{rna.exp<- sqlFetch(dta, "RNADNA") %>% slice(0)}
  157. if (clinics.mod){
  158. ## Importar los datos clínicos de pacientes existentes y generar nueva entrada par los nuevos
  159. upd.clinics<-sqlFetch(conn, "CLINICOS")
  160. umid.new<-sqlFetch(conn, "UMID") %>% filter(NHC %in% nhcs)
  161. upd.clinics<-merge(umid.new,upd.clinics, all.x=T, by="UMID")
  162. upd.clinics$NHC<-as.character(upd.clinics$NHC)
  163. for (i in colnames(upd.clinics)[sapply(upd.clinics, lubridate::is.POSIXct)]){upd.clinics[,i]<-as.Date(upd.clinics[,i])}
  164. }
  165. ## Exportar tablas a la plantilla de entrada para su rellenado
  166. wb <- loadWorkbook(file)
  167. writeData(wb, "NHC", merge(nhc.table, upd.umid, sort = F))
  168. if (samples.mod){writeData(wb,"samples",samples.exp)}
  169. if (clinics.mod){writeData(wb,"CLINICS",upd.clinics)}
  170. writeData(wb,"CNAG",cnag.exp)
  171. writeData(wb,"RNADNA",rna.exp)
  172. saveWorkbook(wb,file,overwrite = TRUE)
  173. }
  174. }
  175. sqlSincBD<-function(conn=dta, filetemp="queryOV.xlsx", sinc.samples=F, sinc.clinics=F, dbtype=NULL){
  176. if(sqlTables(dta) %>% filter(TABLE_NAME == "UMID") %>% nrow > 0){dbtype<-"UM"}
  177. if(sqlTables(dta) %>% filter(TABLE_NAME == "OVID") %>% nrow > 0){dbtype<-"OV"}
  178. ## Añadir código de muestra nueva a la base de datos
  179. if (dbtype == "OV"){
  180. print("DB OV detectada")
  181. nsamples<-sqlFetch(conn, "SAMPLES") %>% nrow
  182. upd.samples<-read.xlsx(filetemp, sheet = "samples", detectDates = T)
  183. if (nrow(upd.samples) > 0){rownames(upd.samples)<-(nsamples+1):(nsamples+nrow(upd.samples)) %>% as.character}
  184. if (sinc.samples & nrow(upd.samples) > 0){
  185. upd.samples$IQ_date<-as.Date(upd.samples$IQ_date)
  186. upd.samples$TIL_date<-as.Date(upd.samples$TIL_date)
  187. ### !! Atención, esto cambia la base de datos:
  188. sqlSave(conn, upd.samples, tablename="SAMPLES", append = T, varTypes = c("IQ_date"="date","TIL_date"="date"))
  189. print("Tabla SAMPLES sincronizada.")
  190. }
  191. }
  192. if (dbtype == "UM"){
  193. print("DB UM detectada")
  194. nsamples<-sqlFetch(conn, "MUESTRAS") %>% nrow
  195. upd.samples<-read.xlsx(filetemp, sheet = "samples", detectDates = T)
  196. if (nrow(upd.samples) > 0){rownames(upd.samples)<-(nsamples+1):(nsamples+nrow(upd.samples)) %>% as.character}
  197. if (sinc.samples & nrow(upd.samples) > 0){
  198. upd.samples$FECHA_RECEPCION<-as.Date(upd.samples$FECHA_RECEPCION)
  199. upd.samples$TIPO<-as.character(upd.samples$TIPO)
  200. upd.samples$OBS<-as.character(upd.samples$OBS)
  201. ### !! Atención, esto cambia la base de datos:
  202. sqlSave(conn, upd.samples, tablename="MUESTRAS", append = T, varTypes = c("FECHA_RECEPCION"="Date"), rownames = F)
  203. print("Tabla MUESTRAS sincronizada.")
  204. }
  205. }
  206. ## Añadir datos clínicos modificados a la base de datos
  207. if (dbtype == "OV"){
  208. upd.clinics<-read.xlsx(filetemp, sheet = "CLINICS",detectDates = T)
  209. ovid.mod<-upd.clinics$OVID[upd.clinics$OVID %in% (sqlFetch(dta, "CLINICS") %>% pull(OVID))]
  210. rnames<-sqlFetch(conn, "CLINICS") %>% filter(OVID %in% ovid.mod) %>% rownames
  211. clinics.mod<-upd.clinics %>% filter(OVID %in% ovid.mod) %>% select(-NHC)
  212. rownames(clinics.mod)<-rnames
  213. ### !! Atención, esto cambia la base de datos:
  214. if (sinc.clinics){
  215. fechas<-colnames(clinics.mod)[grepl("DO|date", colnames(clinics.mod))]
  216. for (i in fechas){
  217. clinics.mod[,i]<-as.Date(clinics.mod[,i])
  218. }
  219. sqlUpdate(conn, clinics.mod,"CLINICS")
  220. print("Tabla CLINICS modificada.")
  221. }
  222. }
  223. if (dbtype == "UM"){
  224. upd.clinics<-read.xlsx(filetemp, sheet = "CLINICS",detectDates = T)
  225. umid.mod<-upd.clinics$UMID[upd.clinics$UMID %in% (sqlFetch(dta, "CLINICOS") %>% pull(UMID))]
  226. rnames<-sqlFetch(conn, "CLINICOS") %>% filter(UMID %in% umid.mod) %>% rownames
  227. clinics.mod<-upd.clinics %>% filter(UMID %in% umid.mod) %>% select(-NHC)
  228. rownames(clinics.mod)<-rnames
  229. ### !! Atención, esto cambia la base de datos:
  230. if (sinc.clinics){
  231. fechas<-colnames(clinics.mod)[grepl("date|MET_DX|DoB", colnames(clinics.mod), ignore.case = T)]
  232. for (i in fechas){
  233. clinics.mod[,i]<-as.Date(clinics.mod[,i])
  234. }
  235. sqlUpdate(conn, clinics.mod,"CLINICOS", index="UMID")
  236. print("Tabla CLINICOS modificada.")
  237. }
  238. }
  239. ## Añadir datos clínicos nuevos a la base de datos
  240. if (dbtype == "OV"){
  241. nsamples.clin<-sqlFetch(conn, "CLINICS") %>% nrow
  242. ovid.new<-upd.clinics$OVID[!upd.clinics$OVID %in% (sqlFetch(conn, "CLINICS") %>% pull(OVID))]
  243. clinics.new<-upd.clinics %>% filter(OVID %in% ovid.new) %>% select(-NHC)
  244. if (length(ovid.new) > 0){rownames(clinics.new)<-(nsamples.clin+1):(nsamples.clin+nrow(clinics.new)) %>% as.character}
  245. ### !! Atención, esto cambia la base de datos:
  246. if (sinc.clinics){
  247. fechas<-colnames(clinics.new)[grepl("DO|date", colnames(clinics.new))]
  248. varTypes<-rep("Date",length(fechas))
  249. names(varTypes)<-fechas
  250. for (i in fechas){
  251. clinics.new[,i]<-as.Date(clinics.new[,i])
  252. }
  253. sqlSave(conn, clinics.new, tablename="CLINICS", append = T, varTypes = varTypes)
  254. print("Tabla CLINICS sincronizada.")
  255. }
  256. }
  257. if (dbtype == "UM"){
  258. nsamples.clin<-sqlFetch(conn, "CLINICOS") %>% nrow
  259. umid.new<-upd.clinics$UMID[!upd.clinics$UMID %in% (sqlFetch(conn, "CLINICOS") %>% pull(UMID))]
  260. clinics.new<-upd.clinics %>% filter(UMID %in% umid.new) %>% select(-NHC)
  261. if (length(umid.new) > 0){rownames(clinics.new)<-(nsamples.clin+1):(nsamples.clin+nrow(clinics.new)) %>% as.character}
  262. ### !! Atención, esto cambia la base de datos:
  263. if (sinc.clinics){
  264. fechas<-colnames(clinics.new)[grepl("date|MET_DX|DoB", colnames(clinics.new), ignore.case = T)]
  265. varTypes<-rep("Date",length(fechas))
  266. names(varTypes)<-fechas
  267. for (i in fechas){
  268. clinics.new[,i]<-as.Date(clinics.new[,i])
  269. }
  270. sqlSave(conn, clinics.new, tablename="CLINICOS", append = T, varTypes = varTypes, rownames = F)
  271. print("Tabla CLINICOS sincronizada.")
  272. }
  273. }
  274. }
  275. sqlMultiSamples<-function(kbl=F, NHC=F, full=F){
  276. query<-sqlQuery(dta, "SELECT U.UMID,U.NHC,M.FECHA_RECEPCION,M.TIPO,M.CODIGO,C.CODIGO,R.CODIGO,R.ESTADO
  277. FROM ((UMID U
  278. LEFT OUTER JOIN MUESTRAS M ON U.UMID=M.UMID)
  279. LEFT OUTER JOIN CNAG C ON M.CODIGO=C.CODIGO)
  280. LEFT OUTER JOIN RNADNA R ON M.CODIGO=R.CODIGO") %>% rename("CNAG"="CODIGO.1","RNADNA"="CODIGO.2")
  281. if (full==F){query<- query %>% filter(!is.na(CODIGO))}
  282. if (NHC==F){query<- query %>% select(-NHC)}
  283. query<-query %>% mutate(
  284. CNAG=case_when(!is.na(CNAG)~"X",TRUE~""),
  285. RNADNA=case_when((!is.na(RNADNA) & (ESTADO=="ENV"))~"X",
  286. (!is.na(RNADNA) & (is.na(ESTADO)))~".",
  287. TRUE~"")
  288. ) %>% select(-ESTADO)
  289. if (kbl==T){
  290. query %>% kableExtra::kbl() %>% kableExtra::kable_styling(full_width = F, bootstrap_options = c("striped"))
  291. }else{
  292. return(query)
  293. }
  294. }