R:通过MLE语法错误拟合具有伽马分布的数据

时间:2018-07-05 19:28:15

标签: r mle exponential-distribution

我的ML估算存在问题:我正在使用的PDF是

我收到此错误消息:

  

NA / Inf替换为最大正值NA / Inf替换为最大值   正值NA / Inf替换为最大正值NA / Inf替换   按最大正值

代码:

x=d.frame    
splice  = function(x,beta,k){
      lambda = exp(beta)
      x_n = quantile(x,k/100)
      dens = (lambda*exp(-lambda*(x-1200000)))/(1-exp(-lambda*(x_n-1200000)))
      log_likelihood =log(sum(dens))
      return(-log_likelihood)
    }

```
MLE 
```{r}
result <- nlm(splice,2000000,hessian=TRUE,x=loss,k=95) 
result
````

0 个答案:

没有答案