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))
You should use a reprex for your example, that way you would know that your code as is is not runable. mledata is not defined.
mledata
You should include meldata or a toy equivalent so that we can run your code and see the results your are getting.
meldata
Here is a pointer to reprex's https://www.tidyverse.org/help/
We really need concrete, run-able examples in a reprex to be able to help you.
If you are having any trouble getting a reprex to work please post a new topic with the issues you are having. You will get a lot of help getting it to work for you.
Thanks