Titration

12 views
Skip to first unread message

Dinko Rekic

unread,
Jun 14, 2017, 8:47:57 AM6/14/17
to mrgsolve, Dinko Rekić
I asked Devin this directly before I discovered the google_group. 

Basically, I want to change the dose for a drug based on the drug's concentration. Here is some code that does not work but illustrates what I want to do. 

Ideally I would like the titration to be dependent on TIME and CONC as well. Any suggestions are appreciated. 

Many thanks and Kind regards..Dinko

##------------------------------------------------------------------------------------------------------------------------------------------------------------------------

code <- '
$PARAM TVCL = 1.3, TVVC=28, TVKA=0.6, WT=70, F1=1

$SET delta=0.1

$CMT GUT CENT

$MAIN

// ... why does not this work? 
 _F(1) = F1;
if(CENT>0.1) _F(1) = F1*0.1;

//if(CENT>0.1 && TIME>48) _F(1) = F1*0.1;


double CLi = exp(log(TVCL) + 0.75*log(WT/70) + ETA(1));
double VCi = exp(log(TVVC) + ETA(2));
double KAi = exp(log(TVKA) + ETA(3));

$OMEGA name="IIV"
0.1 0 0


$ODE
dxdt_GUT = -KAi*GUT;
dxdt_CENT = KAi*GUT - (CLi/VCi)*CENT;

$TABLE
table(CP) = CENT/VCi;
table(ETA1) = ETA(1);
table(ETA2) = ETA(2);

$CAPTURE  ETA(1) ETA(2)  _F(1)
'
##------------------------------------------------------------------------------------------------------------------------------------------------------------------------


Devin Pastoor

unread,
Jun 14, 2017, 8:58:10 AM6/14/17
to Dinko Rekic, mrgsolve, Dinko Rekić
discussing with kyle....

the short answer is 'problem' is in the addl locks the dosing records (this was to try to mimic nonmem's behavior.

please post this to github as we're trying to use github as much as possible over this group so we can track issues and make them more easily discoverable for new people. 

I can post a work around on github


--
MetrumRG Website: http://www.metrumrg.com/opensourcetools.html
---
You received this message because you are subscribed to the Google Groups "mrgsolve" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mrgsolve+u...@metrumrg.com.
To post to this group, send email to mrgs...@metrumrg.com.
Visit this group at https://groups.google.com/a/metrumrg.com/group/mrgsolve/.

Kyle Baron

unread,
Jun 16, 2017, 12:41:10 AM6/16/17
to Devin Pastoor, Dinko Rekic, mrgsolve, Dinko Rekić
Dinko - 

Please do stop by github if you an find the time.  I'd like to work on this with you to get a convenient implementation for dynamic dosing like this.


Thanks,
Kyle

To unsubscribe from this group and stop receiving emails from it, send an email to mrgsolve+unsubscribe@metrumrg.com.

To post to this group, send email to mrgs...@metrumrg.com.
Visit this group at https://groups.google.com/a/metrumrg.com/group/mrgsolve/.

--
MetrumRG Website: http://www.metrumrg.com/opensourcetools.html
---
You received this message because you are subscribed to the Google Groups "mrgsolve" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mrgsolve+unsubscribe@metrumrg.com.

To post to this group, send email to mrgs...@metrumrg.com.
Visit this group at https://groups.google.com/a/metrumrg.com/group/mrgsolve/.



--
Kyle Baron, PharmD, PhD
Senior Scientist
Metrum Research Group

Reply all
Reply to author
Forward
0 new messages