# Para lm
Para los lm es muy fácil. Simplemente en el código de ggplot que estás usando cambias la segunda línea para agregar un mapeado de color a la variable `Trat`.

    dataT1 %>%   
      ggplot(aes(x=Dia, y=IAF, color=Trat)) +
      geom_point() +
      geom_smooth(formula = y ~ x + I(x^2)-1, 
                  method = "lm",
                  se=FALSE) +
          labs(x="Tiempo (dds)",  
          y="IAF (m2)", 
          color="Tratamiento") + 
      theme_minimal()

[![introducir la descripción de la imagen aquí][1]][1]

El resultado son tres curvas de ajuste con diferentes colores. Como lo mapeas directamente en la llamada a la función lo hereda el resto, por eso colorea los puntos de `geom_point()` y hace que `geom_smooth()` ajuste tres modelos, uno para cada nivel del factor `Trat`, y grafique las curvas con colores diferentes entre sí, pero los mismos que usa `geom_point()`. Ya sale todo etiquetado. 

# Para `nlsLM`
Para los `nlsLM` se puede con un método similar, te recomiendo visitar [esta respuesta][2] que trata exactamente ese problema. Si puedes leer inglés ahí está la solución. 

Alternativamente podés generar los puntos que ajusta el modelo y plotearlos. Va el código, probado para nls. Hay bastante uso de iteración sobre listas con `purrr::map()` y familiares. Se me complicaría explicar en esta respuesta todo lo que hace, aunque traté de aclararlo en los muchos comentarios. 
Si conoces `lapply()` es muy similar, sólo cambia ligeramente la sintaxis y es más estricto con el tipo de datos de salida. 

    library(minpack.nls)
    library(broom)        #para `augment()`
    
    MST1 %>%              #Es el df que está en el enlace de DropBox
      split(.$Trat) %>%   #Armo una lista con 3 data.frame, uno por nivel de Trat
    #map es similar a lapply, pasa una función a cada elemento de una lista y regresa una lista.
    #Uso el . (punto) como referencia a los datos porque estoy usando tuberías %>% 
    # El símbolo ~ indica a map que la siguiente expresión es una función, me ahorro declarar una función anónima como lo haría con lapply.
      map( ~nlsLM(MS ~ a + ((b - a)/(1 + exp(-c * (Dia - d)))),
                  start = list(a = min(.$MS),
                               b = max(.$MS),
                               c = 1, 
                               d = median(.$Dia)),
                  trace = TRUE, 
                  algorithm = "port", 
                  data=.)) %>% 
    #augment() crea un data.frame al "aumentado" con los datos originales, la predicción del modelo y el error
      map_df(~augment(.), .id="Trat") %>%  #.id es un argumento para map_df(), indica el nombre de la columna clave, necesario para pasar de lista a data.frame  
      as.tibble() %>%                      
      ggplot(aes(x=Dia, y=MS, color=Trat)) +   #Acá empieza el gráfico propiamente dicho, lo anterior fue prepararle los datos. 
      geom_point() +
      geom_line(aes(y=.fitted))    #a x lo hereda, mapeo un nuevo y para la línea. 
 [![introducir la descripción de la imagen aquí][3]][3]
PD: si estás muy perdido con `ggplot2` te paso un enlace a una [guía básica][4] que hice para mis alumnos. 


  [1]: https://i.sstatic.net/vOxTB.png
  [2]: https://stackoverflow.com/questions/45971661/how-to-use-method-nlslm-in-packages-minpack-lm-in-geom-smooth
  [3]: https://i.sstatic.net/USbkd.png
  [4]: https://martinpaladino.github.io/graficos_con_r_y_ggplot.html