Skip to main content
Agregado de librerías, aclaración al método para nlsLM
Origen Enlace
mpaladino
  • 6k
  • 9
  • 19
#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()

introducir la descripción de la imagen aquí

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.

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

introducir la descripción de la imagen aquí

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

introducir la descripción de la imagen aquí

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.

se añadieron 89 caracteres en el cuerpo
Origen Enlace
mpaladino
  • 6k
  • 9
  • 19
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()
library(minpack.lm)
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.tibbleas_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í PD

PD: si estás muy perdido con ggplot2 te paso un enlace a una guía básica que hice para mis alumnos.

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()
library(minpack.lm)
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í PD: si estás muy perdido con ggplot2 te paso un enlace a una guía básica que hice para mis alumnos.

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()
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í

PD: si estás muy perdido con ggplot2 te paso un enlace a una guía básica que hice para mis alumnos.

se eliminó 1 carácter en el cuerpo
Origen Enlace
mpaladino
  • 6k
  • 9
  • 19
library(minpack.nlslm)
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. 
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. 
library(minpack.lm)
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. 
se añadieron 1443 caracteres en el cuerpo
Origen Enlace
mpaladino
  • 6k
  • 9
  • 19
Loading
Origen Enlace
mpaladino
  • 6k
  • 9
  • 19
Loading