Đọc dữ liệu mNUTRIC stroke vào R

library(readxl)
mNUTRIC <- read_excel("D:/1. Dr Hung 108/Chuyên môn/5. My Journal/Đang thực hiện/Dinh duong/NC BS HUNG/257 tong nhap.xlsx")
View(mNUTRIC)

Thông tin dữ liệu vào mNUTRIC

Có bao nhiêu biến số (variable) và quan sát (observation)

Liệt kê 10 quan sát đầu tiên của dữ liệu

Liệt kê 10 quan sát cuối cùng của dữ liệu

dim(mNUTRIC)
## [1] 174  27
head(mNUTRIC,10)
## # A tibble: 10 × 27
##        id phone   age gender diagnosis  copd    af heartd hypertension diabeties
##     <dbl> <chr> <dbl>  <dbl>     <dbl> <dbl> <dbl>  <dbl>        <dbl>     <dbl>
##  1 2.23e7 0919…    68      1         1     2     2      2            1         2
##  2 2.23e7 0912…    78      2         1     2     2      2            1         2
##  3 2.23e7 0914…    74      2         2     2     2      2            1         2
##  4 2.23e7 0977…    72      1         2     2     2      2            1         2
##  5 2.23e7 0375…    61      1         2     2     2      2            2         2
##  6 2.23e7 0977…    53      2         1     2     2      2            1         1
##  7 2.23e7 0399…    57      1         2     2     2      2            1         1
##  8 2.23e7 0982…    64      2         2     2     2      2            1         2
##  9 2.23e7 0325…    74      1         2     2     2      2            1         2
## 10 2.23e7 0983…    71      2         2     2     2      2            1         2
## # ℹ 17 more variables: gcs <dbl>, sofa <dbl>, apache <dbl>, mnutric <dbl>,
## #   wbc <dbl>, neut <dbl>, lym <dbl>, plt <dbl>, mv <dbl>, timemv <dbl>,
## #   sonde <dbl>, timesonde <dbl>, hap <dbl>, icuday <dbl>, hosday <dbl>,
## #   deadhos <dbl>, dead30 <dbl>
tail(mNUTRIC,10)
## # A tibble: 10 × 27
##        id phone   age gender diagnosis  copd    af heartd hypertension diabeties
##     <dbl> <chr> <dbl>  <dbl>     <dbl> <dbl> <dbl>  <dbl>        <dbl>     <dbl>
##  1 2.39e7 0375…    41      1         3     2     2      2            2         2
##  2 2.39e7 0.90…    72      1         2     2     2      2            2         2
##  3 2.39e7 0912…    49      2         3     2     2      2            2         2
##  4 2.39e7 0394…    60      1         1     2     2      2            1         2
##  5 2.40e7 9.04…    75      1         1     2     2      2            1         1
##  6 2.40e7 0982…    79      1         2     2     2      2            2         2
##  7 2.40e7 0979…    83      2         3     2     2      2            2         2
##  8 2.40e7 .097…    47      1         2     2     2      2            2         2
##  9 2.40e7 0915…    81      2         2     2     2      2            2         2
## 10 2.40e7 9.62…    51      1         1     2     2      2            1         1
## # ℹ 17 more variables: gcs <dbl>, sofa <dbl>, apache <dbl>, mnutric <dbl>,
## #   wbc <dbl>, neut <dbl>, lym <dbl>, plt <dbl>, mv <dbl>, timemv <dbl>,
## #   sonde <dbl>, timesonde <dbl>, hap <dbl>, icuday <dbl>, hosday <dbl>,
## #   deadhos <dbl>, dead30 <dbl>

Tóm tắt dữ liệu mNUTRUC bằng hàm summary

