basic.GOF4 {nonmem2R}

Basic 4- and 6-panel GOF

Description

basic.GOF4 provides a 4-panel GOF plot showing
i) Observations (DV) vs population predictions (PRED)
ii) Observations (DV) vs individual predictions (IPRED)
iii) sqrt(abs(CWRES)) vs individual predictions (IPRED), and
iV) CWRES vs TIME or other columns as set by idv

basic.GOF6 provides 2 additional GOF's showing histogram and qqnorm GOF's for CWRES.

Caption is added as reurned by get.caption

Usage

basic.GOF4(
  data,
  residual = "CWRES",
  idv1 = "IPRED",
  idv2 = "TIME",
  title = "",
  color = "",
  log.scale = FALSE,
  global.ggplot.options = NULL,
  refline = TRUE,
  add.loess = TRUE,
  control = GOF.control()
)

basic.GOF6(
  data,
  residual = "CWRES",
  idv1 = "IPRED",
  idv2 = "TIME",
  title = "",
  color = "",
  log.scale = FALSE,
  global.ggplot.options = NULL,
  refline = TRUE,
  add.loess = TRUE,
  bins = NA,
  control = GOF.control()
)

Arguments

data

data.frame to plot

residual

column name for residuals, default="CWRES"

idv1

independent variable for plot of sqrt of absolute CWRES

idv2

independent variable for plot of CWRES

title

title

color

data columns to set different colors in plot, interpreted as factor

log.scale

use log scale for DV, IPRED and PRED (TRUE) or normal scale (FALSE)

global.ggplot.options

ggplot option added to each ggplot before plotting/returning object

refline

add reference line (default =TRUE) or not (FALSE)

add.loess

add loess smoother to plot (TRUE), or not (FALSE)

control

an optional list of control settings. See GOF.control for the names of the settable control values and their effect.

bins

number of histogram bins, if NA, bins are set using Sturges' formula

Examples

# Get path to the example files included in nonmem2R package
file1 <- system.file("extdata", "sdtab999", package = "nonmem2R")
sdtab<-read.table(file=file1,skip=1,header=TRUE)
set.script.name("MyScript.R")
## Example 4 panel basic GOF
basic.GOF4(subset(sdtab,DV>0),idv2="TAPD")
## Example 6 panel basic GOF
basic.GOF6(subset(sdtab,DV>0),idv1="TAPD",idv2="PRED")