apparently modifed flu stats for npa: check!

flu_stats_unpaired_npa.R

@@ -35,8 +35,12 @@ lf = npa_adults_lf[npa_adults_lf$flustat == 1,]
 table(lf$timepoint)
 lf$timepoint = paste0("t", lf$timepoint)
 
-my_sample_type = "npa"
-
+if (table(lf$flustat) == table(npa_adults_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 at each timepoint b/w groups: wilcoxon UNpaired analysis 
 # with correction