New R function refdist() to generate reference distribution for POTH for a given network structure
Rename R function cumul() to kbest()