RE: Help on a cell transit model
From: "Urien" s.urien@stcloud-huguenin.org
Subject: RE: [NMusers] Help on a cell transit model
Date: Wed, 1 Mar 2006 11:12:58 +0100
Peter
I have already used this model successfully
I suggest to try the following modification :
write the $ERROR block after the $DES
then specify that if CMT3 IPRED=A3
IF (CMT .EQ. 3 .AND. TIME .GT. 0) IPRED = A3
Saik
Sak URIEN, MD, PhD
Directeur de Recherche, INSERM
Service de Pharmacologie
Centre Ren Huguenin
35 rue Dailly
F92210 St Cloud
Tel. 01 4711 1684
Fax 01 4711 1617
email : s.urien@stcloud-huguenin.org
$DES
CP = A(1)/V1
DRUG = SLOP*CP
DADT(1) = -(CL + Q)*A(1)/V1 + Q*A(2)/V2
DADT(2) = Q*A(1)/V1 - Q*A(2)/V2
DADT(3) = KP*A(3) - KP*A(3) ; STEM
DADT(4) = KP*A(3) - KP*A(4) ; TRANSIT 1
DADT(5) = KP*A(4) - KP*A(5) ; TRANSIT 2
DADT(6) = KP*A(5) - KP*A(6) ; TRANSIT 3
DADT(7) = KP*A(6) - KP*(1 + DRUG)*A(7) ; CELL COUNTS
A3 = A(3)
A4 = A(4)
A5 = A(5)
A6 = A(6)
A7 = A(7)
$ERROR
PKY=0
PDY=0
IF (F .GT. 0 .AND. CMT .EQ. 1 .AND. TIME .GT. 0) IPRED = LOG(F)
IF (CMT .EQ. 3 .AND. TIME .GT. 0) IPRED = A3
IF (CMT .EQ. 1 .AND. TIME .GT. 0) PKY = IPRED + EPS(1)
IF (CMT .EQ. 3 .AND. TIME .GT. 0) PDY = IPRED*EXP(EPS(2)) + EPS(3)
Y = PKY + PDY