In $DES block I this expression DSIZE = LOG(BSL*(EXP(-DEC*TIME) + EXP(GRO*TIME) - 1)+1)
I would use T rather than TIME You may also try tol=12 nsig=4 sign=12 Leonid > On Apr 22, 2016, at 1:37 PM, Mai, Tu [MED] <tm...@medicine.bsd.uchicago.edu> > wrote: > > 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 then > got the message error=136 instead of error=134. But minimization still > unsuccessful. Can you please help pointing out how I can solve this problem? > I really appreciate it! > > Thank you > > > > 1) Use the estimates for THETA that are ~10% deviated from the estimates > provided 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 which > it is addressed and may contain information that is privileged and > confidential. > 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, > please > notify the sender and destroy all copies of the transmittal. > > Thank you > University of Chicago Medicine and Biological Sciences > ********************************************************************************