# NOT RUN {
callback <- function(fromwhere, opt, dataset, control, ...) {
# call the standard callback to print a diagnostic line
mphcrm.callback(fromwhere, opt, dataset, control, ...)
# print the distribution and two coefficients
if(fromwhere == 'full') {
print(round(mphdist(opt),6))
print(summary(opt)$coefs[c('job.alpha','job.x1'),])
}
}
# }
Run the code above in your browser using DataLab