estoy tratando de ajustar una curva sigmoidal a varios conjuntos de datos, por lo cual estoy usando este script: plotPoints(MS ~ Dia, data = dataT1)
MST1 <- nls(y ~ a + ((b - a)/(1 + exp(-c * (x - d)))),
data = dataT1,
start = list(a = min(dataT1$y),
b = max(dataT1$y),
c = 1, d = median(dataT1$x)),
trace = TRUE, algorithm = "port")
overview(MST1)
plotfit(MST1, smooth = TRUE)
El problema que estoy teniendo es que con varios grupos de datos el script me funciona pero con otros no, dandome el siguiente error:
Error in nls(Nabs ~ a + ((b - a)/(1 + exp(-c * (Dia - d)))), data = dataT1NPK, :
Convergence failure: false convergence (8)
y
Error in nls(Kabs ~ a + ((b - a)/(1 + exp(-c * (Dia - d)))), data = dataT1NPK, :
Convergence failure: iteration limit reached without convergence (10)
¿A qué se podría deber esto? Les dejo uno de los database que me da error: https://www.dropbox.com/s/f8srl24xmd7jw0v/Duduas_database.xlsx?dl=0
Muchas gracias!