求助，用R求伽马参数极大似然估计

轻灵旋风

Error in optim(c(-1.2, 1), f, dat = dat) :
  function cannot be evaluated at initial parameters
In addition: Warning message:
In log(dat) : NaNs produced