From a5a73e8d76faff5983399386de571431e8824636 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Sat, 13 Sep 2025 09:08:43 +0200 Subject: Refactor for code readability --- R/saem.R | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) (limited to 'R') diff --git a/R/saem.R b/R/saem.R index 45a7f8e9..2a9c2fa5 100644 --- a/R/saem.R +++ b/R/saem.R @@ -871,12 +871,12 @@ parms.saem.mmkin <- function(object, ci = FALSE, covariates = NULL, ...) # estimates for the corresponding covariate values passed # in 'covariates'. covariate_centers <- object$covariate_centers - centered_covariates <- covariates + parm_covariates <- covariates if (!is.na(covariate_centers[1])) { # If covariates were centered by the saem function, we need to center # user specified covariates as well for (i in names(covariates)) { - centered_covariates[i] <- covariates[i] - object$covariate_centers[i] + parm_covariates[i] <- covariates[i] - object$covariate_centers[i] } } est_for_cov <- matrix(NA, @@ -890,7 +890,7 @@ parms.saem.mmkin <- function(object, ci = FALSE, covariates = NULL, ...) beta_degparm_name <- paste0("beta_", covariate, "(", deg_parm_name, ")") est_for_cov[deg_parm_name, ] <- estimate[deg_parm_name] + - centered_covariates[[covariate]] * estimate[beta_degparm_name] + parm_covariates[[covariate]] * estimate[beta_degparm_name] } else { est_for_cov[deg_parm_name, ] <- estimate[deg_parm_name] } -- cgit v1.2.3