boxplot.fnc

Copy, Paste and Adapts

boxplot.fnc(dat, dv=3)

boxplot.fnc(dat, dv=3, which.factor='zone')

boxplot.fnc(dat, dv='rt', identification=T)

Objetive

Permite valorar la distribución de una variable cuantitativa mediante el diagrama de cajas. La gráfica muestra un rectángulo que contiene en su interior al 50% centrado de la distribución con sus lados delimitados por el valor correspondiente al percentil 25 y 75 respectivamente. La mediana se representa mediante una línea gruesa de color negro. Los "bigotes" delimitan aproximandamente 1.5 la distancia intercuatílica (P75-P25). Asimismo toda observación que quede fuera de la superficie delimitada por estos marcadores (bigotes) serán marcados con un círculo.

Boxplot

Para demostrar la facilidad de uso de esta función, generaremos unos datos artificiales para 200 observaciones en la variable medida para dos grupos (A y B).

dat=data.frame(measure=c(rnorm(100,100,15),rnorm(100,90,22)),

group=factor(rep(c('A','B'),each=100)),

gender=factor(rep(c('H','M'),each=50)))

measure group gender

1 111.5089 A H

2 109.8613 A H

3 106.4138 A H

4 117.2250 A H

5 117.4536 A H

6 106.7223 A H

1. boxplot.fnc(dat, dv='measure')

2. boxplot.fnc(dat, dv='measure', which.factor='group')

3. boxplot.fnc(dat, dv='measure', which.factor='group:gender')

Observa en la tercera llamada a la función que hemos incluido la "interacción" (grupo:genero) de los dos factores en el argumento que.factor. Esto provoca que se cree una sola gráfica con los diagramas de cajas correspondientes a las cuatro condiciones del cruce de ambos factores.

OUTLIERS IDENTIFICATION

Ahora provocaremos que los casos 50 y 80 presenten valores bajos y claramente distantes de la media del grupo A.

dat[50,1]=45

dat[80,1]=58

Veamos como "saber" el número del caso o registro al que pertenece un valor extremo en un diagrama de cajas. Lo conseguiremos incluyendo el argumento identification=T.

boxplot.fnc(dat, dv='measure', which.factor='group',

identification=T)

Cuando incluimos el argumento identification=T, una vez creado el gráfico aparecerá una cruz que al situarla encima de los casos extremos y al hacer clic con el ratón revelará el número de registro al que pertenece. Dado que cada muestra se genera al azar, tus resultados gráficos no serán idénticos a los aquí presentados.