Creación de Matrices

Matriz y

#Generación de la matriz y
matriz_y<- matrix(data = c(30,20,36,24,40),
                  nrow = 5,
                  ncol = 1,byrow =TRUE)
colnames(matriz_y)<- c("Y")
print(matriz_y)
##       Y
## [1,] 30
## [2,] 20
## [3,] 36
## [4,] 24
## [5,] 40

Matriz X

#Generación de la 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

Productos de matrices Cálculo 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