summary(mNUTRIC)
##        id              phone                age            gender     
##  Min.   :22271433   Length:174         Min.   :21.00   Min.   :1.000  
##  1st Qu.:22558753   Class :character   1st Qu.:55.00   1st Qu.:1.000  
##  Median :23249395   Mode  :character   Median :64.00   Median :1.000  
##  Mean   :23134788                      Mean   :62.46   Mean   :1.385  
##  3rd Qu.:23658004                      3rd Qu.:73.00   3rd Qu.:2.000  
##  Max.   :24002433                      Max.   :94.00   Max.   :2.000  
##    diagnosis          copd             af            heartd     
##  Min.   :1.000   Min.   :1.000   Min.   :1.000   Min.   :1.000  
##  1st Qu.:1.000   1st Qu.:2.000   1st Qu.:2.000   1st Qu.:2.000  
##  Median :2.000   Median :2.000   Median :2.000   Median :2.000  
##  Mean   :1.828   Mean   :1.983   Mean   :1.966   Mean   :1.948  
##  3rd Qu.:2.000   3rd Qu.:2.000   3rd Qu.:2.000   3rd Qu.:2.000  
##  Max.   :3.000   Max.   :2.000   Max.   :2.000   Max.   :2.000  
##   hypertension     diabeties          gcs             sofa       
##  Min.   :1.000   Min.   :1.000   Min.   : 3.00   Min.   : 0.000  
##  1st Qu.:1.000   1st Qu.:2.000   1st Qu.: 8.00   1st Qu.: 1.000  
##  Median :1.000   Median :2.000   Median :11.50   Median : 4.000  
##  Mean   :1.316   Mean   :1.833   Mean   :11.13   Mean   : 3.592  
##  3rd Qu.:2.000   3rd Qu.:2.000   3rd Qu.:15.00   3rd Qu.: 6.000  
##  Max.   :2.000   Max.   :2.000   Max.   :15.00   Max.   :11.000  
##      apache         mnutric           wbc             neut       
##  Min.   : 1.00   Min.   :0.000   Min.   : 3.76   Min.   : 1.890  
##  1st Qu.: 7.00   1st Qu.:2.000   1st Qu.: 8.72   1st Qu.: 5.992  
##  Median :10.00   Median :2.000   Median :11.84   Median : 9.340  
##  Mean   :10.02   Mean   :2.241   Mean   :12.83   Mean   :10.232  
##  3rd Qu.:13.00   3rd Qu.:2.000   3rd Qu.:15.73   3rd Qu.:13.670  
##  Max.   :25.00   Max.   :6.000   Max.   :44.27   Max.   :40.740  
##       lym             plt              mv            timemv      
##  Min.   :0.380   Min.   :  1.1   Min.   :1.000   Min.   : 0.000  
##  1st Qu.:1.030   1st Qu.:206.2   1st Qu.:1.000   1st Qu.: 0.000  
##  Median :1.420   Median :250.5   Median :1.000   Median : 2.000  
##  Mean   :1.587   Mean   :258.4   Mean   :1.408   Mean   : 4.621  
##  3rd Qu.:1.972   3rd Qu.:304.8   3rd Qu.:2.000   3rd Qu.: 8.000  
##  Max.   :7.240   Max.   :604.0   Max.   :2.000   Max.   :29.000  
##      sonde        timesonde          hap            icuday      
##  Min.   :0.00   Min.   : 0.00   Min.   :1.000   Min.   : 0.000  
##  1st Qu.:1.00   1st Qu.: 0.00   1st Qu.:2.000   1st Qu.: 0.250  
##  Median :1.00   Median :10.00   Median :2.000   Median : 7.000  
##  Mean   :1.31   Mean   :10.74   Mean   :1.799   Mean   : 9.356  
##  3rd Qu.:2.00   3rd Qu.:18.00   3rd Qu.:2.000   3rd Qu.:16.000  
##  Max.   :2.00   Max.   :50.00   Max.   :2.000   Max.   :37.000  
##      hosday         deadhos          dead30     
##  Min.   : 2.00   Min.   :1.000   Min.   :1.000  
##  1st Qu.:10.00   1st Qu.:2.000   1st Qu.:2.000  
##  Median :15.00   Median :2.000   Median :2.000  
##  Mean   :15.35   Mean   :1.908   Mean   :1.845  
##  3rd Qu.:19.00   3rd Qu.:2.000   3rd Qu.:2.000  
##  Max.   :50.00   Max.   :2.000   Max.   :2.000

Biên tập dữ liệu bằng gói phân tích “tidyverse”

Mã hoá biến gender (F/M) thành biến sex với giá trị 0/1 (0= M; 1= F)

Cách đơn giản:

mNUTRIC$sex[mNUTRIC$gender == "1"] = 0
## Warning: Unknown or uninitialised column: `sex`.
mNUTRIC$sex[mNUTRIC$gender == "2"] = 1
table(mNUTRIC$sex)
## 
##   0   1 
## 107  67

Mô tả đặc điểm mẫu nghiên cứu:

