CREACION DE MATRICES

Matriz y

## GENERACION DE LA MATRIZ Y
matriz_y<-matrix(data =  c(20, 30, 36, 24,40 ),
                    nrow = 5,
                    ncol = 1, byrow = TRUE)
colnames(matriz_y)<-("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 ))

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
print(matriz_Xx)
cte  x1  x2

cte 5 25 50 x1 25 141 262 x2 50 262 510

sigma_matriz<-function(matriz)(t(matriz)%*%matriz)
matriz_xx2<-sigma_matriz(matriz_x)
print(matriz_xx2)
##     cte  x1  x2
## cte   5  25  50
## x1   25 141 262
## x2   50 262 510