npc.coverage: Function to plot the coverage of the Numerical
Predictive Check
Description
This function takes the output from the npc command in Perl
Speaks NONMEM (PsN) and makes a coverage plot. A coverage plot for
the NPC looks at different prediction intervals (PIs) for each data
point and calculates the total number of data points in the data set
lying outside of these
PIs. The plot shows the relative amount of data points outside of
their PI compared to the expected amount at that PI. In addition a
confidence interval around these values are computed based on the
simulated data.
Usage
npc.coverage(npc.info="npc_results.csv",
main = "Default",
main.sub=NULL,
main.sub.cex=0.85,
...)
Arguments
npc.info
The results file from the npc command in
PsN. for example npc_results.csv, or if the file is in a
separate directory ./npc_dir1/npc_results.csv.
main
A string giving the plot title or NULL if
none. "Default" creates a default title.
main.sub
Used for names above each plot when using multiple
plots. Should be a vector c("Group 1","Group 2")
library(xpose4)
npc.coverage()
## to read files in a directory different than the current working directory npc.file <- "./another_directory/npc_results.csv"npc.coverage(npc.info=npc.file)