note q nao entro na discussao da validade do modelo.... mas a questao eh a selecao dos valores iniciais... b blah = simple[complete.cases(simple),] nls(log(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=1.68,b=-24.599,c=-0.026), data=blah) curve(-10.48+6.79*exp(.003915*x), 90, 140) with(blah, points(Gest, log(PESO.FETO)))