Time to event analysis

From: Gaurav Bajaj Date: July 28, 2011 technical Source: mail-archive.com
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
Jul 28, 2011 Gaurav Bajaj Time to event analysis
Jul 28, 2011 Nick Holford Re: Time to event analysis