- Contemporary messages sorted: [ by date ] [ by thread ] [ by subject ] [ by author ] [ by messages with attachments ]

From: HUI, Ka Ho <matthew.hui_at_link.cuhk.edu.hk>

Date: Thu, 28 Jul 2016 08:20:51 +0000

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=_at_ 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

Received on Thu Jul 28 2016 - 04:20:51 EDT

Date: Thu, 28 Jul 2016 08:20:51 +0000

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=_at_ 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

Received on Thu Jul 28 2016 - 04:20:51 EDT