NMsim_NWPRI: Simulate with parameter variability using the NONMEM NWPRI...

In NMsim: Seamless 'Nonmem' Simulation Platform

Simulate with parameter variability using the NONMEM NWPRI subroutine

Description

Modify control stream for simulation with uncertainty using inverse-Wishart distribution for OMEGA and SIGMA parameters

This function does not run any simulations. To simulate, using this method, see 'NMsim()'. See examples.

Usage

NMsim_NWPRI(file.sim, file.mod, data.sim, PLEV = 0.999)

Arguments

file.sim	The path to the control stream to be edited. This function overwrites the contents of the file pointed to by file.sim.
file.mod	Path to the path to the original input control stream provided as 'file.mod' to 'NMsim()'.
data.sim	Included for compatibility with 'NMsim()'. Not used.
PLEV	Used in $PRIOR NWPRI PLEV=0.999. This is a NONMEM argument to the NWPRI subroutine. When PLEV < 1, a value of THETA will actually be obtained using a truncated multivariate normal distribution, i.e. from an ellipsoidal region R1 over which only a fraction of mass of the normal occurs. This fraction is given by PLEV.

Details

Simulate with parameter uncertainty. THETA parameters are sampled from a multivariate normal distribution while OMEGA and SIGMA are simulated from the inverse-Wishart distribution. Correlations of OMEGA and SIGMA parameters will only be applied within modeled "blocks".

Value

Path to simulation control stream

Author(s)

Brian Reilly, Philip Delff

References

inverse-Wishart degrees of freedom calculation for OMEGA and SIGMA: NONMEM tutorial part II, supplement 1, part C.

See Also

NMsim_VarCov

Examples

## Not run: 
simres <- NMsim(file.path,method.sim=NMsim_WPRI,typical=TRUE,subproblems=500)

## End(Not run)

NMsim documentation built on Nov. 2, 2024, 9:06 a.m.