…{r} # Ngay 3 : Mo hinh hoi quy tuyen tinh setwd(“E:/cuong_dataR/test”) ob = read.csv(“Obesity data.csv”, header = T) attach(ob) head(ob) names(ob) # Chon bien lien quan dat = ob[, c(“gender”, “age”, “height”, “weight”, “bmi”, “pcfat”)] library(GGally); library(ggplot2); library(visreg); library(tabe1); library(BMA) table1(~age + height + weight + bmi + pcfat | gender, data=ob) ggpairs(dat, mapping = aes(color = gender)) # Task 2 : Mo hinh hoi quy tuyen tinh don gian # Khac biet giua nam va nu m = lm(pcfat ~ gender, data=ob) summary(m) visreg(m) # Khac biet giua nam va nu hieu chinh cho BMI m = lm(pcfat ~ gender + bmi, data=ob) summary(m) # Tuong quan voi BMI m1 = lm(pcfat ~ bmi, data=ob) m2 = lm(pcfat ~ bmi + I(bmi^2), data=ob) m3 = lm(pcfat ~ bmi + I(bmi^2) + I(bmi^3), data=ob) anova(m1, m2, m3) par(mfrow=c(1,3)) visreg(m1) visreg(m2) visreg(m3) …