Nonmem Software
The software consists of three parts. The NONMEM® program The NONMEM® program itself is a very general (non-interactive) model analysis program that can be used to. Pharmacokinetic Software. Integrates with NONMEM and other existing software to expedite population modeling and analysis. For Windows systems and requires.
Background Started in the early 1980s with the development of the NONMEM (acronym based on “NON-linear Mixed-Effects Modeling”) software, “population analysis” has proven to be extremely useful within pharmacometrics, both in the development of new drugs and the improvement of therapy with approved drugs. Development of the NONMEM software continues, and although over time several other modeling software tools have become available, NONMEM is still regarded as the gold standard within the pharmacometric community: a recent survey identified NONMEM (together with PsN) as the most frequently used software tool by far.
Modeling and simulation (M&S) in clinical pharmacology, and the use of NONMEM in particular, however, has a steep learning curve for most starting researchers. This is partially because of the fact that NONMEM is invoked from the command line, and models are implemented using a Fortran-derived syntax (NM-TRAN).
In addition, at its core, NONMEM performs only model estimation (or simulation), and the implementation of essential diagnostic tools such as bootstrap analyses and the creation of goodness-of-fit plots are left to the modeler. Therefore, alongside the development of NONMEM, many third-party tools have been developed that facilitate the use of NONMEM by providing tools for organization and automation. In this tutorial, we will demonstrate the use of some of the most widely used auxiliary software tools: PsN, Xpose, and Pirana. All three tools are released under an open-source license and are freely available (except for the commercial use of Pirana, for which a commercial license is required). Separately, each tool offers useful functionalities, but there is a synergistic benefit as well when used together. The aim of this tutorial is to show how these three tools provide a comprehensive workbench for M&S. This will be performed by showing examples of the most often encountered steps in pharmacokinetic (PK) model development, i.e., a covariate modeling procedure and model evaluation using residual- and simulation-based diagnostics.
However, the aim of the tutorial is not to provide guidance on how to perform a population PK analysis but strictly on these software tools. A detailed overview of important aspects in population PK analyses has recently been presented in this journal, and we refer the reader to that article for guidance. The current tutorial is structured in three parts: first, a brief introduction to these software tools is given explaining their basic purpose. In the main part of the tutorial, we show how a typical PK model-building analysis is performed with these tools and NONMEM. At the end of the tutorial, several interesting additional features are highlighted for each specific tool. We will assume the reader is already somewhat familiar with NONMEM, although we have made efforts to present and discuss the tools in a general manner wherever possible. All of the presented models, data sets, output, and diagnostic plots are available in the Supplementary Table S1 online.
Xpose Xpose is a tool for plotting and analyzing NONMEM output, developed as a module for the R software ( open-source, S-based statistical software). The tools in Xpose can be used from the R command line or from a text-based menu system in R.
Xpose was first released in 1998 and was initially developed for S-Plus. The current version (main version number 4) is, however, released exclusively for R and builds on the lattice module for plotting.
Both PsN and Xpose are developed at the Uppsala University and are released under an open-source license (GNU v2). Pirana Pirana is a graphical user interface for NONMEM, PsN, and Xpose/R. It has functionality for model management, model execution, output generation, interpretation of results, and includes many other tools. Development of Pirana started in 2007 at the Netherlands Cancer Institute/Slotervaart Hospital (Amsterdam, The Netherlands) and is currently continued by Pirana Software & Consulting BV.
Pirana is released under an open-source license (Creative Commons) for academic users as well as a commercial license. Tutorial In this tutorial, file and folder names are shown in italic, Pirana actions are shown red-italic, while commands, arguments, NONMEM syntax, and screen output are shown in fixed-width font. We will step through an example model-building exercise, with the intent of showing how to create, manage, and evaluate PK models for a given data set.
The PK data set used in this tutorial ( pktab1) is available in the Supplementary Table S1 online and contains plasma concentration data obtained from a simulated clinical trial of a novel i.v. Drug, performed in 50 patients, measured at 8 time points. All model files that are mentioned in this article are also available online and can be used as reference. Make sure that all software is installed correctly, and NONMEM runs can be started from both PsN and Pirana (visit the respective websites for installation instructions) and that the Xpose4 package is installed in R.
Create a folder for this analysis somewhere on your hard-drive, and put pktab1 in this folder. Browse to this folder in Pirana and save it as the project “TutorialPSP” (button 4 in ). As a starting point for NONMEM models, it is often easiest to use the PK model wizard in Pirana or start from one of the models available in the model library.
Start the wizard dialog window in Pirana ( Tools → Wizards), choose the PK model wizard, and explore the options. However, for this tutorial, we have already provided the first model ( run1.mod): copy the model run1.mod from the Supplementary Table S1 online to the folder you created. This model should now be visible in the Pirana main model list when you direct Pirana to the right folder (if not, refresh Pirana; button 5 in ).
If you run from a model created by the wizard, make sure that the columns in the data set match up exactly with the records specified in the $INPUT record in the model file. For the provided run1.mod, this has already been done. First runs If we would invoke the native NONMEM run script, we would run, e.g., nmfe72 run1.mod run1.lst, but Pirana automates this in a dialog window: ( select model → right click → NONMEM → nmfe). Several options are available to configure how and where (local or on a cluster) to run the model, or whether to submit it to a job scheduler. If you choose to run your models in Pirana through nmfe, you need to configure a NONMEM installation in Pirana ( Settings → NONMEM). Select the quick-find option to scan your hard-drive for common installation locations for NONMEM.

