creaccion de matrices en R

Matriz Y

#generando matriz y
Matriz_Y<- matrix(data = c(20,30,36,24,40),
                  nrow = 5,
                  ncol = 1, byrow = TRUE)
colnames(Matriz_Y)<-c("y")
print(Matriz_Y)
##       y
## [1,] 20
## [2,] 30
## [3,] 36
## [4,] 24
## [5,] 40

Matriz X

#generando matriz X
Matriz_X<- cbind(rep(1,5), matrix(data = c(4,10,3,8,6,11,4,9,8,12),
                                   nrow = 5,
                                   ncol = 2,byrow = TRUE))
colnames(Matriz_X)<-c("Cte", "x1", "x2")
print(Matriz_X)
##      Cte x1 x2
## [1,]   1  4 10
## [2,]   1  3  8
## [3,]   1  6 11
## [4,]   1  4  9
## [5,]   1  8 12

Producto de Matrices Calculo de X’x(Sigma Matriz)

Matriz_XX<-t(Matriz_X)%*%(Matriz_X)
#Mostrando
print(Matriz_XX)
##     Cte  x1  x2
## Cte   5  25  50
## x1   25 141 262
## x2   50 262 510
Sigma_matriz<- function(matriz){t(Matriz_X)%*%(Matriz_X)}
Matriz_XX2<-Sigma_matriz(Matriz_X)
#Mostrando
print(Matriz_XX2)
##     Cte  x1  x2
## Cte   5  25  50
## x1   25 141 262
## x2   50 262 510