CREACION DE MATRICES

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

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))

print(matriz_x)
##      [,1] [,2] [,3]
## [1,]    1    4   10
## [2,]    1    3    8
## [3,]    1    6   11
## [4,]    1    4    9
## [5,]    1    8   12

producto de matrices

matriz_xx<- t(matriz_x) %*% matriz_x

print(matriz_xx)
##      [,1] [,2] [,3]
## [1,]    5   25   50
## [2,]   25  141  262
## [3,]   50  262  510
sigma_matiz<- function(matriz)(t(matriz)%*%matriz)
matriz_xx2<-sigma_matiz(matriz_x)
print(matriz_xx2)
##      [,1] [,2] [,3]
## [1,]    5   25   50
## [2,]   25  141  262
## [3,]   50  262  510