diff --git a/boxplots_stat_ob_paper.R b/boxplots_stat_ob_paper.R
index 08f073b..98884c7 100644
--- a/boxplots_stat_ob_paper.R
+++ b/boxplots_stat_ob_paper.R
@@ -19,6 +19,47 @@ getwd()
 source("boxplot_stat_function.R")
 source("plot_data_ob_paper.R")
 
+
+
+# numbers
+table(vl_lf_plot_df$mosaic)
+length(unique(vl_lf_plot_df$mosaic))
+table(vl_lf_plot_df$obesity, vl_lf_plot_df$timepoint)
+table(vl_lf_plot_df$obesity, vl_lf_plot_df$timepoint, vl_lf_plot_df$mediator)
+
+length(unique(npa_lf_plot_df$mosaic))
+table(npa_lf_plot_df$obesity, npa_lf_plot_df$timepoint, npa_lf_plot_df$mediator)
+table(npa_lf_plot_df$mediator, npa_lf_plot_df$obesity,  npa_lf_plot_df$timepoint)
+
+ob_grp_npa = npa_wf[npa_wf$obesity==1,]
+nob_grp_npa = npa_wf[npa_wf$obesity==0,]
+
+wilcox.test(ob_grp_npa$il1_npa1,
+            nob_grp_npa$il1_npa1, paired = F)
+wilcox.test(ob_grp_npa$il1_npa2,
+            nob_grp_npa$il1_npa2, paired = F)
+
+wilcox.test(ob_grp_npa$il6_npa1,
+            nob_grp_npa$il6_npa1, paired = F)
+wilcox.test(ob_grp_npa$il6_npa2,
+            nob_grp_npa$il6_npa2, paired = F)
+
+wilcox.test(ob_grp_npa$il8_npa1,
+            nob_grp_npa$il8_npa1, paired = F)
+wilcox.test(ob_grp_npa$il8_npa2,
+            nob_grp_npa$il8_npa2, paired = F)
+
+wilcox.test(ob_grp_npa$ip10_npa1,
+            nob_grp_npa$ip10_npa1, paired = F)
+wilcox.test(ob_grp_npa$ip10_npa2,
+            nob_grp_npa$ip10_npa2, paired = F)
+
+length(unique(serum_lf_plot_df$mosaic))
+table(serum_lf_plot_df$obesity, serum_lf_plot_df$timepoint, serum_lf_plot_df$mediator)
+table(serum_lf_plot_df$mediator, serum_lf_plot_df$obesity,  serum_lf_plot_df$timepoint)
+
+
+
 #=============
 # Output
 #=============