1. generando matriz X’X
options(scipen = 9999999)
matriz_xx<-matrix(data = c(25,4586,2018,
                           4586,1030398,364545,
                           2018,364545,204312),nrow = 3,ncol = 3,
                  byrow = TRUE)
print(matriz_xx)
##      [,1]    [,2]   [,3]
## [1,]   25    4586   2018
## [2,] 4586 1030398 364545
## [3,] 2018  364545 204312

2.Generando matriz X’Y

matriz_xy<-matrix(data = c(55331,12524626,4374490),nrow = 3,ncol = 1,
                  byrow = TRUE)
print(matriz_xy)
##          [,1]
## [1,]    55331
## [2,] 12524626
## [3,]  4374490

3.calculando la matriz inversa de X’X

inv_matriz_xx<-solve(matriz_xx)
print(inv_matriz_xx)
##              [,1]             [,2]             [,3]
## [1,]  0.397982654 -0.0010321198463 -0.0020893284104
## [2,] -0.001032120  0.0000053085599  0.0000007224679
## [3,] -0.002089328  0.0000007224679  0.0000242418099

4.Calculando el estimador de Parametros β

beta<-inv_matriz_xx%*%matriz_xy
print(beta)
##             [,1]
## [1,] -45.8830775
## [2,]  12.5399333
## [3,]  -0.5104348