fixed #60

DESCRIPTION

@@ -1,6 +1,6 @@
 Package: SimDesign
 Title: Structure for Organizing Monte Carlo Simulation Designs
-Version: 2.18.3
+Version: 2.18.4
 Authors@R: c(person("Phil", "Chalmers", email = "[email protected]", role = c("aut", "cre"),
                     comment = c(ORCID="0000-0001-5332-2810")),
              person("Matthew", "Sigal", role = c("ctb")),

NEWS.md

@@ -1,5 +1,10 @@
 # NEWS file for SimDesign
 
+## Changes in SimDesign 2.19
+
+- Multiple parameter `bias()` returned incorrect behavior/errors
+  for `relative` and `abs_relative` types (see #60)
+
 ## Changes in SimDesign 2.18
 
 - Objects built by `createDesign()` gain `[]` and `rbind()` S3 functions for 

R/summary_functions.R

@@ -80,6 +80,8 @@
 #' # different parameter associated with each column
 #' mat <- cbind(M1=rnorm(1000, 2, sd = 0.25), M2 = rnorm(1000, 3, sd = .25))
 #' bias(mat, parameter = c(2,3))
+#' bias(mat, parameter = c(2,3), type='relative')
+#' bias(mat, parameter = c(2,3), type='standardized')
 #'
 #' # same, but with data.frame
 #' df <- data.frame(M1=rnorm(100, 2, sd = 0.5), M2 = rnorm(100, 2, sd = 1))
@@ -128,8 +130,8 @@ bias <- function(estimate, parameter = NULL, type = 'bias', abs = FALSE,
     if(!equal_len)
         stopifnot(ncol(estimate) == length(parameter))
     diff <- t(t(estimate) - parameter)
-    ret <- if(type == 'relative') colMeans(diff / parameter)
-        else if(type == 'abs_relative') colMeans(diff / abs(parameter))
+    ret <- if(type == 'relative') rowMeans(t(diff) / parameter)
+        else if(type == 'abs_relative') rowMeans(t(diff) / abs(parameter))
         else if(type == 'standardized') colMeans(diff) / colSDs(estimate)
         else colMeans(diff)
     if(abs) ret <- abs(ret)

SimDesign.Rproj

@@ -1,4 +1,5 @@
 Version: 1.0
+ProjectId: fe769c1e-d77f-4cb0-8cc9-a5d4a1fc8243
 
 RestoreWorkspace: No
 SaveWorkspace: No