modelling of events

[fkrauer]

Hi

I have just started using diffeqr for R and I would love to use it for my (rather large) disease modelling project. I will have an ODE model with a substantial number of compartments and I have to trigger changes the values of states or parameters at specific time points (externally forced, corresponding to an intervention). In plain R and deSolve() I would do this with an event function. Is there an equivalent for diffeqr? And if yes, is example code available?

Bw, Fabienne

[ChrisRackauckas]

DifferentialEquations.jl has the callback functions defined here: https://diffeq.sciml.ai/stable/features/callback_functions/ . I haven't tried using them from R though. My guess is that they might have a few issues since they rely on mutation and R uses call-by-copy (https://stackoverflow.com/questions/2603184/can-you-pass-by-reference-in-r). So I think we'll need to expose an alternative non-mutating form for R users. Until then though, the solution would be to eval a callback function.