pan-genome.plot.r

[Murukarthick]

Hi,
I've tried to create the pan-genome saturation plot with Arbidopsis dataset that you provided. I get the below error when I run the fit of decay. Can you pls give some suggestions to fix this ?

model_1 <- nls(yvalues ~ Aexp(Bxvalues)+C,start=list(A=800,B=-0.3,C=24000))
Error in numericDeriv(form[[3L]], names(ind), env) :
Missing value or an infinity produced when evaluating the model

Regards,
Muru