R/Bursts.R

In nunesmatt/binhf: Haar-Fisz Functions for Binomial Data

`Bursts` <-
function(t){

I1 <- exp(-(abs(t-0.2)/0.01)^1.2)*(t<=0.2) + exp(-(abs(t-0.2)/0.03)^1.2)*(t>0.2)
I2 <- exp(-(abs(t-0.3)/0.01)^1.2)*(t<=0.3) + exp(-(abs(t-0.3)/0.03)^1.2)*(t>0.3)
I3 <- exp(-(abs(t-0.4)/0.01)^1.2)*(t<=0.4) + exp(-(abs(t-0.4)/0.03)^1.2)*(t>0.4)

y<- (4*I1+3*I2+4.5*I3)

y<-y/max(y)

y

}