schneebergerlab / AMPRIL-genomes

scripts for the project of seven thaliana genomes assembly
38 stars 18 forks source link

pan-genome.plot.r #9

Open Murukarthick opened 3 years ago

Murukarthick commented 3 years ago

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