Re: repetitive parameter estimation
Kehua,
You dont describe how you got the 250,000 simulated patients but the simplest way to achieve what you want is to create a data file template with 50 patients then use NONMEM to simulate then estimate with 1000 subproblems.
e.g.
$PROB simest
$DATA my50patients.dat
$INPUT ID TIME etc
$SIM (12345) NSUB=1000
$EST METHOD=COND etc
$PK etc
The only problem with this method is the poor error trapping in NONMEM. If you have one run out the 1000 that falls over then all the remaining runs will be omitted.
Wings for NONMEM has a workaround for this very long standing and apparently incurable bug with its nmgosim command.
Nick
Quoted reply history
On 21/03/2011 8:55 p.m., kehua wu wrote:
> Hi,
>
> We have a data file with 250,000 simulated patients, and we want to estimate a parameter in NONMEM using a different set of 50 patients at a time (5,000 replicates, but not bootstrapping). Is there an easy way that anybody knows to automatically repeat the parameter estimation while substituting the new data each time?
>
> Thanks,
>
> Kehua
--
Nick Holford, Professor Clinical Pharmacology
Dept Pharmacology& Clinical Pharmacology
University of Auckland,85 Park Rd,Private Bag 92019,Auckland,New Zealand
tel:+64(9)923-6730 fax:+64(9)373-7090 mobile:+64(21)46 23 53
email: [email protected]
http://www.fmhs.auckland.ac.nz/sms/pharmacology/holford