################################################################################
### R BASICS WORKSHOP                                                        ###
### PRESENTACIÓN 9: CONTROL DE FLUJO                                         ### 
###                                                                          ###
### Center for Conservation and Sustainable Development                      ###
### Missouri Botanical Garden                                                ###
### Sitio en la red: rbasicsworkshop.weebly.com                              ### 
################################################################################

### INTRODUCCION ###############################################################

# En R hay una serie de elementos que permiten controlar el flujo del código.
# Hay 4 tipos principales maneras de controlar flujo:

# 1. Bucles (loops) *for*
# 2. Bucles (loops) *while*
# 3. Condicionales 
# 4. Rupturas (breaks)

# Esta presentación se enfoca en *for*, *while* e *if*

## IMPORTANTE: para ayuda con control de flujo:
help(Control)


#########################
## Un ejemplo práctico ##
#########################

# El objetivo del ejemplo es crear un mapa que muestra la distribución de 
# especies de árboles a través de una parcela de bosque. Cada individuo se 
# representa con un punto, y cada especie con un color diferente.

# Abra una versión de los datos de la Parcela De Dinúmicas Forestales de Tyson. 
# Este conjunto de datos contiene información sobre la identidad y distribución 
# de individuos de diferentes especies árboles a través de un área de 25 hectáreas
# en el Centro de  Investigación Tyson cerca de St. Louis, MO. (Los datos han sido 
# aleatorizados. Gracias a Jonathan A. Myers por los datos).

tyson <- read.table(file=file.choose(), header=TRUE, sep="\t")

# Compruebe que los datos est?n bien abiertos
dim(tyson)
head(tyson)


# Ahora, vamos a hacer un mapa de la distribución de las especies de árboles en 
# la parcela de bosque

# OPCIÓN 1: SIN BUCLES #

sp.list <- unique(tyson$spcode)
colors <- adjustcolor(rainbow(length(sp.list)), alpha.f=0.8)

