diff --git a/R/mc_sensitivity.R b/R/mc_sensitivity.R index 169cc185e31b6cd63b0a30ff8ed6262568b302ae..64e2de5eae5725c2b261328e87bc0d1a9938824f 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))