pnbd.EstimateParameters error on big dataset

[nipnipj]

When using 20k rows

params <- pnbd.EstimateParameters(pnbd_rdf, hardie = TRUE)
Error in optimx.check(par, optcfg$ufn, optcfg$ugr, optcfg$uhess, lower,  : 
  Function returned is infinite or NA (non-computable)

When using 1k rows

Warning messages:
1: In max(logpar) : no non-missing arguments to max; returning -Inf
2: In min(logpar) : no non-missing arguments to min; returning Inf