[R] using LSODA in R
Ben Bolker
bolker at zoo.ufl.edu
Wed Nov 5 14:43:53 CET 2003
Try returning list(c(Rprime,Cprime,Pprime),NULL) -- the first element in
the returned list should be a numeric *vector* of the derivatives.
Ben
On Tue, 4 Nov 2003, Ivan Kautter wrote:
> R help list subscribers,
>
> I am a new user of R. I am attempting to use R to explore a set of
> equations specifying the dynamics of a three trophic level food chain. I
> have put together this code for the function that is to be evaluted by
> LSODA. My equations Rprime, Cprime, and Pprime are meant to describe the
> actual equation of the derivative. When I run LSODA, I do not get the
> output that these equations should be giving. Can someone tell me if I have
> set this function up correctly to use with LSODA when the user is specifying
> the equation of the derivative or offer some advice for using LSODA in R?
> An example of how to code for user specified differential equations would be
> great.
>
> function(times,y,p)
> {
> Rprime <-
> (R*(1-R))-((xc*yc*C*R)/(R+R0))-((w*xp*ypr*P*R)/(R02+((1-w)*C)+(w*R)))
> Cprime <-
> (-1*(xc*C)*(1-(yc*R)/(R+R0)))-(((1-w)*xp*ypc*P*C)/((w*R)+((1-w)*C)+C0))
> Pprime <-
> (-1*P)-(((1-w)*xp*ypc*C*P)/((w*R)+((1-w)*C)+C0))+((w*xp*ypr*P*R)/((w*R)+((1-w)*C)+R02))
> list(c(Rprime, Cprime, Pprime))
> }
>
> The above is the function yprime which the documentation for the odesolve
> says that I may specify.
>
> Thanks for any help that anyone can provide.
>
> Ivan Kautter
>
> _________________________________________________________________
> Compare high-speed Internet plans, starting at $26.95.
> https://broadband.msn.com (Prices may vary by service area.)
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
>
--
620B Bartram Hall bolker at zoo.ufl.edu
Zoology Department, University of Florida http://www.zoo.ufl.edu/bolker
Box 118525 (ph) 352-392-5697
Gainesville, FL 32611-8525 (fax) 352-392-3704
More information about the R-help
mailing list