R/prediction.average.R

In Auburngrads/SMRD: Statistical Methods For Reliability Data

prediction.average <-
function (timevec, 
          mu, 
          sigma, 
          distribution) 
{
    zout <- POSTPR(as.double(timevec), 
                   as.integer(length(timevec)), 
                   as.double(mu), 
                   as.double(sigma),
                   as.integer(length(mu)),
                   as.integer(numdist(distribution)), 
                   pdf = double(length(timevec)), 
                   cdf = double(length(timevec)))
    
    return(list(x = timevec, 
                post.pred.den = zout$pdf, 
                post.pred.cdf = zout$cdf))
}