Example 2: We follow [1]} and [2]} to consider the Tobit DGP as
\(y_i=(\mu _i+c\epsilon _i)I(\mu _i+c\epsilon _i>0)=({\bf x}_i^{\rm T}{\beta }+c\epsilon _i)I({\bf x}_i^{\rm T}{\beta }+c\epsilon _i>0),\quad i=1,2,\ldots ,n,\) 
