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
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
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
#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
Preferencias de tabaco:
mean(smoking$type[smoking$smoke=="Yes"]=="Packets")
## [1] 0.7054632
mean(smoking$type[smoking$smoke=="Yes"]=="Hand-Rolled")
## [1] 0.1710214
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.