RE: problem with simulation

From: "kai wu" Subject: RE: [NMusers] problem with simulation Date: Tue, November 30, 2004 3:31 pm Sorry, I mistakenly post part of the code for a one-comp. PK. For the $DES block, it should be: $DES DADT(1)=-KA*A(1) DADT(2)=KA*A(1)+K21*A(3)-K12*A(2)-K10*A(2) DADT(3)=K12*A(2)-K21*A(3) CP=A(2)/S2 EFFC=1-CP*0.1/(EC50+CP*0.1) DADT(4)=RCORT*EFFC-KE*A(4) Cort=A(4)/S4 KTC=30000000 KALB=5000 QTC=0.0000007 QALB=0.00055 MW=362.47*(10**6) CortM=Cort/MW B1=KTC*QTC+KALB*QALB+1-CortM*KTC A1=KALB*QALB*KTC+KTC CORTFM=(-B1+SQRT(B1**2+4*A1*CortM))/(2*A1) Cortf=CORTFM*MW UPP=CP*0.1+Cortf*RRBC/RRBE DOW=EC50CL+CP*0.1+Cortf*RRBC/RRBE EFFL=1-UPP/DOW DADT(5)=KIN*EFFL-KOUT*A(5) To initialize the PD compartment, I tried two methods, one is similar to Nick suggested, I set up like: #ID TIME AMT DV CMT EVID OCC 1 0 115 4 1 1 1 0 100 5 1 1 1 0 0 115 4 0 1 1 0 0 100 5 0 1 The other method is to start actual dosing at time 24 instead of time 0, which give 24 hours for system to achieve sd st. Both method simulate the same concentrations. _______________________________________________________