From 76d5304bcc612da673edf6f50290d5fa823e0976 Mon Sep 17 00:00:00 2001 From: Tobias Benz Date: Thu, 20 Sep 2018 13:10:53 +0200 Subject: [PATCH] Added ROPE to summary output --- R/analyze.stanreg.R | 18 +++++++++++++++--- 1 file changed, 15 insertions(+), 3 deletions(-) diff --git a/R/analyze.stanreg.R b/R/analyze.stanreg.R index c189d2c..f139e56 100644 --- a/R/analyze.stanreg.R +++ b/R/analyze.stanreg.R @@ -190,6 +190,7 @@ analyze.stanreg <- function(x, CI=90, index="overlap", ROPE_bounds=NULL, effsize Median_std = values$effects[[varname]]$std_median, MAD_std = values$effects[[varname]]$std_mad, MPE = values$effects[[varname]]$MPE, + ROPE = values$effects[[varname]]$ROPE, Overlap = values$effects[[varname]]$overlap ) ) @@ -198,6 +199,12 @@ analyze.stanreg <- function(x, CI=90, index="overlap", ROPE_bounds=NULL, effsize if (effsize == FALSE) { summary <- select_(summary, "-Median_std", "-MAD_std") } + + if (index == "ROPE") { + summary <- select_(summary, "-Overlap") + } else { + summary <- select_(summary, "-ROPE") + } # Text -------------------------------------------------------------------- # ------------------------------------------------------------------------- @@ -388,6 +395,7 @@ analyze.stanreg <- function(x, CI=90, index="overlap", ROPE_bounds=NULL, effsize values$MPE <- NA values$MPE_values <- NA values$overlap <- NA + values$ROPE <- NA values$adjusted_r_squared <- R2.adj # Text @@ -465,9 +473,10 @@ analyze.stanreg <- function(x, CI=90, index="overlap", ROPE_bounds=NULL, effsize values$MPE <- NA values$MPE_values <- NA values$overlap <- NA - - - + values$ROPE <- NA + + + # Text values$text <- paste0( " The intercept is at ", @@ -554,6 +563,9 @@ analyze.stanreg <- function(x, CI=90, index="overlap", ROPE_bounds=NULL, effsize rope <- rope(posterior, bounds=ROPE_bounds) values$ROPE_decision <- rope$rope_decision values$ROPE <- rope$rope_probability + }else{ + values$ROPE <- NA + values$ROPE_decision <- NA } if(index == "overlap"){