R/000.R

In nlmixr: Nonlinear Mixed Effects Models in Population PK/PD

#' GOF plots for nlme-based mixed-effect models
#'
#' Generates basic goodness-of-fit plots for nlme-based mixed-effect models
#'
#' @param fit nlme fit object
#' @param ... optional additional arguments
#' @return nothing, displays plots
#' @export
nlme_gof <- function(fit, ...) {
  RxODE::rxReq("lattice")
  df <- getData(fit)
  df <- rbind(
    cbind(df[, c("ID", "TIME", "DV")], grp = 0),
    cbind(df[, c("ID", "TIME")], DV = fit$fitted[, 1], grp = 1),
    cbind(df[, c("ID", "TIME")], DV = fit$fitted[, 2], grp = 2)
  )
  p <- lattice::xyplot(DV ~ TIME | ID, df, group = grp, type = "b", lwd = c(NA, 1, 1), pch = c(1, NA, NA), col = lattice::trellis.par.get("superpose.line")$col[c(1, 1, 2)], ...)
  p
}