# *** Definición del turno de máxima producción continua" ***
# Basado en Bravo et al. 2015. Análisis de datos selvícolas con R.

# ** Definiciones **
# Incremento medio anual (m3/ha.año): vol. / edad
# Incremento periódico anual (m3/ha.año): (vol(t2) - vol(t1)) / (t2 - t1)

# PASO 1. Generamos los datos
# NA: no está disponible (not available)

edad <- seq(from=20, to=70, by=5)
increMe <- c(1.7, 2, 2.3, 2.5, 2.6, 2.7, 2.7, 2.7, 2.7, 2.6, NA)

volT <- edad * increMe

increPe <- c(NA, diff(volT)/5)

#PASO 2. Realizamos los gráficos
# Gráfico de rendimiento
plot(edad, volT,
type = "b",
col = "blue", 
main = "Gráfico de volumen total (rendimiento o producción)",
sub = "Betula alba/Galicia",
xlab = "Edad(años)",
ylab = "Volumen(m^3/ha)")

# Creamos un gráfico con dos curvas en el tiempo
# Incremento periódico (IPA)
plot(edad, increPe,
type = "b",
col = "black", 
main = "Turno de máxima producción continua",
sub = "Betula alba/Galicia",
ylim = c(0, 5),
xlab = "Edad (años)",
ylab = "Crec.(m^3/ha.año)")

# Crecimiento medio (IMA)
lines(edad, increMe, type="b", col="blue")

# Aumento del 10% del crecimiento medio (IMA + 10%)
lines(edad, increMe*1.10, type="b", col="green")

# Disminución del 10% del crecimiento medio (IMA - 10%)
lines(edad, increMe*0.90, type="b", col="red")