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
.
#Librerías necesarios para los dos bloques de código.
library(tidyverse) #Para map(), map_df(), as_tibble(), ggplot()
library(minpack.lm) #Para nlsLM()
library(broom) #Para augment()
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()
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 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.
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.
Atención: el código para NLM no hace exactamente lo mismo que el código para lm. Solamente plotea los puntos en los que hay una estimación del valor de y, geom_line()
se encarga de unirlos y crear la línea. Por eso no se ve como una curva y son visibles los saltos. Se puede corregir usando predict()
con más data points de x, de ese modo se obtienen más puntos predichos de y en la gráfica se asemeja más a una curva. Aunque en sentido matemático estricto siguen siendo puntos discretos al tener suficientes puntos y usar la función de suavizado de geom_line()
se ve como si fuera continua.
PD: si estás muy perdido con ggplot2
te paso un enlace a una guía básica que hice para mis alumnos.