Having difficulty fitting on 0 inflated binomal data set

[alexhallam]

I was just hoping you would have some pointers. I can't seem to get this to work for 0 inflated binomial data.

Here is a reproducible example

library(tidyverse)
library(rmetalog)
surg <- read_delim("https://hastie.su.domains/CASI_files/DATA/nodes.txt")
surg <- surg |>
  mutate(avg = x/n)
fit <- metalog(x=surg$avg, term_limit = 3, term_lower_bound = 3,
               boundedness = "b", bounds = c(0,1))

And I just get the error

Error in lpSolve::lp("min", f.obj, f.con, f.dir, f.rhs) : 
  NA/NaN/Inf in foreign function call (arg 5)

Any pointers would be very much appreciated. Thanks!

This is the PDF I was expecting to see

This is what I am getting for the PDF when I do get an output.