tiff(filename="TysonForestMap_1.tif", width=10, height=10, units="in", res=300)

  # Crea una gráfica vacía
  plot(tyson$gx, tyson$gy, xlab="x", ylab="y", asp=1, type="n", cex.lab=1.5,
       cex.axis=1.5)

 # Grafica los árboles de la especie "ulmrub"
  points(tyson$gx[tyson$spcode=="ulmrub"], tyson$gy[tyson$spcode=="ulmrub"],
         pch=16, col=colors[1])

  # Grafica los ?árboles de la especies "lonmaa"
  points(tyson$gx[tyson$spcode=="lonmaa"], tyson$gy[tyson$spcode=="lonmaa"],
         pch=16, col=colors[2])

  # etc.
  points(tyson$gx[tyson$spcode=="asitri"], tyson$gy[tyson$spcode=="asitri"],
         pch=16, col=colors[3])

  points(tyson$gx[tyson$spcode=="acerub"], tyson$gy[tyson$spcode=="acerub"],
         pch=16, col=colors[4])

  points(tyson$gx[tyson$spcode=="diovir"], tyson$gy[tyson$spcode=="diovir"],
         pch=16, col=colors[5])

  points(tyson$gx[tyson$spcode=="linben"], tyson$gy[tyson$spcode=="linben"],
         pch=16, col=colors[6])

  points(tyson$gx[tyson$spcode=="cercan"], tyson$gy[tyson$spcode=="cercan"],
         pch=16, col=colors[7])

  points(tyson$gx[tyson$spcode=="celocc"], tyson$gy[tyson$spcode=="celocc"],
         pch=16, col=colors[8])

  points(tyson$gx[tyson$spcode=="unknown"], tyson$gy[tyson$spcode=="unknown"],
         pch=16, col=colors[9])

  points(tyson$gx[tyson$spcode=="plaocc"], tyson$gy[tyson$spcode=="plaocc"],
         pch=16, col=colors[10])

  points(tyson$gx[tyson$spcode=="pruser"], tyson$gy[tyson$spcode=="pruser"],
         pch=16, col=colors[11])

  points(tyson$gx[tyson$spcode=="ostvir"], tyson$gy[tyson$spcode=="ostvir"],
         pch=16, col=colors[12])

  points(tyson$gx[tyson$spcode=="carcor"], tyson$gy[tyson$spcode=="carcor"],
         pch=16, col=colors[13])

  points(tyson$gx[tyson$spcode=="tilame"], tyson$gy[tyson$spcode=="tilame"],
         pch=16, col=colors[14])

  points(tyson$gx[tyson$spcode=="prusp"], tyson$gy[tyson$spcode=="prusp"],
         pch=16, col=colors[15])

  points(tyson$gx[tyson$spcode=="celten"], tyson$gy[tyson$spcode=="celten"],
         pch=16, col=colors[16])

  points(tyson$gx[tyson$spcode=="quesp"], tyson$gy[tyson$spcode=="quesp"],
         pch=16, col=colors[17])

  points(tyson$gx[tyson$spcode=="morrub"], tyson$gy[tyson$spcode=="morrub"],
         pch=16, col=colors[18])

  points(tyson$gx[tyson$spcode=="malioe"], tyson$gy[tyson$spcode=="malioe"],
         pch=16, col=colors[19])

  points(tyson$gx[tyson$spcode=="cordru"], tyson$gy[tyson$spcode=="cordru"],
         pch=16, col=colors[20])

  points(tyson$gx[tyson$spcode=="aceneg"], tyson$gy[tyson$spcode=="aceneg"],
         pch=16, col=colors[21])

  points(tyson$gx[tyson$spcode=="quemue"], tyson$gy[tyson$spcode=="quemue"],
         pch=16, col=colors[22])

  points(tyson$gx[tyson$spcode=="junvir"], tyson$gy[tyson$spcode=="junvir"],
         pch=16, col=colors[23])

  points(tyson$gx[tyson$spcode=="carsp"], tyson$gy[tyson$spcode=="carsp"],
         pch=16, col=colors[24])

  points(tyson$gx[tyson$spcode=="quealb"], tyson$gy[tyson$spcode=="quealb"],
         pch=16, col=colors[25])

  points(tyson$gx[tyson$spcode=="querub"], tyson$gy[tyson$spcode=="querub"],
         pch=16, col=colors[26])

  points(tyson$gx[tyson$spcode=="quevel"], tyson$gy[tyson$spcode=="quevel"],
         pch=16, col=colors[27])

  points(tyson$gx[tyson$spcode=="corflo"], tyson$gy[tyson$spcode=="corflo"],
         pch=16, col=colors[28])

  points(tyson$gx[tyson$spcode=="quemar"], tyson$gy[tyson$spcode=="quemar"],
         pch=16, col=colors[29])

  points(tyson$gx[tyson$spcode=="cartex"], tyson$gy[tyson$spcode=="cartex"],
         pch=16, col=colors[30])

  points(tyson$gx[tyson$spcode=="amearb"], tyson$gy[tyson$spcode=="amearb"],
         pch=16, col=colors[31])

  points(tyson$gx[tyson$spcode=="cargla"], tyson$gy[tyson$spcode=="cargla"],
         pch=16, col=colors[32])

  points(tyson$gx[tyson$spcode=="carova"], tyson$gy[tyson$spcode=="carova"],
         pch=16, col=colors[33])

  points(tyson$gx[tyson$spcode=="fraame"], tyson$gy[tyson$spcode=="fraame"],
         pch=16, col=colors[34])

  points(tyson$gx[tyson$spcode=="cartom"], tyson$gy[tyson$spcode=="cartom"],
         pch=16, col=colors[35])

  points(tyson$gx[tyson$spcode=="queste"], tyson$gy[tyson$spcode=="queste"],
         pch=16, col=colors[36])

  points(tyson$gx[tyson$spcode=="pruame"], tyson$gy[tyson$spcode=="pruame"],
         pch=16, col=colors[37])

  points(tyson$gx[tyson$spcode=="jugnig"], tyson$gy[tyson$spcode=="jugnig"],
         pch=16, col=colors[38])

  points(tyson$gx[tyson$spcode=="sasalb"], tyson$gy[tyson$spcode=="sasalb"],
         pch=16, col=colors[39])

  points(tyson$gx[tyson$spcode=="ailalt"], tyson$gy[tyson$spcode=="ailalt"],
         pch=16, col=colors[40])

  points(tyson$gx[tyson$spcode=="crasp"], tyson$gy[tyson$spcode=="crasp"],
         pch=16, col=colors[41])

  points(tyson$gx[tyson$spcode=="vibruf"], tyson$gy[tyson$spcode=="vibruf"],
         pch=16, col=colors[42])

  points(tyson$gx[tyson$spcode=="fracar"], tyson$gy[tyson$spcode=="fracar"],
         pch=16, col=colors[43])

  points(tyson$gx[tyson$spcode=="sidlan"], tyson$gy[tyson$spcode=="sidlan"],
         pch=16, col=colors[44])

  points(tyson$gx[tyson$spcode=="gletri"], tyson$gy[tyson$spcode=="gletri"],
         pch=16, col=colors[45])

  points(tyson$gx[tyson$spcode=="acesac"], tyson$gy[tyson$spcode=="acesac"],
         pch=16, col=colors[46])

