RE: simulation question

From: Jakob Ribbing Date: February 06, 2009 technical Source: mail-archive.com
Ethan, Sebastian is right that a non-parametric bootstrap may be suitable for determining the uncertainty in the population parameters. However, I got the impression that you wanted to investigate possible study designs on how informative they are for a future model-based analysis? If you would like to do simulation based on your current best guess of the model parameters (i.e. the point estimates) and you would like to do model based analysis of the future study in isolation, then PsN has another program which is highly efficient and would automatically provide you with summary statistics, without any programming required. This program is called sse for stochastic simulation and estimation: http://psn.sourceforge.net/PDF_docs/sse_userguide.pdf sse would also allow you to investigate the performance of alternative models, e.g. simulation with a two-compartment model and estimation with a one-compartment model, to see if CL can be estimated with good precision and low bias, even from the sparse data of a future study. You need to generate the data sets with the study designs yourself, but the rest is really slick. If you are planning to analyse the new study in conjunction with the currently-available data (i.e. a pooled analysis) there may be some clever way of tweaking sse into evaluating this, but I could not say exactly how to best achieve that. (maybe someone in Uppsala has a suggestion in that case) Best regards Jakob
Quoted reply history
________________________________ From: [email protected] [mailto:[email protected]] On Behalf Of Sebastian Ueckert Sent: 06 February 2009 09:51 To: Ethan Wu Cc: [email protected] Subject: Re: [NMusers] simulation question Dear Ethan, the simplest solution to solve your problem would be to use the bootstrap command of PSN ( http://psn.sourceforge.net/). With PSN installed you would simply do: bootstrap final_model.mod -samples=200 PSN would take care of unsuccessful runs and provide a nice summary of the individual estimates. Best regards Sebastian On Thu, Feb 5, 2009 at 11:18 PM, Ethan Wu <[email protected]> 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
Feb 05, 2009 Ethan Wu simulation question
Feb 05, 2009 Leonid Gibiansky Re: simulation question
Feb 06, 2009 Ethan Wu simulation question
Feb 06, 2009 Leonid Gibiansky Re: simulation question
Feb 06, 2009 Sebastian Ueckert Re: simulation question
Feb 06, 2009 Elodie Plan RE: simulation question
Feb 06, 2009 Jakob Ribbing RE: simulation question
Feb 06, 2009 Martin Bergstrand RE: simulation question
Feb 06, 2009 Marc Gastonguay Re: simulation question
Feb 06, 2009 Nick Holford Re: simulation question
Feb 06, 2009 Ekaterina Gibiansky RE: simulation question
Feb 06, 2009 Rob ter Heine Re: simulation question
Feb 06, 2009 Serge Guzy RE: simulation question