From 12134a16a10daa4bad0cb0ce056ce49e24714c01 Mon Sep 17 00:00:00 2001 From: Christian Tischer Date: Thu, 14 Mar 2019 17:28:05 +0100 Subject: [PATCH] Add Bonferroni p-value adjustment --- functions.r | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/functions.r b/functions.r index d607370..d36e798 100644 --- a/functions.r +++ b/functions.r @@ -652,14 +652,16 @@ applyQC <- function(df, dfQC){ ##################################################### -################ From Tischi (below) ################ +################ From HTM-Explorer (below) ########## ##################################################### # # Normalization # -htmNormalization <- function(data, measurements, col_Experiment, transformation, gradient_correction, normalisation, negcontrol, col_QC, col_Well, col_Treatment, num_WellX, num_WellY) { +htmNormalization <- function( + data, measurements, col_Experiment, transformation, gradient_correction, normalisation, +negcontrol, col_QC, col_Well, col_Treatment, num_WellX = 0, num_WellY = 0) { echo("*") echo("* Data normalization") @@ -1140,6 +1142,7 @@ htmTreatmentSummary <- function(data, measurements, col_Experiment, col_Treatmen t_test__estimate = rep(NA,numEntries), t_test__p_value = rep(NA,numEntries), + t_test__p_value_adjusted = rep(NA,numEntries), t_test__signCode = rep(NA,numEntries),