No output table with YLO
Dear all,
I just started using the M2 method to fit a two-compartment model to a
left-censored dataset. However, no table, as defined by $TABLE, could be output.
If I remove just the line “YLO=BMT”, leaving all other coding unchanged, the
table can be generated.
I have failed to figure out why. Any advice?
Thanks,
Matthew
↓The code looks like the followings:
$INPUT C ID DOSE_N ID_IDC DOSE_IDC SEX AGE AGED HT WT BMI BSA DX AGET K_SCHW
DOSE_BY_BSA DOSE_MG AMT=DOSE_UM TIME_INF RATE_MG RATE=RATE_UM TIME DV=CP_UM
CP_LN_UM CP_MG SCR_UM SCR_MG CRCL ALT MDV EVID
$DATA ..\DATASET_3.csv IGNORE=@ IGNORE=(DX.EQ.1)
$SUBROUTINE ADVAN3 TRANS4
$PK
;SETTING BMT
IF (NEWIND.EQ.0) THEN
BMT=0.1
ENDIF
;CLEARANCE
TVCL=THETA(1)
CL=TVCL*EXP(ETA(1))
EXP_ETA_CL=EXP(ETA(1))
;CENTRAL VOLUME
TVV1=THETA(2)*(BSA/0.735)**THETA(6)
V1=TVV1*EXP(ETA(2))
EXP_ETA_V1=EXP(ETA(2))
;INTERCOMPARTMENTAL CLEARANCE
TVQ=THETA(3)
Q=TVQ*EXP(ETA(3))
EXP_ETA_Q=EXP(ETA(3))
;PERIPHERAL VOLUME
TVV2=THETA(4)
V2=TVV2*EXP(ETA(4))
EXP_ETA_V2=EXP(ETA(4))
;REQUIRED CONSTANT
S1=V1
$THETA
;INITIAL ESTIMATES
$ERROR
YLO=BMT
W=THETA(5)
Y=F*EXP(W*EPS(1))
IPRED=F
FLAG=0
IF (DV.EQ.0) FLAG=1
EPSI=LOG(FLAG+DV/IPRED)/W
$OMEGA
0.09 ;OMG_CL
0 FIXED ;OMG_V1
0 FIXED ;OMG_Q
0 FIXED ;OMG_V2
$SIGMA
1 FIXED ;SIG_CP
$EST MAXEVAL=9999 PRINT=1 METHOD=1 INTERACTION LAPLACIAN NOABORT
$COV PRINT=E
$TABLE ID DOSE_N ID_IDC DOSE_IDC SEX AGE AGED HT WT BMI BSA AGET DOSE_BY_BSA
DOSE_UM TIME_INF RATE_UM TIME CP_UM CP_LN_UM CRCL ALT MDV EVID CL EXP_ETA_CL V1
EXP_ETA_V1 Q EXP_ETA_Q V2 EXP_ETA_V2 EPSI NOAPPEND NOPRINT ONEHEADER
FILE=OUTPUT.txt