Creación de Matrices

Matriz Y

matriz_D<-matrix(data = c(20,30,36,24,40),
                 nrow = 5,
                 ncol = 1,byrow = TRUE)
colnames(matriz_D)<-c("Y")
print(matriz_D)
##       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_W<-t(matriz_X)%*%matriz_X
print(Matriz_W)
##     cte  x1  x2
## cte   5  25  50
## x1   25 141 262
## x2   50 262 510
Sigma_Matriz<-function(matriz)(t(matriz)%*%matriz)
Matriz_Z<-Sigma_Matriz(matriz_X)
print(Matriz_Z)
##     cte  x1  x2
## cte   5  25  50
## x1   25 141 262
## x2   50 262 510