Gráficos múltiples de barras de error
Tuesday, June 8th, 2010Aquí os dejo el código para hacer un gráfico múltiple de barras de error.
Son las estimaciones con sus desviaciones típicas de dos parámetros de la ecuación de crecimiento de von Bertalanffy para erizos de mar del submareal.
Tuve que hacerlo por partes para que los plots de L8 y K me salieran con el mismo ancho.
Aquí teneis los datos originales por si quereis hacer la prueba:
l8_k_estimates.txt
Y aquí va el código:
figure4 <- read.csv("l8_k_estimates.txt",header=T, sep=" ") pdf("figure4b.pdf", width=8, family="Times") op <- par(mfrow=c(1,3), mai=c(0.5,0,0.5,0.2)) # Se añaden los labels plot.new() axis(2,at=seq(0.01,0.99, length.out=19), labels=figure4$name, tick= F, las=2, pos=c(1.1,1)) # Se añaden los estimadores de l8 dotchart(figure4$l8.value, pch=20, xlab="L8", color="grey20", xlim=c(min(figure4$l8.value)-max(figure4$l8.std),max(figure4$l8.value)+max(figure4$l8.std))) segments(figure4$l8.value,1:19,figure4$l8.value+figure4$l8.std,1:19,col="grey20") segments(figure4$l8.value,1:19,figure4$l8.value-figure4$l8.std,1:19,col="grey20") # Se añaden los estimadores de k dotchart(figure4$k.value, pch=20, xlab="K", color="grey20", xlim=c(min(figure4$k.value)-max(figure4$k.std),max(figure4$k.value)+max(figure4$k.std))) segments(figure4$k.value,1:19,figure4$k.value+figure4$k.std,1:19,col="grey20") segments(figure4$k.value,1:19,figure4$k.value-figure4$k.std,1:19,col="grey20") par(op) dev.off()