Ejemplo 15

Se realiza un estudio acerca de los efectos del tabaquismo sobre los patrones de sueño. La medición que se observa es el tiempo, en minutos, que toma quedar dormido. Se obtienen los siguientes datos:

Fumadores:

       69.3 56.0 22.1 47.6

       53.2 48.1 52.7 34.4
       
       60.2 43.8 23.2 13.8

No fumadores:

          28.6 25.1 26.4 34.9

          29.8 28.4 38.5 30.2
          
          30.6 31.8 41.6 21.1
          
          36.0 37.9 13.9
######### CARGAMOS LOS DATOS ############

#### FUMADORES
x1 <- c(69.3, 56.0, 22.1, 47.6,
        53.2, 48.1, 52.7, 34.4,
        60.2, 43.8, 23.2, 13.8)

#### NO FUMADORES
x2 <- c(28.6, 25.1, 26.4, 34.9,
        29.8, 28.4, 38.5, 30.2,
        30.6, 31.8, 41.6, 21.1,
        36.0, 37.9, 13.9)

Convertimos a data.frame

g1 <- as.data.frame(x1)
g2 <- as.data.frame(x2)

Generamos los histogramas

par(mfrow = c(1, 2))
 hist(g1$x1,plot=TRUE, freq = FALSE,main = "Fumadores")
dx1 <- density(g1$x1)
lines(dx1$x, dx1$y, col = 2, lwd = 4)

 hist(g2$x2,plot=TRUE, freq = FALSE, main = "No Fumadores")
dx2 <- density(g2$x2)
lines(dx2$x, dx2$y, col = 2, lwd = 4)

Tablas de frecuencias con datos agrupados

Tabla1
##   Lower Upper Main Frequency Percentage CF   CPF
## 1    10    20   15         1        8.3  1   8.3
## 2    20    30   25         2       16.7  3  25.0
## 3    30    40   35         1        8.3  4  33.3
## 4    40    50   45         3       25.0  7  58.3
## 5    50    60   55         3       25.0 10  83.3
## 6    60    70   65         2       16.7 12 100.0
Tabla2
##   Lower Upper Main Frequency Percentage CF   CPF
## 1    10    15 12.5         1        6.7  1   6.7
## 2    15    20 17.5         0        0.0  1   6.7
## 3    20    25 22.5         1        6.7  2  13.3
## 4    25    30 27.5         5       33.3  7  46.7
## 5    30    35 32.5         4       26.7 11  73.3
## 6    35    40 37.5         3       20.0 14  93.3
## 7    40    45 42.5         1        6.7 15 100.0

Medidas de tendencia central

##### FUMADORES
mean(g1$x1)
## [1] 43.7
sd(g1$x1)
## [1] 16.92776
median(g1$x1)
## [1] 47.85
##### NO FUMADORES
mean(g2$x2)
## [1] 30.32
sd(g2$x2)
## [1] 7.127833
median(g2$x2)
## [1] 30.2

Coeficiente de variación

####### FUMADORES
cv1 <- sd(g1$x1)/mean(g1$x1)
cv1
## [1] 0.387363
###### NO FUMADORES
cv2 <- sd(g2$x2)/mean(g2$x2)
cv2
## [1] 0.2350868

Índice de asimetría de Pearson

### Fumadores
As1 <- (3*(mean(g1$x1)-median(g1$x1)))/(sd(g1$x1))
As1
## [1] -0.7354782
### No Fumadores
As2 <- (3*(mean(g2$x2)-median(g2$x2)))/(sd(g2$x2))
As2
## [1] 0.05050624