ah, found it, added lf data number check to flu stats

This commit is contained in:
Tanushree Tunstall 2020-10-26 20:18:26 +00:00
parent 6da0bd2e00
commit a50cb6d331
2 changed files with 14 additions and 0 deletions

View file

@ -35,6 +35,13 @@ lf = sam_adults_lf[sam_adults_lf$flustat == 1,]
table(lf$timepoint)
length(unique(lf$mosaic))
lf$timepoint = paste0("t", lf$timepoint)
if (table(lf$flustat) == table(sam_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

View file

@ -35,6 +35,13 @@ lf = serum_adults_lf[serum_adults_lf$flustat == 1,]
table(lf$timepoint)
lf$timepoint = paste0("t", lf$timepoint)
if (table(lf$flustat) == table(serum_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
#######################################################################