Histogram and/or QQ-norm GOF for ETA's

Provides histogram and/or QQ-norm GOF for all ETA's included in the input data.frame with caption added as reurned by get.caption

basic.eta.GOF(
  data,
  title = "",
  global.ggplot.options = NULL,
  type = c("both", "qqnorm", "hist"),
  ETA.subset = NULL,
  refline = TRUE,
  drop.fixed = TRUE,
  id.column = "ID",
  standardize = TRUE,
  bins = NA,
  control = GOF.control()
)

eta.hist.GOF(
  data,
  title = "",
  drop.fixed = TRUE,
  refline = TRUE,
  id.column = "ID",
  standardize = TRUE,
  bins = NA,
  control = GOF.control()
)

eta.qqnorm.GOF(
  data,
  title = "",
  drop.fixed = TRUE,
  refline = TRUE,
  id.column = "ID",
  standardize = TRUE,
  control = GOF.control()
)

Arguments

data

data.frame to plot

title

title

global.ggplot.options

ggplot option added to each ggplot before plotting/returning object

type

do only qq-norm plots (qqnorm), only histograms (hist), or do both (both)

ETA.subset

index for subset of ETA's to plot. If NULL (default) all ETA's are plotted.

refline

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

drop.fixed

drop ETA's that are fixed, i.e. any ETA with the same value in all subjects

id.column

column name that indicate subject identifier

standardize

Standardize all ETA's (e.g. scale to unit variance)

bins

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

control

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

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")
basic.eta.GOF(sdtab)