R Scripts/temporal.R


1
2
3
4
5
6
7
8
9
10

temporal = function(tu,tv,alpha){
  dt = tv - tu - 1
  if (sum(dt<0)==0) {
    return (alpha*exp(-alpha*dt))
#     return ((0.01/alpha)*(1+dt/alpha)^(-1.01))
  } else {
    return (NA)
  }
}