From 2f7090139ef14f04ac101106f4a3fd62b0d86ca3 Mon Sep 17 00:00:00 2001
From: Tanushree Tunstall <tanu@tunstall.in>
Date: Sat, 19 Nov 2022 16:46:34 +0000
Subject: [PATCH] just added minor checks for numbers and stats

---
 boxplots_stat_ob_paper.R | 41 ++++++++++++++++++++++++++++++++++++++++
 1 file changed, 41 insertions(+)

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
 #=============