Creación de Matrices

Matriz Y

##Generador de Matriz
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

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 Cálculo de X’X (Sigma matriz)

Multiplicación matricial debe ser “%*%”

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
f_matriz <- function(matriz)(t(matriz)%*% matriz)
matriz_x2x<-f_matriz(matriz_x)
print(matriz_x2x)
##     Cte  X1  X2
## Cte   5  25  50
## X1   25 141 262
## X2   50 262 510