Nonmem Software Pharmacokinetics
If NONMEM is not found there yet, specify the location yourself. In this tutorial, we will, however, not use nmfe, but only PsN to run models.
There are multiple benefits of using PsN's execute over the regular nmfe command, the main ones being that models are run in separate folders, runs can be restarted automatically upon crashes and unsuccessful minimizations, and initial estimates can be tweaked. Select the model and right click → PsN → execute. Pirana will show a similar dialog window as shown for nmfe, but now the command line for PsN's execute tool is shown ( ).
For our first run, we will use the default command: execute run1.mod. After clicking the run button, if PsN and Pirana are set up correctly, a console window will open with the following message. PsN dialog window. Starting 1 NONMEM executions. 1 in parallel. S:1., which indicates that one model estimation has been started.
The actual NONMEM process will run in a subfolder created by PsN. By default, this will be /modelfitdir1, but a custom folder name can be specified using the -dir argument. If the argument -modeldirname is specified, PsN will automatically choose an informative folder name, in this case /run1.mod.dir.1. Check that inside the folder created by PsN, a subfolder is created in which the NONMEM process actually runs ( /NMrun1). In Pirana, PsN folders are hidden by default, as they may quickly clutter the model overview. To unhide them again, select either “PsN folders” or “All folders” from the folder selector (button 7 in ).
Many additional arguments can be supplied to PsN commands to customize how PsN runs the NONMEM model, e.g.: execute run1.mod -dir=test1 –retries=3 -parafile=mpi6.pnm –nmversion=nm72, which will run model run1.mod in the folder /test1. PsN is also instructed now to run using NONMEM version 7.2, parallelized using MPI (specified in configuration file mpi6.pnm) and to retry three times if the model does not converge in the initial run, each time changing the initial estimates of the parameters.
Hiring staff members is yet another interesting addition that gives you the possibility to improve your club's rating by upgrading your employees skills to higher levels. Another aspect that influenced Career mode was the introduction of Team Chemistry that dictated how well your squad plays together as a result of happiness, morale, and the number of matches played together. Free download games fifa 2006 for pc.
A list of all available arguments including help information is available from the PsN dialog window in Pirana ( ), which is similar to using the –h and –help arguments on the command line. In the PsN configuration file ( psn.conf), a list of default arguments can be supplied as well, so commonly used arguments do not have to be repeated on the command line. Results evaluation After estimation has been completed, refresh the Pirana main overview again. Basic results are now shown for the run, such as the objective function value (OFV), whether minimization was successful (S), and whether the covariance step was successful (C). If the run was successful, PsN will automatically copy back the results file ( run1.lst) and the produced output tables ( sdtab1, patab1) to the parent folder, which are then parsed by Pirana upon refreshing.
Select the model again, and now select the Estimates tab from the right section of the Pirana window (11 in ). The right section will then show the parameter estimates for this run. A more detailed overview that will, e.g., show full random effect blocks can be opened from here.
Create a run summary in Word format (alternatives are HTML, plain text, or LaTeX), from the menu → Results → Run Reports → Word. Diagnostic plots are a main decision tool in pharmacometric model development, in which a distinction can be made between residual-based diagnostics (goodness-of-fit plots based on residuals plotted vs. Time or model predictions) and simulation-based diagnostics (e.g., VPC). Obviously, no strict rule should be applied here, but in our experience, residual-based diagnostics are more useful in the first stages of model development, whereas simulation-based diagnostics are more useful in the latter stages of model development. Depending on what part of the model is diagnosed, different diagnostic plots will be useful, and most likely, no single plot will suffice to evaluate model fit.
Suggestions for diagnostic plots for specific model parts are shown in. Note that this list is not exhaustive, and in addition, more specific plots may be required for appropriate diagnosis. Suggested diagnostic plots for PK model development For the current analysis, we will first have a look at some basic, residual-based goodness-of-fit plots. A general model diagnostic is the plot of weighted residuals vs. Model predictions. For models run with the conditional estimation method (FOCE), it is most appropriate to use conditional weighted residuals (CWRES). Use the Xpose graphical user interface (GUI) within Pirana to create these plots ( ): select model → right click → Xpose → Xpose GUI.
Add the plot cwres.vs.pred and cwres.vs.idv to the list of plots, and possibly some additional goodness-of-fit plots (a list of useful Xpose plots is given in ). In these plots, you should observe an obvious pattern in the residuals. It may not be obvious immediately what is the cause of the misspecification, but a look at the individual plots ( ind.plots) probably gives more insight: the model seems to be missing the bi-phasic nature of the observed data. Therefore, in the next step, we will add distribution of drug into a peripheral compartment to our PK model, but first, we will briefly discuss how to keep track of model development.
DataInspector: a quick and flexible method for creating diagnostic plots based on NONMEM output, or checking input data sets. Select the DataInspector for run1.mod ( right click → Model → Open DataInspector). Pirana will ask you now which input- or output file to open, choose the file sdtab1. In the DataInspector window, the variables on the x and y axes can be easily changed to any variable in the data set, and filters can be applied on patient ID, or on any other variable. From this window, R-code can be generated that will recreate the same graph in R. R-script library: Pirana comes bundled with a library of diagnostic R scripts, which can be run automatically from Pirana. Select a model and select one of the R scripts from the Scripts Tab, e.g., Basic GOF → DV vs IPRED, and then select Run script from the buttons above the list.
The requested plot will be created and opened in a pdf document. Created plots will be listed in the “Reports” tab (10 in ). The standard scripts bundled with Pirana can be extended and customized easily, but that is beyond the scope of this tutorial. Please refer to the manual and default scripts for guidance.
Xpose menu system in R: The Xpose menu system can be started from Pirana ( right click → Xpose → Start Xpose menu in R), or manually from an R session by invoking: library(xpose4) xpose4 Xpose will then ask for the run number, which is used to determine which tables to import. Note that to use Xpose, NONMEM needs to be instructed to output tables that adhere to a specific naming and structure ( Supplementary Table S1 online, or the Xpose Web site). If you did not use the Pirana model wizard to create the tables, add them manually in the NONMEM model file, using: $TABLE ID TIME IPRED IWRES EVID MDV NOPRINT ONEHEADER FILE=sdtab1 From the Xpose menu system, goodness-of-fit plots can be created, e.g., browse to “4: Goodness-of-fit plots” and then “2: Basic goodness-of-fit plots.”. Model management in Pirana In M&S projects in both academia and industry, it is essential to keep track of the model development history. Pirana and PsN offer a tool to aid the modeler with this through the run record.
The run record entails a standardized way of keeping track of model meta-information such as the “parent” model, a description, and/or a label for the model, and other information about the model components. We consider it good modeling practice to create a new model file for every change that is made to a model. Duplicate the first model in Pirana ( right click on run1.mod actions → File actions → duplicate).
In the dialog that will open up, run2.mod is suggested as new model file name, and we can select which model is the reference model for this model ( run1.mod), and optionally update the initial estimates with the parameter estimates from the reference model. The model description and other meta-information are stored in the model file itself using the run record syntax (also see ). Pirana can show the main model overview as a list or as a tree (button 6 in ) and can export a copy of the run record to various formats including comma separated (csv) files, HTML files, and Word documents ( Results → Run records), or create an interactive visualization of the model development tree or “visual run record” such as that shown in ( Results → Run records → Visual run record). Visual run record example. Green indicates a significant improvement in fit, red indicates a worsening, yellow indicates approximate identical fit (ΔOFV.
Covariate model The data set pktab1 contains three continuous covariates: weight (WT), creatinine clearance (CRCL), AGE, and two binary covariates: SEX and study center (CENT). Because the number of covariates is low (5), one might choose to test these covariates manually. However, in the next step of this tutorial, we will perform stepwise covariate modeling (scm) in PsN, to demonstrate this automated method. In the first phase of the scm, covariates are added to the base model in a stepwise manner based on statistical significance (decrease in OFV). In a second “backward elimination” phase, covariates are then removed from the final model obtained in the first phase, if removal of the covariate does not result in significantly worse fit. The backward step is typically done at a higher significance level.
Like all PsN tools, the scm can be run from the command line or from Pirana, in a similar way as done before for execute. However, the scm tool requires that a configuration file is created first. Several examples of configuration files are available on the PsN Web site, but here we will create one using the wizard in Pirana ( Tools → Wizards → scm configuration file). Create a configuration file, with the filename psp.scm, in which all covariates are tested on both CL and V1, at a significance level of 0.1 (forward step) and 0.05 (backward step) and restrict the analysis to linear relationships (set the relationships argument to “1,2”). Check that the correct covariates are included; the ones shown in the wizard are only placeholders. If you are not sure on the correct syntax of the scm file, please have a look at the provided annotated configuration file ( psp.scm in the Supplementary Data online).
Before starting the scm, duplicate run3.mod to run4.mod, and remove any $TABLE records that are present. Start the scm from the PsN dialog window using this command: scm psp.scm –model=run4.mod After the scm has completed, PsN will compile an output file called scmlog.txt in the scm folder. Open this file and interpret the information, it holds the results for all tests performed in the scm and shows how the final covariate model was constructed.
You should find that four covariate relationships are significant: CRCL on CL, WT on CL, WT on V, and SEX on V. This will be our “final model” from the covariate model-building step. If you did not find these results, compare your results with those provided in the Supplementary Data online. More elaborate covariate model-building techniques based on the scm are available in PsN as well, such as the cross-validated scm (xvscm), and the bootstrap of the scm (bootscm). These tools provide more details on the predictive ability of the model, and on selection bias and expected covariate model size, respectively. They can also be run on linearized versions of the model, to reduce the computational workload. Other covariate modeling tools that are available in PsN include the lasso and the fixed random effects model (FREM), whereas in Xpose, the gam and bootstrap of the gam ( bootgam) are implemented.
All of these methods have certain advantages and drawbacks, which are however beyond the scope of this article. Final model If the scm completed successfully, you will find a folder called finalmodels inside the scm folder created by PsN. Inside that folder, locate the final model with included covariate relations and copy that to the main model folder, renaming it to run5.mod. We will use this model to perform some more elaborate diagnostic evaluation, i.e., we will get bootstrap estimates for our parameters and create a VPC. Therefore, first open the PsN bootstrap dialog window ( Right click → PsN → Model diagnostic → bootstrap). A required argument to the bootstrap command is the –samples argument, which specifies how many bootstrap samples should be drawn. More samples will make the bootstrap estimates and percentile estimates more precise; however, at the expense of increased computation times.

The recommended number of samples for an analysis will depend on what statistic is desired by the modeler. For obtaining bootstrap means of the parameter estimates, 100 samples may be enough, but to obtain accurate 5th and 95th percentiles, 1,000 or more are required to obtain precise estimates.
After the bootstrap is finished (on a single core, a bootstrap of 200 samples may easily take an hour for this model), the results from the bootstrap are available in the bootstrap folder, in the files rawresultsrun5.csv and bootstrapresults.csv. Open these files and interpret the output: are the parameters estimated with adequate precision? In Pirana, select the bootstrap folder created by PsN and run the R script to create plots of parameter distributions ( Scripts tab → PsN → Bootstrap results), which facilitates evaluation and interpretation of these bootstrap results. Earlier in this tutorial, we presented how basic goodness-of-fit plots can be created in Pirana and Xpose. Now, we will use all three tools to create a VPC of the final model ( run5.mod). A VPC is an internal validation tool, i.e., it shows how well the model predicts the data on which the model was conditioned.
It visualizes the median model prediction, the variability between patients and within patients (5th and 95th percentiles), and the uncertainty around the predicted percentiles. Before creating a VPC, we would like to stress a few points. First, in VPCs, it is important to show the uncertainty around the model-predicted percentiles, instead of presenting only lines indicating the point estimates of the percentiles. Signs of model improvement are when the confidence areas of the simulated percentiles encompass more of the observation percentiles, but also when the confidence intervals themselves are shrinking (while still encompassing the observation percentiles). Furthermore, because the VPC shows summary statistics (percentiles) of observations and predictions for binned data, the binning method is an important consideration.
To avoid inflation of the uncertainty around the summary statistics, each bin should contain an adequate amount of observations. A quick rule of thumb is that there should not be more bins than the number of observations per individual. PsN currently incorporates several binning approaches: manual bin specification (-binarray=x,y,z,etc), binning into an equal number of data points per bin (-binbycount=1 –noofbins=#), or spacing the bins equally over the x-axis (-binbycount=0, -noofbins=#). Of note, binning across a large variability in dose and/or influential covariates can diminish the diagnostic value of a VPC. One approach to overcome this is to stratify the VPC by the covariate, but this is not always possible, e.g., due to limited data per stratum. Moreover, when data arises from studies with adaptive designs (e.g., dose adjustments), VPCs are misleading unless the adaptation protocol is included in the model.
The prediction-corrected VPC offers a solution to these problems while retaining the visual interpretation of the traditional VPC. Prediction-corrected VPCs can be constructed by adding the argument –predcorr to the PsN VPC command. VPCs can be made in Pirana using the PsN dialog window ( right click → PsN → model evaluation → VPC), or from the command line. As discussed above, essential arguments to pass to the VPC are the number of samples and the binning method. Increasing the number of samples will increase the accuracy of the summary statistics for the simulation and their uncertainty interval (but it will not decrease the uncertainty interval itself). Start the VPC tool using vpc –samples=500 –noofbins=8 –binbycount=1 run5.mod The first of the two NONMEM runs that are started will only output a table with the necessary (observation) data for the VPC, it does not perform parameter estimation. The other model runs repeated Monte Carlo simulations of the model, using the same data set design as the original.
After the two NONMEM runs have finished, PsN will process the output, bin the simulated and observed data, calculate the percentiles and confidence intervals for each bin, and export a csv file with the summary statistics. This csv-file can then be processed by Xpose and turned into a VPC-plot, which can be automated from Pirana using the Xpose GUI as described before, or using the R scripts library. Both the approaches will create and open a pdf file with the VPC (see example in ). Using the latter approach, this is done as follows: select the VPC-folder created by PsN, and open the Scripts tab on the right. In the list of R scripts, choose Xpose → VPC → Basiclogy.R. A pdf file with the VPC will be generated and opened from Pirana. As an exercise, try to create the same plot using the Xpose GUI in Pirana as well: select the model in the main overview and Right click → Xpose → Xpose GUI.
Visual predictive check. A visual predictive check for run3.mod. The solid line connects the observed median values per bin, whereas the dashed lines represent the observed 5th and 95th percentile of the observations.
Blue areas indicate the 95% CIs of. Because only limited space is available for this tutorial, we will not demonstrate other diagnostics. We will however create a run record of the model development: in the menu bar, click Results → Run records → Detailed run record (csv). This will create and open a csv file containing run numbers, descriptions, other meta-information, and run results. If the file does not open automatically, please look it up in the “Files” tab on the right and double click it. As a final step in our PK model development, select the final model again and bundle that model file, the associated result files and output files and the VPC folder into a zip file ( Right click → File actions→ Compress to zip-file).
Moreover, create a bundled report of the goodness-of-fit plots for the final model (in the Reports tab, select the desired pdf's and click Bundle into → ). Additional features For all three software tools presented here, we have only scratched the surface of possibilities. The reader is encouraged to explore more advanced functionality using the documentation available on the respective websites. We will highlight a few examples of more advanced features for the three programs below. The simulation and re-estimation (sse) tool can be used to evaluate trial designs, e.g., to evaluate whether model parameters can be estimated with adequate precision under the intended or alternative experimental designs or with alternative models. It is also useful for evaluating fundamental modeling aspects, e.g., to evaluate how model diagnostics perform under given designs. Another interesting tool for design evaluation, specifically for the calculation of study power and number of study subjects required, is Monte Carlo Mapped Power (mcmp).
In this tool, which is also based on simulation and re-estimation, the individual OFVs calculated in NONMEM are exploited to evaluate designs in a more rapid way than using sse. Case deletion diagnostics (cdd) is a tool that can be used in the final stages of model development to investigate whether model fit depends more heavily on particular strata of the data set. The most important PsN functions are listed in. Overview of commonly used PsN tools Xpose. In this tutorial, we showed how to create goodness-of-fit plots from the Xpose menu system, and from the Xpose interface in Pirana (VPC). However, the most efficient way to use Xpose, especially when doing repetitive jobs, is to use scripts to create the desired goodness-of-fit plots. A list of useful Xpose functions is given in, whereas a more complete overview of functions and example scripts (“bestiarium”) is available on the Xpose website.
The default plots in Xpose can be extended and modified in various ways. First, most Xpose functions take arguments that alter the implementation of the plot. The looks of the plots can also be changed by directly passing lattice arguments to the Xpose function. Second, the input variables can be altered easily to allow creation of nondefault plots.
For example, to create a plot of a covariate METAB vs. IDV instead of DV vs. IDV, we can “trick” Xpose into using METAB as DV while maintaining the plot characteristics: change.xvardef(xpdb, “dv”) xplot print(xplot) Finally, since Xpose is an open-source R module, it is possible to modify the Xpose functions directly, if further modifications are required. A model translator is included that can translate any NONMEM model (written in an NM-TRAN ADVAN subroutine) to differential equation code for R (using deSolve library), Berkeley Madonna (software dedicated to ODE simulations), or Matlab/PopED. Note that many functions in Pirana can be extended and customized by the user, such as the R-scripts library and the wizards. Pirana can also be used for modeling on remote clusters (through ssh-connections) and has native support for SGE, Torque, and Condor job schedulers.
NONMEM stands for NONlinear Mixed Effects Modeling. NONMEM is a computer program that is implemented in Fortran90/95. It solves pharmaceutical statistical problems in which within subject and between subjects variability is taken into account when fitting a pharmacokinetic and/or pharmacodynamic (PK/PD) model to data. The appropriate statistical analysis using the appropriate model helps pharmaceutical companies determine appropriate dosing strategies for their products, and increases their understanding of drug mechanisms and interactions. NONMEM software was originally developed by Lewis Sheiner and Stuart Beal and the NONMEM Project Group at the University of California, and has been used for over 30 years for population analysis by many pharmaceutical companies and the PK/PD modeling community. Its continued development and improvement by ICON Development Solutions assures pharmaceutical companies that they may continue to use the analysis tool with which they are familiar for present day pharmaceutical development. PK/PD are data typically collected from clinical studies of pharmaceutics agents, involving the administration of a drug to individuals and the subsequent observation of drug, metabolite, and/or biomarker levels (most often in the blood, plasma, or urine), as well as clinical outcome measures.
Proper modeling of these data involves accounting for both unexplainable inter- and intra-subject effects (random effects), as well as measured concomitant effects (fixed effects). NONMEM allows this mixed effect modeling. Such modeling is especially useful when there are only a few measurements from each individual sampled in the population, or when the data collection design varies considerably between these individuals. However, NONMEM is a general program which can be used to fit models to a wide variety of data. Date Version Sponsor/Medium/ Language Components Added Documentation Added (NONMEM Users Guide Series) Description 1980 I UCSF Fortran IV source code IBM specific binaries Magnetic tape NONMEM (LS, SB): FO Part I Users Basic (SB) Part II Users Supplemental (SB) FO: Simultaneously discern variability of measured levels of drug or drug response (residual variance) within a subject (intra-subject effects), and variability of PK/PD parameters between subjects (inter-subject effects). Method could determine how population PK/PD parameters related to patient characteristics, even when there were few data points per subject 1984 II Fortran 77 source code All platforms (AB) PREDPP I (SB, LS, AB): ADVAN1-6 Part VI PREDPP (SB) PREDPP: library of routines for specific PK models, designed by Stuart Beal and implemented by Alison Boeckmann Reduced the mathematical barrier which users had to overcome.
1989 III PREDPP II (SB, LS, AB): ADVAN7-9 NM-TRAN I (AB) Part IV NM-TRAN (SB) Part V Introductory (AB,LB) NM-TRAN: Data preprocessor and model translator, made it much easier for PK/PD models to be developed 1992 IV Floppy disks NONMEM (SB): FOCE/Laplace Interaction Mixture models PREDPP III: ALAG (AB), ADVAN10 (SB) NM-TRAN II (AB) Part III Installation (AB) Instructions for Installation on DOS (WB) Part VII Conditional Estimation Methods (SB) FOCE: FO was fast, but very approximate. Sometimes inaccurate assessments occurred if residual error and/or inter-subject variability were large.