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