From: Ethan Wu <*ethan.wu75*>

Date: Fri, 19 Dec 2008 09:34:39 -0800 (PST)

Dear all, I am trying to run a PK/PD model with PD being binary data. Aft=

er looking into previous post I got my control stream running. My question =

is that how to assess the model fit? Can I still use PRED, RES and WRES? At=

tached is the code I found and used after a bit modification. ########=

################################################## From: "Piotrovskij, Vl=

adimir [JanBe]" <VPIOTROV

e to solve some of your problems in NONMEM. However, the best way is to a=

pply generalized linear regression using one of the statistical packages.=

response: fixed effect of schedule $DATA nmd.ssc $INPUT ID AUC SCHD DV=

HD3 = 0 SCHD4 = 0 IF (SCHD.EQ.1) SCHD1 = 1 IF (SCHD.EQ.2) SCHD2=

= 1 IF (SCHD.EQ.3) SCHD3 = 1 IF (SCHD.EQ.4) SCHD4 = 1 SLOPE =

= THETA(1) E50 = SCHD1*THETA(2)+SCHD2*THETA(3)+SCHD3*THETA(4)+SCHD4*T=

HETA(5) INT = -LOG(E50) * SLOPE LOGIT = INT + SLOPE * LOG(AUC) + ET=

A(1) A=EXP(LOGIT) P=A/(1+A) IF (DV.EQ.1) Y=P IF (DV.EQ.0) Y=

=1-P $THETA (2 5 7); 1 SIGM (0 30 50); 2 E50 SCHD=1 (20 50 70)=

; 3 E50 SCHD=2 (40 70 90); 4 E50 SCHD=3 (60 100 200); 5 E50 SCHD=4=

Received on Fri Dec 19 2008 - 12:34:39 EST

Date: Fri, 19 Dec 2008 09:34:39 -0800 (PST)

Dear all, I am trying to run a PK/PD model with PD being binary data. Aft=

er looking into previous post I got my control stream running. My question =

is that how to assess the model fit? Can I still use PRED, RES and WRES? At=

tached is the code I found and used after a bit modification. ########=

################################################## From: "Piotrovskij, Vl=

adimir [JanBe]" <VPIOTROV

e to solve some of your problems in NONMEM. However, the best way is to a=

pply generalized linear regression using one of the statistical packages.=

response: fixed effect of schedule $DATA nmd.ssc $INPUT ID AUC SCHD DV=

HD3 = 0 SCHD4 = 0 IF (SCHD.EQ.1) SCHD1 = 1 IF (SCHD.EQ.2) SCHD2=

= 1 IF (SCHD.EQ.3) SCHD3 = 1 IF (SCHD.EQ.4) SCHD4 = 1 SLOPE =

= THETA(1) E50 = SCHD1*THETA(2)+SCHD2*THETA(3)+SCHD3*THETA(4)+SCHD4*T=

HETA(5) INT = -LOG(E50) * SLOPE LOGIT = INT + SLOPE * LOG(AUC) + ET=

A(1) A=EXP(LOGIT) P=A/(1+A) IF (DV.EQ.1) Y=P IF (DV.EQ.0) Y=

=1-P $THETA (2 5 7); 1 SIGM (0 30 50); 2 E50 SCHD=1 (20 50 70)=

; 3 E50 SCHD=2 (40 70 90); 4 E50 SCHD=3 (60 100 200); 5 E50 SCHD=4=

Received on Fri Dec 19 2008 - 12:34:39 EST