将方差Gamma分布拟合到数据

Question

x<-c(1291,849,238,140,118,108,87,70,63,58,50,47,21,21,19)
library(VarianceGamma)
init<-c(0,0.5,0,0.5)
vgFit(x, freq = NULL, breaks = NULL, paramStart = init, startMethod = "Nelder-Mead", startValues = "SL", method = "Nelder-Mead", hessian = FALSE, plots = TRUE)

我想要拟合数据，但它显示错误。

Error in optim(paramStart, llfunc, NULL, method = "Nelder-Mead", hessian = hessian,  : 

  function cannot be evaluated at initial parameters.

我不知道错误意味着什么，因为我使用了可以在VG分发pdf中使用的初始值。   我上网了，有人说它可能有零点，不能在r。

中最小化