library(table1)
## 
## Attaching package: 'table1'
## The following objects are masked from 'package:base':
## 
##     units, units<-
table1(~ age + sex + diagnosis + copd + af + heartd + hypertension + diabeties + gcs + sofa + apache + mnutric + wbc + neut + lym + plt + mv + timemv + sonde + timesonde + hap + icuday + hosday + dead30, data = mNUTRIC)
Overall
(N=174)
age
Mean (SD) 62.5 (14.6)
Median [Min, Max] 64.0 [21.0, 94.0]
sex
Mean (SD) 0.385 (0.488)
Median [Min, Max] 0 [0, 1.00]
diagnosis
Mean (SD) 1.83 (0.748)
Median [Min, Max] 2.00 [1.00, 3.00]
copd
Mean (SD) 1.98 (0.131)
Median [Min, Max] 2.00 [1.00, 2.00]
af
Mean (SD) 1.97 (0.183)
Median [Min, Max] 2.00 [1.00, 2.00]
heartd
Mean (SD) 1.95 (0.222)
Median [Min, Max] 2.00 [1.00, 2.00]
hypertension
Mean (SD) 1.32 (0.466)
Median [Min, Max] 1.00 [1.00, 2.00]
diabeties
Mean (SD) 1.83 (0.374)
Median [Min, Max] 2.00 [1.00, 2.00]
gcs
Mean (SD) 11.1 (3.45)
Median [Min, Max] 11.5 [3.00, 15.0]
sofa
Mean (SD) 3.59 (2.54)
Median [Min, Max] 4.00 [0, 11.0]
apache
Mean (SD) 10.0 (4.64)
Median [Min, Max] 10.0 [1.00, 25.0]
mnutric
Mean (SD) 2.24 (1.02)
Median [Min, Max] 2.00 [0, 6.00]
wbc
Mean (SD) 12.8 (5.39)
Median [Min, Max] 11.8 [3.76, 44.3]
neut
Mean (SD) 10.2 (5.45)
Median [Min, Max] 9.34 [1.89, 40.7]
lym
Mean (SD) 1.59 (0.885)
Median [Min, Max] 1.42 [0.380, 7.24]
plt
Mean (SD) 258 (82.1)
Median [Min, Max] 251 [1.10, 604]
mv
Mean (SD) 1.41 (0.493)
Median [Min, Max] 1.00 [1.00, 2.00]
timemv
Mean (SD) 4.62 (5.74)
Median [Min, Max] 2.00 [0, 29.0]
sonde
Mean (SD) 1.31 (0.476)
Median [Min, Max] 1.00 [0, 2.00]
timesonde
Mean (SD) 10.7 (10.5)
Median [Min, Max] 10.0 [0, 50.0]
hap
Mean (SD) 1.80 (0.402)
Median [Min, Max] 2.00 [1.00, 2.00]
icuday
Mean (SD) 9.36 (9.06)
Median [Min, Max] 7.00 [0, 37.0]
hosday
Mean (SD) 15.4 (7.77)
Median [Min, Max] 15.0 [2.00, 50.0]
dead30
Mean (SD) 1.84 (0.363)
Median [Min, Max] 2.00 [1.00, 2.00]

Mô tả đặc điểm mẫu nghiên cứu theo hai nhom: nhom song va nhom tu vong 30 ngay

as.factor là báo cho phần mềm hiều biến hypertension và diabetes biến định tính

