Gráficos de Historia y Densidad en R

Puedes crear histogramas con la función hist(x</em>) donde x es un vector numérico de valores a representar. La opción freq=FALSE traza densidades de probabilidad en lugar de frecuencias. La opción breaks= controla el número de bins.

# Simple Histogram
hist(mtcars$mpg)

# Colored Histogram with Different Number of Bins
hist(mtcars$mpg, breaks=12, col="red")

# Add a Normal Curve (Thanks to Peter Dalgaard)
x <- mtcars$mpg
h<-hist(x, breaks=10, col="red", xlab="Miles Per Gallon",
   main="Histogram with Normal Curve")
xfit<-seq(min(x),max(x),length=40)
yfit<-dnorm(xfit,mean=mean(x),sd=sd(x))
yfit <- yfit*diff(h$mids[1:2])*length(x)

lines(xfit, yfit, col="blue", lwd=2)

Los histogramas pueden ser un mal método para determinar la forma de una distribución, porque se ven muy afectados por el número de intervalos utilizados.

Gráficos de densidad del núcleo

Los gráficos de densidad de núcleos suelen ser una forma mucho más eficaz de ver la distribución de una variable. Crea el gráfico utilizando gráfico(densidad(x)) donde x es un vector numérico.

# Kernel Density Plot
d <- density(mtcars$mpg) # returns the density data
plot(d) # plots the results

# Filled Density Plot
d <- density(mtcars$mpg)
plot(d, main="Kernel Density of Miles Per Gallon")
polygon(d, col="red", border="blue")

Comparación de Grupos Densidad del Núcleo VIA

La función sm.density.compare( ) del paquete sm</a ></strong > te permite superponer los gráficos de densidad kernal de dos o más grupos. El formato es sm.densidad.comparar(x , factor) donde x es un vector numérico y factor es la variable de agrupación.

# Compare MPG distributions for cars with
#
4,6, or 8 cylinders
library(sm)
attach(mtcars)

# create value labels
cyl.f <- factor(cyl, levels= c(4,6,8),
  labels = c("4 cylinder", "6 cylinder", "8 cylinder"))

# plot densities
sm.density.compare(mpg, cyl, xlab="Miles Per Gallon")
title(main="MPG Distribution by Car Cylinders")

# add legend via mouse click
colfill<-c(2:(2+length(levels(cyl.f))))
legend(locator(1), levels(cyl.f), fill=colfill)