NONMEM Users Network Archive

Hosted by Cognigen

Time to event analysis using NONMEM

From: Anubha_Gupta
Date: Tue, 9 Jun 2009 16:23:39 +0100



Dear NMUsers

I am trying to analyse time to single event data using=
 NONMEM VI. I am
getting the following termination message. If some body co=
uld point out
what I am doing wrong or may be elaborate on the err=
or message that would
be very helpful.

Thanks in advance

Anubha


0PROGRAM TERMINATED BY OBJ
 ERROR IN CLIK WITH INDIVIDUAL 1 ID=0=
.40000000E+01
 CONDITIONAL LIKELIHOOD FOR DATA FROM OBS RECORD 1
 ESTIMATED TO BE NEGLIGIBLE
 MESSAGE ISSUED FROM ESTIMATION STEP
 AT INITIAL OBJ. FUNCTION EVALUATION


Data set

                         =
                         =
           
    C | STUD | =
 I | TIME | CS | DOSE
                | =
  | D | | | =
       
  --------------+-----------+----+--------+--------+--------
    . | 101 =
| 4 | 111 | 0 | 0.4 =

  --------------+-----------+----+--------+--------+--------
    . | 101 =
| 5 | 51 | 0 | 0.4 =
 
  --------------+-----------+----+--------+--------+--------
    . | 101 =
| 6 | 31 | 0 | 0.8 =
 
  --------------+-----------+----+--------+--------+--------
    C Missing | 101 | 7 | =
 . | . | 0.8
    DV | =
 | | | | =
      
  --------------+-----------+----+--------+--------+--------
    . | 101 =
| 8 | 59 | 0 | 1.6 =
 
  --------------+-----------+----+--------+--------+--------
    . | 101 =
| 9 | 477 | 1 | 3.2 =

  --------------+-----------+----+--------+--------+--------
    . | 101 =
| 1 | 278 | 0 | 3.2 =

                | =
  | 1 | | | =
       
  --------------+-----------+----+--------+--------+--------
    . | 101 =
| 1 | 58 | 0 | 6.4 =
 
                | =
  | 2 | | | =
       
                         =
                         =
           


CS=1 is the censored event
TIME is in days


Model file:

$PROB Time To Event data

$INPUT C STUD ID TIME=DV CS DOSE

$DATA

IGNORE=C

$SUBR ADVAN=6 TOL=6
$MODEL COMP=(RISK)

$PK

  BASE = THETA(1)*EXP(ETA(1)) ;the ETA is a placeholde=
r here

$DES
  DADT(1)=BASE ;hazard

$ERROR


  HZ = A(1) ;haza=
rd
  SUR = EXP(-HZ) ;survival pro=
bability
  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


THIS EMAIL AND ANY ATTACHED FILES ARE CONFIDENTIAL.
If you are not the intended recipient you are notified =
that any disclosure, reproduction,copying, distribution, or =
action taken in reliance on the contents of this info=
rmation is strictly prohibited. If you have received =
this transmission in error please notify the sender immedi=
ately and then delete this email. Email transmission ca=
nnot be guaranteed to be secure or error free as infor=
mation could be intercepted, corrupted, lost, destroyed,arrive=
 late or incomplete, or contain viruses. The company/sende=
r accepts no liability for any changes made to this =
email during transmission or any damage caused by any =
virus transmitted by this email. Any views or opinions =
expressed in this email are solely those of the author =
and do not necessarily represent those of the company. T=
he company/sender accepts no liability for the content =
of this email, or for the consequences of any actions =
taken on the basis of the information provided, unless t=
hat information is subsequently confirmed in writing.
Received on Tue Jun 09 2009 - 11:23:39 EDT

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request@iconplc.com.

Once subscribed, you may contribute to the discussion by emailing: nmusers@globomaxnm.com.