From c0699d16931159b8d83848b9cdcb919126de95e0 Mon Sep 17 00:00:00 2001
From: wbonat <wbonat@gmail.com>
Date: Tue, 10 Nov 2015 16:12:30 +0100
Subject: [PATCH] Tests sensitivity matrix
---
R/mc_sensitivity.R | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/R/mc_sensitivity.R b/R/mc_sensitivity.R
index 169cc18..64e2de5 100644
--- a/R/mc_sensitivity.R
+++ b/R/mc_sensitivity.R
@@ -13,9 +13,9 @@ mc_sensitivity <- function(product) {
Sensitivity1 <- matrix(0, n_par, n_par)
for (i in 1:n_par) {
for (j in 1:n_par) {
- # Sensitivity_temp[i,j] <- -sum(diag(product[[i]]%*%product[[j]]))
+ #Sensitivity_temp[i,j] <- -sum(diag(product[[i]]%*%product[[j]]))
Sensitivity[i, j] <- -sum(t(product[[i]]) * product[[j]])
- # Sensitivity1[i,j] <- -sum(product[[i]]*product[[j]])
+ #Sensitivity1[i,j] <- -sum(product[[i]]*product[[j]])
}
}
# print(forceSymmetric(Sensitivity)) print(forceSymmetric(Sensitivity_temp)) print(forceSymmetric(Sensitivity1))
--
GitLab