A short addendum, resulting from an off-list discussion:
The reason why Colleen's code failed was raising a negative base to a
fractional exponent in the third state equation for certain sets of
parameters, esp. fractional values of beta.
Old versions of odesolve broke down, and recent versio
Hi Colleen,
this error was not uncommon and is usually a sign of a numerically
problematic or wrongly implemented model. Please use package deSolve,
the successor of odesolve, that is more robust and has also a bunch
alternative solvers for difficult cases.
I tested your code with deSolve (o
Dear list -
Does anyone have any ideas / comments about why I am receiving the following
warning when I run lsoda:
1: lsoda-- at t (=r1), too much accuracy requested in: lsoda(start, times,
model, parms)
2: for precision of machine.. see tolsf (=r2) in: lsoda(start, times,
model,
3 matches
Mail list logo
