Re: 오늘의 마지막 문제. 다중회귀 이론 수업 맨 아래쪽에 나오는 연립방정식 문제를...

[19기_이동근]

# 필요한 패키지 로드
library(MASS)

# 다중 회귀 분석을 위한 reg 함수 정의
reg <- function(y, x) { 
  x <- as.matrix(x)
  x <- cbind(Intercept = 1, x)
  b <- ginv(t(x) %*% x) %*% t(x) %*% y   
  colnames(b) <- "estimate"
  print(b)
}

# 주어진 연립일차방정식의 행렬 A와 벡터 b
A <- matrix(c(1, 3, 2, 4, -1, 1), nrow = 2, byrow = TRUE)
A
b <- matrix(c(5, 4), nrow = 2, byrow = TRUE)
b

# 함수를 사용하여 해를 구함
reg(b, A)

###################### cbind 줄 제외

# 필요한 패키지 로드
library(MASS)

# 다중 회귀 분석을 위한 reg 함수 정의
reg <- function(y, x) { 
  x <- as.matrix(x)
  b <- ginv(t(x) %*% x) %*% t(x) %*% y   
  colnames(b) <- "estimate"
  print(b)
}

# 주어진 연립일차방정식의 행렬 A와 벡터 b
A <- matrix(c(1, 3, 2, 4, -1, 1), nrow = 2, byrow = TRUE)
A
b <- c(5, 4)
b

# 함수를 사용하여 해를 구함
reg(b, A)