Re: [R] deSolve question

2013-06-20 Thread Chris Campbell
, Chippenham, Wiltshire , SN14 OGB UK -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Andras Farkas Sent: 19 June 2013 00:46 To: r-help@r-project.org; r-sig-dynamic-mod...@r-project.org Subject: [R] deSolve question Dear All

[R] deSolve question

2013-06-18 Thread Andras Farkas
Dear All wonder if you could provide some insights on the following: currently I have this code which produces the expected results: require(deSolve) pars - list(k = 0.08,v=15) intimes - c(0,0.5,12) input - c(800,0,0) forc - approxfun(intimes, input, method=constant, rule=2) derivs -

Re: [R] deSolve question

2009-06-12 Thread Thomas Petzoldt
insun nam wrote: Dear All, I like to simulate a physiologically based pharmacokinetics model using R but am having a problem with the daspk routine. The same problem has been implemented in Berkeley madonna and Winbugs so that I know that it is working. However, with daspk it is not, and the

Re: [R] deSolve question

2009-06-12 Thread Karline
In-Sun, It is very simple. You define your state variables in the following order: y - c(Agi = 0,Alu = 0, Abr = 0, Ah = 0, Ali = 0, Ak = 0, Am = 0, Ask = 0, Aad = 0, Apa = 0, Asp = 0, Aar = 0, Ave = 0) and your rates of change in another order: dy = c(dAar, dAve, dAlu, dAli, dAbr,

[R] deSolve question

2009-06-11 Thread insun nam
Dear All, I like to simulate a physiologically based pharmacokinetics model using R but am having a problem with the daspk routine. The same problem has been implemented in Berkeley madonna and Winbugs so that I know that it is working. However, with daspk it is not, and the numbers are

Re: [R] deSolve question

2009-06-11 Thread spencerg
Dear In-Sun: I have not seen a reply, so I will offer some suggestions, hoping I can help without understanding all the details. 1. Have you run that code with options(warn=2)? It produced over 50 warnings for me, and options(warn=2) will convert the warnings to errors,