This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.

Try executing this chunk by clicking the Run button within the chunk or by placing your cursor inside it and pressing Ctrl+Shift+Enter.

#data ##data list

library(gtsummary)
airquality
NA

#analysis ##univariate

```{r)




<!-- rnb-text-end -->


<!-- rnb-chunk-begin -->


<!-- rnb-output-begin eyJkYXRhIjoiXG48IS0tIHJuYi1zb3VyY2UtYmVnaW4gZXlKa1lYUmhJam9pWUdCZ2NseHViR2xpY21GeWVTaGtjR3g1Y2lsY2JtRnBjbkYxWVd4cGRIa2dKVDRsSUZ4dUlDQjBZbXhmYzNWdGJXRnllU2dwWEc1Z1lHQWlmUT09IC0tPlxuXG5gYGByXG5saWJyYXJ5KGRwbHlyKVxuYWlycXVhbGl0eSAlPiUgXG4gIHRibF9zdW1tYXJ5KClcbmBgYFxuXG48IS0tIHJuYi1zb3VyY2UtZW5kIC0tPlxuIn0= -->


<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxubGlicmFyeShkcGx5cilcbmFpcnF1YWxpdHkgJT4lIFxuICB0Ymxfc3VtbWFyeSgpXG5gYGAifQ== -->

```r
library(dplyr)
airquality %>% 
  tbl_summary()
Characteristic N = 1531
Ozone 32 (18, 64)
    Unknown 37
Solar.R 205 (115, 259)
    Unknown 7
Wind 9.7 (7.4, 11.5)
Temp 79 (72, 85)
Month
    5 31 (20%)
    6 30 (20%)
    7 31 (20%)
    8 31 (20%)
    9 30 (20%)
Day 16 (8, 23)
1 Median (Q1, Q3); n (%)
NA

#selected variable

head(airquality)
NA
NA
NA
NA
airquality %>% 
tbl_summary( 
by = "Month") 
Characteristic 5
N = 31
1
6
N = 30
1
7
N = 31
1
8
N = 31
1
9
N = 30
1
Ozone 18 (11, 32) 23 (20, 37) 60 (35, 80) 52 (28, 84) 23 (16, 36)
    Unknown 5 21 5 5 1
Solar.R 194 (66, 290) 189 (127, 273) 253 (175, 274) 198 (99, 233) 192 (112, 236)
    Unknown 4 0 0 3 0
Wind 11.5 (8.6, 14.3) 9.7 (8.0, 11.5) 8.6 (6.9, 10.9) 8.6 (6.3, 11.5) 10.3 (7.4, 12.6)
Temp 66 (59, 69) 78 (76, 83) 84 (81, 86) 82 (79, 89) 76 (71, 81)
Day 16 (8, 24) 16 (8, 23) 16 (8, 24) 16 (8, 24) 16 (8, 23)
1 Median (Q1, Q3)

#p value

airquality %>% 
tbl_summary( by = "Month")%>%
  add_p()
Error in airquality %>% tbl_summary(by = "Month") %>% add_p() : 
  could not find function "%>%"
airquality %>%
  tbl_summary(
    by = "Month")
Error in airquality %>% tbl_summary(by = "Month") : 
  could not find function "%>%"
LS0tDQp0aXRsZTogIlIgTm90ZWJvb2siDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KDQpUaGlzIGlzIGFuIFtSIE1hcmtkb3duXShodHRwOi8vcm1hcmtkb3duLnJzdHVkaW8uY29tKSBOb3RlYm9vay4gV2hlbiB5b3UgZXhlY3V0ZSBjb2RlIHdpdGhpbiB0aGUgbm90ZWJvb2ssIHRoZSByZXN1bHRzIGFwcGVhciBiZW5lYXRoIHRoZSBjb2RlLiANCg0KVHJ5IGV4ZWN1dGluZyB0aGlzIGNodW5rIGJ5IGNsaWNraW5nIHRoZSAqUnVuKiBidXR0b24gd2l0aGluIHRoZSBjaHVuayBvciBieSBwbGFjaW5nIHlvdXIgY3Vyc29yIGluc2lkZSBpdCBhbmQgcHJlc3NpbmcgKkN0cmwrU2hpZnQrRW50ZXIqLiANCg0KYGBge3J9DQoNCmBgYA0KI2RhdGENCiMjZGF0YSBsaXN0DQpgYGB7cn0NCmxpYnJhcnkoZ3RzdW1tYXJ5KQ0KYWlycXVhbGl0eQ0KDQpgYGANCiNhbmFseXNpcw0KIyN1bml2YXJpYXRlDQoNCmBgYHtyKQ0KDQogIA0KYGBgDQoNCg0KYGBge3J9DQpsaWJyYXJ5KGRwbHlyKQ0KYWlycXVhbGl0eSAlPiUgDQogIHRibF9zdW1tYXJ5KCkNCiAgDQpgYGANCiNzZWxlY3RlZCB2YXJpYWJsZQ0KYGBge3J9DQpoZWFkKGFpcnF1YWxpdHkpDQoNCg0KDQoNCmBgYA0KYGBge3J9DQphaXJxdWFsaXR5ICU+JSANCnRibF9zdW1tYXJ5KCANCmJ5ID0gIk1vbnRoIikgDQpgYGANCiNwIHZhbHVlDQoNCmBgYHtyfQ0KYWlycXVhbGl0eSAlPiUgDQp0Ymxfc3VtbWFyeSggYnkgPSAiTW9udGgiKSU+JQ0KICBhZGRfcCgpDQoNCg0KYGBgDQpgYGB7cn0NCmFpcnF1YWxpdHkgJT4lDQogIHRibF9zdW1tYXJ5KA0KICAgIGJ5ID0gIk1vbnRoIikNCiAgICgpIHN0YXRpc3RpYyA9IGxpc3QoDQogICAgICBhbGxfY29udGludW91cygpIH4gInttZWFufSAoe3NkfSkiLA0KICAgICAgYWxsX2NhdGVnb3JpY2FsKCkgfiAie259IC8ge059ICh7cH0lKSINCiAgICApDQogICkNCg0KYGBgDQoNCg0KYGBge3J9DQpgYGANCg0KDQpgYGB7cn0NCmBgYA0KDQo=