search

ArtPoon / kamphir

Phylogenetic inference using a tree-shape kernel in an Approximate Bayesian Computation framework
BSD 3-Clause "New" or "Revised" License
6 stars 2 forks source link

rcolgem: ODE error #21

Closed ArtPoon closed 9 years ago

ArtPoon commented 9 years ago 
DLSODE-  At T (=R1) and step size H (=R2), the    
      corrector convergence failed repeatedly     
      or with ABS(H) = HMIN   
In above message, R1 = 0, R2 = 7.62939e-06

DLSODE-  At start of problem, too much accuracy   
      requested for precision of machine..  See TOLSF (=R1) 
In above message, R1 = nan

Error in lsode(y, times, func, parms, mf = 10, ...) : 
  illegal input detected before taking any integration steps - see written message
Calls: simulate.binary.dated.tree.fgy ... lapply -> FUN -> .solve.Q.A.L -> ode -> lsode -> .Call
In addition: Warning messages:
1: In lsode(y, times, func, parms, mf = 10, ...) :
  repeated convergence test failures on a step, but integration was successful - inaccurate Jacobian matrix?
2: In lsode(y, times, func, parms, mf = 10, ...) :
  Returning early. Results are accurate, as far as they go
Execution halted
ArtPoon commented 9 years ago

Probably associated with integration method adams under a particular parameter range. beta was at a very low value in this case.

ArtPoon commented 9 years ago

Switching to rk4 seems to address this. Closing.