Merge branch 'hotfix/0.7.1' into devel

DESCRIPTION

@@ -1,6 +1,6 @@
 Package: heemod
 Title: Models for Health Economic Evaluation
-Version: 0.7.0
+Version: 0.7.1
 Authors@R: c(
     person("Antoine", "Filipovic-Pierucci", email = "[email protected]", role = c("aut", "cre")),
     person("Kevin", "Zarca", email = "[email protected]", role = "aut"),

NEWS.md

@@ -1,3 +1,9 @@
+# heemod 0.7.1
+
+## Bugfixes
+
+  * Fixed an error resulting in incorrect covariance analysis results when relations between values and parameters were negative.
+
 # heemod 0.7.0
 
 ## Breaking changes

R/covariance.R

@@ -30,14 +30,16 @@ compute_cov <- function(psa) {
     ))
     res <- stats::lm(form, data = data_trans)
 
-    tot <- sum(stats::coef(res)[-1])
+    val <- abs(stats::coef(res)[-1])
+    tot <- sum(val)
     r2 <- summary(res)$r.squared
     if (r2 < .99) {
       warning(sprintf(
-        "Only %.0f%% of variance explained, results may be inaccurate."
+        "Only %.0f%% of variance explained, results may be inaccurate.",
+        r2 * 100
       ))
     }
-    as.data.frame(as.list(stats::coef(res)[-1] / tot * r2))
+    as.data.frame(as.list(val / tot * r2))
   }
 
   psa$psa %>% 

inst/examples/example_update.R

@@ -44,7 +44,7 @@ res <- run_model(
 
 # generating table with new parameter sets
 new_tab <- data.frame(
-  age_init = 40:80
+  age_init = 40:45
 )
 
 # with run_model result
@@ -55,7 +55,7 @@ summary(ndt)
 # using weights
 
 new_tab2 <- data.frame(
-  age_init = 40:80,
+  age_init = 40:45,
   .weights = runif(41)
 )
 ndt2 <- update(res, newdata = new_tab2)