Estoy tratando de ajustar una curva de acumulación de nutrientes y materia seca de maíz en R. Se que la curva que describe es sigmoidal y que la acumulación de materia seca o de los nutrientes depende del tiempo. Mi primera duda es como ajusto y gráfico esto en R y luego, teniendo varias curvas de distintos tratamientos, como comparo las curvas entre si. Paso algo de lo que escribi en R y algunos datos en un enlace:
dataT1 <- subset(C, Trat == "T1")
dataT3 <- subset(C, Trat == "T3")
dataT7 <- subset(C, Trat == "T7")
plotPoints(MS ~ Dia, data = dataT1)
MST1 <- nls(MS ~ a + ((b - a)/(1 + exp(-c * (Dia - d)))), data = dataT1, start = list(a = min(MS),
b = max(Dia), c = 1, d = round(median(MS))), trace = TRUE)
Al escribir esto me salta el error Error in typeof(x) : objeto 'MS' no encontrado
Database: https://www.dropbox.com/s/kxt2t6sfs52ilq3/database_dudas.xlsx?dl=0
Muchas gracias!