model
{
for(i in 1:N){
sright[i]<-(p)+(1-p)
*pow(1+pow(t[i]/alpha,beta),-1)-k[i]
sinterval[i]<-(1-p)*(pow(1+pow(t[i]/alpha,beta),-1)-pow(1+pow((t[i]+k[i])/alpha,beta),-1))
L[i]<-pow(sinterval[i],d[i])*pow(sright[i],1-d[i])
logL[i]<-log(L[i])
zeros[i]<-0
zeros[i]~dloglik(logL[i])
}
beta~dgamma(1,1)
alpha~dgamma(1,1)
p~dbeta(1,1)
}
list(t=c(5,5,5,5,6,6,6,6,6,6,7,7,7,7,7,7,7,7,18,18),d=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0),k=c(1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,0,0),N=20)