Inference for cub models: A program in R