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

added TASK at the beginning of the scripts

Browse source
This commit is contained in:
Tanushree Tunstall 2020-11-02 16:36:52 +00:00
parent 8033c1785f
commit a4d305d44b
3 changed files with 37 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

10
boxplot_stat_all.R
View file

@ -3,6 +3,16 @@ getwd()
setwd("~/git/mosaic_2020/")
getwd()
###############################################################
# TASK: boxplots with stats for NPA, SAM and SERUM
# for each mediator, at each timepoint by group
# TODO: check boxplot_stat_function.R
# output plots even if stats fail due to data containing only LLODs
# however for a given mediator, if there is error for ANY timepoint
# stats will not be added for the entire plot! To be handled in
# next iteration!
###############################################################
#=============
# Input
#=============

17
boxplot_stat_function.R
View file

@ -3,7 +3,22 @@ getwd()
setwd("~/git/mosaic_2020/")
getwd()
############################################################
# TASK: boxplots with stats
# TASK: Boxplots with stats
# for each mediator, at each timepoint by group
# Handles some error: i.e will continue plotting even if there
# is stats error. However, needs fixing see below point 1):
# FIXME:
# 1) For a given mediator, if there is error for ANY timepoint
# stats will not be added for the entire plot!
# 2) make the function more generic to take the following inputs
# df, column, group
# this way you can directly pass the log_calculated value column
# To be handled in
# next iteration!
############################################################
is.error <- function(x) inherits(x, "try-error")

12
boxplot_stat_log.R
View file

@ -3,11 +3,21 @@ getwd()
setwd("~/git/mosaic_2020/")
getwd()
###############################################################
# TASK: boxplots with stats for NPA, SAM and SERUM but with log scale
# for each mediator, at each timepoint, by group
# TODO: check boxplot_stat_function.R
# output plots even if stats fail due to data containing only LLODs
# however for a given mediator, if there is error for ANY timepoint
# stats will not be added for the entire plot! To be handled in
# next iteration!
###############################################################
#=============
# Input
#=============
source("boxplot_stat_function.R")
source("boxplot_function.R")
#source("boxplot_function.R") for sam
source("plot_data.R")
#=============