R o Stata I: Graficando frecuencias relativas acumulativas.

Graficar las frecuencias relativas acumulativas.

Haciendo uso de la variable intervalos generada en unos de los ejemplos anteriores podemos graficar un scatter o gráfico de puntos que represente las frecuencias relativas acumuladas. Intervalos iría en el eje horizontal y las frecuencias relativas acumuladas en el eje vertical. Haría falta generar una variable de frecuencias relativas acumuladas.

Datos:
3.45, 2.76, 4.06, 2.65, 6.43, 4.70, 3.24, 4.12, 3.25, 2.67, 3.65, 4.56, 3.45, 4.53, 4.56, 2.43, 3.65, 2.48, 5.46, 2.54, 1.36, 4.30, 2.54, 3.54, 1.23, 3.43, 5.43, 2.76, 3.14, 5.36

*Ordenar de menor a mayor la muestra

sort variable6

*Generar frecuencias absolutas

egen b=group(intervalos)

*Generar n o total de muestra.

egen n = count(intervalos)

*Dejar valores como únicos de las variables intervalos y frecuencias.

collapse (count) b, by(intervalos n)

*Generar las frecuencias relativas

gen frecuencia_rela = b/n

*Calcular las frecuencias relativas acumuladas

gen frecuencias_rel_acum=sum(frecuencias_relativas)

Ya se tienen las dos variables necesitadas. Los límites de clases especificados en la variable intervalos y las frecuencias relativas acumuladas especificada en su correspondiente variable. Para graficar aplicamos.

scatter frecuencis_rel_acum intervalos, connect(1) mlabel(frecuencias_rel_acum)

 

Notará que para el comando scatter primero se coloca el eje vertical y luego el eje horizontal. También se agregó otra opción al gráfico como son las etiquetas con la opción mlabel.

Graficando las frecuencias relativas acumuladas en R.

variable6 <- c(3.45, 2.76, 4.06, 2.65, 6.43, 4.70, 3.24, 4.12, 3.25, 2.67, 3.65, 4.56, 3.45, 4.53, 4.56, 2.43, 3.65, 2.48, 5.46, 2.54, 1.36, 4.30, 2.54, 3.54, 1.23, 3.43, 5.43, 2.76, 3.14, 5.36)

intervalos <- cut(variable6, breaks=6)

frecuencias <- as.data.frame(table(intervalos))

Calular las frecuencias relativas

frecuencias_absolutas <- frecuencias$Freq

frecuencias_relativas<-frecuencias_absolutas/length(variable6)

Calcular las frecuencias relativas acumulativas

frecuencias_relativas_acumuladas<-cumsum(frecuencias_relativas)

Obtener los límites de los intervalos

Primero, obtener la longitud de los intervalos, se quieren 6 intervalos:

diff(range(variable6))/6

Tenemos que la longitud es 0.8666667

Segundo, construir los límites:

limites<-seq(min(variable6)+0.8666667,max(variable6)+0.8666667, by=0.8666667)

Hemos creado las dos variables necesarias a graficar: limites y frecuencias_relativas_acumuladas. Una diferencia entre Stata y R es que R permite ubicar en la sintaxis la variable para el eje horizontal de primero.

plot(limites, frecuencias_relativas_acumuladas, col="red", type="b")

text(limites, frecuencias_relativas_acumuladas, frecuencias_relativas_acumuladas, pos=4,  col="green")

Nota: si tienes un método diferente o quieres compartir cómo se hace esto en otros programas no dude de compartirlos a morales.economia@gmail.com

Deybi Morales León

morales.economia@gmail.com