# Have left endpoints of intervals in a vector
# ``lefts'' and right endpoints in ``rights''and
# set ``N'' equal to the common length.
# Partition support (0,infty) into sets B
grid <- c( seq(0,100,by=1), 1000 )
ngrid <- length(grid)
grt <- grid[2:ngrid]; glt <- grid[1:(ngrid-1)]
# Exponential prior guess and prior sample size
gg <- exp( -glt/(365/2) ) - exp( -grt/(365/2) )
alpha <- 1
# Identify partition elements that may contain
# each survival time
less <- function(x,y){ return( x