I'm interested to find mle estimates of lindley distribution in R. But the nlm() command returns me the initial guess. Here is the code

```
library(stats)
library(stats4)
x1<-mledata$x
ll<-function(theta){
sum(-(1+x1)-(1-theta[1])*log(1-(1+theta[2]+theta[2]*x1)/1+theta[2])-theta[2]*x1)
}
nlm(ll,theta<-c(1,1))
```