[NMusers] Re: Rounding error 134 or obj. func. is infinite (error=136)

From: Karovic, Sanja [MED] <skarovic_at_medicine.bsd.uchicago.edu>
Date: Sat, 23 Apr 2016 02:52:08 +0000

Please discard this email. My friend already posted this question.








Dear NMusers,

Our model keeps having the rounding error (error=134). I tried to fix it =
using the following methods or their combinations, however, sometimes I the=
n got the message error=136 instead of error=134. But minimization stil=
l unsuccessful. Can you please help pointing out how I can solve this probl=
em? I really appreciate it!



Thank you



1) Use the estimates for THETA that are ~10% deviated from the estimates pr=
ovided after scm run.
2) Set TOL=9, NSIG=2, SIGL=6

I attached here the code:

$PROBLEM IRM code

$INPUT ID TIME EX DV DVID ECOG AGE HT WT RACE SITE DIAG=DROP HEM

            ALK ALB LAC TES STUDY DRUG PRIORDOC

$DATA survival.csv

$SUBROUTINE ADVAN13 TOL=9

$MODEL NCOMP=1 COMP(HAZ)

$PK



;;; DECECOG-DEFINITION START

IF(ECOG.EQ.1) DECECOG = 1 ; Most common

IF(ECOG.EQ.0) DECECOG = ( 1 + THETA(27))

IF(ECOG.EQ.2) DECECOG = ( 1 + THETA(28))

;;; DECECOG-DEFINITION END





;;; GROPRIORDOC-DEFINITION START

IF(PRIORDOC.EQ.0) GROPRIORDOC = 1 ; Most common

IF(PRIORDOC.EQ.1) GROPRIORDOC = ( 1 + THETA(26))

;;; GROPRIORDOC-DEFINITION END





;;; GROECOG-DEFINITION START

IF(ECOG.EQ.1) GROECOG = 1 ; Most common

IF(ECOG.EQ.0) GROECOG = ( 1 + THETA(24))

IF(ECOG.EQ.2) GROECOG = ( 1 + THETA(25))

;;; GROECOG-DEFINITION END





;;; GROALK-DEFINITION START

IF(ALK.EQ.-99) THEN

   GROALK = 1

ELSE

   GROALK = ( 1 + THETA(23)*(ALK - 94.00))

ENDIF

;;; GROALK-DEFINITION END



;;; GRO-RELATION START

GROCOV=GROALK*GROECOG*GROPRIORDOC

;;; GRO-RELATION END





;;; DECDRUG-DEFINITION START

IF(DRUG.EQ.3) DECDRUG = 1 ; Most common

IF(DRUG.EQ.2) DECDRUG = ( 1 + THETA(21))

IF(DRUG.EQ.1) DECDRUG = ( 1 + THETA(22))

;;; DECDRUG-DEFINITION END





;;; DECAGE-DEFINITION START

DECAGE = ( 1 + THETA(20)*(AGE - 64.00))

;;; DECAGE-DEFINITION END



;;; DEC-RELATION START

DECCOV=DECAGE*DECDRUG*DECECOG

;;; DEC-RELATION END





;;; BSLHEM-DEFINITION START

IF(HEM.LE.124.00) BSLHEM = ( 1 + THETA(18)*(HEM - 124.00))

IF(HEM.GT.124.00) BSLHEM = ( 1 + THETA(19)*(HEM - 124.00))

IF(HEM.EQ.-99) BSLHEM = 1

;;; BSLHEM-DEFINITION END





;;; BSLECOG-DEFINITION START

IF(ECOG.EQ.1) BSLECOG = 1 ; Most common

IF(ECOG.EQ.0) BSLECOG = ( 1 + THETA(16))

IF(ECOG.EQ.2) BSLECOG = ( 1 + THETA(17))

;;; BSLECOG-DEFINITION END





;;; BSLALK-DEFINITION START

IF(ALK.EQ.-99) THEN

   BSLALK = 1

ELSE

   BSLALK = ( 1 + THETA(15)*(ALK - 94.00))

ENDIF

;;; BSLALK-DEFINITION END



;;; BSL-RELATION START

BSLCOV=BSLALK*BSLECOG*BSLHEM

;;; BSL-RELATION END





;;; BSHZPRIORDOC-DEFINITION START

IF(PRIORDOC.EQ.0) BSHZPRIORDOC = 1 ; Most common

IF(PRIORDOC.EQ.1) BSHZPRIORDOC = ( 1 + THETA(14))

;;; BSHZPRIORDOC-DEFINITION END





;;; BSHZALK-DEFINITION START

IF(ALK.EQ.-99) THEN

   BSHZALK = 1

ELSE

   BSHZALK = ( 1 + THETA(13)*(ALK - 94.00))

ENDIF

;;; BSHZALK-DEFINITION END



;;; BSHZ-RELATION START

BSHZCOV=BSHZALK*BSHZPRIORDOC

;;; BSHZ-RELATION END





;;; BETAHEM-DEFINITION START

IF(HEM.EQ.-99) THEN

   BETAHEM = 1

ELSE

   BETAHEM = ( 1 + THETA(12)*(HEM - 124.00))

ENDIF

;;; BETAHEM-DEFINITION END





;;; BETAECOG-DEFINITION START

IF(ECOG.EQ.1) BETAECOG = 1 ; Most common

IF(ECOG.EQ.0) BETAECOG = ( 1 + THETA(10))

IF(ECOG.EQ.2) BETAECOG = ( 1 + THETA(11))

;;; BETAECOG-DEFINITION END





;;; BETADRUG-DEFINITION START

IF(DRUG.EQ.3) BETADRUG = 1 ; Most common

