Função Gnls

Olá, então to com o seguinte problema no R:

Tenho este modelo abaixo:

Gompertz=nls(Ht~SSgompertz(Dap,b0,b1,b2))
summary(Gompertz)

Com os seguintes resultados para os seus parâmetros

Formula: Ht ~ SSgompertz(Dap, b0, b1, b2)

Parameters:
Estimate Std. Error t value Pr(>|t|)
b0 16.35440 0.98953 16.53 <2e-16 ***
b1 1.67231 0.67153 2.49 0.0133 *
b2 0.81996 0.05814 14.10 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 4.704 on 310 degrees of freedom

Number of iterations to convergence: 7
Achieved convergence tolerance: 9.538e-07

Porem quando incluo a covariavel Idade (I) aparece o seguinte erro:

Gompertz.covariavel=gnls(Ht~SSgompertz(Dap,b0,b1,b2),params=list(b0~I,b1~I,b2~I),start=c(b0=16,0,b1=2,0,b2=1,0))

Erro em gnls(Ht ~ SSgompertz(Dap, b0, b1, b2), params = list(b0 ~ I, :
Step halving factor reduced below minimum in NLS step

Alguem poderia me dar algumas dicas