I am using the package of Icens but for the function Bisect I do not know what ndir is it just says the direction of pvec, but what kind of values should I use ? thanks for reply pls give me an example.. Bisect(tA, pvec, ndir, Meps, tolbis=1e-07) Arguments: tA: The transpose of the clique matrix. pvec: The current estimate of the probability vector. ndir: The direction to explore. Meps: Machine epsilon, elements of 'pvec' that are less than this are assumed to be zero. tolbis: The tolerance used to determine if the algorithm has converged.