Rのサンプル airquality を使用
dat <- na.omit(airquality)
head(dat)

fit <- lm(Solar.R ~ Wind+Temp, data=dat)
summary(fit)

下記の係数の推定値を求めてみましょう
summary(fit)$coef[1:3]

\(\bf b=
\left(
\begin{array}{ccc}
-49.8135134 \\
0.6478038 \\
2.9331301
\end{array}
\right)
\)
データフレームを行列に変換して変数行列 \(\bf X\) を作成
(切片用の列を追加するのを忘れずに)
dat$b <- c(rep(1, dim(dat)[1]))
dat2 <- dat[ ,c("b", "Wind", "Temp")]
mat <- as.matrix(dat2)
\(\bf b = (X^TX)^{-1}X^Ty\)より
solve(
t(mat) %*% mat) %*%
t(mat) %*%
dat$Solar.R

ダメ出し 間違い、分かりにくい部分などのご意見をお待ちします