IF(DRUG.EQ.2) BETADRUG = ( 1 + THETA(8))

IF(DRUG.EQ.1) BETADRUG = ( 1 + THETA(9))

;;; BETADRUG-DEFINITION END





;;; BETAAGE-DEFINITION START

BETAAGE = ( 1 + THETA(7)*(AGE - 64.00))

;;; BETAAGE-DEFINITION END



;;; BETA-RELATION START

BETACOV=BETAAGE*BETADRUG*BETAECOG*BETAHEM

;;; BETA-RELATION END







IF (NEWIND.LE.1) THEN

   SRVZ=1 ; Survivor function at TIME=0

ENDIF



;-----------BASELINE PSA--------------



TVBSL = THETA(1)



TVBSL = BSLCOV*TVBSL

BSL = TVBSL*EXP(ETA(1))



;-----------PSA PARAMETERS--------------



TVGRO = THETA(2)



TVGRO = GROCOV*TVGRO

TVDEC = THETA(3)



TVDEC = DECCOV*TVDEC



GRO = TVGRO*EXP(ETA(2))

DEC = TVDEC*EXP(ETA(3))



;-----------SURVIVAl MODEL PARAMETERS--------------



TVBSHZ = THETA(5) ; Baseline Hazard



TVBSHZ = BSHZCOV*TVBSHZ

  BSHZ = TVBSHZ

TVBETA = THETA(6) ; Parameter relating dropout hazard to PSA estimate



TVBETA = BETACOV*TVBETA

  BETA = TVBETA



$DES

  DSIZE = LOG(BSL*(EXP(-DEC*TIME) + EXP(GRO*TIME) - 1)+1)

  TEMP = BETA*DSIZE

  DADT(1) = BSHZ*EXP(TEMP)



$ERROR



IF (DVID.EQ.1) THEN

 F_FLAG=0 ; If continuous type data



 SIZE = LOG(BSL*(EXP(-DEC*TIME) + EXP(GRO*TIME) - 1)+1)

 IPRED = SIZE



   ADD = THETA(4)

SDSIZE = SQRT(ADD*ADD)

     Y = SIZE + SDSIZE*EPS(1)



     W = Y

  IRES = DV-IPRED

 IWRES = IRES/(W)

ENDIF



;-----------HAZARD PREDICTION--------------



CMHZ = A(1)

SRVT = EXP(-CMHZ)



IF (DVID.EQ.2.AND.DV.EQ.0) THEN

  F_FLAG=1

  Y = SRVT

ENDIF



IF (DVID.EQ.2.AND.DV.EQ.1) THEN

  F_FLAG=1

  Y = SRVZ-SRVT

ENDIF





$THETA (50,85,150) ; POP_BSL

 (0,0.0004,0.1) ; POP_GRO

 (0,0.02,0.1) ; POP_DEC

 (0,0.3,0.35) ; ADD_SD

 (0,0.0009,0.0014) ; BSHZ

 (0.1,0.3,0.7) ; BETA

$THETA (-0.038,-0.003,0.016) ; BETAAGE1

$THETA (-1,0.3,5) ; BETADRUG1

 (-1,0.1,5) ; BETADRUG2

$THETA (-1,-0.2,5) ; BETAECOG1

 (-1,0.3,5) ; BETAECOG2

$THETA (-0.023,-0.01,0.016) ; BETAHEM1

$THETA (0,0.001,0.010) ; BSHZALK1

$THETA (-1,-0.4,5) ; BSHZPRIORDOC1

$THETA (0,0.002,0.010) ; BSLALK1

$THETA (-1,-0.3,5) ; BSLECOG1

 (-1,0.2,5) ; BSLECOG2

$THETA (-1000000,-0.05,0.016) ; BSLHEM1

 (-0.023,-0.0002,100) ; BSLHEM2

$THETA (-0.038,-0.0009,0.016) ; DECAGE1

$THETA (-1,-0.6,5) ; DECDRUG1

 (-1,-0.8,5) ; DECDRUG2

$THETA (0,98E-006,0.010) ; GROALK1

$THETA (-1,-0.02,5) ; GROECOG1

 (-1,0.5,5) ; GROECOG2

$THETA (-1,5,7) ; GROPRIORDOC1

$THETA (-1,0.0004,5) ; DECECOG1

 (-1,-0.2,5) ; DECECOG2

$OMEGA 2.5 ; PPV_BSL

 2 ; PPV_GRO

 1.5 ; PPV_DEC

;-------------------------

$SIGMA 1. FIX

;-------------------------

$ESTIMATION METHOD=COND INTER NUMERICAL LAPLACE SLOW MAX=9999 NSIG=2

            SIGL=6 NOABORT PRINT=5

$COVARIANCE

$TABLE ID TIME IPRED BSL GRO DEC STUDY DRUG WRES CWRES ETA1 ETA2

            ETA3 ONEHEADER NOPRINT

            FILE=C:\Users\fulldata.tab

$TABLE ID TIME PRED RES WRES ONEHEADER NOPRINT

            FILE=C:\Users\original.tab










***************************************************************************=
*****
This e-mail is intended only for the use of the individual or entity to whi=
ch
it is addressed and may contain information that is privileged and confiden=
tial.
If the reader of this e-mail message is not the intended recipient, you are=
 
hereby notified that any dissemination, distribution or copying of this
communication is prohibited. If you have received this e-mail in error, ple=
ase
notify the sender and destroy all copies of the transmittal.

Thank you
University of Chicago Medicine and Biological Sciences
***************************************************************************=
*****
Received on Fri Apr 22 2016 - 22:52:08 EDT

This archive was generated by hypermail 2.3.0 : Fri Sep 27 2019 - 16:50:09 EDT