hoiquy <- function(data, x, y) {
  model <- lm(formula = paste(y, "~", x), data = data)
  summary(model)
  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))
 # Tìm dộ tin cậy 95%
  confint(model) 
  re <- resid(model)
  # Giả thiết 1: Sai số ngẫu nhiên có phân phối chuẩn
  shapiro.test(re)
  # Giả thiết 2: Kỳ vọng của sai số ngẫu nhiên tại mỗi giá trị bằng 0
  t.test(re, mu = 0)
  # Giả thiết 3: Phương sai của sai số ngẫu nhiên không đổi
  ncvTest(model)
  
}
LS0tDQp0aXRsZTogIlVudGl0bGVkIg0KYXV0aG9yOiAia2nhu4NtIHRyYSINCmRhdGU6ICIyMDIzLTA3LTI1Ig0Kb3V0cHV0OiANCiAgaHRtbF9kb2N1bWVudDoNCiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlDQogICAgY29kZV9mb2xkaW5nOiBoaWRlDQogICAgaGlnaGxpZ2h0OiBweWdtZW50cw0KICAgICMgbnVtYmVyX3NlY3Rpb25zOiB5ZXMNCiAgICB0aGVtZTogImZsYXRseSINCiAgICB0b2M6IFRSVUUNCiAgICB0b2NfZmxvYXQ6IFRSVUUNCi0tLQ0KDQpgYGB7ciBzZXR1cCwgaW5jbHVkZT1GQUxTRX0NCmtuaXRyOjpvcHRzX2NodW5rJHNldChlY2hvID0gVFJVRSkNCmBgYA0KDQoNCmBgYHtyfQ0KaG9pcXV5IDwtIGZ1bmN0aW9uKGRhdGEsIHgsIHkpIHsNCiAgbW9kZWwgPC0gbG0oZm9ybXVsYSA9IHBhc3RlKHksICJ+IiwgeCksIGRhdGEgPSBkYXRhKQ0KICBzdW1tYXJ5KG1vZGVsKQ0KICByZXR1cm4oc3VtbWFyeShtb2RlbCkpIA0KICAjIEtp4buDbSDEkeG7i25oIHTGsMahbmcgcXVhbg0KICBjb3IudGVzdChkYXRhLCB4KQ0KICAjIE3hu5l0IHPhu5EgxJHhu5MgdGjhu4sgdHJvbmcgcGjDom4gdMOtY2ggdMawxqFuZyBxdWFuDQogIHJhcU1hdHJpeCA8LSBjb3IoZGF0YSAlPiUgc2VsZWN0KHksIHgpKQ0KICMgVMOsbSBk4buZIHRpbiBj4bqteSA5NSUNCiAgY29uZmludChtb2RlbCkgDQogIHJlIDwtIHJlc2lkKG1vZGVsKQ0KICAjIEdp4bqjIHRoaeG6v3QgMTogU2FpIHPhu5Egbmfhuqt1IG5oacOqbiBjw7MgcGjDom4gcGjhu5FpIGNodeG6qW4NCiAgc2hhcGlyby50ZXN0KHJlKQ0KICAjIEdp4bqjIHRoaeG6v3QgMjogS+G7syB24buNbmcgY+G7p2Egc2FpIHPhu5Egbmfhuqt1IG5oacOqbiB04bqhaSBt4buXaSBnacOhIHRy4buLIGLhurFuZyAwDQogIHQudGVzdChyZSwgbXUgPSAwKQ0KICAjIEdp4bqjIHRoaeG6v3QgMzogUGjGsMahbmcgc2FpIGPhu6dhIHNhaSBz4buRIG5n4bqrdSBuaGnDqm4ga2jDtG5nIMSR4buVaQ0KICBuY3ZUZXN0KG1vZGVsKQ0KICANCn0NCmBgYA==