koefisien_regresi <- function(x1, x2, y) {
# Membuat vektor x0 yang berisi 1 sepanjang panjang x1
x0 <- rep(1, length(x1))
# Menggabungkan x0, x1, dan x2 menjadi matriks X
X <- as.matrix(cbind(x0, x1, x2))
# Mengubah vektor y menjadi matriks
y <- as.matrix(y)
# Menghitung koefisien regresi
beta_duga <- solve(t(X) %*% X) %*% t(X) %*% y
return(beta_duga)
}
# Contoh data
x1 <- c(10, 20, 30, 40, 40, 50, 50, 60, 60)
x2 <- c(15, 25, 35, 45, 55, 65, 75, 85, 95)
y <- c(18.3, 23, 25.6, 33, 40.7, 44, 54.5, 55, 60.8)
koefisien <- koefisien_regresi(x1, x2, y)
koefisien
## [,1]
## x0 10.2136150
## x1 -0.2416901
## x2 0.7070423