Re: [R] Solution to differential equation

2010-12-17 Thread Ravi Varadhan
.@jhmi.edu -Original Message- From: dave fournier [mailto:da...@otter-rsch.com] Sent: Friday, December 17, 2010 7:02 AM To: Ravi Varadhan Cc: r-help@r-project.org Subject: Re: [R] Solution to differential equation Ravi Varadhan wrote: Because the numerical solution is more flexible. In the exampl

Re: [R] Solution to differential equation

2010-12-17 Thread dave fournier
ave fournier Sent: Friday, December 17, 2010 11:23 AM To: r-help@r-project.org Subject: Re: [R] Solution to differential equation It is not very difficult to integrate this DE numerically. For parameter estimation it is a good idea for stability to use a semi-implicit formulation. The idea is d

Re: [R] Solution to differential equation

2010-12-17 Thread Ravi Varadhan
AM To: r-help@r-project.org Subject: Re: [R] Solution to differential equation It is not very difficult to integrate this DE numerically. For parameter estimation it is a good idea for stability to use a semi-implicit formulation. The idea is described here. http://otter-rsch.com/a

Re: [R] Solution to differential equation

2010-12-17 Thread dave fournier
It is not very difficult to integrate this DE numerically. For parameter estimation it is a good idea for stability to use a semi-implicit formulation. The idea is described here. http://otter-rsch.com/admodel/cc4.html __ R-help@r-project.org mail

Re: [R] Solution to differential equation

2010-12-17 Thread Ben Bolker
Ben Bolker gmail.com> writes: > Mike Marchywka hotmail.com> writes: [snip] > The gsl package has this function, apparently -- it agrees with > Mathematica/Wolfram Alpha's Hypergeometric2F1 for a single set of > inputs (2,3,4,0.5), although apparently the algorithm that GSL has > only conv

Re: [R] Solution to differential equation

2010-12-17 Thread Ravi Varadhan
oject.org [mailto:r-help-boun...@r-project.org] On Behalf Of Mike Marchywka Sent: Friday, December 17, 2010 6:02 AM To: r-help@r-project.org; msamt...@gmail.com Subject: Re: [R] Solution to differential equation sorry, wanted to CC list hit wrong button no caffeine > > From: rva

Re: [R] Solution to differential equation

2010-12-17 Thread Ben Bolker
Mike Marchywka hotmail.com> writes: [snip] > > did you see my earlier post with link to wolfram integrator? Where i also > > requested anyone wanting to get rid of a copy of G&R Integral Tables to > > contact me off list since a dog really did eat mine? I think it came up > > with "F" or hype

Re: [R] Solution to differential equation

2010-12-17 Thread Mike Marchywka
sorry, wanted to CC list hit wrong button no caffeine > > From: rvarad...@jhmi.edu > > To: rvarad...@jhmi.edu > > Date: Thu, 16 Dec 2010 22:37:17 -0500 > > CC: r-help@r-project.org; msamt...@gmail.com > > Subject: Re: [R] Solution to differential equation >

Re: [R] Solution to differential equation

2010-12-16 Thread Ravi Varadhan
School of Medicine Johns Hopkins University Ph. (410) 502-2619 email: rvarad...@jhmi.edu - Original Message - From: Ravi Varadhan Date: Thursday, December 16, 2010 4:11 pm Subject: RE: [R] Solution to differential equation To: 'Ravi Varadhan' , 'Scionforbai' ,

Re: [R] Solution to differential equation

2010-12-16 Thread Ravi Varadhan
-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Scionforbai Sent: Wednesday, December 15, 2010 12:28 PM To: mahesh samtani Cc: r-help@r-project.org Subject: Re: [R] Solution to differential equation > I am trying to find the analytical solution to this differential equation >

Re: [R] Solution to differential equation

2010-12-15 Thread Ravi Varadhan
PM To: mahesh samtani Cc: r-help@r-project.org Subject: Re: [R] Solution to differential equation > I am trying to find the analytical solution to this differential equation > > dR/dt = k1*(R^k2)*(1-(R/Rmax)); R(0) = Ro > If there is an analytial solution to this differential equation th

Re: [R] Solution to differential equation

2010-12-15 Thread Scionforbai
> I am trying to find the analytical solution to this differential equation > > dR/dt = k1*(R^k2)*(1-(R/Rmax)); R(0) = Ro > If there is an analytial solution to this differential equation then it It is a polynomial function of R, so just develop the expression and when you get the two terms in R

Re: [R] Solution to differential equation

2010-12-15 Thread Mike Marchywka
> Date: Wed, 15 Dec 2010 11:46:40 -0500 > From: msamt...@gmail.com > To: r-help@r-project.org > Subject: [R] Solution to differential equation > > Hello, > I am trying to find the analytical solution to this differential equation

Re: [R] Solution to differential equation

2010-12-15 Thread Jeff Newmiller
ODEs don't get much easier than this... integration by parts is overkill. See the posting guide re homework. "mahesh samtani" wrote: >Hello, >I am trying to find the analytical solution to this differential >equation > >dR/dt = k1*(R^k2)*(1-(R/Rmax)); R(0) = Ro > >k1 and k2 are parameters th

[R] Solution to differential equation

2010-12-15 Thread mahesh samtani
Hello, I am trying to find the analytical solution to this differential equation dR/dt = k1*(R^k2)*(1-(R/Rmax)); R(0) = Ro k1 and k2 are parameters that need to fitted, while Ro and Rmax are the baseline and max value (which can be fitted or fixed). The response (R) increases initially at an exp

Re: [R] Solution to differential equation for nls function

2010-06-30 Thread Ravi Varadhan
t;) lines(t, y2, col=2) lines(t, y3, col=3) lines(t, y4, col=4) Hope this helps, Ravi. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of mahesh samtani Sent: Wednesday, June 30, 2010 10:44 AM To: r-help@r-project.org Subject: [R] Solution

[R] Solution to differential equation for nls function

2010-06-30 Thread mahesh samtani
Hello, I am trying to find the analytical solution to this differential equation dR/dt = k1*R (1-(R/Rmax)^k2); R(0) = Ro k1, k2, and Rmax are parameters that need to fitted, while Ro is the baseline value (which can be fitted or fixed). The response (R) increases initially at an exponential rate