just added minor checks for numbers and stats

This commit is contained in:
Tanushree Tunstall 2022-11-19 16:46:34 +00:00
parent 7011dc75fa
commit 2f7090139e

View file

@ -19,6 +19,47 @@ getwd()
source("boxplot_stat_function.R") source("boxplot_stat_function.R")
source("plot_data_ob_paper.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 # Output
#============= #=============