From 702ff4e327814dd8f5362c6a2e95fe0f4b350af2 Mon Sep 17 00:00:00 2001
From: Tanushree Tunstall <tanu@tunstall.in>
Date: Wed, 18 Nov 2020 16:11:35 +0000
Subject: [PATCH] renamed reg_cols_extraction.R to colnames_clinical_meds.R

---
 ...s_extraction.R => colnames_clinical_meds.R |  0
 data_extraction_formatting_clinical.R         |  2 +-
 data_extraction_formatting_regression.R       |  2 +-
 flu_stats_unpaired_clinical.R                 | 30 ++-----------------
 4 files changed, 5 insertions(+), 29 deletions(-)
 rename reg_cols_extraction.R => colnames_clinical_meds.R (100%)

diff --git a/reg_cols_extraction.R b/colnames_clinical_meds.R
similarity index 100%
rename from reg_cols_extraction.R
rename to colnames_clinical_meds.R
diff --git a/data_extraction_formatting_clinical.R b/data_extraction_formatting_clinical.R
index 8ba6f7a..b6ff2c3 100644
--- a/data_extraction_formatting_clinical.R
+++ b/data_extraction_formatting_clinical.R
@@ -10,7 +10,7 @@ getwd()
 # Input: source data
 #====================
 source("read_data.R")
-source("reg_cols_extraction.R")
+source("colnames_clinical_meds.R")
 ########################################################################
 # quick sanity checks
 table(adult_df$ia_exac_copd==1 & adult_df$asthma == 1) # check this is 4 
diff --git a/data_extraction_formatting_regression.R b/data_extraction_formatting_regression.R
index 676d671..14dc360 100644
--- a/data_extraction_formatting_regression.R
+++ b/data_extraction_formatting_regression.R
@@ -10,7 +10,7 @@ getwd()
 # Input: source data for clinical
 #====================
 source("data_extraction_formatting_clinical.R")
-#source("reg_cols_extraction.R")
+#source("colnames_clinical_meds.R")
 #=======================================
 # Data for mediator to include in regression
 #=======================================
diff --git a/flu_stats_unpaired_clinical.R b/flu_stats_unpaired_clinical.R
index 989cc01..b17dd2f 100755
--- a/flu_stats_unpaired_clinical.R
+++ b/flu_stats_unpaired_clinical.R
@@ -20,31 +20,7 @@ source("data_extraction_formatting_clinical.R")
 outfile_name_clinical = paste0("flu_stats_clinical_unpaired.csv")
 outfile_clinical_unpaired = paste0(outdir_stats, outfile_name_clinical)
 outfile_clinical_unpaired 
-#===============================
-# data assignment for stats
-#================================
-wf = npa_wf[npa_wf$flustat == 1,]
-lf = npa_lf[npa_lf$flustat == 1,]
-lf$timepoint = paste0("t", lf$timepoint)
-lf = lf[!lf$mediator == "vitd",]
-########################################################################
-# clear variables
-rm(sam_lf, sam_wf
-, serum_lf, serum_wf) 
-rm(colnames_sam_df, expected_rows_sam_lf
-, colnames_serum_df, expected_rows_serum_lf)
 
-rm(pivot_cols)
-
-# sanity checks
-table(lf$timepoint)
-
-#if (table(lf$flustat) == table(npa_lf$flustat)[[2]]){
-#  cat("Analysing Flu positive patients for:", my_sample_type)
-#}else{
-#  cat("FAIL: problem with subsetting data for:", my_sample_type)
-#  quit()
-#}
 ########################################################################
 # Unpaired stats for clinical data b/w groups: wilcoxon UNpaired analysis 
 # No correction required
@@ -73,9 +49,9 @@ gathercols <- c("age", "los", "onset2final", "onsfindeath", "onset_2_initial", "
 
 clinical_lf = gather_(clinical_df_numerical, keycol, valuecol, gathercols)
 
-#==============
-# unpaired: t1
-#==============
+#====================
+# unpaired: clinical
+#====================
 sum(is.na(clinical_lf$value))
 
 foo = clinical_lf[which(is.na(clinical_lf$value)),]