MATRIZ Y ## CREACIÃ’N DE MATRICES

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