dev.off()

# ¿Cómo se debería modificar el código anterior para cambiar el tamaño de los 
# símbolos que representan a cada árbol?



# OPCIÓN 2: SI POR FAVOR, NECESITO UN BUCLE! #

sp.list <- unique(tyson$spcode)
colors <- adjustcolor(rainbow(length(sp.list)), alpha.f=0.8)

tiff(filename="TysonForestMap_2.tif", width=10, height=10, units="in", res=300)

  plot(tyson$gx, tyson$gy, xlab="x", ylab="y", asp=1, type="n", cex.lab=1.5,
       cex.axis=1.5)
  
  for(i in 1:length(sp.list))
    points(tyson$gx[tyson$spcode==sp.list[i]], tyson$gy[tyson$spcode==sp.list[i]],
           pch=16, col=colors[i])

dev.off()

# ?C?mo se debería modificar el código anterior (con el bucle) para cambiar el 
# tamaño de los símbolos que representan a cada árbol?


################################################################################
### 1. BUCLES *for* ############################################################
################################################################################

# Un bucle permite repetir un pedazo de código múltiples veces sin tener que 
# repetirlo.

# *for* es la manera mas común de construir bucles. Este tipo de bucle repite 
# un pedazo de código un numero determinado de veces.


# La estructura general de un bucle *for* es la siguiente:
#
# for(i in vector)
# {
#   code
# }


# Esto quiere decir aproximadamente:
#
# para cada valor que i toma del vector repetir
# {
#   este código
# }


#####################
## Ejemplo fácil 1 ##
#####################

for(i in 1:10)
{
    print(i)
}

# Esto se traduce como:
#
# Crear una secuencia del 1 al 10
# 
# Para cada valor *i* de la secuencia hacer lo siguiente:
# imprimir el valor de i
#


#####################
## Ejemplo fácil 2 ##
#####################

letters

for(i in letters)
{
    print(i)
}


#####################
## Ejemplo fácil 3 ##
#####################

length(letters)

result <- 0

for(i in letters)
{
	print(i)
	result <- result + 1
}

result


#####################
## Ejemplo fácil 4 ##
#####################

v <- c(1,3,5,2,4)

result <- 0

for(i in 1:length(v))
{
	print( c(i, v[i]) )

	result <- result + v[i]
}

result


#####################
## Ejemplo fácil 5 ##
#####################

# La secuencia de Fibonacci es una secuencia famosa en matemáticas. Los primeros 
# dos elementos son 1 y 1. Los elementos posteriores se definen como la suma de 
# los dos inmediatamente anteriores. Por ejemplo, el tercer elemento es 2 
# (1 + 1), el cuarto elemento es 3 (2 + 1), y así sucesivamente. En este ejemplo, 
# vamos a calcular los primeros 'n' números en la secuencia de Fibonacci.

# Esto crea una variable que determina la longitud de la secuencia de Fibonacci:
n <- 25

# A menudo es útil para crear un objeto con valores NA que almacenará los
# valores creados en cada iteración de un bucle. En este caso, creamos un
# vector de NAs de longitud 'n':
fibonacci <- rep(NA, times=n)

# Comprobamos el contenido de 'fibonacci':
fibonacci

# Por definición, los primeros dos elementos de la secuencia son 1 y 1:
fibonacci[1] <- 1
fibonacci[2] <- 1

# Este bucle calculará los elementos 3 a 'n' de la secuencia:
for(i in 3:n)
{
  # El elemento 'i' se calcula como la suma de los elementos'i-1' e 'i-2'
  fibonacci[i] <- fibonacci[i-1] + fibonacci[i-2]
}

# Ahora podemos ver el resultado del bucle - la secuencia de Fibonacci de 1 al 
# elemento "n":
fibonacci