random_start: Number of random starting points, default is 0.
maxit: Maximum number of iteration.
Returns
A list contains: theta is the estimate causal effect, se is standard error of estimated theta, l is negative log-likelihood, r_est is estimated r vector.