The goal of pmplots is to create exploratory and diagnostic plots commonly used in pharmacometrics.

Installation

You can install the released version of pmplots from MPN with:

mpn <- "https://mpn.metworx.com/snapshots/stable/2021-06-20"
install.packages("pmplots", repos = mpn)

This installs from a specific, dated snapshot. Please check https://mpn.metworx.com/docs/snapshots/ for the most recent snapshot date.

You can install the development version from GitHub with:

# install.packages("devtools")
devtools::install_github("metrumresearchgroup/pmplots")

Examples

library(pmplots)
library(dplyr)

data <- pmplots_data_obs()

The default pmplots behavior is to expect names following NONMEM convention, for example DV refers to observed values, PRED refers to population predicted values. With that convention, we can easily make observed versus predicted plots

dv_pred(data, yname = "MyDrug (ng/ml)")

dv_pred(data, loglog = TRUE)

NPDE plots

npde_time(data)

npde_pred(data)

QQ plots

cwres_q(data)

Exploratory plots

id <- pmplots_data_id()
cont_cat(id, x = "STUDYc", y = c("WT", "SCR", "AAG")) %>% pm_grid()

pairs_plot(id, y = c("WT//Weight (kg)", "SCR//Creat (mg/dL)", "BMI//BMI (kg/m2)"))

Some faceted versions

wrap_cont_time(data, y = c("RES","WRES", "CWRES", "NPDE"), ncol = 2)

Function and workflow overview

A functions and features (code + output) listing is available at pmplots_complete.md.

Suggested packages

Consider installing the cowplot package to help arranging plots on a page. Also, consider installing latex2exp to allow you to use latex in axis titles.