Zonación caracoles

Zonación en caracoles Playa Tamarindo

Datos desde csv

zona <- read.csv("zonacion.csv")
head(zona)

Estadísticas descriptivas para gráfica de barras

library(Rmisc)

## Loading required package: lattice

## Loading required package: plyr

zonaLang <- summarySE(zona, measurevar="Lang", groupvars="cuadrado")
zonaLang

zonaLlin <- summarySE(zona, measurevar="Llin", groupvars="cuadrado")
zonaLlin

zonaNpel <- summarySE(zona, measurevar="Npel", groupvars="cuadrado")
zonaNpel

zonaNver <- summarySE(zona, measurevar="Nver", groupvars="cuadrado")
zonaNver

Gráfica de zonación

library(ggplot2)
# Error bars represent standard error of the mean
# Lang
ggplot(zonaLang, aes(x=cuadrado, y=Lang)) +
  geom_bar(position=position_dodge(), stat="identity", fill="cornflowerblue",
           width=0.5) +
  geom_errorbar(aes(ymin=Lang-se, ymax=Lang+se),
                  width=.2,  # Width of the error bars
                  position=position_dodge(.9)) +
  scale_x_discrete("Gradiente splash-sumergido") +
  scale_y_continuous("Frecuencia Relativa")

# Llin
ggplot(zonaLlin, aes(x=cuadrado, y=Llin)) +
  geom_bar(position=position_dodge(), stat="identity", fill="green",
           width=0.5) +
  geom_errorbar(aes(ymin=Llin-se, ymax=Llin+se),
                  width=.2,  # Width of the error bars
                  position=position_dodge(.9)) +
  scale_x_discrete("Gradiente splash-sumergido") +
  scale_y_continuous("Frecuencia Relativa")

# Npel
ggplot(zonaNpel, aes(x=cuadrado, y=Npel)) +
  geom_bar(position=position_dodge(), stat="identity", fill="gray",
           width=0.5) +
  geom_errorbar(aes(ymin=Npel-se, ymax=Npel+se),
                  width=.2,  # Width of the error bars
                  position=position_dodge(.9)) +
  scale_x_discrete("Gradiente splash-sumergido") +
  scale_y_continuous("Frecuencia Relativa")

# Nver
ggplot(zonaNver, aes(x=cuadrado, y=Nver)) +
  geom_bar(position=position_dodge(), stat="identity", fill="red",
           width=0.5) +
  geom_errorbar(aes(ymin=Nver-se, ymax=Nver+se),
                  width=.2,  # Width of the error bars
                  position=position_dodge(.9)) +
  scale_x_discrete("Gradiente splash-sumergido") +
  scale_y_continuous("Frecuencia Relativa")