[NMusers] No output table with YLO

[HUI, Ka Ho]

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