-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathrun_aldex2.R
30 lines (26 loc) · 1005 Bytes
/
run_aldex2.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
library(ALDEx2)
run_aldex <- function(counts, cond, mc.samples, test){
#' Run all modules of ALDEx2 tool
#'
#'
#' @param counts The dataframe containing counts with samples as columns
#' @param cond A list with len(num_of_samples) with conditions
#' @param mc.samples Number of mc.samples
#' @param test Type of statistical test
#
results <- aldex(counts, cond, mc.samples=mc.samples, test=test, effect=TRUE,
include.sample.summary=FALSE, denom="all", verbose=FALSE, paired.test=FALSE)
return(results)
}
get_clr <- function(counts, cond, mc.samples,instance){
#' Run clr module of ALDEx2 tool
#'
#'
#' @param counts The dataframe containing counts with samples as columns
#' @param cond A list with len(num_of_samples) with conditions
#' @param mc.samples Number of mc.samples
#
results_clr <- aldex.clr(counts, cond, mc.samples=mc.samples,, denom="all", verbose=FALSE)
results <- getMonteCarloSample(results_clr,instance)
return(results)
}