From f13d1212e95a705c0212bbdf9fe7ee0f28430876 Mon Sep 17 00:00:00 2001
From: Tanushree Tunstall <tanu@tunstall.in>
Date: Mon, 26 Oct 2020 20:15:44 +0000
Subject: [PATCH] apparently modifed flu stats for npa: check!

---
 flu_stats_unpaired_npa.R | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

diff --git a/flu_stats_unpaired_npa.R b/flu_stats_unpaired_npa.R
index d439f10..c9448d0 100755
--- a/flu_stats_unpaired_npa.R
+++ b/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