diff --git a/flu_stats_unpaired_npa.R b/flu_stats_unpaired_npa.R
index bcc0c81..ff04e8e 100755
--- a/flu_stats_unpaired_npa.R
+++ b/flu_stats_unpaired_npa.R
@@ -61,7 +61,7 @@ my_adjust_method = "BH"
 lf_t1 = lf[lf$timepoint == "t1",]
 sum(is.na(lf_t1$value))
 
-#foo = lf_t1[which(is.na(lf_t1$value)),]
+foo = lf_t1[which(is.na(lf_t1$value)),]
 #ci = which(is.na(lf_t1$value))
 #lf_t1_comp = lf_t1[-ci,]
 
diff --git a/flu_stats_unpaired_sam.R b/flu_stats_unpaired_sam.R
index 67be452..a3b9ff1 100755
--- a/flu_stats_unpaired_sam.R
+++ b/flu_stats_unpaired_sam.R
@@ -42,7 +42,6 @@ rm(pivot_cols)
 # sanity checks
 table(lf$timepoint)
 length(unique(lf$mosaic))
-lf$timepoint = paste0("t", lf$timepoint)
 
 if (table(lf$flustat) == table(sam_lf$flustat)[[2]]){
   cat("Analysing Flu positive patients for:", my_sample_type)
@@ -64,7 +63,7 @@ my_adjust_method = "BH"
 lf_t1 = lf[lf$timepoint == "t1",]
 sum(is.na(lf_t1$value))
 
-#foo = lf_t1[which(is.na(lf_t1$value)),]
+foo = lf_t1[which(is.na(lf_t1$value)),]
 #ci = which(is.na(lf_t1$value))
 #lf_t1_comp = lf_t1[-ci,]
 
diff --git a/flu_stats_unpaired_serum.R b/flu_stats_unpaired_serum.R
index 42c9687..9fa2033 100755
--- a/flu_stats_unpaired_serum.R
+++ b/flu_stats_unpaired_serum.R
@@ -39,7 +39,6 @@ rm(pivot_cols)
 
 # sanity checks
 table(lf$timepoint)
-lf$timepoint = paste0("t", lf$timepoint)
 
 if (table(lf$flustat) == table(serum_lf$flustat)[[2]]){
   cat("Analysing Flu positive patients for:", my_sample_type)