R code for grenander est