diff --git a/flu_stats_unpaired_npa.R b/flu_stats_unpaired_npa.R
index f9ac727..d439f10 100755
--- a/flu_stats_unpaired_npa.R
+++ b/flu_stats_unpaired_npa.R
@@ -20,10 +20,12 @@ rm(colnames_sam_df, expected_rows_sam_lf
 
 rm(pivot_cols)
 
+my_sample_type = "npa"
 #=============
 # Output: unpaired analysis of time for npa
 #=============
-stats_time_unpaired_flu_npa = paste0(outdir_stats, "flu_stats_time_unpaired_npa.csv")
+outfile_name = paste0("flu_stats_time_unpaired_", my_sample_type, ".csv")
+flu_stats_time_unpaired = paste0(outdir_stats, outfile_name)
 
 #%%========================================================
 # data assignment for stats
@@ -316,5 +318,5 @@ colnames(combined_unpaired_stats_f)
 #******************
 # write output file
 #******************
-cat("UNpaired stats for groups will be:", stats_time_unpaired_flu_npa)
-write.csv(combined_unpaired_stats_f, stats_time_unpaired_flu_npa, row.names = FALSE)
+cat("UNpaired stats for groups will be:", flu_stats_time_unpaired)
+write.csv(combined_unpaired_stats_f, flu_stats_time_unpaired, row.names = FALSE)
diff --git a/flu_stats_unpaired_sam.R b/flu_stats_unpaired_sam.R
index 4af022e..1ed0664 100755
--- a/flu_stats_unpaired_sam.R
+++ b/flu_stats_unpaired_sam.R
@@ -20,10 +20,13 @@ rm(colnames_npa_df, expected_rows_npa_lf
 
 rm(pivot_cols)
 
+my_sample_type = "sam"
 #=============
 # Output: unpaired analysis of time for sam
 #=============
-stats_time_unpaired_flu_sam = paste0(outdir_stats, "flu_stats_time_unpaired_sam.csv")
+outfile_name = paste0("flu_stats_time_unpaired_", my_sample_type, ".csv")
+flu_stats_time_unpaired = paste0(outdir_stats, outfile_name)
+
 #%%========================================================
 # data assignment for stats
 wf = sam_df_adults_clean[sam_df_adults_clean$flustat == 1,]
@@ -32,8 +35,6 @@ lf = sam_adults_lf[sam_adults_lf$flustat == 1,]
 table(lf$timepoint)
 length(unique(lf$mosaic))
 lf$timepoint = paste0("t", lf$timepoint)
-my_sample_type = "sam"
-
 ########################################################################
 # Unpaired stats at each timepoint b/w groups: wilcoxon UNpaired analysis 
 # with correction
@@ -317,5 +318,5 @@ colnames(combined_unpaired_stats_f)
 #******************
 # write output file
 #******************
-cat("UNpaired stats for groups will be:", stats_time_unpaired_flu_sam)
-write.csv(combined_unpaired_stats_f, stats_time_unpaired_flu_sam, row.names = FALSE)
+cat("UNpaired stats for groups will be:", flu_stats_time_unpaired)
+write.csv(combined_unpaired_stats_f, flu_stats_time_unpaired, row.names = FALSE)
diff --git a/flu_stats_unpaired_serum.R b/flu_stats_unpaired_serum.R
index 296a4c2..29e2f6a 100755
--- a/flu_stats_unpaired_serum.R
+++ b/flu_stats_unpaired_serum.R
@@ -20,10 +20,13 @@ rm(colnames_sam_df, expected_rows_sam_lf
 
 rm(pivot_cols)
 
+my_sample_type = "serum"
 #=============
 # Output: unpaired analysis of time for serum
 #=============
-stats_time_unpaired_flu_serum = paste0(outdir_stats, "flu_stats_time_unpaired_serum.csv")
+outfile_name = paste0("flu_stats_time_unpaired_", my_sample_type, ".csv")
+flu_stats_time_unpaired = paste0(outdir_stats, outfile_name)
+
 #%%========================================================
 # data assignment for stats
 wf = serum_df_adults_clean[serum_df_adults_clean$flustat == 1,]
@@ -32,8 +35,6 @@ lf = serum_adults_lf[serum_adults_lf$flustat == 1,]
 table(lf$timepoint)
 lf$timepoint = paste0("t", lf$timepoint)
 
-my_sample_type = "serum"
-
 ########################################################################
 # Unpaired stats at each timepoint b/w groups: wilcoxon UNpaired analysis with correction
 #######################################################################
@@ -315,5 +316,5 @@ colnames(combined_unpaired_stats_f)
 #******************
 # write output file
 #******************
-cat("UNpaired stats for groups will be:", stats_time_unpaired_flu_serum)
-write.csv(combined_unpaired_stats_f, stats_time_unpaired_flu_serum, row.names = FALSE)
+cat("UNpaired stats for groups will be:", flu_stats_time_unpaired)
+write.csv(combined_unpaired_stats_f, flu_stats_time_unpaired, row.names = FALSE)