Hello NONMEM Users,
As SAEM looks like a fast alternative to some other methods, I tried it many
times as a stratup method. It worked well except for one imortant occasion
when NONMEM crashed without any written warnings. A pop-up message was saying
"nonmem.exe has stopped working." The *.out file looks as follows:
GENERAL NONLINEAR KINETICS MODEL USING LSODA (ADVAN13)
0MODEL SUBROUTINE USER-SUPPLIED - ID NO. 9999
0MAXIMUM NO. OF BASIC PK PARAMETERS: 14
0COMPARTMENT ATTRIBUTES
COMPT. NO. FUNCTION INITIAL ON/OFF DOSE DEFAULT DEFAULT
STATUS ALLOWED ALLOWED FOR DOSE FOR OBS.
1 C1 ON YES YES YES YES
2 C2 ON YES YES NO NO
3 C3 ON YES YES NO NO
4 C4 ON YES YES NO NO
5 C5 ON YES YES NO NO
6 C6 ON YES YES NO NO
7 C7 ON YES YES NO NO
8 C8 ON YES YES NO NO
9 C9 ON YES YES NO NO
10 C10 ON YES YES NO NO
11 C11 ON YES YES NO NO
12 OUTPUT OFF YES NO NO NO
0NRD VALUE(S) FROM SUBROUTINE TOL: 5
1
ADDITIONAL PK PARAMETERS - ASSIGNMENT OF ROWS IN GG
COMPT. NO. INDICES
SCALE BIOAVAIL. ZERO-ORDER ZERO-ORDER ABSORB
FRACTION RATE DURATION LAG
1 * 17 * * *
2 * 18 * * *
3 15 * * * *
4 * * * * *
5 * * * * *
6 * * * * *
7 16 * * * *
8 * * * * *
9 * * * * *
10 * * * * *
11 * * * * *
12 * - - - -
- PARAMETER IS NOT ALLOWED FOR THIS MODEL
* PARAMETER IS NOT SUPPLIED BY PK SUBROUTINE;
WILL DEFAULT TO ONE IF APPLICABLE
0DATA ITEM INDICES USED BY PRED ARE:
EVENT ID DATA ITEM IS DATA ITEM NO.: 10
TIME DATA ITEM IS DATA ITEM NO.: 4
DOSE AMOUNT DATA ITEM IS DATA ITEM NO.: 5
DOSE RATE DATA ITEM IS DATA ITEM NO.: 7
COMPT. NO. DATA ITEM IS DATA ITEM NO.: 11
PRED. COMPT. DATA ITEM IS DATA ITEM NO.: 12
0PK SUBROUTINE CALLED WITH EVERY EVENT RECORD.
PK SUBROUTINE NOT CALLED AT NONEVENT (ADDITIONAL OR LAGGED) DOSE TIMES.
0PK SUBROUTINE INDICATES THAT COMPARTMENT AMOUNTS ARE INITIALIZED.
0ERROR SUBROUTINE CALLED WITH EVERY EVENT RECORD.
0DES SUBROUTINE USES FULL STORAGE MODE.
1
#TBLN: 1
#METH: Stochastic Approximation Expectation-Maximization
EM/BAYES SETUP
THETAS THAT ARE MU MODELED:
1 2 3 4 5 6
THETAS THAT ARE SIGMA-LIKE:
MONITORING OF SEARCH:
Stop Time:
Sat 11/17/2012
05:25 AM
When I tried Bayes, the program did not crash, but the method did not converge
well.
Is there any way to fix the SAEM run? What can cause this crash?
Thanks,
Pavel