Dear All,
I like to use nlmeode but so far having a trouble with setting up a function
part. It does not allow me to use Cliver etc in the DiffEq and it does not
allow me to have Vabs in the ObsEq. Would you kindly let me know what to do?
Best wishes,
In-Sun
FunODE <- list(
DiffEq=list(
dAliverdt = ~ -Qh *Cinlet - Qh * Cliver/Kp - fb*CLint * Cliver/Kp,
dAinletdt = ~ Qh * Csys + Vabs - Qh * Cinlet,
dAsysdt = ~ Qh*Cliver/Kp - Qh*Csys - CLnh*Csys),
ObsEq=list(
Cliver = ~ Aliver/Vliver,
Cinlet = ~ Ainliet/Vinlet,
Csys = ~ Asys/Vsys,
Vabs = ~ka*Dose*Fa*exp(-ka*times)),
States=c("Aliver","Ainlet","Asys"),
Parms=c( "Dose","Fa","fb","ka","CLint","CLnh","Kp","Qh", "Vliver", "Vsys",
"Vinlet", "Vabs"),
Init=list(0,0,0))
--
Dr In-Sun Nam Knutsson
Research Associate
The Centre for Applied Pharmacokinetic Research (CAPKR)
School of Pharmacy and Pharmaceutical Sciences
University of Manchester
Stopford Building
Oxford Road
Manchester
U.K.
Phone: +44 161 275 2355
Email: in-sunnam.knuts...@manchester.ac.uk
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
