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==