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

Comentarios

Envíanos un mensaje

Nombre

Correo electrónico *

Mensaje *