NONMEM Users Network Archive

Hosted by Cognigen

assess model for logistic regression

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

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request@iconplc.com.

Once subscribed, you may contribute to the discussion by emailing: nmusers@globomaxnm.com.