No sé si esto es exactamente lo que estás buscando, pero el gráfico resultante está hecho con los datos que cargaste y la especificación del modelo lineal.
Librerías
library(tidyverse) #Para la estructura de datos tribble y el gráfico
Capturo los datos:
tribble(
~Trat, ~Bloque, ~Dia, ~IAF,
"T1", "A", 1, 0,
"T1", "B", 1, 0,
"T1", "C", 1, 0,
"T1", "A", 44, 0.88,
"T1", "B", 44, 0.75,
"T1", "C", 44, 0.80,
"T1", "A", 53, 2.16,
"T1", "B", 53, 1.63,
"T1", "C", 53, 1.63,
"T1", "A", 67, 1.81,
"T1", "B", 67, 1.90,
"T1", "C", 67, 1.68,
"T1", "A", 82, 2.10,
"T1", "B", 82, 1.29,
"T1", "C", 82, 1.16,
"T1", "A", 97, 2.01,
"T1", "B", 97, 1.46,
"T1", "C", 97, 1.16) -> dataT1IAF
Ajusto el modelo
En sentido estricto no es necesario, ggplot lo hará de nuevo.
IAFT1 <- lm(IAF ~ Dia + I(Dia^2)-1, data = dataT1IAF)
Genero el gráfico con ggplot:
Voy a usar dos elementos geométricos de esta librería: geom_point()
para los puntos x y y geom_smooth()
para graficar la curva de ajuste del modelo. Por defecto grafica tb la banda del intervalo de confianza. A geom_smooth()
hay que especificarle la fórmula del modelo en su idioma (i.e. usando los nombres de variable que toma ggplot y que especificamos con argumento aes()
) y el tipo de modelo que estamos ajustando (lm, glm, loess, etc.).
dataT1IAF %>% #Datos
ggplot(aes(x=Dia, y=IAF)) + #Mapeo los ejes x y y, en adelante así los llamaré.
geom_point() + #Genero los puntos, uno para cada par de x y y
geom_smooth(formula = y ~ x + I(x^2)-1, #Espefico el modelo: uso las variables internas de ggplot.
method = "lm", #Por defecto ajusta un loess, aclaro que quiero un lineal.
se=TRUE) + #Graficar el intervalo de confianza. FALSE sólo la curva de ajuste.
labs(title="Crecimiento en MS de T1", #Etiqueta de título: string encomillado.
x="Tiempo (dds)",
y="KgMS/ha",
caption="Elaboración propia con datos de \nla n antecedida por una barra invertida\n marca un corte de línea.") +
theme_minimal() #Un tema simple para los gráficos.
labs()
con las anotaciones y señalé una opción para desactivar la banda de ic. Si esta es la respuesta que estabas buscando podés aceptarla.