Đưa dữ liệu vào Rstudio

library(readxl)
lar1 <- read_excel("~/Desktop/R - R studio/Chạy mẫu lần 1/lar1.xlsx")
View(lar1)
attach(lar1)
head(lar1)
## # A tibble: 6 × 57
##     stt   sex   age  ngay   tha   dtd   btm   dqn   bpm  bthm   bgm    ut   nkp
##   <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1     2     1  69.6     6     0     0     0     0     0     0     0     1     1
## 2     3     1  68.6    20     1     1     0     0     0     0     0     0     0
## 3     4     1  71.2    12     1     1     1     0     0     0     0     0     0
## 4     5     0  65.6    13     1     1     0     0     0     0     0     0     0
## 5     7     1  72.7     8     0     0     0     0     0     0     0     0     1
## 6     8     0  86.6     1     0     0     0     1     0     0     0     0     1
## # ℹ 44 more variables: nktn <dbl>, nkth <dbl>, nkdmm <dbl>, khac <dbl>,
## #   gcs <dbl>, hatb <dbl>, hatt <dbl>, hattr <dbl>, mach <dbl>, nhietdo <dbl>,
## #   nhiptho <dbl>, sp02 <dbl>, fi02 <dbl>, PFr <dbl>, caymau <dbl>, rltk <dbl>,
## #   rlhh <dbl>, rltm <dbl>, rlthan <dbl>, rlgan <dbl>, rldongmau <dbl>,
## #   suydacoquan <dbl>, thongkhich <dbl>, locmau <dbl>, bc <dbl>, hct <dbl>,
## #   tc <dbl>, bilitp <dbl>, lactate <dbl>, albumin <dbl>, lar <dbl>, pct <dbl>,
## #   na <dbl>, k <dbl>, ast <dbl>, alt <dbl>, ure <dbl>, crea <dbl>, …
head(lar1,10)
## # A tibble: 10 × 57
##      stt   sex   age  ngay   tha   dtd   btm   dqn   bpm  bthm   bgm    ut   nkp
##    <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
##  1     2     1  69.6     6     0     0     0     0     0     0     0     1     1
##  2     3     1  68.6    20     1     1     0     0     0     0     0     0     0
##  3     4     1  71.2    12     1     1     1     0     0     0     0     0     0
##  4     5     0  65.6    13     1     1     0     0     0     0     0     0     0
##  5     7     1  72.7     8     0     0     0     0     0     0     0     0     1
##  6     8     0  86.6     1     0     0     0     1     0     0     0     0     1
##  7     9     0  63.2     1     0     0     0     1     0     0     0     1     1
##  8    10     1  66.7    11     1     0     1     0     0     0     0     0     1
##  9    11     0  84.6    17     1     0     1     1     0     1     0     0     1
## 10    13     0  45.4     7     0     0     0     0     0     0     0     0     0
## # ℹ 44 more variables: nktn <dbl>, nkth <dbl>, nkdmm <dbl>, khac <dbl>,
## #   gcs <dbl>, hatb <dbl>, hatt <dbl>, hattr <dbl>, mach <dbl>, nhietdo <dbl>,
## #   nhiptho <dbl>, sp02 <dbl>, fi02 <dbl>, PFr <dbl>, caymau <dbl>, rltk <dbl>,
## #   rlhh <dbl>, rltm <dbl>, rlthan <dbl>, rlgan <dbl>, rldongmau <dbl>,
## #   suydacoquan <dbl>, thongkhich <dbl>, locmau <dbl>, bc <dbl>, hct <dbl>,
## #   tc <dbl>, bilitp <dbl>, lactate <dbl>, albumin <dbl>, lar <dbl>, pct <dbl>,
## #   na <dbl>, k <dbl>, ast <dbl>, alt <dbl>, ure <dbl>, crea <dbl>, …

Các test kiểm định phân bố chuẩn

1. Vẽ biểu đồ Histogram và QQ plot

A. Histogram

library(ggplot2)
hist.exam = ggplot(data = lar1,aes(albumin))+
  geom_histogram(aes(y=after_stat(density)),color="black",fill="white", bins=30)+
  stat_function(fun = dnorm,
                args = list(mean = mean(lar1$albumin,na.rm = TRUE),
                            sd = sd(lar1$albumin,na.rm = TRUE)),
                color ='blue',linewidth = 1)
hist.exam

Giải thích: Phải tạo cái package là ggplot2 –> gọi tên package là library(ggplot2) Đặt tên cho biểu đồ hist.exam Dữ liệu là lar1, cái aes(y=..density..) là phối hợp biểu đồ density plot Color là viền cột. fill là màu trong ô Stat_function là để vẽ đường liên tục

B. QQ plot

library(ggpubr)
QQplot = ggqqplot(lar1, x = "albumin")
QQplot

2. Sử dụng kiểm định Shapiro-Wilk và Sử dụng kiểm định Anderson-Darling

2.1. Sử dụng kiểm định Shapiro-Wil

shapiro.test(lar1$albumin)
## 
##  Shapiro-Wilk normality test
## 
## data:  lar1$albumin
## W = 0.97703, p-value = 0.07236

2.2. Sử dụng kiểm định Anderson-Darling

library(nortest)
ad.test(lar1$albumin)
## 
##  Anderson-Darling normality test
## 
## data:  lar1$albumin
## A = 0.27772, p-value = 0.6451