# 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