library(datasets) 
library(dplyr) 
library(gtsummary)
data(airquality)
head(airquality)
summary(airquality)
     Ozone           Solar.R           Wind       
 Min.   :  1.00   Min.   :  7.0   Min.   : 1.700  
 1st Qu.: 18.00   1st Qu.:115.8   1st Qu.: 7.400  
 Median : 31.50   Median :205.0   Median : 9.700  
 Mean   : 42.13   Mean   :185.9   Mean   : 9.958  
 3rd Qu.: 63.25   3rd Qu.:258.8   3rd Qu.:11.500  
 Max.   :168.00   Max.   :334.0   Max.   :20.700  
 NA's   :37       NA's   :7                       
      Temp           Month            Day      
 Min.   :56.00   Min.   :5.000   Min.   : 1.0  
 1st Qu.:72.00   1st Qu.:6.000   1st Qu.: 8.0  
 Median :79.00   Median :7.000   Median :16.0  
 Mean   :77.88   Mean   :6.993   Mean   :15.8  
 3rd Qu.:85.00   3rd Qu.:8.000   3rd Qu.:23.0  
 Max.   :97.00   Max.   :9.000   Max.   :31.0  
                                               
summary_table <- airquality %>% tbl_summary(by =Wind)
Error in tbl_summary(., by = Wind) : 
  could not find function "tbl_summary"
print(summary_table)
Error: object 'summary_table' not found

#Univariate Analysis

table(airquality$Month)

 5  6  7  8  9 
31 30 31 31 30 
table(airquality$Day)

 1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 
 5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5 
22 23 24 25 26 27 28 29 30 31 
 5  5  5  5  5  5  5  5  5  3 
hist(airquality$Ozone, main = "Ozone Distribution", xlab = "Ozone", col = "lightblue", border = "black")

barplot(table(airquality$Ozone), main = "Ozone Distribution", xlab = "Ozone", ylab = "Frequency", col = "lightgreen")

airquality %>%
  tbl_summary(
    by = Ozone)
Error in tbl_summary(., by = Ozone) : 
  could not find function "tbl_summary"
airquality %>%
  tbl_summary(
    by = day,
    statistic = list(
      all_continuous() ~ "{mean} ({sd})",
      all_categorical() ~ "{n} / {N} ({p}%)"
    )
  )
Error in tbl_summary(., by = day, statistic = list(all_continuous() ~  : 
  could not find function "tbl_summary"
LS0tDQp0aXRsZTogIlIgTm90ZWJvb2siDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KDQpgYGB7cn0NCmxpYnJhcnkoZGF0YXNldHMpIA0KbGlicmFyeShkcGx5cikgDQpgYGANCg0KDQpgYGB7cn0NCmxpYnJhcnkoZ3RzdW1tYXJ5KQ0KYGBgDQoNCg0KYGBge3J9DQpkYXRhKGFpcnF1YWxpdHkpDQpgYGANCg0KDQpgYGB7cn0NCmhlYWQoYWlycXVhbGl0eSkNCmBgYA0KDQoNCmBgYHtyfQ0Kc3VtbWFyeShhaXJxdWFsaXR5KQ0KYGBgDQpgYGB7cn0NCnN1bW1hcnlfdGFibGUgPC0gYWlycXVhbGl0eSAlPiUgdGJsX3N1bW1hcnkoYnkgPVdpbmQpDQpgYGANCg0KDQpgYGB7cn0NCnByaW50KHN1bW1hcnlfdGFibGUpDQpgYGANCiNVbml2YXJpYXRlIEFuYWx5c2lzDQoNCmBgYHtyfQ0KdGFibGUoYWlycXVhbGl0eSRNb250aCkNCmBgYA0KDQoNCmBgYHtyfQ0KdGFibGUoYWlycXVhbGl0eSREYXkpDQpgYGANCg0KYGBge3J9DQpoaXN0KGFpcnF1YWxpdHkkT3pvbmUsIG1haW4gPSAiT3pvbmUgRGlzdHJpYnV0aW9uIiwgeGxhYiA9ICJPem9uZSIsIGNvbCA9ICJsaWdodGJsdWUiLCBib3JkZXIgPSAiYmxhY2siKQ0KYGBgDQoNCg0KYGBge3J9DQpiYXJwbG90KHRhYmxlKGFpcnF1YWxpdHkkT3pvbmUpLCBtYWluID0gIk96b25lIERpc3RyaWJ1dGlvbiIsIHhsYWIgPSAiT3pvbmUiLCB5bGFiID0gIkZyZXF1ZW5jeSIsIGNvbCA9ICJsaWdodGdyZWVuIikNCmBgYA0KDQoNCmBgYHtyfQ0KDQpgYGANCg0KDQpgYGB7cn0NCmFpcnF1YWxpdHkgJT4lDQogIHRibF9zdW1tYXJ5KA0KICAgIGJ5ID0gT3pvbmUpDQpgYGANCg0KDQpgYGB7cn0NCmFpcnF1YWxpdHkgJT4lDQogIHRibF9zdW1tYXJ5KA0KICAgIGJ5ID0gZGF5LA0KICAgIHN0YXRpc3RpYyA9IGxpc3QoDQogICAgICBhbGxfY29udGludW91cygpIH4gInttZWFufSAoe3NkfSkiLA0KICAgICAgYWxsX2NhdGVnb3JpY2FsKCkgfiAie259IC8ge059ICh7cH0lKSINCiAgICApDQogICkNCmBgYA0KDQoNCmBgYHtyfQ0KYGBgDQoNCg==