Bài tập lượng giá 27/10

MỞ DATA

library(readxl)
datacnls <- read_excel("C:/PTSL YHDP16/datacnls.xls")
View(datacnls)

TẠO BIẾN SỐ MỚI

attach(datacnls)

Biến nhóm tuổi mẹ

datacnls$ntme[datacnls$tuoime <25] = "<25"
datacnls$ntme[datacnls$tuoime >=25 & datacnls$tuoime <30] = "25-29"
datacnls$ntme[datacnls$tuoime >=30 & datacnls$tuoime <35] = "30-34"
datacnls$ntme[datacnls$tuoime >=35] = ">=35"

Biến nhóm khám thai

datacnls$nkt[datacnls$solankt <8] = "<8 lần"
datacnls$nkt[datacnls$solankt >=8 & datacnls$solankt <=10] = "8-10 lần"
datacnls$nkt[datacnls$solankt >10] = ">10 lần"

Biến sinh non

datacnls$sinhnon[datacnls$tuoithai<37]<- "Co"
datacnls$sinhnon[datacnls$tuoithai>=37] <- "Khong"

Biến nhẹ cân

datacnls$nhecan[datacnls$cnls<2500]<- "Co"
datacnls$nhecan[datacnls$cnls>=2500] <- "Khong"

BẢNG 1

table1(~ ntme + luongoi + cha + nkt, data = datacnls)

Bảng 2

table1(~ gioi + cnls + sinhnon + nhecan, data = datacnls)

Bảng 3

Tính cân nặng lúc sinh trung bình theo nhóm tuổi

datacnls %>% select(cnls, ntme) %>% group_by(ntme) %>% summarise_all(mean)

datacnls %>% select(cnls, ntme) %>% group_by(ntme) %>% summarise_all(sd)

Kiểm định anova được sử dụng do phân phối chuẩn

a <- lm(cnls~ntme)
anova(a)
summary(a)

Tính cân nặng lúc sinh trung bình theo lượng ối bất thường

datacnls %>% select(cnls, luongoi) %>% group_by(luongoi) %>% summarise_all(mean)
datacnls %>% select(cnls, luongoi) %>% group_by(luongoi) %>% summarise_all(sd)

Kiểm định t.test không bắt cặp

t.test(cnls~luongoi)

Tính cân nặng lúc sinh trung bình theo cao huyết áp

datacnls %>% select(cnls, cha) %>% group_by(cha) %>% summarise_all(mean)

datacnls %>% select(cnls, cha) %>% group_by(cha) %>% summarise_all(sd)

Kiểm định t.test không bắt cặp

t.test(cnls~cha)

Tính cân nặng lúc sinh trung bình theo nhóm khám thai

datacnls %>% select(cnls, nkt) %>% group_by(nkt) %>% summarise_all(mean)

datacnls %>% select(cnls, nkt) %>% group_by(nkt) %>% summarise_all(sd)

pt anova

b <- lm(cnls~nkt)
anova(b)
summary(b)