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
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
See the documentation for $DATA. Set up multiple problems within a single
control file. Use the RECORDS= option to select the number of data records
for 50 patients, use the NOREWIND option so that for each problem, the data
for the next 50 patients.
_____
Quoted reply history
From: [email protected] [mailto:[email protected]] On
Behalf Of kehua wu
Sent: Monday, March 21, 2011 2:55 PM
To: [email protected]
Subject: [NMusers] repetitive parameter estimation
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