#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()
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.