Scribd
Upload
8 views

Copy of File Code BTL

Uploaded by

maiquockhoi426
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
8 views

Copy of File Code BTL

Uploaded by

maiquockhoi426
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 2

cereal_data <- read.csv("C:/Users/PC/Downloads/cereal.

csv")
head(cereal_data,15)

library(dplyr)#Tai package "dplyr"


new_data <- cereal_data %>% select(-c("name","mfr","shelf","weight","cups")) # loc bien
head(new_data,25)

library(questionr) #Tai package "questionr" cho lenh freg.na


freq.na(new_data) #Lenh dung de kiem tra tong va ti le cua du lieu khuyet

new_data[new_data == -1] <- 0

num_columns <- sapply(new_data, is.numeric)


num_data <- new_data[, num_columns]

summary_data <- sapply(num_data, function(x) {


c(
Mean = mean(x), #trung binh
SD = sd(x), #Do lech chuan
Min = min(x),#Gia tri nho nhat
Max = max(x), #Gia tri lon nhat
PV1 = quantile(x, 0.25), #phan vi 1
Median = median(x), #Trung vi
PV3 = quantile(x, 0.75) #Phan vi 3

)
})
t(summary_data)

result <- table(new_data$type)


names(result) <- paste("type:", names(result))
print(result)

hist(new_data$rating,main="Consumer
Ratings",xlab="Ratings",ylab="freq",col="lightgreen",label=T,ylim = c(0,25),xlim = c(0,100))
boxplot(rating~type,new_data,main="Boxplot of Rating for
Type",col=c("yellow","lightblue"),ylim= c(20,100))

par(mfrow=c(1,3))
plot(rating~calories,new_data,main="rating & calories",pch=17)
plot(rating~protein,new_data,main="rating & protein",pch=17)
plot(rating~fat,new_data,main="rating & fat",pch=17)
plot(rating~sodium,new_data,main="rating & sodium",pch=17)
plot(rating~fiber,new_data,main="rating & fiber",pch=17)
plot(rating~carbo,new_data,main="rating & carbo",pch=17)
plot(rating~sugars,new_data,main="rating & sugars",pch=17)
plot(rating~potass,new_data,main="rating & potass",pch=17)
plot(rating~vitamins,new_data,main="rating & vitamins",pch=17)

par(mfrow=c(1,1))
library(corrplot)
corrplot(cor(new_data[,2:11]),method="number")

model_1<-lm(rating~type+calories+protein+fat+sodium+fiber+carbo+sugars+potass+
vitamins, data=new_data)
summary(model_1)

model_2<-
lm(rating~calories+protein+fat+sodium+fiber+carbo+sugars+potass+vitamins,data=new_data)
summary(model_2)

par(mfrow=c(2,2))
plot(model_2)

library(car)
dw_test <- durbinWatsonTest(model_2)
print(dw_test)

library(car)
vif(model_2)

You might also like