diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2020-11-13 10:27:00 +0100 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2020-11-13 10:27:00 +0100 |
| commit | 8d21c9966305627f9ea8599503fb3709ca916f9b (patch) | |
| tree | f17b0a995759860d153385c7d8f38fc1e7ac656c /R/plot.mixed.mmkin.R | |
| parent | f6fc5a050442c60d74a1943fa1181b8338322f30 (diff) | |
Fixed plotting of non-analytical solutions in saem
Diffstat (limited to 'R/plot.mixed.mmkin.R')
| -rw-r--r-- | R/plot.mixed.mmkin.R | 8 |
1 files changed, 3 insertions, 5 deletions
diff --git a/R/plot.mixed.mmkin.R b/R/plot.mixed.mmkin.R index 0d6b15c0..d8a1c2ac 100644 --- a/R/plot.mixed.mmkin.R +++ b/R/plot.mixed.mmkin.R @@ -75,8 +75,8 @@ plot.mixed.mmkin <- function(x, rownames(degparms_optim) <- ds_names degparms_optim_names <- setdiff(names(fit_1$par), names(fit_1$errparms)) colnames(degparms_optim) <- degparms_optim_names - residual_type = ifelse(standardized, "iwres", "ires") - residuals <- x$so@results@predictions[[residual_type]] + residual_type = ifelse(standardized, "residual", "standardized") + residuals <- x$data[[residual_type]] degparms_pop <- x$so@results@fixed.effects names(degparms_pop) <- degparms_optim_names } @@ -95,9 +95,7 @@ plot.mixed.mmkin <- function(x, odeini_names <- grep("_0$", degparms_all_names, value = TRUE) odeparms_names <- setdiff(degparms_all_names, odeini_names) - residual_type = ifelse(standardized, "iwres", "ires") - - observed <- cbind(x$data, + observed <- cbind(x$data[c("ds", "name", "time", "value")], residual = residuals) solution_type = fit_1$solution_type |