Re: simulation question
Hi Ethan,
If think you get the message "message issued from estimation step" you still
have the $ESTIMATION block in your control stream, although you're using
SIMONLY option. Try removing the $ESTIMATION block and try again. Perhaps you
could post your control stream as well?
Sincerely,
Rob
_______________________________
Rob ter Heine, MSc, PharmD
Department of Pharmacology, Slotervaart Hospital
Amsterdam, The Netherlands
E: [email protected]
T: +31-20-5124737
>>> Ethan Wu <[email protected]> 02/06/09 2:28 >>>
Hi Robm
the simulatio step runs fine as I checked by using SIMONLY
the termination says
"0PROGRAM TERMINATED BY OBJ
ERROR IN NCONTR WITH INDIVIDUAL 16 ID=0.16000000E+02
NUMERICAL HESSIAN OF OBJ. FUNC. FOR COMPUTING CONDITIONAL ESTIMATE
IS NON POSITIVE DEFINITE
MESSAGE ISSUED FROM ESTIMATION STEP
AT 0TH ITERATION, UPON EVALUATION OF GRADIENT WITH RESPECT TO OMEGA
0TRY SUBSTANTIALLY INCREASING "INITIAL ESTIMATE OF SIGMA"
OR DECREASING "INITIAL ESTIMATE OF OMEGA""
Quoted reply history
________________________________
From: Rob ter Heine <[email protected]>
To: [email protected]; [email protected]
Sent: Friday, February 6, 2009 5:54:33 AM
Subject: Betr.: Re: [NMusers] simulation question
Dear Ethan,
Have you tried removing the $ESTIMATION step and replace all values with your
final estimations? What does your result file say after termination?
Sincerely
Rob
>>> Sebastian Ueckert <[email protected]> 02/06/09 10:51 >>>
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
>
>
>