Re: [R] ODE's in R

2008-01-23 Thread Spencer Graves
Markku Karhunen wrote: > Thank you all. > > We must think about implementing these packages. In the meantime, I > should clarify my question: Is there any evidence that doing the dumb > for loop discretisation is any more dangerous in R, than in any other > language? Apparently not? I know of

Re: [R] ODE's in R

2008-01-23 Thread Markku Karhunen
Thank you all. We must think about implementing these packages. In the meantime, I should clarify my question: Is there any evidence that doing the dumb for loop discretisation is any more dangerous in R, than in any other language? Apparently not? Best, Markku Karhunen > have you looked at ls

Re: [R] ODE's in R

2008-01-22 Thread Spencer Graves
have you looked at lsoda{odesolve}? have you looked at the scripts\CSTR subdirectory in the fda package? it includes an example worked in both R and Matlab with slightly better answers in R but with a much longer compute time. sg The fda package Peter Dalgaard wrote: > Markku Karhunen wrote:

Re: [R] ODE's in R

2008-01-22 Thread Peter Dalgaard
Markku Karhunen wrote: > Thanks, Dr. Maechler. > >> No, there's no such track. >> [ Matlab users coming to R may produce wrong R code >> by using 0:n-1 instead of 0:(n-1) ; but I don't assume this >> would be the case ] >> >> >> > Been there, done that! > >> MK> We use just a

Re: [R] ODE's in R

2008-01-22 Thread Markku Karhunen
Thanks, Dr. Maechler. > No, there's no such track. > [ Matlab users coming to R may produce wrong R code > by using 0:n-1 instead of 0:(n-1) ; but I don't assume this > would be the case ] > > Been there, done that! > MK> We use just a simple discretisation written in a for loop >

Re: [R] ODE's in R

2008-01-22 Thread Martin Maechler
> "MK" == Markku Karhunen <[EMAIL PROTECTED]> > on Tue, 22 Jan 2008 10:57:25 +0200 writes: MK> Dear all, MK> Has anybody tried numerical solving of ODE's and/or transport equations MK> in R? (Don't ask how we ended up in using R for this job, in the first MK> place!)

[R] ODE's in R

2008-01-22 Thread Markku Karhunen
Dear all, Has anybody tried numerical solving of ODE's and/or transport equations in R? (Don't ask how we ended up in using R for this job, in the first place!) More precisely, does anybody know any technical issue that could make the work insecure in the sense of propagation of errors? Is the