0.1 INTRODUCCIÓN

Este trabajo se centra en analizar una base de datos exhaustiva para explorar el panorama del tabaquismo en Inglaterra.

library(readr)
## Warning: package 'readr' was built under R version 4.2.3
smoking <- read_csv("C:/Users/usuario/Downloads/smoking.csv")
## Rows: 1691 Columns: 12
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): gender, marital_status, highest_qualification, nationality, ethnici...
## dbl (3): age, amt_weekends, amt_weekdays
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
#porcentaje de fumadores

summary(smoking$smoke=="Yes")
##    Mode   FALSE    TRUE 
## logical    1270     421
length(smoking$smoke)
## [1] 1691
(421/1691)*100
## [1] 24.89651

Haciendo un summary en en la base de datos obtenemos que las personas que fuman son 421 y las no lo hacen son 1270 con lo que se con concluye que el 24.89% de las personas recogidas en este banco de datos fuman. ## Tabaquismo según genero

#hombres
summary(smoking$gender[smoking$smoke=="Yes"]=="Male")
##    Mode   FALSE    TRUE 
## logical     234     187
#¿De los hombres cuantos fuman mas de 15 cigarros?
table(smoking$gender[smoking$amt_weekdays>15]=="Male")
## 
## FALSE  TRUE 
##    63    76
#mujeres
summary(smoking$gender[smoking$smoke=="Yes"]=="Female")
##    Mode   FALSE    TRUE 
## logical     187     234
#¿De las mujeres cuantas fuman mas de 15 cigarros?
table(smoking$gender[smoking$amt_weekdays>15]=="Female")
## 
## FALSE  TRUE 
##    76    63


0.2 Rangos de edad

summary(smoking$smoke[smoking$age<=18]=="Yes")
##    Mode   FALSE    TRUE 
## logical      34      16
#Número de personas menores de 18 años que fuman
length(smoking$smoke[smoking$age<=18]=="Yes")
## [1] 50
#Personas menores de 40 años que fuman
table(smoking$smoke[smoking$age>18]=="Yes")-table(smoking$smoke[smoking$age>=40]=="Yes")
## 
## FALSE  TRUE 
##   357   181
summary(smoking$smoke[smoking$age>=40]=="Yes")
##    Mode   FALSE    TRUE 
## logical     879     224
#Número de personas mayores de 40 años que fuman 
length(smoking$smoke[smoking$age>=40]=="Yes")
## [1] 1103


0.3 Tabaquismo según estado civil

summary(smoking$smoke[smoking$marital_status=="Single"]=="Yes")
##    Mode   FALSE    TRUE 
## logical     269     158
length(smoking$smoke[smoking$marital_status=="Single"]=="Yes")
## [1] 427
summary(smoking$smoke[smoking$marital_status=="Married"]=="Yes")
##    Mode   FALSE    TRUE 
## logical     669     143
length(smoking$smoke[smoking$marital_status=="Married"]=="Yes")
## [1] 812
table(smoking$marital_status[smoking$smoke=="Yes"]=="Divorced")+table(smoking$marital_status[smoking$smoke=="Yes"]=="Separated")
## 
## FALSE  TRUE 
##   762    80
summary(smoking$smoke[smoking$marital_status=="Widowed"]=="Yes")
##    Mode   FALSE    TRUE 
## logical     183      40
length(smoking$smoke[smoking$marital_status=="Widowed"]=="Yes")
## [1] 223


0.4 Cantidad de cigarros fumados durante el fin de semana

#menos de 10 cigarros
summary(smoking$smoke=="Yes")
##    Mode   FALSE    TRUE 
## logical    1270     421
table(smoking$smoke[smoking$amt_weekends<=10]=="Yes")
## 
## TRUE 
##  156
#entre 10 y 20 cigarros
summary(smoking$smoke=="Yes")
##    Mode   FALSE    TRUE 
## logical    1270     421
table(smoking$smoke[smoking$amt_weekends>10]=="Yes")-table(smoking$smoke[smoking$amt_weekends>=20]=="Yes")
## 
## TRUE 
##   77
#mas de 20 cigarros
summary(smoking$smoke=="Yes")
##    Mode   FALSE    TRUE 
## logical    1270     421
table(smoking$smoke[smoking$amt_weekends>=20]=="Yes")
## 
## TRUE 
##  188


0.5 Media de tipo de tabaco

Preferencias de tabaco:

mean(smoking$type[smoking$smoke=="Yes"]=="Packets")
## [1] 0.7054632
mean(smoking$type[smoking$smoke=="Yes"]=="Hand-Rolled")
## [1] 0.1710214


0.6 EJERCICIO PROPUESTO

Utilizando el archivo Smoking.csv.La media de cigarros fumados durante el fin de semana puede ser 18 con una confianza del 90%

t.test(smoking$amt_weekends[smoking$smoke=="Yes"], mu =18, conf.level = 0.9)
## 
##  One Sample t-test
## 
## data:  smoking$amt_weekends[smoking$smoke == "Yes"]
## t = -3.2958, df = 420, p-value = 0.001065
## alternative hypothesis: true mean is not equal to 18
## 90 percent confidence interval:
##  15.61610 17.20575
## sample estimates:
## mean of x 
##  16.41093

Como se puede comprobar por el p-valor obtenido, se acepta la hipótesis alternativa.