Objetivo

Modelar el desarrollo de área foliar del cultivo de la Caña de Azúcar (Saccharum officinarum) usando la temperatura de Neiva (Huila).

Procedimiento

Para la realización de los cálculos se tuvo en cuenta una base de datos, las cuales contenían entre otras variables la temperatura máxima y mínima, importantes para poder cumplir el objetivo.

Pasos implementados

  1. Inicialmente se hizo la lectura del dataframe y la actualización de los encabezados:
getwd()
## [1] "C:/Users/TEMP"
setwd("C:/Users/Carlos Mojocó/Downloads/Modelacion/Rstudio/Ejercicio_1")
datosNeiva = read.table("Neiva_Huila.csv", header = TRUE, sep = ",")
names(datosNeiva)[2:7]=c("Día/Mes","Mes/Año","R. Solar", "Tmax (°C)", "Tmin (°C)", "Precipitación (mm/día)")
  1. Una vez hecha la lectura de los datos se creó un nuevo dataframe que contenía los datos desde el 30 de marzo de 2010 hasta el 27 de junio de 2010:
nuevos.datosNeiva = datosNeiva[89:178,]
  1. Luego de separar la data se inició con los cálculos necesarios para conocer el desarrollo de área foliar del cultivo de la Caña de Azúcar, en el primer cálculo se creó una nueva columna en el nuevo dataframe para agregar los grados día o tiempo térmico (DTU):
nuevos.datosNeiva$DTU = ((nuevos.datosNeiva$`Tmax (°C)`+nuevos.datosNeiva$`Tmin (°C)`)/2)-10
  1. Seguidamente se calculó el desarrollo diario de los nudos de las hojas (INODE). En esta ecuación es importante tener en cuenta que el PHYL para el cultivo de la Caña de Azúcar está entre 100 y 120 °C, pero por consideraciones del ejercició se definió que el PHLY era de 100 °C.
nuevos.datosNeiva$INODE = (nuevos.datosNeiva$DTU/100)
  1. Una vez obtenido el cálculo de desarrollo diario de los nudos de las hojas (INODE), se calculó el número total de nudos en la planta (MSNN) cada día:
nuevos.datosNeiva$MSNN = nuevos.datosNeiva$INODE[1]
i = 2

while (i !=91){
  nuevos.datosNeiva$MSNN[i] = nuevos.datosNeiva$MSNN[i-1]+nuevos.datosNeiva$INODE[i]
  i=i+1
}
  1. Cálculo del área foliar de la planta (PLA) en m^2 por planta:
a = 0.0235
b = 0.166

nuevos.datosNeiva$PLA = round(a*exp(b*nuevos.datosNeiva$MSNN),3)
  1. índice de área foliar (LAI) En el caso de este ejercicio se tuvo en cuenta 10 plantas / m^2 como densidad del cultivo.
PDEN = 10

nuevos.datosNeiva$LAI = round(nuevos.datosNeiva$PLA*PDEN,3)

Resultados

DT::datatable(nuevos.datosNeiva)
plot(nuevos.datosNeiva$MSNN, nuevos.datosNeiva$PLA, col="dark green", main="Número de nudos Vs. Área foliar en el Nieva - Huila", xlab= "Número de nudos", ylab="Área foliar")

plot(nuevos.datosNeiva$MSNN, col="dark green", main="Formación acumulada de nodulos en Neiva - Huila", xlab= "Días desde el 30/03/2010 hasta el 27/06/2010", ylab="MSNN")

Conclusiones