我在R中使用deSolve来解决diff eq问题。在diff eq函数中,我需要使用一个参数,该参数取决于自上次事件发生以来的时间。例如,像这样:
SAMPLE <- function(t, A, parms) {
K1 <- K10*exp(-TSLE*0.1);
dA1 = -K1*A[1] ;
list(dA1)
}
TSLE是自上次事件以来的时间。所有的想法都受到赞赏。
修改:添加了评论部分的详细信息
diff eq的一个参数是时间相关的。在事件发生之前,该值是常量。事件发生后,参数会以某种速率衰减。所以,让我们假设一个事件发生在t = 10。从时间0到9.999,参数的值为K10。从10开始,参数以一定的速率常数0.1衰减,直到最终达到零。因此,我需要能够计算自事件发生以来在diff eq中的时间。功能