NonParametric engine

Applicable to Gaussian or user-defined likelihood data. The nonparametric engine is intended as a post-processor after one of the parametric engines has been run. It makes no assumptions regarding the random effects distribution, conceptually modeling the distribution as a discrete distribution on an arbitrarily fine grid in random effects space. It can be used, for example, to detect bimodality in a parameter such as a clearance. In the nonparametric log likelihood function, the parameters to be fit are the probabilities associated with each grid point in the random effects space.

If the grid is very fine, there can be an enormous number of these probabilities. However, mathematically it can be shown that at the maximum likelihood distribution, almost all of the probabilities are zero, which can be used to simplify the computation. The optimal nonparametric distribution takes the form of a discrete distribution on at most N support points, that is, at most N of the probabilities are non-zero, regardless of how many starting grid points were used.

An iteration for the nonparametric engine involves:

Selection of a set of candidate support points, which usually includes all of the support points with non-zero probability from the previous iteration plus generation of some additional candidates that are likely to improve the likelihood of the nonparametric distribution.
Computation of the optimal probabilities on the candidate support points.

The user specifies the number of iterations to apply. On the first iteration, the support points are set at the optimal post-hoc estimates from the initial parametric run. Any fixed effects associated with the residual error model or covariate models are frozen to the values from the parametric run for all of the iterations. A specially designed convex primal dual optimization engine then computes optimal probabilities on these support points. The results of the first iteration are in principle the same as the results of the nonparametric algorithm in NONMEM. However, NONMEM cannot perform any additional iterations, and the final NONMEM support points are fixed at the parametric POSTHOC values, which can be highly suboptimal.

If subsequent iterations are desired, Phoenix first discards any current iteration support points with a probability of zero, and then introduces additional candidate support points and the primal dual algorithm is reapplied to compute a new discrete distribution, which in general will include at least some of the new candidate support points. From iteration to iteration, the likelihood improves monotonically, and the support points migrate to optimal positions. The Phoenix algorithm has the capability of optimizing both probabilities and support point positions using multiple iterations. The NONMEM nonparametric algorithm can only perform a single pass that optimizes probabilities on support points fixed at POSTHOC estimates from a preceding parametric run.

The primary raw result of the Phoenix nonparametric algorithm is the optimal discrete distribution of ETAs in terms of support points and associated probabilities. The means, any covariances, and marginal distributions of each ETA of this distribution are reported. In addition to the optimal population distribution, each individual has a discrete posterior distribution from which a mean ETA value can be computed. Tables of nonparametric ETA means are produced, as are covariate vs. nonparametric ETA mean plots, which can be used to screen for potential covariate relationships.

Note:If the Model engine gives an exception, it is a general exception caused by a bad fit to data. Should an exception occur, try reconsidering the engine, initial parameters estimates, and number of compartments.