Post

1 follower Follow
0
Avatar

MLE estimation for lindley distribution

Hello I'm interested to find out mle estimates for lindley distribution for my wind speed data. 

 

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))
```

Harsh