diff options
author | Johannes Ranke <johannes.ranke@jrwb.de> | 2025-09-13 09:08:43 +0200 |
---|---|---|
committer | Johannes Ranke <johannes.ranke@jrwb.de> | 2025-09-13 09:08:43 +0200 |
commit | a5a73e8d76faff5983399386de571431e8824636 (patch) | |
tree | 6b9c749096e5eebf39eb377eddbeba2280766fde /R/saem.R | |
parent | cbfd3aa2b4cf1481b0466f6398a0c2ef8d75e546 (diff) |
Diffstat (limited to 'R/saem.R')
-rw-r--r-- | R/saem.R | 6 |
1 files changed, 3 insertions, 3 deletions
@@ -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] } |