tanu/mosaic_2020
1
0
Fork 0
Code Issues Pull requests Releases Wiki Activity

saving work

Browse source
This commit is contained in:
Tanushree Tunstall 2020-12-08 12:56:28 +00:00
parent 74c4dd878c
commit 9acc36c8de
4 changed files with 17 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

2
plot_data_ns.R
View file

@ -55,4 +55,4 @@ lf_fp_serum$obesity = as.factor(lf_fp_serum$obesity)
head(lf_fp_sam$value[lf_fp_serum$mediator == "vitd"])
########################################################################
rm(npa_lf, npa_wf, sam_lf, sam_wf, serum_lf, serum_wf, fp_adults_ics)
rm(npa_lf, npa_wf, sam_lf, sam_wf, serum_lf, serum_wf)

2
plot_data_s.R
View file

@ -55,4 +55,4 @@ lf_fp_serum$obesity = as.factor(lf_fp_serum$obesity)
head(lf_fp_sam$value[lf_fp_serum$mediator == "vitd"])
########################################################################
rm(npa_lf, npa_wf, sam_lf, sam_wf, serum_lf, serum_wf, fp_adults_ics)
rm(npa_lf, npa_wf, sam_lf, sam_wf, serum_lf, serum_wf)

6
read_data.R
View file

@ -23,6 +23,12 @@ ifelse(!dir.exists(outdir_stats_na), dir.create(outdir_stats_na), FALSE)
outdir_stats_ns = paste0(maindir, "output/stats/non_severe/")
ifelse(!dir.exists(outdir_stats_ns), dir.create(outdir_stats_ns), FALSE)
outdir_stats_s = paste0(maindir, "output/stats/severe/")
ifelse(!dir.exists(outdir_stats_s), dir.create(outdir_stats_s), FALSE)
outdir_stats_a = paste0(maindir, "output/stats/asthma/")
ifelse(!dir.exists(outdir_stats_a), dir.create(outdir_stats_a), FALSE)
outdir_plots = paste0(maindir, "output/plots/")
ifelse(!dir.exists(outdir_plots), dir.create(outdir_plots), FALSE)
########################################################################

10
stats_check.R
View file

@ -15,7 +15,15 @@ wf_serum = wf_fp_serum
# Numbers
table(fp_adults_ics$asthma, fp_adults_ics$obesity)
table(fp_adults_ics$obesity)
table(fp_adults_ics$obesity); table(fp_adults_ics$ia_cxr)
table(fp_adults_ics$obesity, fp_adults_ics$ia_cxr)
fisher.test(table(fp_adults_ics$obesity, fp_adults_ics$ia_cxr))
table(fp_adults_ics$obesity, fp_adults_ics$com_noasthma)
fisher.test(table(fp_adults_ics$obesity, fp_adults_ics$com_noasthma))
chisq.test(table(fp_adults_ics$obesity, fp_adults_ics$com_noasthma))
#=========