Gráfico de dispersión con barras de error

Una alumna me preguntó el otro día cómo se podían añadir con R unas barras de error horizontales a los puntos de un diagrama de dispersión. He escrito una función muy simple para hacerlo y la dejo aquí por si acaso a alguien más le viene bien. El argumento error corresponde al radio de las barras.


PlotError <- function(x, y, error = 0){
 x.inf <- x - error
 x.sup <- x + error
 plot(x,y, pch=16, bty='l')
 segments(x.inf, y, x.sup, y, col='blue')
}

# Ejemplo
n <- 20
x <- rnorm(n)
y <- x + rnorm(n, sd=0.5)

PlotError(x, y) # igual que plot, con algunas opciones gráficas diferentes
PlotError(x, y, error = 0.15)

El resultado de los ejemplos debe ser parecido a los gráficos siguientes:

ejemplo-barras

Anuncios
Esta entrada fue publicada en estadística, R y etiquetada . Guarda el enlace permanente.

Responder

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s