Re: simulation question
You can use S+ (or R that has an advantage of being free) to control nonmem, but you may need a bit of programming. One possible scheme is to prepare template control stream (with $SIM SUBPROB=1 option) and then run the loop
for(i in 1:NumerOfIterations) {
# prepare control file from a template
# (need to change random seed in the $SIM record)
system(NonmemCOmmandToRunNonmem)
}
The control stream should include INFN routine or INFN block that outputs parameter estimates to some file (adds them to the file, not over-writes it). R (or S+) working directory should contain all the files needed to run nonmem (in R, use setwd(NonmemDirName) to do it)
Leonid
--------------------------------------
Leonid Gibiansky, Ph.D.
President, QuantPharm LLC
web: www.quantpharm.com
e-mail: LGibiansky at quantpharm.com
tel: (301) 767 5566
Ethan Wu wrote:
> Dear users, I am trying to compare several specific PK/PD study designs by: -- run 200 simulations with the final model (develope from original dataset)
>
> -- fit the final model to the 200 simulated dataset
> To achieve above, I used $SIM SUBPROB=200 option
>
> however, nonmem would completely stop after running into estimation problem at one specific simulation/estimation cycles, for some designs it stop even before 10th iterations.
>
> Is there anyway nonmem could continue go on?
> Or, does someone know alternative way to achieve the goal?
> thanks