add cumprob.R