ESTADÍSTICA DESCRIPTIVA CON R ANÁLISIS DE DISTRIBUCIÓN DE FRECUENCIAS

ESTADÍSTICA DESCRIPTIVA CON R

ANÁLISIS DE DISTRIBUCIÓN DE FRECUENCIAS

VARIABLE CUANTITATIVA

Una vez ejecutado el programa aparece la consola de trabajo. El prompt de R es el signo > el cual indica que el programa está listo para recibir comandos.

Primero se importa desde EXCEL la base de datos que contiene la variable en estudio. Este archivo se llama ALEACIÓN y se ha guardado con extensión .csv.

En la consola de trabajo de R introducimos la secuencia de comandos:

> data<-read.delim("clipboard",header=T,dec=".",check.names=T)

Seleccionamos la tabla en EXCEL y copiamos (CTRL+C). En la consola de R damos enter

> data

Aparece la base de datos conteniendo 58 valores de resistencia a la compresión

> attach(data)

Los datos ya se pueden utilizar

Análisis de la variable resistencia

Se halla el valor mínimo y máximo de la distribución:

> min(resistencia)

[1] 66.4

> max(resistencia)

[1] 75.3

Primero se construye la tabla de frecuencias utilizando la función table. Se especifica el valor mínimo, máximo y el ancho de los intervalos. La secuencia de comandos es la siguiente:

> frecuencia<-table(cut(resistencia,breaks=seq(66,76,1)))

> frelativa<-frecuencia/58

> facumulada<-cumsum(frecuencia)

> facurelativa<-facumulada/58  

> data.frame(frecuencia,frelativa,facumulada,facurelativa)

Aparece la tabla de distribución de frecuencias. Para incluir la tabla en un informe, se guardan los resultados en un archivo con formato .doc mediante

Archivo ------- Guardar en Archivo --------- Resultados.doc ---------- Guardar

Posteriormente se abre este archivo utilizando un procesador de texto como Microsoft Word y mediante copiar y pegar se anexa la tabla en el informe.

Una forma alternativa consiste en seleccionar la tabla, copiar y pegar en un archivo de Word utilizando como fuente la letra Courier con tamaño 10

Para dibujar el Histograma se utiliza la función hist:

> > hist(resistencia,xlab="resistencia",main="Histograma",breaks=seq(66,76,1))

Aparece el histograma en la ventana gráfica. Puede guardarse con

Archivo --------- Guardar como --------- Jpeg

O copiarse directamente mediante: CTRL+C (Copiar como un Bitmap) y pegar en un documento de texto.

El argumento col permite dar color a las celdas. Los códigos de los colores son los siguientes: blanco (0,8), negro (1,9), rojo (2,10), verde (3,11), azul (4,12), magenta (5,13), violeta (6,14), amarillo (7,15).

> > hist(resistencia,xlab="resistencia",main="Histograma",breaks=seq(66,76,1),col="red")

La instrucción

> hist(resistencia)

Permite obtener el histograma que por defecto construye R con 10 intervalos

Puede indicarse un número determinado de intervalos, por ejemplo 8:

> > hist(resistencia,xlab="resistencia",main="Histograma",breaks=8,col=2)

R realiza un ajuste del número de intervalos apropiado para el Histograma, en este caso dibuja 10.

Si se desea obtener el histograma de frecuencias relativas:

> > hist(resistencia,xlab="resistencia", main="Histograma de frecuencias relativas", breaks=seq(66,76,1),prob=T)

Finalmente se dibuja la curva de densidad:

> lines(density(resistencia))

EJERCICIOS PROPUESTOS

1. Los datos de la Tabla corresponden a un muestreo de ruido ambiental del nivel de presión sonora (LP) medida en decibeles (dB) en diferentes estaciones de la ciudad de Cali durante el día (LPD) y la noche (LPN).

...