hoiquy <- function(data, x, y) {
model <- lm(formula = paste(y_col, "~", x_col), data = data)
return(summary(model))
# Kiểm định tương quan
cor.test(data, x)
# Một số đồ thị trong phân tích tương quan
raqMatrix <- cor(data %>% select(y, x))
}