r - 通过 R 中的“泊松”包生成非齐次泊松过程事件

Question

我正在尝试使用nhpp.sim()R 中“泊松”包中的函数，但是，我得到的返回值超出了 24 小时的范围。函数方法是nhpp.sim(rate, num.events, prob.func, num.sims = 1, t0 = 0, prepend.t0 = T) 据我了解，rate是事件时间之间的平均值，num.events是事件的数量，prob.func是强度函数（事件在那个时间段的概率）。我正在尝试根据强度函数在 24 小时内模拟 10 个事件

intensity  <- function(t) {    
  if (t >= 0 && t < 1) {lambda = 0.0100456621004566}   
  else if (t >= 1 && t < 2) {lambda = 0.00650684931506849}    
  else if (t >= 2 && t < 3) {lambda = 0.00707762557077626}        
.

. 

.  
else if (t >= 23 && t < 24) {lambda = 0.0219178082191781}    
  else if (t >= 24) {lambda = 0}

  return(lambda)    
}

然后我通过以下方式调用该函数：nhpp.sim(0.4316, 10, prob.func=intensity,1,0,24) 并且我得到的事件时间超出了 24 小时的限制，这是为什么呢？