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
