Time to event analysis
Dear All,
I am trying to run an example of time to single event data from NONMEM
archive
( http://www.mail-archive.com/[email www.mail-archive.com%2fnmusers%40globomaxnm.com%2fmsg01743.html">http://redir.aspx/?C=652afed2db90470ba9ed0ca22fecfac9&URL=http%3a%2f%2fwww.mail-archive.com%2fnmusers%40globomaxnm.com%2fmsg01743.html),
using NONMEM 7.2. It gives an error in control statement: $EST "laplace
required with likelihood, Pop. data, unless omega=0". I am guessing that
it might be a NONMEM version issue. Since I am new to NONMEM, it will be
good to have your suggestions.
*Dataset:*
C STUD ID TIME=DV CS
DOSE
MDV . 101 4 0 1
0.4 1 . 101 4 111 1
0.4 0 . 101 5 0 1
0.4 1
. 101 5 51 1
0.4 0 . 101 6 0 1
0.8 1
. 101 6 31 1
0.8 0 . 101 7 0 1
0.8 1
. 101 7 . 50
1
0.8 0 . 101 8 0 0 1.6 1
. 101 8 59 1 1.6 0 . 101 9 0 0 3.2 1
. 101 9 477 1
3.2 0 . 101 11 0 1
3.2 1
. 101 11 278 1 3.2 0 . 101 12 0 1
6.4 1
. 101 12 58 1
6.4 0
CS=0 is the censored event
Time is in days
---------------------------------------------------------------------------------------------
*Model:*
$PROB Time To Event data
$INPUT C STUD ID TIME=DV CS DOSE MDV
$DATA hz1.csv IGNORE=C
$SUBR ADVAN=6 TOL=6
$MODEL COMP=(RISK)
$PK
BASE = THETA(1)*EXP(ETA(1)) ;the ETA is a placeholder here
$DES DADT(1)=BASE ;hazard
$ERROR
HZ = A(1) ;hazard
SUR = EXP(-HZ) ;survival probability
DENS=HZ*SUR
Y = (1-CS)*DENS + CS*SUR
$THETA (0,.025) ; BASE $OMEGA 0 FIX
$ESTIM MAXEVAL=9990 METHOD = COND LAPLACE LIKE PRINT=1 MSFO=msfb01
$COV PRINT=E
----
Thanks,
Gaurav