Medicina Basada en la Evidencia

Estadística básica – 524 – Anexo 1. Histogramas: procedimiento en RCommander/KMggplot2 y Script Anexo 2. Gráfico de tallo y hojas: procedimiento en RCommander y Script Figura 1 . En RCommander, desde KMggplot2 acceder a Histogram - Variable: Talla -Horizontal axis label: Talla en cm - Vertical axis label: Frecuencia absoluta - No of bins: Sturges - Options: Heat map - Axis scaling: Frequency counts - Colour pattern: Paired R script Figura 1 : require(“ggplot2”) .df <- na.omit(data.frame(x = Funda_graf$Talla)) .nbins <- pretty(range(.df$x), n = nclass.Sturges(.df$x), min.n = 1) .plot <- ggplot(data = .df, aes(x = x, y = ..count..)) + geom_histogram(aes(fill = ..count..), breaks = .nbins) + scale_y_continuous(expand = c(0.01, 0)) + scale_fill_gradient(low = RColorBrewer::brewer.pal(3, “Paired”)[2], high = RColorBrewer::brewer.pal(3, “Paired”)[1]) + xlab(“Talla en cm”) + ylab(“Frecuencia absoluta”) + theme_bw(base_size = 12, base_family = “sans”) + theme(legend.position = “none”) print(.plot) rm(.df, .nbins, .plot) Figura 2 . En RCommander, desde KMggplot2 acceder a Histogram - Variable: Talla -Facet variable in rows: Sexo - Horizontal axis label: Talla en cm - Vertical axis label: Frecuencia absoluta - No of bins: Sturges - Options: Heat map - Axis scaling: Frequency counts - Colour pattern: Paired R script Figura 2: require(“ggplot2”) .df <- na.omit(data.frame(x = Funda_graf$Talla, s = Funda_graf$Sexo)) .nbins <- pretty(range(.df$x), n = nclass.Sturges(.df$x), min.n = 1) .plot <- ggplot(data = .df, aes(x = x, y = ..count..)) + geom_histogram(aes(fill = ..count..), breaks = .nbins) + scale_y_continuous(expand = c(0.01, 0)) + scale_fill_gradient(low = RColorBrewer::brewer.pal(3, “Paired”)[2], high = RColorBrewer::brewer.pal(3, “Paired”)[1]) + facet_wrap( ~ s) + xlab(“Talla en cm”) + ylab(“Frecuencia absoluta”) + theme_bw(base_size = 12, base_family = “sans”) + theme(panel.spacing = unit(0.3, “lines”), legend.position = “none”) print(.plot) rm(.df, .nbins, .plot) Figura 3 . En RCommander accederemos a Gráficas-Gráfica de tallos y hojas-Datos: Talla -Opciones: Partes del tallo: Automático-Estilo de división de los tallos: Tukey R script Figura 3 : library(aplpack, pos=25) with(Funda_graf, stem.leaf(Talla, unit=1, na.rm =TRUE))