table1(~ age + as.factor(sex) + as.factor(diagnosis) + as.factor(copd) + as.factor(af) + as.factor(heartd) + as.factor(hypertension) + as.factor(diabeties) + gcs + sofa + apache + mnutric + wbc + neut + lym + plt + as.factor(mv) + timemv + as.factor(sonde) + timesonde + as.factor(hap) + icuday + hosday | dead30, data = mNUTRIC)
## Warning in table1.formula(~age + as.factor(sex) + as.factor(diagnosis) + :
## Terms to the right of '|' in formula 'x' define table columns and are expected
## to be factors with meaningful labels.
1
(N=27)		 2
(N=147)		 Overall
(N=174)
age
Mean (SD) 65.7 (10.4) 61.9 (15.2) 62.5 (14.6)
Median [Min, Max] 66.0 [48.0, 94.0] 63.0 [21.0, 92.0] 64.0 [21.0, 94.0]
as.factor(sex)
0 16 (59.3%) 91 (61.9%) 107 (61.5%)
1 11 (40.7%) 56 (38.1%) 67 (38.5%)
as.factor(diagnosis)
1 10 (37.0%) 56 (38.1%) 66 (37.9%)
2 9 (33.3%) 63 (42.9%) 72 (41.4%)
3 8 (29.6%) 28 (19.0%) 36 (20.7%)
as.factor(copd)
1 1 (3.7%) 2 (1.4%) 3 (1.7%)
2 26 (96.3%) 145 (98.6%) 171 (98.3%)
as.factor(af)
1 2 (7.4%) 4 (2.7%) 6 (3.4%)
2 25 (92.6%) 143 (97.3%) 168 (96.6%)
as.factor(heartd)
1 1 (3.7%) 8 (5.4%) 9 (5.2%)
2 26 (96.3%) 139 (94.6%) 165 (94.8%)
as.factor(hypertension)
1 21 (77.8%) 98 (66.7%) 119 (68.4%)
2 6 (22.2%) 49 (33.3%) 55 (31.6%)
as.factor(diabeties)
1 3 (11.1%) 26 (17.7%) 29 (16.7%)
2 24 (88.9%) 121 (82.3%) 145 (83.3%)
gcs
Mean (SD) 7.63 (2.63) 11.8 (3.19) 11.1 (3.45)
Median [Min, Max] 7.00 [3.00, 13.0] 13.0 [5.00, 15.0] 11.5 [3.00, 15.0]
sofa
Mean (SD) 5.93 (2.02) 3.16 (2.39) 3.59 (2.54)
Median [Min, Max] 6.00 [2.00, 11.0] 3.00 [0, 10.0] 4.00 [0, 11.0]
apache
Mean (SD) 13.3 (3.94) 9.43 (4.52) 10.0 (4.64)
Median [Min, Max] 12.0 [8.00, 25.0] 8.00 [1.00, 25.0] 10.0 [1.00, 25.0]
mnutric
Mean (SD) 2.96 (1.43) 2.11 (0.869) 2.24 (1.02)
Median [Min, Max] 2.00 [1.00, 6.00] 2.00 [0, 6.00] 2.00 [0, 6.00]
wbc
Mean (SD) 14.9 (4.41) 12.4 (5.48) 12.8 (5.39)
Median [Min, Max] 15.5 [8.75, 23.8] 11.4 [3.76, 44.3] 11.8 [3.76, 44.3]
neut
Mean (SD) 12.6 (4.19) 9.80 (5.56) 10.2 (5.45)
Median [Min, Max] 12.1 [5.78, 21.1] 8.58 [1.89, 40.7] 9.34 [1.89, 40.7]
lym
Mean (SD) 1.27 (0.640) 1.65 (0.913) 1.59 (0.885)
Median [Min, Max] 1.16 [0.440, 3.16] 1.51 [0.380, 7.24] 1.42 [0.380, 7.24]
plt
Mean (SD) 260 (79.6) 258 (82.8) 258 (82.1)
Median [Min, Max] 255 [68.0, 418] 248 [1.10, 604] 251 [1.10, 604]
as.factor(mv)
1 26 (96.3%) 77 (52.4%) 103 (59.2%)
2 1 (3.7%) 70 (47.6%) 71 (40.8%)
timemv
Mean (SD) 8.15 (6.29) 3.97 (5.41) 4.62 (5.74)
Median [Min, Max] 7.00 [0, 24.0] 1.00 [0, 29.0] 2.00 [0, 29.0]
as.factor(sonde)
0 0 (0%) 1 (0.7%) 1 (0.6%)
1 27 (100%) 91 (61.9%) 118 (67.8%)
2 0 (0%) 55 (37.4%) 55 (31.6%)
timesonde
Mean (SD) 11.6 (9.89) 10.6 (10.7) 10.7 (10.5)
Median [Min, Max] 8.00 [2.00, 40.0] 10.0 [0, 50.0] 10.0 [0, 50.0]
as.factor(hap)
1 9 (33.3%) 26 (17.7%) 35 (20.1%)
2 18 (66.7%) 121 (82.3%) 139 (79.9%)
icuday
Mean (SD) 10.4 (8.78) 9.17 (9.12) 9.36 (9.06)
Median [Min, Max] 7.00 [2.00, 31.0] 7.00 [0, 37.0] 7.00 [0, 37.0]
hosday
Mean (SD) 11.6 (9.89) 16.0 (7.14) 15.4 (7.77)
Median [Min, Max] 8.00 [2.00, 40.0] 15.0 [3.00, 50.0] 15.0 [2.00, 50.0]