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
>
>