# NOT RUN {
cbs <- cdnow.sample()$cbs # load CDNow summary data
params <- bgcnbd.EstimateParameters(cbs)
palive <- bgcnbd.PAlive(params, cbs$x, cbs$t.x, cbs$T.cal)
mean(palive) # Estimated share of customers to be still alive
# }
Run the code above in your browser using DataLab